Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn energy stable finite element scheme for the three-component Cahn-Hilliard-type model for macromolecular microsphere composite hydrogels
243
0
0.0
(
0
)
Added by
Maoqin Yuan
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this article, we present and analyze a finite element numerical scheme for a three-component macromolecular microsphere composite (MMC) hydrogel model, which takes the form of a ternary Cahn-Hilliard-type equation with Flory-Huggins-deGennes energy potential. The numerical approach is based on a convex-concave decomposition of the energy functional in multi-phase space, in which the logarithmic and the nonlinear surface diffusion terms are treated implicitly, while the concave expansive linear terms are explicitly updated. A mass lumped finite element spatial approximation is applied, to ensure the positivity of the phase variables. In turn, a positivity-preserving property can be theoretically justified for the proposed fully discrete numerical scheme. In addition, unconditional energy stability is established as well, which comes from the convexity analysis. Several numerical simulations are carried out to verify the accuracy and positivity-preserving property of the proposed scheme.
rate research
Read More
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.
We present and analyze an unconditionally energy stable and convergent finite difference scheme for the Functionalized Cahn-Hilliard equation. One key difficulty associated with the energy stability is based on the fact that one nonlinear energy functional term in the expansion appears as non-convex, non-concave. To overcome this subtle difficulty, we add two auxiliary terms to make the combined term convex, which in turns yields a convex-concave decomposition of the physical energy. As a result, an application of the convex splitting methodology assures both the unique solvability and the unconditional energy stability of the proposed numerical scheme. To deal with a 4-Laplacian solver in an $H^{-1}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm to solve the corresponding nonlinear systems. In addition, a global in time $H_{rm per}^2$ stability of the numerical scheme is established at a theoretical level, which in turn ensures the full order convergence analysis of the scheme. A few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme.
In this paper, we construct and analyze a uniquely solvable, positivity preserving and unconditionally energy stable finite-difference scheme for the periodic three-component Macromolecular Microsphere Composite (MMC) hydrogels system, a ternary Cahn-Hilliard system with a Flory-Huggins-deGennes free energy potential. The proposed scheme is based on a convex-concave decomposition of the given energy functional with two variables, and the centered difference method is adopted in space. We provide a theoretical justification that this numerical scheme has a pair of unique solutions, such that the positivity is always preserved for all the singular terms, i.e., not only two phase variables are always between $0$ and $1$, but also the sum of two phase variables is between $0$ and $1$, at a point-wise level. In addition, we use the local Newton approximation and multigrid method to solve this nonlinear numerical scheme, and various numerical results are presented, including the numerical convergence test, positivity-preserving property test, energy dissipation and mass conservation properties.
A thermodynamically consistent phase-field model is introduced for simulating motion and shape transformation of vesicles under flow conditions. In particular, a general slip boundary condition is used to describe the interaction between vesicles and the wall of the fluid domain. A second-order accurate in both space and time C0 finite element method is proposed to solve the model governing equations. Various numerical tests confirm the convergence, energy stability, and conservation of mass and surface area of cells of the proposed scheme. Vesicles with different mechanical properties are also used to explain the pathological risk for patients with sickle cell disease.
We present a second-order-in-time finite difference scheme for the Cahn-Hilliard-Hele-Shaw equations. This numerical method is uniquely solvable and unconditionally energy stable. At each time step, this scheme leads to a system of nonlinear equations that can be efficiently solved by a nonlinear multigrid solver. Owing to the energy stability, we derive an $ell^2 (0,T; H_h^3)$ stability of the numerical scheme. To overcome the difficulty associated with the convection term $ abla cdot (phi boldsymbol{u})$, we perform an $ell^infty (0,T; H_h^1)$ error estimate instead of the classical $ell^infty (0,T; ell^2)$ one to obtain the optimal rate convergence analysis. In addition, various numerical simulations are carried out, which demonstrate the accuracy and efficiency of the proposed numerical scheme.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
