A level set-based solver for two-phase incompressible flows: Extension to magnetic fluids

Paria Makaremi-Esfarjani,Andrew J. Higgins,Alireza Najafi-Yazdi

DOI: https://doi.org/10.1080/10618562.2024.2334073

2024-06-03

Abstract:Development of a two-phase incompressible solver for magnetic flows in the magnetostatic case is presented. The proposed numerical toolkit couples the Navier-Stokes equations of hydrodynamics with Maxwell's equations of electromagnetism to model the behaviour of magnetic flows in the presence of a magnetic field. To this end, a rigorous implementation of a second-order two-phase solver for incompressible nonmagnetic flows is introduced first. This solver is implemented in the finite-difference framework, where a fifth-order conservative level set method is employed to capture the evolution of the interface, along with an incompressible solver based on the projection scheme to model the fluids. The solver demonstrates excellent performance even with high density ratios across the interface (Atwood number $\approx 1$), while effectively preserving the mass conservation property. Subsequently, the numerical discretization of Maxwell's equations under the magnetostatic assumption is described in detail, utilizing the vector potential formulation. The primary second-order solver for two-phase flows is extended to the case of magnetic flows, by incorporating the Lorentz force into the momentum equation, accounting for high magnetic permeability ratios across the interface. The implemented solver is then utilized for examining the deformation of ferrofluid droplets in both quiescent and shear flow regimes across various susceptibility values of the droplets. The results suggest that increasing the susceptibility value of the ferrofluid droplet can affect its deformation and rotation in low capillary regimes. In higher capillary flows, increasing the magnetic permeability jump across the interface can further lead to droplet breakup as well. The effect of this property is also investigated for the Rayleigh-Taylor instability growth in magnetic fluids.

Fluid Dynamics,Computational Physics

What problem does this paper attempt to address?

The problems that this paper attempts to solve mainly focus on developing a numerical tool for simulating the two - phase incompressible flow of magnetohydrodynamics in the presence of a magnetic field. Specifically, the paper aims to:

1. **Develop an efficient two - phase incompressible flow solver**: First, the paper introduces a second - order two - phase incompressible flow solver based on a finite - difference framework. This solver uses a fifth - order conservative level - set method to capture the evolution of the interface and adopts the projection method to simulate the incompressible nature of the fluid. The solver can perform excellently in the case of a high density ratio (Atwood number is about 1) while effectively maintaining mass conservation.

2. **Extend the solver to handle magnetohydrodynamics**: Based on the basic solver, the paper further describes how to numerically discretize Maxwell's equations under the magnetostatic assumption, using the vector potential formula. By incorporating the Lorentz force into the momentum equation, the solver is extended to be able to handle magnetohydrodynamics, especially when the ratio of magnetic permeabilities on both sides of the interface is high.

3. **Verify the performance of the solver**: The paper evaluates the performance of the developed solver through several test cases, including the deformation behavior of magnetohydrodynamic droplets in static and shear flow fields, and the growth of Rayleigh - Taylor instability in magnetohydrodynamics. The research results show that increasing the magnetic susceptibility of magnetohydrodynamic droplets will affect their deformation and rotation behaviors, especially in low - capillary - number flow fields. In high - capillary - number flow fields, increasing the ratio of magnetic permeabilities on both sides of the interface may also lead to droplet breakage.

4. **Explore the influence of different magnetic susceptibilities on the deformation of magnetohydrodynamic droplets**: One of the main contributions of the paper is to study the deformation behavior of magnetohydrodynamic droplets in shear flow fields under different magnetic susceptibility values. Although previous studies have explored the influence of multiple factors on droplet deformation, such as the viscosity ratio between the droplet and the surrounding medium, the strength and direction of the applied magnetic field, the study of deformation behavior under different magnetic susceptibility values is still an important research direction.

### Formula Summary

- **Definition of the level - set function**:
\[
|\phi(x, t)| = |x - x_\Gamma|
\]
where \(x_\Gamma\) represents the distance from point \(x\) to the nearest point on the interface.

- **Transport equation of level - set motion**:
\[
\frac{\partial \phi}{\partial t}+\mathbf{u}\cdot\nabla\phi = 0
\]

- **Definition of the conservative level - set method**:
\[
\psi(x, t_0)=\frac{1}{2}\left(\tanh\left(\frac{\phi(x, t_0)}{2\epsilon}\right)+1\right)
\]
where \(\epsilon\) is a parameter representing the interface thickness.

- **Conservative form of the level - set transport equation**:
\[
\frac{\partial \psi}{\partial t}+\nabla\cdot(\mathbf{u}\psi)=0
\]

- **Level - set re - initialization equation**:
\[
\frac{\partial \psi}{\partial \tau}+\nabla\cdot(\psi(1 - \psi)\mathbf{n})=\nabla\cdot(\epsilon(\nabla\psi\cdot\mathbf{n})\mathbf{n})
\]
where \(\tau\) is the pseudo - time, \(\mathbf{n}\) is the normal vector at the initial time, calculated as \(\mathbf{n}=\frac{\nabla\psi}{|\nabla\psi|}\).

These formulas and methods together form the core of the numerical tool proposed in the paper for efficiently and accurately simulating the two - phase incompressible flow of magnetohydrodynamics.