In this paper, we consider numerical approximation of an electrically conductive ferrofluid model, which consists of Navier-Stokes equations, magnetization equation, and magnetic induction equation. To solve this highly coupled, nonlinear, and multiphysics system efficiently, we develop a decoupled, linear, second-order in time, and unconditionally energy stable finite element scheme. We incorporate several distinct numerical techniques, including reformulations of the equations and a scalar auxiliary variable to handle the coupled nonlinear terms,a symmetric implicit-explicit treatment for the symmetric positive definite nonlinearity, and stable finite element approximations. We also prove that the numerical scheme is provably uniquely solvable and unconditionally energy stable rigorously. A series of numerical examples are presented to illustrate the accuracy and performance of our scheme.
翻译:本文研究一种导电铁磁流体模型的数值逼近方法,该模型由Navier-Stokes方程、磁化方程和磁感应方程构成。为高效求解这一强耦合、非线性、多物理场系统,我们提出了一种解耦、线性、时间二阶精度且无条件能量稳定的有限元格式。我们融合了多种独特的数值技术,包括:通过方程重构与标量辅助变量处理耦合非线性项;对对称正定非线性项采用对称隐式-显式处理;以及稳定的有限元逼近。我们严格证明了该数值格式具有可证明的唯一可解性与无条件能量稳定性。通过一系列数值算例验证了所提格式的精度与性能。