Abstract:In this research, we introduce and investigate an approximation method that preserves the structural integrity of the non-isothermal Cahn-Hilliard-Navier-Stokes system. Our approach extends a previously proposed technique [1], which utilizes conforming (inf-sup stable) finite elements in space, coupled with implicit time discretization employing convex-concave splitting. Expanding upon this method, we incorporate the unstable P1|P1 pair for the Navier-Stokes contributions, integrating Brezzi-Pitkäranta stabilization. Additionally, we improve the enforcement of incompressibility conditions through grad div stabilization. While these techniques are well-established for Navier-Stokes equations, it becomes apparent that for non-isothermal models, they introduce additional coupling terms to the equation governing internal energy. To ensure the conservation of total energy and maintain entropy production, these stabilization terms are appropriately integrated into the internal energy equation.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a numerical method for solving the non - isothermal Cahn - Hilliard - Navier - Stokes (CHNST) system. Specifically, the researchers hope to introduce and study an approximation method that can maintain the structural integrity of the system, ensuring that physical properties such as mass conservation, total energy conservation, and entropy production are still retained during the discretization process.
### Specific problem description:
1. **Numerical simulation of non - isothermal two - phase flow models**:
- The CHNST system is widely used to study various phenomena, such as two - phase flow and fluid - phase coupling interactions. These phenomena are of great significance in scientific and industrial fields, for example, in additive manufacturing (3D printing) and ink - jet printing.
- In particular, in the powder bed fusion additive manufacturing (PBF - AM) process, the non - isothermal CHNST model is used to describe coupling processes such as fluid - phase interactions, melting flow dynamics, and heat transfer.
2. **Improvement of numerical methods**:
- Previous methods used the inf - sup stable finite element method for spatial discretization, combined with implicit time discretization and convex - concave splitting techniques. This paper further extends on this basis, introducing the unstable P1|P1 element pair to handle the Navier - Stokes equations, and combining it with the Brezzi - Pitkäranta stabilization method.
- At the same time, the incompressible condition constraint is strengthened through the grad div stabilization method. Although these techniques are well - established in the Navier - Stokes equations, in non - isothermal models, they will introduce additional coupling terms into the internal energy equation to ensure total energy conservation and entropy production.
### Mathematical models and equations:
The partial differential equations of this system are as follows:
\[
\begin{aligned}
& \partial_t \phi+\mathbf{u} \cdot \nabla \phi-\nabla \cdot (L_{11} \nabla \mu + L_{12} \nabla \theta)=0, \\
& \mu = -\gamma \Delta \phi+\frac{\partial \Psi(\phi, \theta)}{\partial \phi}, \\
& \partial_t e+\mathbf{u} \cdot \nabla e-\nabla \cdot (L_{12} \nabla \mu - L_{22} \nabla \theta)-(\eta D \mathbf{u}-\sigma): \nabla \mathbf{u}=0, \\
& \partial_t \mathbf{u}+(\mathbf{u} \cdot \nabla) \mathbf{u}-\nabla \cdot (\eta D \mathbf{u}-pI - \sigma)=0, \\
& \nabla \cdot \mathbf{u}=0.
\end{aligned}
\]
where:
- \(\phi\) represents the conservative phase field,
- \(\mathbf{u}\) represents the flow velocity,
- \(\theta\) represents the inverse temperature,
- \(e \equiv e(\phi, \theta)\) represents the internal energy.
### Main contributions of the research:
1. **Structure preservation**: The proposed method can maintain the structural integrity of the system during the discretization process, especially mass conservation, total energy conservation, and entropy production.
2. **Treatment of unstable finite element pairs**: The Brezzi - Pitkäranta stabilization method is introduced to handle the unstable P1|P1 element pairs in the Navier - Stokes equations.
3. **Theoretical verification**: The effectiveness of the method is verified through appropriate convergence tests, and its application effects in two - dimensional and three - dimensional simulations are demonstrated.
In summary, this research aims at the non - isothermal Cahn - Hilliard - Navier - Stokes system.