Phase-field models are widely used in studying multiphase flow dynamics.Given the complexity and strong nonlinearity,designing accurate,efficient,and stable numerical algorithms to solve these models has been an activ...Phase-field models are widely used in studying multiphase flow dynamics.Given the complexity and strong nonlinearity,designing accurate,efficient,and stable numerical algorithms to solve these models has been an active research field for decades.This paper proposes a novel numerical scheme to solve a highly cited and used phase field hydrodynamic model for simulating ternary phase fluid flows.The main novelty is the introduction of a supplementary variable to reformulate the original problem into a constrained optimization problem.This reformulation leads to several advantages for our proposed numerical algorithms compared with many existing numerical techniques for solving this model.First,the developed schemes allow more straightforward calculations for the hydrodynamic phase-field models by solving a few decoupled Helmholtz or Poisson-type systems with a constant precomputable coefficient matrix,remarkably reducing the computational cost.Secondly,the numerical schemes can maintain mass conservation and energy dissipation at the discrete level.Additionally,the developed scheme based on the secondorder backward difference formula respects the original energy dissipation law that differs from many existing schemes,such as the IEQ,SAV,and Lagrange multiplier approaches for which a modified energy dissipation law is respected.Furthermore,rigorous proof of energy stability and practical implementation strategies are provided.We conduct adequate 2D and 3D numerical tests to demonstrate the proposed schemes’accuracy and effectiveness.展开更多
基金partially supported by the National Natural Science Foundation of China(Grant No.12201297)by the LCP Fund for Young Scholar(Grant No.6142A05QN22005)+4 种基金partially supported by the Foundation of Jiangsu Key Laboratory for Numerical Simulation of Large Scale Complex Systems,China(Grant No.202002)by the Fundamental Research Funds for the Central Universities(Grant No.NS2022070)by the Natural Science Foundation of Jiangsu Province(Grant No.BK20220131)by the National Natural Science Foundation of China(Grant Nos.12271252,12071216)partially supported by the National Science Foundation(NSF)USA with DMS-2111479.
文摘Phase-field models are widely used in studying multiphase flow dynamics.Given the complexity and strong nonlinearity,designing accurate,efficient,and stable numerical algorithms to solve these models has been an active research field for decades.This paper proposes a novel numerical scheme to solve a highly cited and used phase field hydrodynamic model for simulating ternary phase fluid flows.The main novelty is the introduction of a supplementary variable to reformulate the original problem into a constrained optimization problem.This reformulation leads to several advantages for our proposed numerical algorithms compared with many existing numerical techniques for solving this model.First,the developed schemes allow more straightforward calculations for the hydrodynamic phase-field models by solving a few decoupled Helmholtz or Poisson-type systems with a constant precomputable coefficient matrix,remarkably reducing the computational cost.Secondly,the numerical schemes can maintain mass conservation and energy dissipation at the discrete level.Additionally,the developed scheme based on the secondorder backward difference formula respects the original energy dissipation law that differs from many existing schemes,such as the IEQ,SAV,and Lagrange multiplier approaches for which a modified energy dissipation law is respected.Furthermore,rigorous proof of energy stability and practical implementation strategies are provided.We conduct adequate 2D and 3D numerical tests to demonstrate the proposed schemes’accuracy and effectiveness.