Second Order, linear and unconditionally energy stable schemes for a hydrodynamic model of Smectic-A Liquid Crystals

In this paper, we consider the numerical approximations for a hydrodynamical model of smectic-A liquid crystals. The model, derived from the variational approach of the modified Oseen-Frank energy, is a highly nonlinear system that couples the incompressible Navier-Stokes equations and a constitutive equation for the layer variable. We develop two linear, second-order time-marching schemes based on the "Invariant Energy Quadratization" method for nonlinear terms in the constitutive equation, the projection method for the Navier-Stokes equations, and some subtle implicit-explicit treatments for the convective and stress terms. Moreover, we prove the well-posedness of the linear system and their unconditionally energy stabilities rigorously. Various numerical experiments are presented to demonstrate the stability and the accuracy of the numerical schemes in simulating the dynamics under shear flow and the magnetic field

A new class of efficient and robust energy stable schemes for gradient flows

We propose a new numerical technique to deal with nonlinear terms in gradient flows. By introducing a scalar auxiliary variable (SAV), we construct efficient and robust energy stable schemes for a large class of gradient flows. The SAV approach is not restricted to specific forms of the nonlinear part of the free energy, and only requires to solve {\it decoupled} linear equations with {\it constant coefficients}. We use this technique to deal with several challenging applications which can not be easily handled by existing approaches, and present convincing numerical results to show that our schemes are not only much more efficient and easy to implement, but can also better capture the physical properties in these models. Based on this SAV approach, we can construct unconditionally second-order energy stable schemes; and we can easily construct even third or fourth order BDF schemes, although not unconditionally stable, which are very robust in practice. In particular, when coupled with an adaptive time stepping strategy, the SAV approach can be extremely efficient and accurate

Numerical approximations of the Cahn-Hilliard and Allen-Cahn Equations with general nonlinear potential using the Invariant Energy Quadratization approach

In this paper, we carry out stability and error analyses for two first-order, semi-discrete time stepping schemes, which are based on the newly developed Invariant Energy Quadratization approach, for solving the well-known Cahn-Hilliard and Allen-Cahn equations with general nonlinear bulk potentials. Some reasonable sufficient conditions about boundedness and continuity of the nonlinear functional are given in order to obtain optimal error estimates. These conditions are naturally satisfied by two commonly used nonlinear potentials including the double-well potential and regularized logarithmic Flory-Huggins potential. The well-posedness, unconditional energy stabilities and optimal error estimates of the numerical schemes are proved rigorously

Decoupled, Energy Stable Scheme for Hydrodynamic Allen-Cahn Phase Field Moving Contact Line Model

In this paper, we present an efficient energy stable scheme to solve a phase field model incorporating contact line condition. Instead of the usually used Cahn-Hilliard type phase equation, we adopt the Allen-Cahn type phase field model with the static contact line boundary condition that coupled with incompressible Navier-Stokes equations with Navier boundary condition. The projection method is used to deal with the Navier-Stokes equa- tions and an auxiliary function is introduced for the non-convex Ginzburg-Landau bulk potential. We show that the scheme is linear, decoupled and energy stable. Moreover, we prove that fully discrete scheme is also energy stable. An efficient finite element spatial discretization method is implemented to verify the accuracy and efficiency of proposed schemes. Numerical results show that the proposed scheme is very efficient and accurat

Efficient schemes with unconditionally energy stability for the anisotropic Cahn-Hilliard Equation using the stabilized-Scalar Augmented Variable (S-SAV) approach

In this paper, we consider numerical approximations for the anisotropic Cahn-Hilliard equation. The main challenge of constructing numerical schemes with unconditional energy stabilities for this model is how to design proper temporal discretizations for the nonlinear terms with the strong anisotropy. We propose two, second order time marching schemes by combining the recently developed SAV approach with the linear stabilization approach, where three linear stabilization terms are added. These terms are shown to be crucial to remove the oscillations caused by the anisotropic coefficients, numerically. The novelty of the proposed schemes is that all nonlinear terms can be treated semi-explicitly, and one only needs to solve three decoupled linear equations with constant coefficients at each time step. We further prove the unconditional energy stabilities rigorously, and present various 2D and 3D numerical simulations to demonstrate the stability and accuracy

Energy Stable Second Order Linear Schemes for the Allen-Cahn Phase-Field Equation

Phase-field model is a powerful mathematical tool to study the dynamics of interface and morphology changes in fluid mechanics and material sciences. However, numerically solving a phase field model for a real problem is a challenge task due to the non-convexity of the bulk energy and the small interface thickness parameter in the equation. In this paper, we propose two stabilized second order semi-implicit linear schemes for the Allen-Cahn phase-field equation based on backward differentiation formula and Crank-Nicolson method, respectively. In both schemes, the nonlinear bulk force is treated explicitly with two second-order stabilization terms, which make the schemes unconditional energy stable and numerically efficient. By using a known result of the spectrum estimate of the linearized Allen-Cahn operator and some regularity estimate of the exact solution, we obtain an optimal second order convergence in time with a prefactor depending on the inverse of the characteristic interface thickness only in some lower polynomial order. Both 2-dimensional and 3-dimensional numerical results are presented to verify the accuracy and efficiency of proposed schemes.Comment: keywords: energy stable, stabilized semi-implicit scheme, second order scheme, error estimate. related work arXiv:1708.09763, arXiv:1710.0360

Energy Stable Multigrid Method for Local and Non-local Hydrodynamic Models for Freezing

In this paper we present a numerical method for hydrodynamic models that arise from time dependent density functional theories of freezing. The models take the form of compressible Navier-Stokes equations whose pressure is determined by the variational derivative of a free energy, which is a functional of the density field. We present unconditionally energy stable and mass conserving implicit finite difference methods for the models. The methods are based on a convex splitting of the free energy and that ensures that a discrete energy is non-increasing for any choice of time and space step. The methods are applicable to a large class of models, including both local and non-local free energy functionals. The theoretical basis for the numerical method is presented in a general context. The method is applied to problems using two specific free energy functionals: one local and one non-local functional. A nonlinear multigrid method is used to solve the numerical method, which is nonlinear at the implicit time step. The non-local functional, which is a convolution operator, is approximated using the Discrete Fourier Transform. Numerical simulations that confirm the stability and accuracy of the numerical method are presented

An energy stable fourth order finite difference scheme for the Cahn-Hilliard equation

In this paper we propose and analyze an energy stable numerical scheme for the Cahn-Hilliard equation, with second order accuracy in time and the fourth order finite difference approximation in space. In particular, the truncation error for the long stencil fourth order finite difference approximation, over a uniform numerical grid with a periodic boundary condition, is analyzed, via the help of discrete Fourier analysis instead of the the standard Taylor expansion. This in turn results in a reduced regularity requirement for the test function. In the temporal approximation, we apply a second order BDF stencil, combined with a second order extrapolation formula applied to the concave diffusion term, as well as a second order artificial Douglas-Dupont regularization term, for the sake of energy stability. As a result, the unique solvability, energy stability are established for the proposed numerical scheme, and an optimal rate convergence analysis is derived in the $\ell^\infty (0,T; \ell^2) \cap \ell^2 (0,T; H_h^2)$ norm. A few numerical experiments are presented, which confirm the robustness and accuracy of the proposed scheme.Comment: 29 pages, 2 figure

Second Order Linear Energy Stable Schemes for Allen-Cahn Equations with Nonlocal Constraints

We present a set of linear, second order, unconditionally energy stable schemes for the Allen-Cahn equation with nonlocal constraints that preserves the total volume of each phase in a binary material system. The energy quadratization strategy is employed to derive the energy stable semi-discrete numerical algorithms in time. Solvability conditions are then established for the linear systems resulting from the semi-discrete, linear schemes. The fully discrete schemes are obtained afterwards by applying second order finite difference methods on cell-centered grids in space. The performance of the schemes are assessed against two benchmark numerical examples, in which dynamics obtained using the volumepreserving Allen-Cahn equations with nonlocal constraints is compared with those obtained using the classical Allen-Cahn as well as the Cahn-Hilliard model, respectively, demonstrating slower dynamics when volume constraints are imposed as well as their usefulness as alternatives to the Cahn-Hilliard equation in describing phase evolutionary dynamics for immiscible material systems while preserving the phase volumes. Some performance enhancing, practical implementation methods for the linear energy stable schemes are discussed in the end

Linearly decoupled energy-stable numerical methods for multi-component two-phase compressible flow

In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are constructed based on the scalar auxiliary variable (SAV) approaches for Helmholtz free energy and the intermediate velocities that are designed to decouple the tight relationship between velocity and molar densities. The intermediate velocities are also involved in the discrete momentum equation to ensure the consistency with the mass balance equations. Moreover, we propose a component-wise SAV approach for a multi-component fluid, which requires solving a sequence of linear, separate mass balance equations. We prove that the methods preserve the unconditional energy-dissipation feature. Numerical results are presented to verify the effectiveness of the proposed methods.Comment: 22 page