A new variational discretization technique for initial value problems bypassing governing equations

Abstract

Motivated by the fact that both the classical and quantum description of nature rest on causality and a variational principle, we develop a novel and highly versatile discretization prescription for classical initial value problems (IVPs). It is based on an optimization (action) functional with doubled degrees of freedom, which is discretized using a single regularized summation-by-parts (SBP) operator. Formulated as optimization task it allows us to obtain classical trajectories without the need to derive an equation of motion. The novel regularization we develop in this context is inspired by the weak imposition of initial data, often deployed in the modern treatment of IVPs and is implemented using affine coordinates. We demonstrate numerically the stability, accuracy and convergence properties of our approach in systems with classical equations of motion featuring both first and second order derivatives in time. onvergence properties of our approach in systems with classical equations of motion featuring both first and second order derivatives in time.

arXiv comment28 pages, 15 figures, accepted in Journal of Computational Physics

References in corpus (10)

• Keldysh Field Theory for Driven Open Quantum Systems
• Review of Summation-by-parts schemes for initial-boundary-value problems
• The classical mechanics of non-conservative systems
• Summation-by-parts operators for correction procedure via reconstruction
• Quantum versus classical statistical dynamics of an ultracold Bose gas
• Geometric Computational Electrodynamics with Variational Integrators and Discrete Differential Forms
• Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems
• Slimplectic Integrators: Variational Integrators for General Nonconservative Systems
• Implicit schemes for real-time lattice gauge theory
• Discrete Vector Calculus and Helmholtz Hodge Decomposition for Classical Finite Difference Summation by Parts Operators

Cited by in corpus (2)

• Variational approach to nonholonomic and inequality-constrained mechanics
• Exact symmetry conservation and automatic mesh refinement in discrete initial boundary value problems