Numerical approximation of port-Hamiltonian systems for hyperbolic or parabolic PDEs with boundary control

07/16/2020
by   Andrea Brugnoli, et al.
0

The present manuscript concerns the design of structure-preserving discretization methods for the solution of systems of boundary controlled Partial Differential Equations (PDEs) thanks to the port- Hamiltonian formalism. We first provide a general structure of infinite-dimensional port-Hamiltonian systems (pHs) for which the Partitioned Finite Element Method (PFEM) straightforwardly applies. The proposed strategy is particularised to abstract multidimensional linear hyperbolic and parabolic systems of PDEs. Then we show that instructional model problems based on the wave equation, Mindlin equation and heat equation fit within this unified framework. Secondly we introduce the ongoing project SCRIMP (Simulation and ContRol of Interactions in Multi-Physics) developed for the numerical simulation of infinite-dimensional pHs. SCRIMP notably relies on the FEniCS open-source computing platform for the finite element spatial discretization. Finally, we illustrate how to solve the considered model problems within this framework by carefully explaining the methodology. As additional support, companion interactive Jupyter notebooks are provided.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
06/14/2019

A Partitioned Finite Element Method for power-preserving discretization of open systems of conservation laws

This paper presents a structure-preserving spatial discretization method...
research
12/14/2022

Port-Hamiltonian Discontinuous Galerkin Finite Element Methods

A port-Hamiltonian (pH) system formulation is a geometrical notion used ...
research
03/18/2022

Convex Optimization-Based Structure-Preserving Filter For Multidimensional Finite Element Simulations

In simulation sciences, it is desirable to capture the real-world proble...
research
05/21/2018

CUDACLAW: A high-performance programmable GPU framework for the solution of hyperbolic PDEs

We present cudaclaw, a CUDA-based high performance data-parallel framewo...
research
07/23/2023

Proximal Galerkin: A structure-preserving finite element method for pointwise bound constraints

The proximal Galerkin finite element method is a high-order, nonlinear n...
research
07/01/2018

Bit Complexity of Computing Solutions for Symmetric Hyperbolic Systems of PDEs with Guaranteed Precision

We establish upper bounds of bit complexity of computing solution operat...

Please sign up or login with your details

Forgot password? Click here to reset