Geometric and energy-aware decomposition of the Navier-Stokes equations: A port-Hamiltonian approach

Federico Califano,Ramy Rashad,Frederic P. Schuller,Stefano Stramigioli

DOI: https://doi.org/10.1063/5.0048359

2021-03-03

Abstract:A port-Hamiltonian model for compressible Newtonian fluid dynamics is presented in entirely coordinate-independent geometric fashion. This is achieved by use of tensor-valued differential forms that allow to describe describe the interconnection of the power preserving structure which underlies the motion of perfect fluids to a dissipative port which encodes Newtonian constitutive relations of shear and bulk stresses. The relevant diffusion and the boundary terms characterizing the Navier-Stokes equations on a general Riemannian manifold arise naturally from the proposed construction.

Fluid Dynamics,Mathematical Physics

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to perform geometric and energy - aware decomposition of the Navier - Stokes equations by introducing the port - Hamiltonian (pH) framework. Specifically, the goals of the paper are:

1. **Extension of the port - Hamiltonian model for ideal fluids**: The paper first reviews the port - Hamiltonian model for ideal fluids (i.e., the incompressible or inviscid Euler equations) and extends it on this basis to Newtonian fluids (Navier - Stokes equations) that include viscous effects. The key to this extension lies in introducing a distributed port to connect the ideal fluid - dynamic system and the purely dissipative system that describes the shear and volumetric stresses of Newtonian fluids.

2. **Geometric and energy - aware decomposition**: The paper proposes a completely coordinate - independent geometric method, using tensor - valued differential forms to describe the energy - exchange subsystems in fluid dynamics. This method is applicable not only to flat spaces but also to Riemannian manifolds of arbitrary dimensions. In this way, the paper can naturally derive the diffusion and boundary terms of the Navier - Stokes equations on general Riemannian manifolds.

3. **Analysis of energy conservation and dissipation**: The paper analyzes in detail the energy balance in the fluid system, especially the energy dissipation due to viscous effects. By introducing new pairs of dual variables (effort and flow), the paper defines the "off - shell" and "on - shell" expressions of power balance and reveals the energy - exchange mechanism of the system through these expressions.

4. **Modular multi - physics system network framework**: The method proposed in the paper is applicable not only to a single - fluid system but can also be used to construct a modular multi - physics system network framework. This framework can maintain energy continuity between different physical domains, thus providing a powerful tool for the simulation and control of multi - physics systems.

### Main contributions

- **Geometric method**: The paper uses tensor - valued differential forms and geometric tools on Riemannian manifolds to provide a completely coordinate - independent representation method of the Navier - Stokes equations.
- **Energy - awareness**: By introducing the port - Hamiltonian framework, the paper can explicitly describe the energy - exchange and dissipation processes in the fluid system.
- **Modularity and extensibility**: The proposed framework is highly modular and extensible and can be easily applied to different types of fluids and other physical systems.

### Formula examples

- **Momentum conservation equation**:
\[
\dot{\nu}+d\left(\frac{1}{2} \star(\nu \wedge \star \nu)\right)+\iota_{v} d \nu =-\frac{d p}{\star \mu}+\frac{\star_{2} d^{\nabla} T}{\star \mu}
\]
where $\nu$ is the velocity 1 - form, $\mu$ is the mass - density n - form, $p$ is the pressure scalar field, and $T$ is the stress - tensor - valued (n - 1)-form.

- **Energy functional**:
\[
H(\nu, \mu)=\int_{M}\left(\frac{1}{2} \star \mu(\nu \wedge \star \nu)+U(\star \mu) \mu\right)
\]
where $U(\star \mu)$ is the internal - energy - density function.

- **Power balance**:
\[
\dot{H}=\int_{M} \delta_{\nu} H \wedge \dot{\nu}+\delta_{\mu} H \wedge \dot{\mu}
\]
where $\delta_{\nu} H$ and $\delta_{\mu} H$ are the variational derivatives of the energy functional $H$ with respect to the state variables $\nu$ and $\mu$, respectively.

Through these formulas, the paper shows how to handle complex fluid - dynamics problems in a geometric framework, especially in non - flat spaces.