Irreversible Port-Hamiltonian Formulations for 1-Dimensional fluid systems

Imagine you are trying to write the ultimate instruction manual for how fluids (like water, air, or oil) move and change temperature. For a long time, scientists had two main ways to look at this:

The "Riverbank" View (Eulerian): You stand still on a bridge and watch the water rush past you. You measure the speed and temperature of the water at that specific spot as it flows by.
The "Surfer" View (Lagrangian): You hop on a surfboard and ride a specific drop of water. You follow that single drop as it moves, changes speed, and heats up.

The Problem:
Most modern engineering tools (called "Port-Hamiltonian Systems") are great at describing things that just sit there and diffuse heat (like a cooling cup of coffee). But they struggle when things are moving fast (convection) and losing energy to friction (viscosity) or heat. When you add movement and friction, the math gets messy, and the "laws of physics" (specifically the First and Second Laws of Thermodynamics) can get lost in the equations.

The Solution (The Paper's Big Idea):
The authors of this paper have built a new, super-flexible "instruction manual" called Irreversible Port-Hamiltonian Systems (IPHS). Think of this as a universal translator that can take the messy, moving, friction-filled reality of fluids and translate it into a clean, structured language that engineers can use to design better controls.

Here is the breakdown using simple analogies:

1. The "Energy Bank" vs. The "Entropy Tax"

In this new framework, the system is viewed through two lenses:

The Hamiltonian (Total Energy): Think of this as the Total Bank Account of the fluid. It includes the money from motion (kinetic energy) and the money from heat (internal energy). The First Law of Thermodynamics says this account can only change if you deposit or withdraw money from the outside (boundary). It never just disappears.
The Entropy (The Tax): Think of entropy as a Tax on Disorder. The Second Law of Thermodynamics says you can never get a tax refund; you can only pay more. Friction (viscosity) and heat flow are like "tax collectors" that constantly take energy from the "Bank Account" and turn it into "Entropy" (waste heat). You can't get it back.

2. The "Traffic Controller" (The Operators)

The paper's main trick is how it handles the "traffic" of the fluid.

In the old "Riverbank" view (Eulerian), the fluid is constantly swapping places. A drop of hot water moves into a cold spot, and a drop of cold water moves into a hot spot.
The authors realized that to make the math work, they had to tweak the "traffic lights" (the differential operators). They added special "convective terms" that act like a moving walkway at an airport. This walkway carries the energy and entropy along with the fluid particles, ensuring that the math correctly accounts for the fact that the fluid is moving, not just sitting there.

3. The "Boundary Doors" (Control)

Imagine the fluid is in a pipe with doors at the ends (the boundaries).

The Old Way: You could only push or pull the fluid.
The New Way (IPHS): The authors designed a special set of "Doors and Handles" (Boundary Port Variables). These handles allow you to control the fluid by managing:
- Pressure and Velocity: Pushing the fluid.
- Temperature and Heat: Heating or cooling the fluid.
- Entropy: Managing the "waste."

The magic of their new formula is that no matter how you turn these handles, the system guarantees that:

Energy is conserved (you don't create energy out of thin air).
Entropy always increases or stays the same (you can't un-mix the coffee and milk).

4. Why Does This Matter?

Why do we care about this math?

Better Design: Engineers use these models to design things like jet engines, chemical reactors, and nuclear power plants.
Stability: Because the model respects the laws of physics by design, the control systems built on top of it are naturally stable. They won't accidentally blow up the engine because the math "knows" that energy can't be created.
Versatility: This paper bridges the gap between the "Riverbank" view and the "Surfer" view, showing that you can use the powerful new control tools on both types of fluid descriptions.

The Takeaway

Think of this paper as upgrading the operating system for fluid dynamics. It takes the chaotic, real-world behavior of moving, heating, and friction-filled fluids and organizes it into a neat, rule-abiding structure. This allows engineers to write "apps" (control algorithms) that can steer these complex systems safely and efficiently, ensuring they obey the fundamental laws of the universe.

Here is a detailed technical summary of the paper "Irreversible Port-Hamiltonian Formulations for 1-Dimensional fluid systems."

1. Problem Statement

Classical Port-Hamiltonian Systems (PHS) are powerful frameworks for modeling and controlling physical systems, but they traditionally struggle with irreversible thermodynamic systems where dynamics depend strongly on temperature and entropy production (e.g., viscous fluids, heat conduction). While the Irreversible Port-Hamiltonian Systems (IPHS) framework was previously developed for diffusion-driven systems and extended to 1D spatial domains, it was not suitable for convective systems (fluid flows).

The specific challenge addressed in this paper is the modeling of non-isentropic compressible fluids with viscous dissipation within the Eulerian description (fixed spatial coordinates).

The Conflict: In the Eulerian frame, fluid particles move past a fixed point. Standard IPHS formulations for diffusion assume local equilibrium tracking material particles (Lagrangian), making the direct application of existing IPHS operators to convective terms inconsistent.
The Goal: To extend the IPHS framework to encompass convective transport in the Eulerian description while strictly adhering to the First and Second Laws of Thermodynamics.

2. Methodology

The authors employ a systematic derivation process involving three main stages:

A. Governing Equations Derivation

The paper first derives the governing equations for 1D non-isentropic compressible fluids in two reference frames:

Eulerian Coordinates: Using Reynolds transport theorem and divergence theorem, the authors derive continuity, momentum, and entropy equations. Key variables include density ( $\rho$ ), velocity ( $v$ ), pressure ( $p$ ), and specific entropy ( $s$ ). Viscous stress ( $\tau$ ) and heat flux ( $q$ ) are modeled via Newtonian and Fourier laws, respectively.
Lagrangian Coordinates: The authors transform the equations to a material coordinate system ( $\eta$ ) to track fluid particles, deriving equivalent equations for specific volume ( $\phi = 1/\rho$ ), velocity, and entropy.

B. IPHS Formulation

The authors reformulate these governing equations into the IPHS structure:
$\dot{x} = (J(x) - R(x)) \frac{\delta H}{\delta x}$
Where $H$ is the total energy (Hamiltonian), and the system matrix is split into a skew-symmetric interconnection part ( $J$ ) and a symmetric positive semi-definite dissipation part ( $R$ ).

State Variables: Defined as energy variables (e.g., $\rho, v, s$ in Eulerian; $\phi, v, s$ in Lagrangian).
Co-energy Variables: Derived as functional derivatives of the Hamiltonian (e.g., enthalpy, momentum, temperature).
Operator Modification: The core methodological innovation is the modification of the underlying differential operators. The authors introduce specific operators ( $A, B, C$ $A, B, C$ in Eulerian; $A_1, B_1, C_1$ $A_{1}, B_{1}, C_{1}$ in Lagrangian) that handle the convective terms ( $v \partial_\zeta$ $v \partial_{ζ}$ ) and dissipative terms (viscosity $\mu$ $μ$ , thermal conductivity $k$ $k$ ).
- They prove that the dissipative operators are formal adjoints of each other (e.g., $B = A^*$ ) and that the entropy production operator is skew-adjoint, ensuring the mathematical structure holds.

C. Boundary Control Parametrization

The paper extends the definition of Boundary Controlled IPHS (BC-IPHS).

They define a general structure for infinite-dimensional systems with convective flows and irreversible processes.
They introduce a specific parametrization of boundary port variables (inputs $u$ and outputs $y$ ) using modified effort variables.
This parametrization is designed to ensure that the power balance at the boundaries satisfies the First Law (energy conservation) and that the internal entropy production satisfies the Second Law (non-negative entropy generation).

3. Key Contributions

Extension of IPHS to Convective Flows: The paper successfully bridges the gap between IPHS theory and Eulerian fluid dynamics. It demonstrates that convective transport can be consistently modeled within the IPHS framework by modifying the differential operators, rather than switching entirely to a Lagrangian framework.
Dual Formulations (Eulerian & Lagrangian): The authors provide rigorous IPHS formulations for the same physical system in both coordinate systems. This validates the framework's robustness and highlights the structural differences in the operators required for each frame.
Thermodynamic Consistency: The proposed formulations are proven to strictly satisfy:
- First Law: The rate of change of total energy equals the power supplied through boundaries (when internal dissipation is accounted for).
- Second Law: The rate of change of total entropy is the sum of boundary entropy flux and non-negative internal entropy production (due to viscosity and heat conduction).
Generalized Boundary Port Definition: A new parametrization for boundary inputs/outputs is derived (Definition 4.2) that accommodates the coupling between material and thermal domains caused by convection.

4. Results

Eulerian Formulation: The system is represented by a matrix equation involving skew-symmetric operators for convection and symmetric operators for diffusion/viscosity. The boundary inputs/outputs are identified as combinations of enthalpy, viscous stress, and entropy flux.
Lagrangian Formulation: A parallel formulation is derived where the convective terms vanish (as the frame moves with the fluid), but the operators are adjusted for the material density $\rho_0$ .
Verification: Through Lemmas 4 and 5 and Propositions 4.3 and 4.4, the authors mathematically prove that:
- $\frac{dH}{dt} = y^T u$ (Energy Balance).
- $\frac{dS}{dt} = \int \sigma_s d\zeta + y_s^T u_s \geq 0$ (Entropy Balance), where $\sigma_s$ is the internal entropy production rate.
Operator Properties: The paper rigorously proves the adjoint relationships between the convective-dissipative operators (e.g., Lemma 1 and Lemma 3), which is essential for the stability and passivity of the system.

5. Significance

Control Design: By casting complex fluid-thermal systems into the IPHS framework, the paper opens the door for energy-shaping control and entropy assignment strategies. This is crucial for stabilizing systems like chemical reactors, nuclear processes, or active material flows where temperature and entropy are critical.
Unified Framework: It unifies the modeling of diffusion (previously handled by IPHS) and convection (previously difficult) under a single thermodynamic framework.
Physical Interpretation: The explicit separation of reversible (interconnection) and irreversible (dissipation) mechanisms in the model structure provides deep physical insight, making it easier to analyze energy flows and entropy generation in engineering applications.
Future Applications: The methodology serves as a foundation for modeling more complex 1D and potentially multi-dimensional fluid systems, including those with chemical reactions or phase changes, within a passivity-based control context.

In summary, this paper successfully generalizes the IPHS framework to handle convective, viscous, and thermal fluid dynamics in the Eulerian description, providing a mathematically rigorous and thermodynamically consistent tool for the analysis and control of irreversible distributed parameter systems.