Abstract:In this paper we introduce discrete gradient methods to discretize irreversible port-Hamiltonian systems showing that the main qualitative properties of the continuous system are preserved using this kind discretizations methods.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to design a discrete - gradient method for Irreversible Port - Hamiltonian Systems (IPHS) to ensure that the main mass characteristics of the system are preserved during the numerical discretization process. Specifically, the paper focuses on how to discretize the irreversible port - Hamiltonian system through the discrete - gradient method, thereby ensuring that the properties of the energy and entropy balance equations are maintained.
### Core problems of the paper
1. **Maintaining energy conservation**: Ensure that the total energy of the system is conserved during the discretization process.
2. **Maintaining the principle of entropy increase**: Ensure that the discretized system can correctly reflect the entropy increase phenomenon in the irreversible process.
3. **Preserving the geometric structure**: Ensure that the discretization method can preserve the geometric structure of the original continuous system, such as the almost - Poisson bracket.
### Specific background
The irreversible port - Hamiltonian system is a special Hamiltonian system, which is suitable for describing open physical systems, especially those related to thermodynamic irreversible processes. This type of system uses two thermodynamic potential functions - total energy and entropy function, and is defined through the almost - Poisson bracket, rather than the traditional symplectic structure or metriplectic structure. Therefore, designing appropriate numerical integration methods is crucial for accurately simulating the dynamic behavior of these systems.
### Application of the discrete - gradient method
The discrete - gradient method is a method for numerically solving ordinary differential equations (ODEs). It constructs a numerical integrator by discretizing the gradient, thereby ensuring that certain key physical quantities (such as energy) remain unchanged during the discretization process. For the irreversible port - Hamiltonian system, the discrete - gradient method is used to ensure:
- **Energy conservation**: The discretized system still satisfies the law of energy conservation.
- **Entropy increase**: The discretized system can correctly reflect the entropy increase phenomenon in the irreversible process.
- **Preserving the geometric structure**: The discretization method can preserve the geometric structure of the original system, such as the almost - Poisson bracket.
### Conclusion
The paper proposes a discretization scheme suitable for the irreversible port - Hamiltonian system. This scheme can simultaneously verify the energy balance equation and the positivity of the entropy production term during the numerical discretization process. Future work will extend this scheme to be applied to a wider range of system types and apply it to more realistic physical systems.
### Summary of mathematical formulas
- **Energy and entropy balance equations**:
\[
\frac{dH}{dt}-y^{T}u = 0
\]
\[
\frac{dS}{dt}-y^{T}\tau=\sum_{i}\hat{\gamma}_{i}(x)\{S, H\}_{J_{i}}^{2}+\sum_{j}\hat{\gamma}_{\text{port}, j}(x, u)\{S_{\text{tot}}, H_{\text{tot}}\}_{J_{\text{port}, gj}}^{2}\geq0
\]
- **Discrete - gradient conditions**:
\[
\bar{\nabla}H(x, x')^{T}(x' - x)=H(x') - H(x)
\]
\[
\bar{\nabla}H(x, x)=\nabla H(x)
\]
These formulas ensure that the discrete - gradient method can maintain the properties of energy conservation and entropy increase during the discretization process.