Structure-preserving space discretization of differential and nonlocal constitutive relations for port-Hamiltonian systems
For researchers in port-Hamiltonian systems and structure-preserving discretization, this work provides a novel method that preserves key physical invariants for a class of systems with differential and nonlocal constitutive relations.
This paper develops a structure-preserving space discretization method for port-Hamiltonian systems with differential and nonlocal constitutive relations, using a Stokes-Lagrange structure and finite element method. The method preserves enstrophy and kinetic energy evolution at both semi-discrete and fully-discrete levels, demonstrated on 1D nanorod, shear beam, and 2D Navier-Stokes equations.
We study the structure-preserving space discretization of port-Hamiltonian (pH) systems defined with differential constitutive relations. Using the concept of Stokes-Lagrange structure to describe these relations, these are reduced to a finite-dimensional Lagrange subspace of a pH system thanks to a structure-preserving Finite Element Method. To illustrate our results, the 1D nanorod case and the shear beam model are considered, which are given by differential and implicit constitutive relations for which a Stokes-Lagrange structure along with boundary energy ports naturally occur. Then, these results are extended to the nonlinear 2D incompressible Navier-Stokes equations written in a vorticity-stream function formulation. It is first recast as a pH system defined with a Stokes-Lagrange structure along with a modulated Stokes-Dirac structure. A careful structure-preserving space discretization is then performed, leading to a finite-dimensional pH system. Theoretical and numerical results show that both enstrophy and kinetic energy evolutions are preserved both at the semi-discrete and fully-discrete levels.