Energy-Based Modeling and Structure-Preserving Discretization of Physical Systems
M. H. M Rashid
TL;DR
The paper develops a unified energy-based framework that generalizes port-Hamiltonian systems to handle constraints and high-index differential-algebraic equations, with a focus on preserving energy balance and dissipation through both modeling and discretization. It introduces structure-preserving regularization to reduce index, proves exponential stability under coercivity, and provides nonlinear structure-preserving discretizations (midpoint and discrete gradient) that maintain the energy-dissipation structure at the discrete level. The theory is demonstrated across diverse applications, including poroelasticity, nonlinear circuits, constrained mechanics, phase-field models, and quantum-thermodynamic memory systems, illustrating robust numerical behavior and physically meaningful long-time dynamics. The work lays a solid foundation for computational physics and engineering, with future directions toward adaptive time stepping, stochastic port-Hamiltonian systems, and model reduction to tackle large-scale multiphysics problems.
Abstract
This paper develops a comprehensive mathematical framework for energy-based modeling of physical systems, with particular emphasis on preserving fundamental structural properties throughout the modeling and discretization process. The approach provides systematic methods for handling challenging system classes including high-index differential-algebraic equations and nonlinear multiphysics problems. Theoretical foundations are established for regularizing constrained systems while maintaining physical consistency, analyzing stability properties, and constructing numerical discretizations that inherit the energy dissipation structure of the continuous models. The versatility and practical utility of the framework are demonstrated through applications across multiple domains including poroelastic media, nonlinear circuits, constrained mechanics, and phase-field models. The results ensure that essential physical properties such as energy balance and dissipation are maintained from the continuous formulation through to numerical implementation, providing robust foundations for computational physics and engineering applications.
