scholarly journals Accelerated Battery Lifetime Simulations Using Adaptive Inter-Cycle Extrapolation Algorithm

Author(s):  
Valentin Sulzer ◽  
Peyman Mohtat ◽  
Sravan Pannala ◽  
Jason Siegel ◽  
Anna Stefanopoulou

Abstract We propose algorithms to speed up physics-based battery lifetime simulations by one to two orders of magnitude compared to the state-of-the-art. First, we propose a reformulation of the Single Particle Model with side reactions to remove algebraic equations and hence reduce stiffness, with 3x speed-up in simulation time (intra-cycle reformulation). Second, we introduce an algorithm that makes use of the difference between the `fast' timescale of battery cycling and the `slow' timescale of battery degradation by adaptively selecting and simulating representative cycles, skipping other cycles, and hence requires fewer cycle simulations to simulate the entire lifetime (adaptive inter-cycle extrapolation). This algorithm is demonstrated with a specific degradation mechanism but can be applied to various models of aging phenomena. In the particular case study considered, simulations of the entire lifetime are performed in under 5 seconds. This opens the possibility for much faster and more accurate model development, testing, and comparison with experimental data.

2021 ◽  
Author(s):  
Valentin Sulzer ◽  
Peyman Mohtat ◽  
Sravan Pannala ◽  
Jason B. Siegel ◽  
Anna G. Stefanopoulou

We propose algorithms to speed up physics-based battery lifetime simulations by one to two orders of magnitude compared to the state-of-the-art. First, we propose a reformulation of the Single Particle Model with side reactions to remove algebraic equations and hence reduce stiffness, with 3x speed-up in simulation time (intra-cycle reformulation). Second, we introduce an algorithm that makes use of the difference between the `fast' timescale of battery cycling and the `slow' timescale of battery degradation by adaptively selecting and simulating representative cycles, skipping other cycles, and hence requires fewer cycle simulations to simulate the entire lifetime (adaptive inter-cycle extrapolation). This algorithm is demonstrated with a specific degradation mechanism but can be applied to various models of aging phenomena. In the particular case study considered, simulations of the entire lifetime are performed in under 5 seconds. This opens the possibility for much faster and more accurate model development, testing, and comparison with experimental data.


2019 ◽  
Vol 26 (3) ◽  
pp. 341-349 ◽  
Author(s):  
Givi Berikelashvili ◽  
Manana Mirianashvili

Abstract A three-level finite difference scheme is studied for the initial-boundary value problem of the generalized Benjamin–Bona–Mahony–Burgers equation. The obtained algebraic equations are linear with respect to the values of the desired function for each new level. The unique solvability and absolute stability of the difference scheme are shown. It is proved that the scheme is convergent with the rate of order {k-1} when the exact solution belongs to the Sobolev space {W_{2}^{k}(Q)} , {1<k\leq 3} .


2016 ◽  
Vol 849 ◽  
pp. 76-83
Author(s):  
Jiří Náprstek ◽  
Cyril Fischer

The exact coincidence of external excitation and basic eigen-frequency of a single degree of freedom (SDOF) nonlinear system produces stationary response with constant amplitude and phase shift. When the excitation frequency differs from the system eigen-frequency, various types of quasi-periodic response occur having a character of a beating process. The period of beating changes from infinity in the resonance point until a couple of excitation periods outside the resonance area. Theabove mentioned phenomena have been identified in many papers including authors’ contributions. Nevertheless, investigation of internal structure of a quasi-period and its dependence on the difference of excitation and eigen-frequency is still missing. Combinations of harmonic balance and small parameter methods are used for qualitative analysis of the system in mono- and multi-harmonic versions. They lead to nonlinear differential and algebraic equations serving as a basis for qualitativeanalytic estimation or numerical description of characteristics of the quasi-periodic system response. Zero, first and second level perturbation techniques are used. Appearance, stability and neighborhood of limit cycles is evaluated. Numerical phases are based on simulation processes and numerical continuation tools. Parametric evaluation and illustrating examples are presented.


2021 ◽  
Vol 28 (2) ◽  
pp. 163-182
Author(s):  
José L. Simancas-García ◽  
Kemel George-González

Shannon’s sampling theorem is one of the most important results of modern signal theory. It describes the reconstruction of any band-limited signal from a finite number of its samples. On the other hand, although less well known, there is the discrete sampling theorem, proved by Cooley while he was working on the development of an algorithm to speed up the calculations of the discrete Fourier transform. Cooley showed that a sampled signal can be resampled by selecting a smaller number of samples, which reduces computational cost. Then it is possible to reconstruct the original sampled signal using a reverse process. In principle, the two theorems are not related. However, in this paper we will show that in the context of Non Standard Mathematical Analysis (NSA) and Hyperreal Numerical System R, the two theorems are equivalent. The difference between them becomes a matter of scale. With the scale changes that the hyperreal number system allows, the discrete variables and functions become continuous, and Shannon’s sampling theorem emerges from the discrete sampling theorem.


Author(s):  
Vladimir P. Gerdt ◽  
Mikhail D. Malykh ◽  
Leonid A. Sevastianov ◽  
Yu Ying

The article considers the midpoint scheme as a finite-difference scheme for a dynamical system of the form ̇ = (). This scheme is remarkable because according to Cooper’s theorem, it preserves all quadratic integrals of motion, moreover, it is the simplest scheme among symplectic Runge-Kutta schemes possessing this property. The properties of approximate solutions were studied in the framework of numerical experiments with linear and nonlinear oscillators, as well as with a system of several coupled oscillators. It is shown that in addition to the conservation of all integrals of motion, approximate solutions inherit the periodicity of motion. At the same time, attention is paid to the discussion of introducing the concept of periodicity of an approximate solution found by the difference scheme. In the case of a nonlinear oscillator, each step requires solving a system of nonlinear algebraic equations. The issues of organizing computations using such schemes are discussed. Comparison with other schemes, including those symmetric with respect to permutation of and .̂


2021 ◽  
Vol 1016 ◽  
pp. 407-412
Author(s):  
Marcin Drajewicz ◽  
Tadeusz Kubaszek ◽  
Barbara Kościelniak ◽  
Marek Góral ◽  
Dawid Dziadosz

The thermal barrier coatings (TBC) are commonly used for protection of jet engine parts. In presented article the influence of chemical composition of TBCs bond coats on isothermal oxidation resistance was analysed. The bond coat was plasma sprayed (APS) using different MCrAlY alloys produced by Oerlikon-Metco: Metco 4451, AMDRY 997, AMDRY 962, AMDRY 365-1, AMDRY 995C. The conducted research showed big difference in coating thickness despite the same spraying parameters. The difference in porosity was not observed. The isothermal oxidation test was conducted during 500h at 1100°C in static laboratory air. The obtained results showed, that degradation mechanism of MCrAlY bond coats regardless of chemical composition is similar. The formation of scale contained aluminium and chromium oxides was observed. The internal oxidation effect was also observed in produced coatings.


2019 ◽  
Vol 97 (Supplement_3) ◽  
pp. 24-24
Author(s):  
Luis O Tedeschi

Abstract The establishment of credibility for a mathematical model’s (MM) predictive ability is an essential component for improving the MM because it stimulates the evolutionary thinking (i.e., the next generation of the model) of mental conceptualizations, assumptions, and boundaries of the MM. Its predictive adequacy is commonly assessed through its ability to precisely or accurately predict observed (real) values. The precision component measures how closely the model predicted values are of each other or whether a defined pattern of predictions exists. The accuracy component, on the other hand, measures how closely the average of the model predicted values are to the actual (true) average. Many statistics exist to determine precision and accuracy of MM such as mean bias, resistant coefficient of determination, coefficient of determination, modeling efficiency, concordance correlation coefficient (CCC), the mean square error of prediction, Kleijnen’s statistic (regression of the difference between predicted and observed on their sum), and Altman and Bland’s limits of agreement statistics among many more. However, for complex models that use skewed data or repeated data in which the data is not independent (e.g., multiple measurements on the same subject), simple statistics may not suffice. For instance, four methods to compute CCC exist (moment, variance components, U-statistics, and generalized estimating equations—GEE), but only the last two methods are resilient to lightly skewed data. Another type of complexity arises when meta-analytical approaches are used at the model development phase or the model evaluation phase. In general, meta-analytical approaches remove errors (i.e., variation) associated with random variables that are believed to be known. Under these circumstances, MM tends to overperform (i.e., they have greater predictive adequacy) and their future performance may be deceitful when trying to forecast at scenarios in which the random variable(s) is(are) indeterminable or unquantifiable.


1999 ◽  
Vol 23 (1B) ◽  
pp. 159-172 ◽  
Author(s):  
J. Phillips ◽  
D. Naylor ◽  
S.J. Harrison ◽  
P.H. Oosthuizen

The present numerical study examines the influence of horizontal louvres (Venetian blinds) on the convective heat transfer from a vertical isothermal surface. A steady, laminar, two dimensional, conjugate conduction/convection solution to this problem has been obtained using the finite element method. Detailed comparisons are made with temperature field and local Nusselt number data obtained using a Mach-Zehnder Interferometer. Also, numerically predicted stream function contours are compared to streamlines obtained from flow visualization experiments. All results are obtained for Pr = 0.7, as the intended application of this study is for air. Results show that as the blind tip-to-plate spacing decreases, the difference between numerical solutions and experimental results increase. This suggests that radiation heat transfer may be a significant factor at smaller blind spacings.


2018 ◽  
Vol 13 (5) ◽  
pp. 873-878
Author(s):  
Noriko N. Ishizaki ◽  
Koji Dairaku ◽  
Genta Ueno ◽  
◽  

A new method was proposed for the probabilistic projection of future climate that introduced quantile mapping to a regression method using a multi-model ensemble (QM_RMME). Results of this method were then compared with those of the traditional regression method (RMME). Six stations in Japan where 100 year observation records were available were used to evaluate the performance of the methods. An initial 50-year period (1901–1950) was used to develop the regression models and the final period (1951–2000) was used for evaluation. Results showed that the estimation errors at the 50th and 90th percentile were smaller for QM_RMME as compared to RMME at most sites. Conversely, when the model development and evaluation periods were limited to 20 years (1901–1920 and 1951–1970, respectively), the 90th percentile error was larger for QM_RMME. This was attributed to quantile mapping resulting in over-fitting of the data during the model development period. Furthermore, the QM_RMME error increased when the difference of observations between the model development and verification periods was large. Therefore, results indicated that the RMME method was more stable for relatively short data verification periods.


2015 ◽  
Vol 15 (21) ◽  
pp. 31247-31286
Author(s):  
L. K. Whalley ◽  
D. Stone ◽  
B. Bandy ◽  
R. Dunmore ◽  
J. F. Hamilton ◽  
...  

Abstract. Near-continuous measurements of OH reactivity in the urban background atmosphere of central London during the summer of 2012 are presented. OH reactivity behaviour is seen to be broadly dependent on airmass origin with the highest reactivity and the most pronounced diurnal profile observed when air had passed over central London to the East, prior to measurement. Averaged over the entire observation period of 26 days, OH reactivity peaked at ~ 27 s−1 in the morning with a minimum of ~ 15 s−1 during the afternoon. A maximum OH reactivity of 116 s−1 was recorded on one day during morning rush hour. A detailed box model using the Master Chemical Mechanism was used to calculate OH reactivity, and was constrained with an extended measurement dataset of volatile organic compounds (VOCs) derived from GC-FID and a two-dimensional GC instrument which included heavier molecular weight (up to C12) aliphatic VOCs, oxygenated VOCs and the biogenic VOCs of α pinene and limonene. Comparison was made between observed OH reactivity and modelled OH reactivity using (i) a standard suite of VOC measurements (C2-C8 hydrocarbons and a small selection of oxygenated VOCs) and (ii) a more comprehensive inventory including species up to C12. Modelled reactivities were lower than those measured (by 33 %) when only the reactivity of the standard VOC suite was considered. The difference between measured and modelled reactivity was improved, to within 15 %, if the reactivity of the higher VOCs (&amp;geq; C9) was also considered, with the reactivity of the biogenic compounds of α pinene and limonene and their oxidation products almost entirely responsible for this improvement. Further improvements in the model's ability to reproduce OH reactivity (to within 6 %) could be achieved if the reactivity and degradation mechanism of unassigned two-dimensional GC peaks were estimated. Neglecting the contribution of the higher VOCs (&amp;geq; C9) (particularly α pinene and limonene) and model-generated intermediates worsened the agreement between modelled and observed OH concentrations (by 41 %) and the magnitude of in situ ozone production calculated from the production of RO2 was significantly lower (60 %). This work highlights that any future ozone abatement strategies should consider the role that biogenic emissions play alongside anthropogenic emissions in influencing London's air quality.


Sign in / Sign up

Export Citation Format

Share Document