Abstract:A Cahn-Hilliard-Allen-Cahn phase-field model coupled with a heat transfer equation, particularly with full non-diagonal mobility matrices, is studied. After reformulating the problem w.r.t. the inverse of temperature, we proposed and analysed a structure-preserving approximation for the semi-discretisation in space and then a fully discrete approximation using conforming finite elements and time-stepping methods. We prove structure-preserving property and discrete stability using relative entropy methods for the semi-discrete and fully discrete case. The theoretical results are illustrated by numerical experiments.
What problem does this paper attempt to address?
This paper attempts to solve the problems of structure preservation and nonlinear stability in the phase - field model under non - isothermal conditions. Specifically, the author studied a coupled Cahn - Hilliard - Allen - Cahn phase - field model and introduced the heat transfer equation. This system is characterized by a fully non - diagonal mobility matrix, which makes the model more complex and closer to the actual situation.
### Main problem description
1. **Elimination of artificial interface effects**: Traditional phase - field models have problems such as violating conservation laws and discontinuity of chemical/thermodynamic potentials at the interface. In order to eliminate these artificial interface effects, the author proposed a variational quantitative phase - field model, which contains cross - coupling terms.
2. **Structure preservation and stability**: The author aims to develop a structure - preserving discretization method to ensure that the numerical solution maintains the physical properties of the system (such as mass, energy conservation and entropy production) during the discretization process in space and time.
3. **Nonlinear stability estimation**: Through the relative entropy method, the author proved the nonlinear stability estimations of the semi - discrete and fully - discrete schemes. This is crucial for verifying the effectiveness and reliability of the numerical method.
### Model introduction
The temperature - dependent Cahn - Hilliard - Allen - Cahn system considered in this paper is:
\[
\frac{\partial \rho}{\partial t}=\nabla \cdot\left(L_{11} \nabla \mu_{\rho}-L_{12} \nabla \frac{1}{T}+L_{13} \mu_{\eta}\right), \quad \mu_{\rho}=-\gamma_{\rho} \Delta \rho+\frac{1}{T} \frac{\partial f}{\partial \rho}
\]
\[
\frac{\partial e(\rho, T, \eta)}{\partial t}=\nabla \cdot\left(L_{12} \nabla \mu_{\rho}-L_{22} \nabla \frac{1}{T}+L_{23} \mu_{\eta}\right)
\]
\[
\frac{\partial \eta}{\partial t}=-L_{13} \cdot \nabla \mu_{\rho}+L_{23} \nabla \frac{1}{T}-L_{33} \mu_{\eta}, \quad \mu_{\eta}=-\gamma_{\eta} \Delta \eta+\frac{1}{T} \frac{\partial f}{\partial \eta}
\]
where $\rho$ is the conservative order parameter, $\eta_i$ is the non - conservative order parameter, $f$ is the free - energy density, $e$ is the internal - energy density, and $s$ is the entropy density.
### Solutions
1. **Variable transformation**: The author introduced new variables $\theta = 1/T$, $\tilde{\mu}_{\rho}=\mu_{\rho}/T$ and $\tilde{\mu}_{\eta}=\mu_{\eta}/T$ to simplify the problem and reveal its variational structure.
2. **Discretization method**: Spatial discretization is carried out by the finite element method, and time discretization is carried out by the implicit time - stepping method. In particular, the author adopted the classical convex - concave splitting method to deal with the driving potential function.
3. **Nonlinear stability estimation**: The nonlinear stability estimations of the semi - discrete and fully - discrete schemes are proved by using the relative entropy method, thus ensuring the stability and convergence of the numerical solution.
### Summary
This paper reformulates the problem by introducing the inverse temperature as the main variable and proposes a structure - preserving discretization method. This method can not only eliminate the artificial interface effects in the traditional phase - field model, but also maintain the physical properties of the system numerically. In addition, the author proved the nonlinear stability estimations by the relative entropy method, providing a theoretical basis for the effectiveness and reliability of the numerical method.