Decoupled, Linear, and Unconditionally Energy Stable Fully Discrete Finite Element Numerical Scheme for a Two-Phase Ferrohydrodynamics Model
We consider in this paper numerical approximations of a phase field model for twophase ferrofluids, which consists of the Navier-Stokes equations, the Cahn-Hilliard equation, the magnetostatic equations, and the magnetic field equation. By combining the projection method for the Navier-Stokes equations and some subtle implicit-explicit treatments for coupled nonlinear terms, we construct a linear, decoupled, fully discrete finite element scheme to solve the highly nonlinear and coupled multiphysics system efficiently. The scheme is provably unconditionally energy stable and leads to a series of decoupled linear equations to solve at each time step. Through numerous numerical examples in simulating benchmark problems such as the Rosensweig instability and droplet deformation, we demonstrate the stability and accuracy of the numerical scheme.
G. D. Zhang et al., "Decoupled, Linear, and Unconditionally Energy Stable Fully Discrete Finite Element Numerical Scheme for a Two-Phase Ferrohydrodynamics Model," SIAM Journal on Scientific Computing, vol. 43, no. 1, pp. B167 - B193, Society for Industrial and Applied Mathematics (SIAM), Jan 2021.
The definitive version is available at https://doi.org/10.1137/19M1288280
Mathematics and Statistics
Center for High Performance Computing Research
Keywords and Phrases
Ferrofluid; Ferrohydrodynamics; Magnetic field; Phase field; Unconditional energy stability
International Standard Serial Number (ISSN)
Article - Journal
© 2021 Society for Industrial and Applied Mathematics (SIAM), All rights reserved.
01 Jan 2021