DAEs that should not be solved
For a Hamiltonian system with an autonomous Hamiltonian , the author is seeking for a numerical method which conserves the Hamiltonian or the total energy. There exist various schemes which stabilize the discret system with respect to the invariant. Here, the invariant is projected and the obtained pure index-2 DAE is discretized by a modified midpoint method. Such a projection midpoint scheme retains reversibility but not behave always well. For the harmonic oscillator equations in the highly oscillatory case, the projected DAEs enforcing preservation of the energy are linearized about a current iterate. Then, the underlying ODE will be unstable for some interval and the DAEs should not be solved. Three concret examples are discussed. For the entire collection see [Zbl 0930.00048].
- scientific article; zbMATH DE number 4108066
- Numerical integration of differential algebraic systems and invariant manifolds
- Comparing stability properties of three methods in DAEs or ODEs with invariants
- Stabilization of DAEs and invariant manifolds
- Stability of numerical methods for solving linear index-3 DAEs
This page was built for publication: DAEs that should not be solved
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4945806)