A Positivity Preserving, Energy Stable Finite Difference Scheme for the Flory-Huggins-Cahn-Hilliard-Navier-Stokes System
In this paper, we propose and analyze a finite difference numerical scheme for the Cahn-Hilliard-Navier-Stokes system, with logarithmic Flory-Huggins energy potential. in the numerical approximation to the singular chemical potential, the logarithmic term and the surface diffusion term are implicitly updated, while an explicit computation is applied to the concave expansive term. Moreover, the convective term in the phase field evolutionary equation is approximated in a semi-implicit manner. Similarly, the fluid momentum equation is computed by a semi-implicit algorithm: implicit treatment for the kinematic diffusion term, explicit update for the pressure gradient, combined with semi-implicit approximations to the fluid convection and the phase field coupled term, respectively. Such a semi-implicit method gives an intermediate velocity field. Subsequently, a Helmholtz projection into the divergence-free vector field yields the velocity vector and the pressure variable at the next time step. This approach decouples the Stokes solver, which in turn drastically improves the numerical efficiency. the positivity-preserving property and the unique solvability of the proposed numerical scheme is theoretically justified, i.e., the phase variable is always between -1 and 1, following the singular nature of the logarithmic term as the phase variable approaches the singular limit values. in addition, an iteration construction technique is applied in the positivity-preserving and unique solvability analysis, motivated by the non-symmetric nature of the fluid convection term. the energy stability of the proposed numerical scheme could be derived by a careful estimate. a few numerical results are presented to validate the robustness of the proposed numerical scheme.
W. Chen et al., "A Positivity Preserving, Energy Stable Finite Difference Scheme for the Flory-Huggins-Cahn-Hilliard-Navier-Stokes System," Journal of Scientific Computing, vol. 92, no. 2, article no. 31, Springer, Aug 2022.
The definitive version is available at https://doi.org/10.1007/s10915-022-01872-1
Mathematics and Statistics
Keywords and Phrases
Cahn-Hilliard-Navier-Stokes System; Energy Stability; Flory Huggins Energy Potential; Monotonicity Analysis; Numerical Accuracy; Positivity Preserving
International Standard Serial Number (ISSN)
Article - Journal
© 2023 Springer, All rights reserved.
01 Aug 2022
Fudan University, Grant None