Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTwo linear, unconditionally stable, second order decoupling methods for the Allen--Cahn--Navier--Stokes phase field model
83
0
0.0
(
0
)
Added by
Ruonan Cao
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Hydrodynamics coupled phase field models have intricate difficulties to solve numerically as they feature high nonlinearity and great complexity in coupling. In this paper, we propose two second order, linear, unconditionally stable decoupling methods based on the Crank--Nicolson leap-frog time discretization for solving the Allen--Cahn--Navier--Stokes (ACNS) phase field model of two-phase incompressible flows. The ACNS system is decoupled via the artificial compression method and a splitting approach by introducing an exponential scalar auxiliary variable. We prove both algorithms are unconditionally long time stable. Numerical examples are provided to verify the convergence rate and unconditional stability.
rate research
Read More
We propose a novel second order in time numerical scheme for Cahn-Hilliard-Navier- Stokes phase field model with matched density. The scheme is based on second order convex-splitting for the Cahn-Hilliard equation and pressure-projection for the Navier-Stokes equation. We show that the scheme is mass-conservative, satisfies a modified energy law and is therefore unconditionally stable. Moreover, we prove that the scheme is uncondition- ally uniquely solvable at each time step by exploring the monotonicity associated with the scheme. Thanks to the weak coupling of the scheme, we design an efficient Picard iteration procedure to further decouple the computation of Cahn-Hilliard equation and Navier-Stokes equation. We implement the scheme by the mixed finite element method. Ample numerical experiments are performed to validate the accuracy and efficiency of the numerical scheme.
A time-fractional Allen-Cahn equation with volume constraint is first proposed by introducing a nonlocal time-dependent Lagrange multiplier. Adaptive linear second-order energy stable schemes are developed for the proposed model by combining invariant energy quadratization and scalar auxiliary variable approaches with the recent L1$^{+}$ formula. The new developed methods are proved to be volume-preserving and unconditionally energy stable on arbitrary nonuniform time meshes. The accelerated algorithm and adaptive time strategy are employed in numerical implement. Numerical results show that the proposed algorithms are computationally efficient in multi-scale simulations, and appropriate for accurately resolving the intrinsically initial singularity of solution and for efficiently capturing the fast dynamics away initial time.
A mathematical model describing the flow of two-phase fluids in a bounded container $Omega$ is considered under the assumption that the phase transition process is influenced by inertial effects. The model couples a variant of the Navier-Stokes system for the velocity $u$ with an Allen-Cahn-type equation for the order parameter $varphi$ relaxed in time in order to introduce inertia. The resulting model is characterized by second-order material derivatives which constitute the main difficulty in the mathematical analysis. Actually, in order to obtain a tractable problem, a viscous relaxation term is included in the phase equation. The mathematical results consist in existence of weak solutions in 3D and, under additional assumptions, existence and uniqueness of strong solutions in 2D. A partial characterization of the long-time behavior of solutions is also given and in particular some issues related to dissipation of energy are discussed.
We present several first-order and second-order numerical schemes for the Cahn-Hilliard equation with discrete unconditional energy stability. These schemes stem from the generalized Positive Auxiliary Variable (gPAV) idea, and require only the solution of linear algebraic systems with a constant coefficient matrix. More importantly, the computational complexity (operation count per time step) of these schemes is approximately a half of those of the gPAV and the scalar auxiliary variable (SAV) methods in previous works. We investigate the stability properties of the proposed schemes to establish stability bounds for the field function and the auxiliary variable, and also provide their error analyses. Numerical experiments are presented to verify the theoretical analyses and also demonstrate the stability of the schemes at large time step sizes.
The Allen-Cahn equation is solved numerically by operator splitting Fourier spectral methods. The basic idea of the operator splitting method is to decompose the original problem into sub-equations and compose the approximate solution of the original equation using the solutions of the subproblems. Unlike the first and the second order methods, each of the heat and the free-energy evolution operators has at least one backward evaluation in higher order methods. We investigate the effect of negative time steps on a general form of third order schemes and suggest three third order methods for better stability and accuracy. Two fourth order methods are also presented. The traveling wave solution and a spinodal decomposition problem are used to demonstrate numerical properties and the order of convergence of the proposed methods.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
