Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPerformance Assessment of Energy-preserving, Adaptive Time-Step Variational Integrators
71
0
0.0
(
0
)
Added by
Harsh Sharma
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A fixed time-step variational integrator cannot preserve momentum, energy, and symplectic form simultaneously for nonintegrable systems. This barrier can be overcome by treating time as a discrete dynamic variable and deriving adaptive time-step variational integrators that conserve the energy in addition to being symplectic and momentum-preserving. Their utility, however, is still an open question due to the numerical difficulties associated with solving the governing discrete equations. In this work, we investigate the numerical performance of energy-preserving, adaptive time-step variational integrators. First, we compare the time adaptation and energy performance of the energy-preserving adaptive algorithm with the adaptive variational integrator for Keplers problem. We also study the effect of variable precision arithmetic on the energy conservation properties. Second, we apply tools from Lagrangian backward error analysis to investigate numerical stability of the energy-preserving adaptive algorithm. Finally, we consider a simple mechanical system example to illustrate our backward stability approach by constructing a modified Lagrangian for the modified equation of an energy-preserving, adaptive time-step variational integrator.
rate research
Read More
The goal of this paper is to develop energy-preserving variational integrators for time-dependent mechanical systems with forcing. We first present the Lagrange-dAlembert principle in the extended Lagrangian mechanics framework and derive the extended forced Euler-Lagrange equations in continuous-time. We then obtain the extended forced discrete Euler-Lagrange equations using the extended discrete mechanics framework and derive adaptive time step variational integrators for time-dependent Lagrangian systems with forcing. We consider two numerical examples to study the numerical performance of energy-preserving variational integrators. First, we consider the example of a nonlinear conservative system to illustrate the advantages of using adaptive time-stepping in variational integrators. We show a trade-off between energy-preserving performance and accurate discrete trajectories while choosing an initial time step. In addition, we demonstrate how the implicit equations become more ill-conditioned as the adaptive time step increases through a condition number analysis. As a second example, we numerically simulate a damped harmonic oscillator using the adaptive time step variational integrator framework. The adaptive time step increases monotonically for the dissipative system leading to unexpected energy behavior.
The paper deals with numerical discretizations of separable nonlinear Hamiltonian systems with additive noise. For such problems, the expected value of the total energy, along the exact solution, drifts linearly with time. We present and analyze a time integrator having the same property for all times. Furthermore, strong and weak convergence of the numerical scheme along with efficient multilevel Monte Carlo estimators are studied. Finally, extensive numerical experiments illustrate the performance of the proposed numerical scheme.
We perform a numerical analysis of a class of randomly perturbed {H}amiltonian systems and {P}oisson systems. For the considered additive noise perturbation of such systems, we show the long time behavior of the energy and quadratic Casimirs for the exact solution. We then propose and analyze a drift-preserving splitting scheme for such problems with the following properties: exact drift preservation of energy and quadratic Casimirs, mean-square order of convergence one, weak order of convergence two. These properties are illustrated with numerical experiments.
A novel class of high-order linearly implicit energy-preserving exponential integrators are proposed for the nonlinear Schrodinger equation. We firstly done that the original equation is reformulated into a new form with a modified quadratic energy by the scalar auxiliary variable approach. The spatial derivatives of the system are then approximated with the standard Fourier pseudo-spectral method. Subsequently, we apply the extrapolation technique to the nonlinear term of the semi-discretized system and a linearized system is obtained. Based on the Lawson transformation, the linearized system is rewritten as an equivalent one and we further apply the symplectic Runge-Kutta method to the resulting system to gain a fully discrete scheme. We show that the proposed scheme can produce numerical solutions along which the modified energy is precisely conserved, as is the case with the analytical solution and is extremely efficient in the sense that only linear equations with constant coefficients need to be solved at every time step. Numerical results are addressed to demonstrate the remarkable superiority of the proposed schemes in comparison with other high-order structure-preserving method.
A time-fractional Allen-Cahn equation with volume constraint is first proposed by introducing a nonlocal time-dependent Lagrange multiplier. Adaptive linear second-order energy stable schemes are developed for the proposed model by combining invariant energy quadratization and scalar auxiliary variable approaches with the recent L1$^{+}$ formula. The new developed methods are proved to be volume-preserving and unconditionally energy stable on arbitrary nonuniform time meshes. The accelerated algorithm and adaptive time strategy are employed in numerical implement. Numerical results show that the proposed algorithms are computationally efficient in multi-scale simulations, and appropriate for accurately resolving the intrinsically initial singularity of solution and for efficiently capturing the fast dynamics away initial time.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
