Authors

Andreas Bartel, Malak Diab, Andreas Frommer, Michael Günther, Nicole Marheineke

Abstract

In the simulation of differential-algebraic equations (DAEs), it is essential to employ numerical schemes that take into account the inherent structure and maintain explicit or hidden algebraic constraints. This paper focuses on operator splitting techniques for coupled systems and aims at preserving the structure in the port-Hamiltonian framework. The study explores two decomposition strategies: one considering the underlying coupled subsystem structure and the other addressing energy-associated properties such as conservation and dissipation. We show that for coupled index-1 DAEs with and without private index-2 variables, the splitting schemes on top of a dimension-reducing decomposition achieve the same convergence rate as in the case of ordinary differential equations. Additionally, we discuss an energy-associated decomposition for linear time-invariant port-Hamiltonian index-1 DAEs and introduce generalized Cayley transforms to uphold energy conservation. The effectiveness of both strategies is evaluated using port-Hamiltonian benchmark examples from electric circuits.

Keywords

DAEs; Port-Hamiltonian systems; Cayley transform; Strang splitting

Citation

  • Journal: Applied Numerical Mathematics
  • Year: 2025
  • Volume: 214
  • Issue:
  • Pages: 28–53
  • Publisher: Elsevier BV
  • DOI: 10.1016/j.apnum.2025.03.004

BibTeX

@article{Bartel_2025,
  title={{Splitting techniques for DAEs with port-Hamiltonian applications}},
  volume={214},
  ISSN={0168-9274},
  DOI={10.1016/j.apnum.2025.03.004},
  journal={Applied Numerical Mathematics},
  publisher={Elsevier BV},
  author={Bartel, Andreas and Diab, Malak and Frommer, Andreas and Günther, Michael and Marheineke, Nicole},
  year={2025},
  pages={28--53}
}

Download the bib file

References