-
Convex neural network synthesis for robustness in the 1-norm
Authors:
Ross Drummond,
Chris Guiver,
Matthew C. Turner
Abstract:
With neural networks being used to control safety-critical systems, they increasingly have to be both accurate (in the sense of matching inputs to outputs) and robust. However, these two properties are often at odds with each other and a trade-off has to be navigated. To address this issue, this paper proposes a method to generate an approximation of a neural network which is certifiably more robu…
▽ More
With neural networks being used to control safety-critical systems, they increasingly have to be both accurate (in the sense of matching inputs to outputs) and robust. However, these two properties are often at odds with each other and a trade-off has to be navigated. To address this issue, this paper proposes a method to generate an approximation of a neural network which is certifiably more robust. Crucially, the method is fully convex and posed as a semi-definite programme. An application to robustifying model predictive control is used to demonstrate the results. The aim of this work is to introduce a method to navigate the neural network robustness/accuracy trade-off.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Mapping back and forth between model predictive control and neural networks
Authors:
Ross Drummond,
Pablo R Baldivieso-Monasterios,
Giorgio Valmorbida
Abstract:
Model predictive control (MPC) for linear systems with quadratic costs and linear constraints is shown to admit an exact representation as an implicit neural network. A method to "unravel" the implicit neural network of MPC into an explicit one is also introduced. As well as building links between model-based and data-driven control, these results emphasize the capability of implicit neural networ…
▽ More
Model predictive control (MPC) for linear systems with quadratic costs and linear constraints is shown to admit an exact representation as an implicit neural network. A method to "unravel" the implicit neural network of MPC into an explicit one is also introduced. As well as building links between model-based and data-driven control, these results emphasize the capability of implicit neural networks for representing solutions of optimisation problems, as such problems are themselves implicitly defined functions.
△ Less
Submitted 18 April, 2024;
originally announced April 2024.
-
Tuning the feedback controller gains is a simple way to improve autonomous driving performance
Authors:
Wenyu Liang,
Pablo R. Baldivieso,
Ross Drummond,
Donghwan Shin
Abstract:
Typical autonomous driving systems are a combination of machine learning algorithms (often involving neural networks) and classical feedback controllers. Whilst significant progress has been made in recent years on the neural network side of these systems, only limited progress has been made on the feedback controller side. Often, the feedback control gains are simply passed from paper to paper wi…
▽ More
Typical autonomous driving systems are a combination of machine learning algorithms (often involving neural networks) and classical feedback controllers. Whilst significant progress has been made in recent years on the neural network side of these systems, only limited progress has been made on the feedback controller side. Often, the feedback control gains are simply passed from paper to paper with little re-tuning taking place, even though the changes to the neural networks can alter the vehicle's closed loop dynamics. The aim of this paper is to highlight the limitations of this approach; it is shown that re-tuning the feedback controller can be a simple way to improve autonomous driving performance. To demonstrate this, the PID gains of the longitudinal controller in the TCP autonomous vehicle algorithm are tuned. This causes the driving score in CARLA to increase from 73.21 to 77.38, with the results averaged over 16 driving scenarios. Moreover, it was observed that the performance benefits were most apparent during challenging driving scenarios, such as during rain or night time, as the tuned controller led to a more assertive driving style. These results demonstrate the value of developing both the neural network and feedback control policies of autonomous driving systems simultaneously, as this can be a simple and methodical way to improve autonomous driving system performance and robustness.
△ Less
Submitted 7 February, 2024;
originally announced February 2024.
-
When are selector control strategies optimal for constrained monotone systems?
Authors:
Hamed Taghavian,
Ross Drummond,
Mikael Johansson
Abstract:
This paper considers optimal control problems defined by a monotone dynamical system, a monotone cost, and monotone constraints. We identify families of such problems for which the optimal solution is bang-ride, i.e., always operates on the constraint boundaries, and prove that the optimal policy switches between a finite number of state feedback controllers. This motivates the use of simpler poli…
▽ More
This paper considers optimal control problems defined by a monotone dynamical system, a monotone cost, and monotone constraints. We identify families of such problems for which the optimal solution is bang-ride, i.e., always operates on the constraint boundaries, and prove that the optimal policy switches between a finite number of state feedback controllers. This motivates the use of simpler policies, such as selector control, that can be designed without perfect models and full state measurements. The approach is successfully applied to several variations of the health-aware fast charging problem for lithium-ion batteries.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
Logarithmically Completely Monotonic Rational Functions
Authors:
Hamed Taghavian,
Ross Drummond,
Mikael Johansson
Abstract:
This paper studies the class of logarithmically completely monotonic (LCM) functions. These functions play an important role in characterising externally positive linear systems which find applications in important control problems such as non-overshooting reference tracking. Conditions are proposed to ensure a rational function is LCM, a result that enables the known space of linear continuous-ti…
▽ More
This paper studies the class of logarithmically completely monotonic (LCM) functions. These functions play an important role in characterising externally positive linear systems which find applications in important control problems such as non-overshooting reference tracking. Conditions are proposed to ensure a rational function is LCM, a result that enables the known space of linear continuous-time externally positive systems to be enlarged and an efficient and optimal pole-placement procedure for the monotonic tracking controller synthesis problem to be developed. The presented conditions are shown to be less conservative than existing approaches whilst being computationally tractable.
△ Less
Submitted 17 February, 2023;
originally announced February 2023.
-
Constrained optimal control of monotone systems with applications to battery fast-charging
Authors:
Ross Drummond,
Nicola E. Courtier,
David A. Howey,
Luis D. Couto,
Chris Guiver
Abstract:
Enabling fast charging for lithium ion batteries is critical to accelerating the green energy transition. As such, there has been significant interest in tailored fast-charging protocols computed from the solutions of constrained optimal control problems. Here, we derive necessity conditions for a fast charging protocol based upon monotone control systems theory.
Enabling fast charging for lithium ion batteries is critical to accelerating the green energy transition. As such, there has been significant interest in tailored fast-charging protocols computed from the solutions of constrained optimal control problems. Here, we derive necessity conditions for a fast charging protocol based upon monotone control systems theory.
△ Less
Submitted 30 October, 2023; v1 submitted 18 January, 2023;
originally announced January 2023.
-
Bounding the difference between model predictive control and neural networks
Authors:
Ross Drummond,
Stephen R. Duncan,
Matthew C. Turner,
Patricia Pauli,
Frank Allgöwer
Abstract:
There is a growing debate on whether the future of feedback control systems will be dominated by data-driven or model-driven approaches. Each of these two approaches has their own complimentary set of advantages and disadvantages, however, only limited attempts have, so far, been developed to bridge the gap between them. To address this issue, this paper introduces a method to bound the worst-case…
▽ More
There is a growing debate on whether the future of feedback control systems will be dominated by data-driven or model-driven approaches. Each of these two approaches has their own complimentary set of advantages and disadvantages, however, only limited attempts have, so far, been developed to bridge the gap between them. To address this issue, this paper introduces a method to bound the worst-case error between feedback control policies based upon model predictive control (MPC) and neural networks (NNs). This result is leveraged into an approach to automatically synthesize MPC policies minimising the worst-case error with respect to a NN. Numerical examples highlight the application of the bounds, with the goal of the paper being to encourage a more quantitative understanding of the relationship between data-driven and model-driven control.
△ Less
Submitted 13 April, 2022;
originally announced April 2022.
-
Regional Stability Analysis of Transitional Fluid Flows
Authors:
Leonardo F. Toso,
Ross Drummond,
Stephen R. Duncan
Abstract:
A method to bound the maximum energy perturbation for which regional stability of transitional fluid flow models can be guaranteed is introduced. The proposed method exploits the fact that the fluid model's nonlinearities are both lossless and locally bounded and uses the axes lengths of the ellipsoids for the trajectory set containment as variables in the stability conditions. Compared to existin…
▽ More
A method to bound the maximum energy perturbation for which regional stability of transitional fluid flow models can be guaranteed is introduced. The proposed method exploits the fact that the fluid model's nonlinearities are both lossless and locally bounded and uses the axes lengths of the ellipsoids for the trajectory set containment as variables in the stability conditions. Compared to existing approaches, the proposed method leads to an average increase in the maximum allowable energy perturbation of 29% for the Waleffe-KimHamilton (WKH) shear flow model and of 38% for the 9-state reduced model of Couette flow.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
A mixture model for determining SARS-Cov-2 variant composition in pooled samples
Authors:
Renan Valieris,
Rodrigo Drummond,
Alexandre Defelicibus,
Emannuel Dias-Neto,
Rafael A. Rosales,
Israel Tojal da Silva
Abstract:
Despite of the fast development of highly effective vaccines to control the current COVID$-$19 pandemic, the unequal distribution and availability of these vaccines worldwide and the number of people infected in the world lead to the continuous emergence of SARS-CoV-2 (Severe Acute Respiratory Syndrome coronavirus 2) variants of concern. It is likely that real-time genomic surveillance will be con…
▽ More
Despite of the fast development of highly effective vaccines to control the current COVID$-$19 pandemic, the unequal distribution and availability of these vaccines worldwide and the number of people infected in the world lead to the continuous emergence of SARS-CoV-2 (Severe Acute Respiratory Syndrome coronavirus 2) variants of concern. It is likely that real-time genomic surveillance will be continuously needed as an unceasing monitoring tool, necessary to follow the spillover of the disease spread and the evolution of the virus. In this context, new genomic variants of SARS-CoV-2 that may emerge as a response to selective pressure, including variants refractory to current vaccines, makes genomic surveillance programs tools of utmost importance. Here propose a statistical model for the estimation of the relative frequencies of SARS-CoV-2 variants in pooled samples. This model is built by considering a previously defined selection of genomic polymorphisms that characterize SARS-CoV-2 variants. The methods described here support both raw sequencing reads for polymorphisms-based markers calling and predefined markers in the VCF format. Results obtained by using simulated data show that our method is quite effective in recovering the correct variant proportions. Further, results obtained by considering longitudinal data from wastewater samples of two locations in Switzerland agree well with those describing the epidemiological evolution of COVID-19 variants in clinical samples of these locations. Our results show that the described method can be a valuable tool for tracking the proportions of SARS-CoV-2 variants.
△ Less
Submitted 3 October, 2021;
originally announced October 2021.
-
Cell-Level State of Charge Estimation for Battery Packs Under Minimal Sensing
Authors:
Dong Zhang,
Luis D. Couto,
Ross Drummond,
Shashank Sripad,
Venkatasubramanian Viswanathan
Abstract:
This manuscript presents an algorithm for individual Lithium-ion (Li-ion) battery cell state of charge (SOC) estimation in a large-scale battery pack under minimal sensing, where only pack-level voltage and current are measured. For battery packs consisting of up to thousands of cells in electric vehicle or stationary energy storage applications, it is desirable to estimate individual cell SOCs wi…
▽ More
This manuscript presents an algorithm for individual Lithium-ion (Li-ion) battery cell state of charge (SOC) estimation in a large-scale battery pack under minimal sensing, where only pack-level voltage and current are measured. For battery packs consisting of up to thousands of cells in electric vehicle or stationary energy storage applications, it is desirable to estimate individual cell SOCs without cell local measurements in order to reduce sensing costs. Mathematically, pure series connected cells yield dynamics given by ordinary differential equations under classical full voltage sensing. In contrast, parallel--series connected battery packs are evidently more challenging because the dynamics are governed by a nonlinear differential--algebraic equations (DAE) system. The majority of the conventional studies on SOC estimation for battery packs benefit from idealizing the pack as a lumped single cell which ultimately lose track of cell-level conditions and are blind to potential risks of cell-level over-charge and over-discharge. This work explicitly models a battery pack with high fidelity cell-by-cell resolution based on the interconnection of single cell models, and examines the observability of cell-level state with only pack-level measurements. A DAE-based state observer with linear output error injection is formulated, where the individual cell SOC and current can be reconstructed from minimal number of pack sensing. The mathematically guaranteed asymptotic convergence of differential and algebraic state estimates is established by considering local Lipschitz continuity property of system nonlinearities. Simulation results for Graphite/NMC cells illustrate convergence for cell SOCs, currents, and voltages.
△ Less
Submitted 16 September, 2021;
originally announced September 2021.
-
Ariel: Enabling planetary science across light-years
Authors:
Giovanna Tinetti,
Paul Eccleston,
Carole Haswell,
Pierre-Olivier Lagage,
Jérémy Leconte,
Theresa Lüftinger,
Giusi Micela,
Michel Min,
Göran Pilbratt,
Ludovic Puig,
Mark Swain,
Leonardo Testi,
Diego Turrini,
Bart Vandenbussche,
Maria Rosa Zapatero Osorio,
Anna Aret,
Jean-Philippe Beaulieu,
Lars Buchhave,
Martin Ferus,
Matt Griffin,
Manuel Guedel,
Paul Hartogh,
Pedro Machado,
Giuseppe Malaguti,
Enric Pallé
, et al. (293 additional authors not shown)
Abstract:
Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths.…
▽ More
Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution.
△ Less
Submitted 10 April, 2021;
originally announced April 2021.
-
Reduced-Order Neural Network Synthesis with Robustness Guarantees
Authors:
Ross Drummond,
Mathew C. Turner,
Stephen R. Duncan
Abstract:
In the wake of the explosive growth in smartphones and cyberphysical systems, there has been an accelerating shift in how data is generated away from centralised data towards on-device generated data. In response, machine learning algorithms are being adapted to run locally on board, potentially hardware limited, devices to improve user privacy, reduce latency and be more energy efficient. However…
▽ More
In the wake of the explosive growth in smartphones and cyberphysical systems, there has been an accelerating shift in how data is generated away from centralised data towards on-device generated data. In response, machine learning algorithms are being adapted to run locally on board, potentially hardware limited, devices to improve user privacy, reduce latency and be more energy efficient. However, our understanding of how these device orientated algorithms behave and should be trained is still fairly limited. To address this issue, a method to automatically synthesize reduced-order neural networks (having fewer neurons) approximating the input/output mapping of a larger one is introduced. The reduced-order neural network's weights and biases are generated from a convex semi-definite programme that minimises the worst-case approximation error with respect to the larger network. Worst-case bounds for this approximation error are obtained and the approach can be applied to a wide variety of neural networks architectures. What differentiates the proposed approach to existing methods for generating small neural networks, e.g. pruning, is the inclusion of the worst-case approximation error directly within the training cost function, which should add robustness. Numerical examples highlight the potential of the proposed approach. The overriding goal of this paper is to generalise recent results in the robustness analysis of neural networks to a robust synthesis problem for their weights and biases.
△ Less
Submitted 7 December, 2021; v1 submitted 18 February, 2021;
originally announced February 2021.
-
Structural Identifiability of a Pseudo-2D Li-ion Battery Electrochemical Model
Authors:
Ross Drummond,
Stephen R. Duncan
Abstract:
Growing demand for fast charging and optimised battery designs is fuelling significant interest in electrochemical models of Li-ion batteries. However, estimating parameter values for these models remains a major challenge. In this paper, a structural identifiability analysis was applied to a pseudo-2D Li-ion electrochemical battery model that can be considered as a linearised and decoupled form o…
▽ More
Growing demand for fast charging and optimised battery designs is fuelling significant interest in electrochemical models of Li-ion batteries. However, estimating parameter values for these models remains a major challenge. In this paper, a structural identifiability analysis was applied to a pseudo-2D Li-ion electrochemical battery model that can be considered as a linearised and decoupled form of the benchmark Doyle-Fuller-Newman model. From an inspection of the impedance function, it was shown that this model is uniquely parametrised by 21 parameters, being combinations of the electrochemical parameters like the conductivities and diffusion coefficients. The well-posedness of the parameter estimation problem with these parameters was then established. This result could lead to more realistic predictions about the internal state of the battery by identifying the parameter set that can be uniquely identified from the data.
△ Less
Submitted 3 December, 2020;
originally announced December 2020.
-
Robust error bounds for quantised and pruned neural networks
Authors:
Jiaqi Li,
Ross Drummond,
Stephen R. Duncan
Abstract:
With the rise of smartphones and the internet-of-things, data is increasingly getting generated at the edge on local, personal devices. For privacy, latency and energy saving reasons, this shift is causing machine learning algorithms to move towards decentralisation with the data and algorithms stored, and even trained, locally on devices. The device hardware becomes the main bottleneck for model…
▽ More
With the rise of smartphones and the internet-of-things, data is increasingly getting generated at the edge on local, personal devices. For privacy, latency and energy saving reasons, this shift is causing machine learning algorithms to move towards decentralisation with the data and algorithms stored, and even trained, locally on devices. The device hardware becomes the main bottleneck for model capability in this set-up, creating a need for slimmed down, more efficient neural networks. Neural network pruning and quantisation are two methods that have been developed for this, with both approaches demonstrating impressive results in reducing the computational cost without sacrificing significantly on model performance. However, the understanding behind these reduction methods remains underdeveloped. To address this issue, a semi-definite program is introduced to bound the worst-case error caused by pruning or quantising a neural network. The method can be applied to many neural network structures and nonlinear activation functions with the bounds holding robustly for all inputs in specified sets. It is hoped that the computed bounds will provide certainty to the performance of these algorithms when deployed on safety-critical systems.
△ Less
Submitted 27 April, 2021; v1 submitted 30 November, 2020;
originally announced December 2020.
-
Resolving Kirchhoff's laws for state-estimator design of Li-ion battery packs connected in parallel
Authors:
Ross Drummond,
Luis D. Couto,
Dong Zhang
Abstract:
A state-space model for Li-ion battery packs with parallel connected cells is introduced. The key feature of this model is an explicit solution to Kirchhoff's laws for parallel connected packs, which expresses the branch currents directly in terms of the model's states, applied current and cell resistances. This avoids the need to solve these equations numerically. To illustrate the potential of t…
▽ More
A state-space model for Li-ion battery packs with parallel connected cells is introduced. The key feature of this model is an explicit solution to Kirchhoff's laws for parallel connected packs, which expresses the branch currents directly in terms of the model's states, applied current and cell resistances. This avoids the need to solve these equations numerically. To illustrate the potential of the proposed model for pack-level control and estimation, a state-estimator is introduced for the nonlinear parallel pack model. By exploiting the system structure seen in the solution to Kirchhoff's laws, algebraic conditions for the observer gains are obtained that guarantee convergence of the estimator's error. Error convergence is demonstrated through an argument based upon Aizerman's conjecture. It is hoped that the insight brought by this model formulation will allow the wealth of results developed for series connected packs to be applied to those with parallel connections.
△ Less
Submitted 30 October, 2020;
originally announced October 2020.
-
Impact of Disturbances on Mixed Traffic Control with Autonomous Vehicles
Authors:
Ross Drummond,
Yang Zheng
Abstract:
This paper investigates the impact of disturbances on controlling an autonomous vehicle to smooth mixed traffic flow in a ring road setup. By exploiting the ring structure of this system, it is shown that velocity perturbations impacting any vehicle on the ring enter an uncontrollable and marginally stable mode defined by the sum of relative vehicle spacings. These disturbances are then integrated…
▽ More
This paper investigates the impact of disturbances on controlling an autonomous vehicle to smooth mixed traffic flow in a ring road setup. By exploiting the ring structure of this system, it is shown that velocity perturbations impacting any vehicle on the ring enter an uncontrollable and marginally stable mode defined by the sum of relative vehicle spacings. These disturbances are then integrated up by the system and cannot be unwound via controlling the autonomous vehicle. In particular, if the velocity disturbances are zero-mean Gaussians, then the traffic flow on the ring will undergo a random walk with the variance growing indefinitely and independently of the control policy applied. In contrast, the impact of acceleration disturbances is benign as these disturbances do no enter the uncontrollable mode, meaning that they can be easily regulated using the autonomous vehicle. Our results support and complement the existing theoretic analysis and field experiments.
△ Less
Submitted 12 August, 2020;
originally announced August 2020.
-
Accelerated Gradient Methods with Memory
Authors:
Ross Drummond,
Stephen Duncan
Abstract:
A set of accelerated first order algorithms with memory are proposed for minimising strongly convex functions. The algorithms are differentiated by their use of the iterate history for the gradient step. The increased convergence rate of the proposed algorithms comes at the cost of robustness, a problem that is resolved by a switching controller based upon adaptive restarting. Several numerical ex…
▽ More
A set of accelerated first order algorithms with memory are proposed for minimising strongly convex functions. The algorithms are differentiated by their use of the iterate history for the gradient step. The increased convergence rate of the proposed algorithms comes at the cost of robustness, a problem that is resolved by a switching controller based upon adaptive restarting. Several numerical examples highlight the benefits of the proposed approach over the fast gradient method. For example, it is shown that these gradient based methods can minimise the Rosenbrock banana function to $7.58 \times 10^{-12}$ in 43 iterations from an initial condition of $(-1,1)$.
△ Less
Submitted 29 August, 2018; v1 submitted 23 May, 2018;
originally announced May 2018.
-
Circuit Synthesis of Electrochemical Supercapacitor Models
Authors:
Ross Drummond,
Shi Zhao,
David A. Howey,
Stephen R. Duncan
Abstract:
This paper is concerned with the synthesis of RC electrical circuits from physics-based supercapacitor models describing conservation and diffusion relationships. The proposed synthesis procedure uses model discretisation, linearisation, balanced model order reduction and passive network synthesis to form the circuits. Circuits with different topologies are synthesized from several physical models…
▽ More
This paper is concerned with the synthesis of RC electrical circuits from physics-based supercapacitor models describing conservation and diffusion relationships. The proposed synthesis procedure uses model discretisation, linearisation, balanced model order reduction and passive network synthesis to form the circuits. Circuits with different topologies are synthesized from several physical models. This work will give greater understanding to the physical interpretation of electrical circuits and will enable the development of more generalised circuits, since the synthesized impedance functions are generated by considering the physics, not from experimental fitting which may ignore certain dynamics.
△ Less
Submitted 30 March, 2016;
originally announced March 2016.
-
The 2003-4 multisite photometric campaign for the Beta Cephei and eclipsing star 16 (EN) Lacertae with an Appendix on 2 Andromedae, the variable comparison star
Authors:
M. Jerzykiewicz,
G. Handler,
J. Daszynska-Daszkiewicz,
A. Pigulski,
E. Poretti,
E. Rodriguez,
P. J. Amado,
Z. Kolaczkowski,
K. Uytterhoeven,
T. N. Dorokhova,
N. I. Dorokhov,
D. Lorenz,
D. Zsuffa,
S. -L. Kim,
P. -O. Bourge,
B. Acke,
J. De Ridder,
T. Verhoelst,
R. Drummond,
A. I. Movchan,
J. -A. Lee,
M. Steslicki,
J. Molenda-Zakowicz,
R. Garrido,
S. -H. Kim
, et al. (4 additional authors not shown)
Abstract:
A multisite photometric campaign for the Beta Cephei and eclipsing variable 16 Lacertae is reported. 749 h of high-quality differential photoelectric Stromgren, Johnson and Geneva time-series photometry were obtained with ten telescopes during 185 nights. After removing the pulsation contribution, an attempt was made to solve the resulting eclipse light curve by means of the computer program EBOP.…
▽ More
A multisite photometric campaign for the Beta Cephei and eclipsing variable 16 Lacertae is reported. 749 h of high-quality differential photoelectric Stromgren, Johnson and Geneva time-series photometry were obtained with ten telescopes during 185 nights. After removing the pulsation contribution, an attempt was made to solve the resulting eclipse light curve by means of the computer program EBOP. Although a unique solution was not obtained, the range of solutions could be constrained by comparing computed positions of the secondary component in the Hertzsprung-Russell diagram with evolutionary tracks.
For three high-amplitude pulsation modes, the uvy and the Geneva UBG amplitude ratios are derived and compared with the theoretical ones for spherical-harmonic degrees l <= 4. The highest degree, l = 4, is shown to be incompatible with the observations. One mode is found to be radial, one is l = 1, while in the remaining case l = 2 or 3.
The present multisite observations are combined with the archival photometry in order to investigate the long-term variation of the amplitudes and phases of the three high-amplitude pulsation modes. The radial mode shows a non-sinusoidal variation on a time-scale of 73 yr. The l = 1 mode is a triplet with unequal frequency spacing, giving rise to two beat-periods, 720.7 d and 29.1 yr. The amplitude and phase of the l = 2 or 3 mode vary on time-scales of 380.5 d and 43 yr.
The light variation of 2 And, one of the comparison stars, is discussed in the Appendix.
△ Less
Submitted 21 August, 2015;
originally announced August 2015.
-
Low-Order Mathematical Modelling of Electric Double Layer Supercapacitors Using Spectral Methods
Authors:
Ross Drummond,
David A. Howey,
Stephen R. Duncan
Abstract:
This work investigates two physics-based models that simulate the non-linear partial differential algebraic equations describing an electric double layer supercapacitor. In one model the linear dependence between electrolyte concentration and conductivity is accounted for, while in the other model it is not. A spectral element method is used to discretise the model equations and it is found that t…
▽ More
This work investigates two physics-based models that simulate the non-linear partial differential algebraic equations describing an electric double layer supercapacitor. In one model the linear dependence between electrolyte concentration and conductivity is accounted for, while in the other model it is not. A spectral element method is used to discretise the model equations and it is found that the error convergence rate with respect to the number of elements is faster compared to a finite difference method. The increased accuracy of the spectral element approach means that, for a similar level of solution accuracy, the model simulation computing time is approximately 50% of that of the finite difference method. This suggests that the spectral element model could be used for control and state estimation purposes. For a typical supercapacitor charging profile, the numerical solutions from both models closely match experimental voltage and current data. However, when the electrolyte is dilute or where there is a long charging time, a noticeable difference between the numerical solutions of the two models is observed. Electrical impedance spectroscopy simulations show that the capacitance of the two models rapidly decreases when the frequency of the perturbation current exceeds an upper threshold.
△ Less
Submitted 28 November, 2014;
originally announced December 2014.
-
Astronomical Sky Quality Near Eureka, in the Canadian High Arctic
Authors:
Eric Steinbring,
William Ward,
James R. Drummond
Abstract:
Nighttime visible-light sky brightness and transparency are reported for the Polar Environment Research Laboratory (PEARL), located on a 610-m high ridge near the Eureka research station, on Ellesmere Island, Canada. Photometry of Polaris obtained in V band with the PEARL All Sky Imager (PASI) over two winters is supported by standard meteorological measurements and visual estimates of sky conditi…
▽ More
Nighttime visible-light sky brightness and transparency are reported for the Polar Environment Research Laboratory (PEARL), located on a 610-m high ridge near the Eureka research station, on Ellesmere Island, Canada. Photometry of Polaris obtained in V band with the PEARL All Sky Imager (PASI) over two winters is supported by standard meteorological measurements and visual estimates of sky conditions from sea level. These data show that during the period of the study, October through March of 2008/09 and 2009/10, the sky near zenith had a mean surface brightness of 19.7 mag/square-arcsec when the sun was more than 12 deg below the horizon, reaching 20.7 mag/square-arcsec during astronomical darkness with no moon. Skies were without thick cloud and potentially usable for astronomy 86% of the time (extinction <2 mag). Up to 68% of the time was spectroscopic (<0.5 mag), attenuated by ice crystals, or clear with stable atmospheric transparency. Those conditions can persist for over 100 hours at a time. Further analysis suggests the sky was entirely free of ice crystals (truly photometric) 48+/-3% of the time at PEARL in winter, and that a higher elevation location nearby may be better.
△ Less
Submitted 13 December, 2011;
originally announced December 2011.
-
A Minimized Mutual Information retrieval for simultaneous atmospheric pressure and temperature
Authors:
Prabhat K. Koner,
James R. Drummond
Abstract:
The primary focus of the Mars Trace Gas Orbiter (TGO) collaboration between NASA and ESA is the detection of the temporal and spatial variation of the atmospheric trace gases using a solar occultation Fourier transform spectrometer. To retrieve any trace gas mixing ratios from these measurements, the atmospheric pressure and temperature have to be known accurately. Thus, a prototype retrieval mode…
▽ More
The primary focus of the Mars Trace Gas Orbiter (TGO) collaboration between NASA and ESA is the detection of the temporal and spatial variation of the atmospheric trace gases using a solar occultation Fourier transform spectrometer. To retrieve any trace gas mixing ratios from these measurements, the atmospheric pressure and temperature have to be known accurately. Thus, a prototype retrieval model for the determination of pressure and temperature from a broadband high resolution infrared Fourier Transform spectrometer experiment with the Sun as a source on board a spacecraft orbiting the planet Mars is presented. It is found that the pressure and temperature can be uniquely solved from remote sensing spectroscopic measurements using a Regularized Total Least Squares method and selected pairs of micro-windows without any a-priori information of the state space parameters and other constraints.
The selection of the pairs of suitable micro-windows is based on the information content analysis. A comparative information content calculation using Bayes theory and a hyperspace formulation are presented to understand the information available in measurement. A method of minimization of mutual information is used to search the suitable micro-windows for a simultaneous pressure and temperature retrieval.
△ Less
Submitted 21 December, 2010;
originally announced December 2010.
-
Towards ensemble asteroseismology of the young open clusters Chi Persei and NGC 6910
Authors:
S. Saesen,
A. Pigulski,
F. Carrier,
G. Michalska,
C. Aerts,
J. De Ridder,
M. Briquet,
G. Handler,
Z. Kolaczkowski,
B. Acke,
E. Bauwens,
P. G. Beck,
Y. Blom,
J. Blommaert,
E. Broeders,
M. Cherix,
G. Davignon,
J. Debosscher,
P. Degroote,
L. Decin,
S. Dehaes,
W. De Meester,
P. Deroo,
M. Desmet,
R. Drummond
, et al. (49 additional authors not shown)
Abstract:
As a result of the variability survey in Chi Persei and NGC6910, the number of Beta Cep stars that are members of these two open clusters is increased to twenty stars, nine in NGC6910 and eleven in Chi Persei. We compare pulsational properties, in particular the frequency spectra, of Beta Cep stars in both clusters and explain the differences in terms of the global parameters of the clusters. We…
▽ More
As a result of the variability survey in Chi Persei and NGC6910, the number of Beta Cep stars that are members of these two open clusters is increased to twenty stars, nine in NGC6910 and eleven in Chi Persei. We compare pulsational properties, in particular the frequency spectra, of Beta Cep stars in both clusters and explain the differences in terms of the global parameters of the clusters. We also indicate that the more complicated pattern of the variability among B type stars in Chi Persei is very likely caused by higher rotational velocities of stars in this cluster. We conclude that the sample of pulsating stars in the two open clusters constitutes a very good starting point for the ensemble asteroseismology of Beta Cep-type stars and maybe also for other B-type pulsators.
△ Less
Submitted 16 April, 2010;
originally announced April 2010.
-
Photometric multi-site campaign on the open cluster NGC 884 I. Detection of the variable stars
Authors:
S. Saesen,
F. Carrier,
A. Pigulski,
C. Aerts,
G. Handler,
A. Narwid,
J. N. Fu,
C. Zhang,
X. J. Jiang,
J. Vanautgaerden,
G. Kopacki,
M. Stȩślicki,
B. Acke,
E. Poretti,
K. Uytterhoeven,
C. Gielen,
R. Østensen,
W. De Meester,
M. D. Reed,
Z. Kołaczkowski,
G. Michalska,
E. Schmidt,
K. Yakut,
A. Leitner,
B. Kalomeni
, et al. (40 additional authors not shown)
Abstract:
CONTEXT: Recent progress in the seismic interpretation of field beta Cep stars has resulted in improvements of the physics in the stellar structure and evolution models of massive stars. Further asteroseismic constraints can be obtained from studying ensembles of stars in a young open cluster, which all have similar age, distance and chemical composition.
AIMS: To improve our comprehension of…
▽ More
CONTEXT: Recent progress in the seismic interpretation of field beta Cep stars has resulted in improvements of the physics in the stellar structure and evolution models of massive stars. Further asteroseismic constraints can be obtained from studying ensembles of stars in a young open cluster, which all have similar age, distance and chemical composition.
AIMS: To improve our comprehension of the beta Cep stars, we studied the young open cluster NGC 884 to discover new B-type pulsators, besides the two known beta Cep stars, and other variable stars.
METHODS: An extensive multi-site campaign was set up to gather accurate CCD photometry time series in four filters (U, B, V, I) of a field of NGC884. Fifteen different instruments collected almost 77500 CCD images in 1286 hours. The images were calibrated and reduced to transform the CCD frames into interpretable differential light curves. Various variability indicators and frequency analyses were applied to detect variable stars in the field. Absolute photometry was taken to deduce some general cluster and stellar properties.
RESULTS: We achieved an accuracy for the brightest stars of 5.7 mmag in V, 6.9 mmag in B, 5.0 mmag in I and 5.3 mmag in U. The noise level in the amplitude spectra is 50 micromag in the V band. Our campaign confirms the previously known pulsators, and we report more than one hundred new multi- and mono-periodic B-, A- and F-type stars. Their interpretation in terms of classical instability domains is not straightforward, pointing to imperfections in theoretical instability computations. In addition, we have discovered six new eclipsing binaries and four candidates as well as other irregular variable stars in the observed field.
△ Less
Submitted 7 January, 2010;
originally announced January 2010.
-
Post-AGB stars with hot circumstellar dust: binarity of the low-amplitude pulsators
Authors:
Hans Van Winckel,
Tom Lloyd Evans,
Maryline Briquet,
Peter De Cat,
Pieter Degroote,
Wim De Meester,
Joris De Ridder,
Pieter Deroo,
Maarten Desmet,
Rachel Drummond,
Laurent Eyer,
Martin A. T. Groenewegen,
Katrien Kolenberg,
David Kilkenny,
Djazia Ladjal,
Karolien Lefever,
Thomas Maas,
Fred Marang,
Peter Martinez,
Roy H. Østensen,
Gert Raskin,
Maarten Reyniers,
Pierre Royer,
Sophie Saesen,
Katrien Uytterhoeven
, et al. (6 additional authors not shown)
Abstract:
While the first binary post-AGB stars were serendipitously discovered, the distinct characteristics of their Spectral Energy Distribution (SED) allowed us to launch a more systematic search for binaries. We selected post-AGB objects which show a broad dust excess often starting already at H or K, pointing to the presence of a gravitationally bound dusty disc in the system. We started a very exte…
▽ More
While the first binary post-AGB stars were serendipitously discovered, the distinct characteristics of their Spectral Energy Distribution (SED) allowed us to launch a more systematic search for binaries. We selected post-AGB objects which show a broad dust excess often starting already at H or K, pointing to the presence of a gravitationally bound dusty disc in the system. We started a very extensive multi-wavelength study of those systems and here we report on our radial velocity and photometric monitoring results for six stars of early F type, which are pulsators of small amplitude. To determine the radial velocity of low signal-to-noise time-series, we constructed dedicated auto-correlation masks. The radial velocity variations were subjected to detailed analysis to differentiate between pulsational variability and variability due to orbital motion. Finally orbital minimalisation was performed to constrain the orbital elements. All of the six objects are binaries, with orbital periods ranging from 120 to 1800 days. Five systems have non-circular orbits. The mass functions range from 0.004 to 0.57 solar mass and the companions are likely unevolved objects of (very) low initial mass. We argue that these binaries must have been subject to severe binary interaction when the primary was a cool supergiant. Although the origin of the circumstellar disc is not well understood, the disc is generally believed to be formed during this strong interaction phase. The eccentric orbits of these highly evolved objects remain poorly understood. With the measured orbits and mass functions we conclude that the circumbinary discs seem to have a major impact on the evolution of a significant fraction of binary systems.
△ Less
Submitted 24 June, 2009;
originally announced June 2009.
-
Long-term photometric monitoring with the Mercator telescope. Frequencies and multicolour amplitudes of $γ$ Doradus stars
Authors:
J. Cuypers,
C. Aerts,
P. De Cat,
J. De Ridder,
K. Goossens,
C. Schoenaers,
K. Uytterhoeven,
B. Acke,
G. Davignon,
J. Debosscher,
L. Decin,
W. De Meester,
P. Deroo,
R. Drummond,
K. Kolenberg,
K. Lefever,
G. Raskin,
M. Reyniers,
S. Saesen,
B. Vandenbussche,
R. Van Malderen,
T. Verhoelst,
H. Van Winckel,
C. Waelkens
Abstract:
Gamma Doradus stars are excellent targets for asteroseismology since the gravity modes present in these stars probe the deep stellar interiors. Mode identification will improve the knowledge of these stars considerably. A selected group of Gamma Doradus stars and some candidates were observed with the Mercator telescope to find and/or confirm the periodicities in the light variations and to deri…
▽ More
Gamma Doradus stars are excellent targets for asteroseismology since the gravity modes present in these stars probe the deep stellar interiors. Mode identification will improve the knowledge of these stars considerably. A selected group of Gamma Doradus stars and some candidates were observed with the Mercator telescope to find and/or confirm the periodicities in the light variations and to derive reliable amplitude ratios in different pass bands. A frequency analysis was performed on all new data obtained in the Geneva photometric system. In order to have more reliable and accurate frequencies, the new data were combined with similar data from the literature and with Hipparcos observations. A set of frequencies that minimized the the residuals in a harmonic fit was searched for while allowing means and amplitudes to vary from one observation set to another. Frequencies and amplitudes in the photometric passbands of the Geneva system are given for 21 Gamma Doradus stars. We report the discovery of HD 74504 as a newly found Gamma Doradus star. Our study provides the first extensive multicolour database for the understanding of gravity modes in F-type stars.
△ Less
Submitted 30 March, 2009;
originally announced March 2009.
-
The Corot Book: Chap. V.5/ Extraction of the photometric information : corrections
Authors:
R. Samadi,
F. Fialho,
J. E. S. Costa,
R. Drummond,
L. Pinheiro Da Silva,
F. Baudin,
P. Boumier,
L. Jorda
Abstract:
We present here the set of corrections that will be applied to the raw data of the CoRoT mission. The aim is to correct the data for instrumental and environmental perturbations as well as to optimize the duty-cycle, in order to reach the expected performances of the mission. The main corrections are : the correction of the electromagnetic interferences, suppression of outliers, the background c…
▽ More
We present here the set of corrections that will be applied to the raw data of the CoRoT mission. The aim is to correct the data for instrumental and environmental perturbations as well as to optimize the duty-cycle, in order to reach the expected performances of the mission. The main corrections are : the correction of the electromagnetic interferences, suppression of outliers, the background correction, the jitter correction and the correction of the integration time variations. We focus here on these corrections and emphasize their efficiency.
△ Less
Submitted 14 March, 2007;
originally announced March 2007.
-
Asteroseismology of the Beta Cephei star 12 (DD) Lacertae: photometric observations, pulsational frequency analysis and mode identification
Authors:
G. Handler,
M. Jerzykiewicz,
E. Rodriguez,
K. Uytterhoeven,
P. J. Amado,
T. N. Dorokhova,
N. I. Dorokhov,
E. Poretti,
J. -P. Sareyan,
L. Parrao,
D. Lorenz,
D. Zsuffa,
R. Drummond,
J. Daszynska-Daszkiewicz,
T. Verhoelst,
J. De Ridder,
B. Acke,
P. -O. Bourge,
A. I. Movchan,
R. Garrido,
M. Paparo,
T. Sahin,
V. Antoci,
S. N. Udovichenko,
K. Csorba
, et al. (6 additional authors not shown)
Abstract:
We report a multisite photometric campaign for the Beta Cephei star 12 Lacertae. 750 hours of high-quality differential photoelectric Stromgren, Johnson and Geneva time-series photometry were obtained with 9 telescopes during 190 nights. Our frequency analysis results in the detection of 23 sinusoidal signals in the light curves. Eleven of those correspond to independent pulsation modes, and the…
▽ More
We report a multisite photometric campaign for the Beta Cephei star 12 Lacertae. 750 hours of high-quality differential photoelectric Stromgren, Johnson and Geneva time-series photometry were obtained with 9 telescopes during 190 nights. Our frequency analysis results in the detection of 23 sinusoidal signals in the light curves. Eleven of those correspond to independent pulsation modes, and the remainder are combination frequencies. We find some slow aperiodic variability such as that seemingly present in several Beta Cephei stars. We perform mode identification from our colour photometry, derive the spherical degree l for the five strongest modes unambiguously and provide constraints on l for the weaker modes. We find a mixture of modes of 0 <= l <= 4. In particular, we prove that the previously suspected rotationally split triplet within the modes of 12 Lac consists of modes of different l; their equal frequency splitting must thus be accidental.
One of the periodic signals we detected in the light curves is argued to be a linearly stable mode excited to visible amplitude by nonlinear mode coupling via a 2:1 resonance. We also find a low-frequency signal in the light variations whose physical nature is unclear; it could be a parent or daughter mode resonantly coupled. The remaining combination frequencies are consistent with simple light-curve distortions.
The range of excited pulsation frequencies of 12 Lac may be sufficiently large that it cannot be reproduced by standard models. We suspect that the star has a larger metal abundance in the pulsational driving zone, a hypothesis also capable of explaining the presence of Beta Cephei stars in the LMC.
△ Less
Submitted 11 October, 2005;
originally announced October 2005.