Next Article in Journal
Fuzzy Fractional-Order PD Vibration Control of Uncertain Building Structures
Next Article in Special Issue
Fractional View Analysis of Swift–Hohenberg Equations by an Analytical Method and Some Physical Applications
Previous Article in Journal
Finite Iterative Forecasting Model Based on Fractional Generalized Pareto Motion
Previous Article in Special Issue
Exact Solutions of the Nonlinear Modified Benjamin-Bona-Mahony Equation by an Analytical Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New L2-Gradient Flow-Based Fractional-in-Space Modified Phase-Field Crystal Equation and Its Mass Conservative and Energy Stable Method

Department of Mathematics, Kwangwoon University, Seoul 01897, Korea
Fractal Fract. 2022, 6(9), 472; https://doi.org/10.3390/fractalfract6090472
Submission received: 19 July 2022 / Revised: 7 August 2022 / Accepted: 24 August 2022 / Published: 27 August 2022
(This article belongs to the Special Issue Recent Advances in Computational Physics with Fractional Application)

Abstract

:
In this paper, we introduce a new fractional-in-space modified phase-field crystal equation based on the L 2 -gradient flow approach, where the mass of atoms is conserved by using a nonlocal Lagrange multiplier. To solve the L 2 -gradient flow-based fractional-in-space modified phase-field crystal equation, we present a mass conservative and energy stable method based on the convex splitting idea. Numerical examples together with standard tests in the classical H 1 -gradient flow-based modified phase-field crystal equation are provided to illustrate the applicability of the proposed framework.

1. Introduction

The modified phase-field crystal (MPFC) equation was introduced to describe both diffusive dynamics and elastic interactions [1]:
2 ϕ t 2 + β ϕ t = Δ μ , μ : = δ E δ ϕ = ϕ 3 + ( 1 ϵ ) ϕ + 2 Δ ϕ + Δ 2 ϕ ,
where ϕ is the density field, β > 0 is a constant, μ is the chemical potential, δ δ ϕ denotes the variational derivative, and E is the Swift–Hohenberg energy functional [2]:
E ( ϕ ) : = Ω 1 4 ϕ 4 + 1 ϵ 2 ϕ 2 | ϕ | 2 + 1 2 ( Δ ϕ ) 2 d x ,
where 0 < ϵ < 1 is a constant with physical significance. We assume that ϕ and μ are periodic on Ω . Note that (i) d d t Ω ϕ d x = Ω ϕ t d x = 0 if we employ an initial condition satisfying Ω ϕ t ( x , 0 ) d x = 0 ; (ii) solutions of the MPFC equation dissipate the following energy [3]:
E ˜ ( ϕ ) : = E ( ϕ ) + 1 2 ϕ t H 1 2 ,
where ( · , · ) H 1 denotes the H 1 -inner product with respect to Ω .
Fractional differential equations have been proved to be valuable tools for modeling diffusive processes associated with anomalous diffusion [4,5,6,7,8]. To describe superdiffusion in the MPFC model, we consider a fractional-in-space version of the MPFC equation
2 ϕ t 2 + β ϕ t = Δ μ s , μ s : = δ E s δ ϕ = ϕ 3 + ( 1 ϵ ) ϕ 2 ( Δ ) s 2 ϕ + ( Δ ) s ϕ ,
where ( Δ ) s 2 is the fractional Laplacian ( 1 < s 2 ), and
E s ( ϕ ) : = Ω 1 4 ϕ 4 + 1 ϵ 2 ϕ 2 | ( Δ ) s 4 ϕ | 2 + 1 2 ( ( Δ ) s 2 ϕ ) 2 d x .
Note that the fractional-in-space MPFC equation reduces to the MPFC equation when s = 2 .
Various methods have been proposed to solve the MPFC equation [9,10,11,12,13,14,15]. In [9,12], the authors presented second-order methods based on the convex splitting idea. Galenko et al. [10] introduced a second-order method with a higher-order generalization of the trapezoidal rule. In [13], Li et al. developed second-order methods based on the invariant energy quadratization approach. Li and Shen [14] constructed a second-order method based on the scalar auxiliary variable approach. In [15], Shin et al. presented high-order (up to fourth-order) methods by combining the concept of energy quadratization and the Runge–Kutta method. Note that the MPFC equation needs to discretize the linear sixth-order, nonlinear second-order, and inverse Laplacian terms.
The aim of this paper is to reformulate the fractional-in-space MPFC equation through the variational approach in the L 2 -space, where a nonlocal Lagrange multiplier is added to cancel out the variation of mass without influencing the energy dissipation property. Then, the new L 2 -gradient flow-based fractional-in-space MPFC ( L 2 -fMPFC) equation only needs to discretize the linear fourth-order and nonlinear terms. To solve the L 2 -fMPFC equation with preserving the mass conservation and energy stability, we propose a numerical method by applying a convex splitting to both the L 2 -gradient flow and nonlocal terms. In addition, to have second-order time accuracy, the convex and concave parts are treated using a second-order secant type approach and a second-order extrapolation, respectively. We prove that the method is mass conservative, unconditionally uniquely solvable, and unconditionally stable with respect to the energy of the L 2 -fMPFC equation.
The organization of this paper is as follows. In Section 2, we introduce the L 2 -fMPFC equation and its energy. In Section 3, we develop the second-order method for the L 2 -fMPFC equation and give proofs of its mass conservation, unconditional unique solvability, and unconditional energy stability. In Section 4, we provide numerical examples illustrating the applicability of the proposed framework. Finally, conclusions are drawn in Section 5.

2. L 2 -Gradient Flow-Based Fractional-In-Space Modified Phase-Field Crystal Equation

We introduce the new L 2 -fMPFC equation:
2 ϕ t 2 + β ϕ t = μ s + α ( t ) ,
where α ( t ) = 1 | Ω | Ω μ s d x is the nonlocal Lagrange multiplier. In addition, we introduce the energy:
F ( ϕ ) : = E s ( ϕ ) + 1 2 ϕ t L 2 2 ,
where ( · , · ) L 2 denotes the L 2 -inner product with respect to Ω .
Integrating Equation (6) over Ω and letting Ψ ( t ) = Ω ϕ t d x , we obtain
d Ψ ( t ) d t + β Ψ ( t ) = Ω μ s d x + Ω μ s d x = 0 .
A solution of this equation is Ψ ( t ) = Ψ ( 0 ) e β t . Thus, if we employ an initial condition satisfying Ψ ( 0 ) = 0 , then
Ψ ( t ) = 0 , t 0 .
The energy F ( ϕ ) is nonincreasing in time:
d F d t = μ s , ϕ t L 2 + ϕ t , 2 ϕ t 2 L 2 = ϕ t , β ϕ t + α ( t ) L 2 = β ϕ t L 2 2 + α ( t ) Ψ ( t ) = β ϕ t L 2 2 0 .

3. Mass Conservative and Energy Stable Method

Using a new variable ψ , Equation (6) can be split as
ψ t = μ s + α ( t ) β ψ ,
ϕ t = ψ
and the energy F ( ϕ ) can be redefined as
F ( ϕ , ψ ) : = E s ( ϕ ) + 1 2 ψ L 2 2 .
Then, Equations (8) and (9) can be discretized with second-order time accuracy as follows:
ψ n + 1 ψ n Δ t = μ s , n + 1 2 + α n + 1 2 β ψ n + 1 + ψ n 2 ,
μ s , n + 1 2 = χ ( ϕ n + 1 , ϕ n ) + ( Δ ) s ϕ n + 1 + ϕ n 2 ( Δ ) s 2 ( 3 ϕ n ϕ n 1 ) ,
ϕ n + 1 ϕ n Δ t = ψ n + 1 + ψ n 2 ,
where χ ( ϕ , φ ) = ϕ 2 + φ 2 2 ϕ + φ 2 + ( 1 ϵ ) ϕ + φ 2 , α n + 1 2 = 1 | Ω | μ s , n + 1 2 , 1 L 2 , and ϕ 1 ϕ 0 . Note that we discretize μ s , n + 1 2 and α n + 1 2 based on the following convex splitting of the energy E s ( ϕ ) [16]:
E s ( ϕ ) = E c s ( ϕ ) E e s ( ϕ ) = Ω 1 4 ϕ 4 + 1 ϵ 2 ϕ 2 + 1 2 ( ( Δ ) s 2 ϕ ) 2 d x Ω | ( Δ ) s 4 ϕ | 2 d x .
Theorem 1.
The second-order method (11)–(13) with an initial condition satisfying Ψ(0) = 0 is mass conservative.
Proof. 
From Equation (11), we obtain
( ψ n + 1 ψ n , 1 ) L 2 = Δ t ( μ s , n + 1 2 + α n + 1 2 , 1 ) L 2 β Δ t 2 ( ψ n + 1 + ψ n , 1 ) L 2 = β Δ t 2 ( ψ n + 1 + ψ n , 1 ) L 2 .
This yields the relation
( ψ n + 1 , 1 ) L 2 = 2 β Δ t 2 + β Δ t ( ψ n , 1 ) L 2 .
With an initial condition satisfying Ψ ( 0 ) = 0 , i.e., ( ψ 0 , 1 ) L 2 = 0 , the relation guarantees that ( ψ n + 1 , 1 ) L 2 = 0 for all n 0 . Now, from Equation (13), we see that
( ϕ n + 1 ϕ n , 1 ) L 2 = 0 if and only if ( ψ n + 1 , 1 ) L 2 = 0 ,
and the result follows: ( ϕ n + 1 , 1 ) L 2 = ( ϕ n , 1 ) L 2 . □
Before proving the unique solvability of the second-order method (11)–(13), we simplify the method as follows:
ϕ n + 1 ϕ n Δ t = Δ t 2 + β Δ t μ s , n + 1 2 α n + 1 2 + 2 2 + β Δ t ψ n ,
ψ n + 1 = 2 ϕ n + 1 ϕ n Δ t ψ n .
Theorem 2.
The second-order method (12) and (14)–(15) with an initial condition satisfying Ψ(0) = 0 is unconditionally uniquely solvable.
Proof. 
We consider the following functional for ϕ defined in the constraint space ( ϕ , 1 ) L 2 = ( ϕ n , 1 ) L 2 given by Theorem 1:
G ( ϕ ) = 1 2 Δ t ϕ ϕ n L 2 2 + Δ t 2 + β Δ t H ( ϕ ) 2 2 + β Δ t ( ψ n , ϕ ) L 2 ,
where
H ( ϕ ) = 1 4 ϕ 4 4 + ϕ 3 3 ϕ n + ϕ 2 2 ( ϕ n ) 2 + ϕ ( ϕ n ) 3 , 1 L 2 + 1 ϵ 2 ϕ 2 2 + ϕ ϕ n , 1 L 2 + 1 4 ( Δ ) s 2 ϕ L 2 2 + 1 2 ( Δ ) s ϕ n ( Δ ) s 2 ( 3 ϕ n ϕ n 1 ) , ϕ L 2 .
It may be shown that ϕ n + 1 is the unique minimizer of G ( ϕ ) if and only if it solves, for any φ with ( φ , 1 ) L 2 = 0 ,
d G ( ϕ + η φ ) d η η = 0 = ϕ ϕ n Δ t , φ L 2 + Δ t 2 + β Δ t χ ( ϕ , ϕ n ) + 1 2 ( Δ ) s ϕ , φ L 2 + Δ t 2 + β Δ t 1 2 ( Δ ) s ϕ n ( Δ ) s 2 ( 3 ϕ n ϕ n 1 ) , φ L 2 2 2 + β Δ t ( ψ n , φ ) L 2 = 0 ,
because G ( ϕ ) is strictly convex by
d 2 G ( ϕ + η φ ) d η 2 η = 0 = 1 Δ t φ L 2 2 + Δ t 2 + β Δ t ϕ 2 2 + 1 4 ( ϕ + ϕ n ) 2 + 1 ϵ 2 , φ 2 L 2 + Δ t 2 ( 2 + β Δ t ) ( Δ ) s 2 φ L 2 2 0 .
Since 1 is in the null space of the inner product with φ , the solution of (16) can be rewritten as follows:
ϕ ϕ n Δ t = Δ t 2 + β Δ t χ ( ϕ , ϕ n ) + ( Δ ) s ϕ + ϕ n 2 ( Δ ) s 2 ( 3 ϕ n ϕ n 1 ) + 2 2 + β Δ t ψ n + Λ 1 .
Using the constraint ( ϕ , 1 ) L 2 = ( ϕ n , 1 ) L 2 and ( ψ n , 1 ) L 2 = 0 given by Theorem 1, we have Λ = Δ t ( 2 + β Δ t ) | Ω | χ ( ϕ , ϕ n ) + ( Δ ) s ϕ + ϕ n 2 ( Δ ) s 2 ( 3 ϕ n ϕ n 1 ) , 1 L 2 , which concludes the proof. □
Theorem 3.
The second-order method (12) and (14)–(15) with an initial condition satisfying Ψ(0) = 0 is unconditionally energy stable,
F ˜ ( ϕ n + 1 , ψ n + 1 , ϕ n ) F ˜ ( ϕ n , ψ n , ϕ n 1 ) ,
where
F ˜ ( ϕ , ψ , φ ) : = F ( ϕ , ψ ) + 1 2 ( Δ ) s 4 ( ϕ φ ) L 2 2 .
Proof. 
Adding the identities
ϕ n + 1 ϕ n , χ ( ϕ n + 1 , ϕ n ) + ( Δ ) s ϕ n + 1 + ϕ n 2 L 2 = E s ( ϕ n + 1 ) E s ( ϕ n ) + ( Δ ) s 4 ϕ n + 1 L 2 2 ( Δ ) s 4 ϕ n L 2 2
and
ϕ n + 1 ϕ n , ( Δ ) s 2 ( 3 ϕ n ϕ n 1 ) L 2 + ( Δ ) s 4 ϕ n + 1 L 2 2 ( Δ ) s 4 ϕ n L 2 2 = 1 2 ( Δ ) s 4 ( ϕ n + 1 ϕ n ) L 2 2 1 2 ( Δ ) s 4 ( ϕ n ϕ n 1 ) L 2 2 + 1 2 ( Δ ) s 4 ( ϕ n + 1 2 ϕ n + ϕ n 1 ) L 2 2 ,
we have from Equations (12) and (15)
( Δ t ψ n + 1 2 , μ s , n + 1 2 ) L 2 = ( ϕ n + 1 ϕ n , μ s , n + 1 2 ) L 2 = E s ( ϕ n + 1 ) + 1 2 ( Δ ) s 4 ( ϕ n + 1 ϕ n ) L 2 2 E s ( ϕ n ) 1 2 ( Δ ) s 4 ( ϕ n ϕ n 1 ) L 2 2 + 1 2 ( Δ ) s 4 ( ϕ n + 1 2 ϕ n + ϕ n 1 ) L 2 2 ,
where ψ n + 1 2 = ψ n + 1 + ψ n 2 . In addition, we obtain from Equation (14) and Theorem 1
( ψ n + 1 2 , Δ t μ s , n + 1 2 ) L 2 = ψ n + 1 2 , ψ n + 1 ψ n Δ t α n + 1 2 + β Δ t ψ n + 1 2 L 2 = 1 2 ψ n + 1 L 2 2 1 2 ψ n L 2 2 Δ t α n + 1 2 ( ψ n + 1 2 , 1 ) L 2 + β Δ t ψ n + 1 2 L 2 2 = 1 2 ψ n + 1 L 2 2 1 2 ψ n L 2 2 + β Δ t ψ n + 1 2 L 2 2 .
Then, we have
F ˜ ( ϕ n + 1 , ψ n + 1 , ϕ n ) F ˜ ( ϕ n , ψ n , ϕ n 1 ) = β Δ t ψ n + 1 2 L 2 2 1 2 ( Δ ) s 4 ( ϕ n + 1 2 ϕ n + ϕ n 1 ) L 2 2 0 .

Numerical Implementation

To handle ( ϕ n + 1 ) 3 and ( ϕ n + 1 ) 2 in μ s , n + 1 2 and α n + 1 2 in Equation (14), a Newton-type linearization is used [9,12,17]:
( ϕ n , m + 1 ) 3 ( ϕ n , m ) 3 + 3 ( ϕ n , m ) 2 ( ϕ n , m + 1 ϕ n , m ) , ( ϕ n , m + 1 ) 2 ( ϕ n , m ) 2 + 2 ϕ n , m ( ϕ n , m + 1 ϕ n , m )
for m = 0 , 1 , . We then develop a Newton-type fixed point iteration method for Equation (14) starting with ϕ n , 0 = ϕ n as
I + Δ t 2 2 + β Δ t A n , m + 1 2 ( Δ ) s 1 | Ω | A n , m , · L 2 ϕ n , m + 1 = ϕ n Δ t 2 2 + β Δ t B n , m + 1 2 ( Δ ) s ϕ n ( Δ ) s 2 ( 3 ϕ n ϕ n 1 ) 1 | Ω | B n , m , 1 L 2 + 2 Δ t 2 + β Δ t ψ n ,
where A n , m = 3 ( ϕ n , m ) 2 + 2 ϕ n , m ϕ n + ( ϕ n ) 2 4 + 1 ϵ 2 , A n , m , · L 2 ϕ = A n , m , ϕ L 2 ,
B n , m = 2 ( ϕ n , m ) 3 ( ϕ n , m ) 2 ϕ n + ( ϕ n ) 3 4 + 1 ϵ 2 ϕ n , and we set
ϕ n + 1 = ϕ n , m + 1
if a relative l 2 -norm of the consecutive error ϕ n , m + 1 ϕ n , m 2 ϕ n , m 2 is less than a tolerance t o l .
The system (17) is solved using the biconjugate gradient method with the preconditioner P = I + Δ t 2 2 ( 2 + β Δ t ) ( Δ ) s . The stopping criterion for the biconjugate gradient iteration is that the relative residual norm is less than t o l . Then, we employ the Fourier spectral method [15,18,19,20,21,22] for the spatial discretization, which gives a full diagonal representation of the fractional operator.

4. Numerical Experiments

4.1. Accuracy Test

We estimate the accuracy of the proposed method with initial conditions [17]
ϕ ( x , y , 0 ) = 0.02 cos 2 π ( x + 10 ) 32 sin 2 π ( y + 3 ) 32 0.01 sin 2 π x 8 sin 2 π ( y 6 ) 8 0.02 cos π ( x 12 ) 16 sin π ( y 1 ) 16 + 0.07 , ψ ( x , y , 0 ) = 0
on Ω = [ 0 , 32 ] × [ 0 , 32 ] . We set ϵ = 0.025 , β = 1 , Δ x = Δ y = 1 3 , and t o l = 10 8 Δ t , and we compute ϕ ( x , y , t ) for 0 < t 10 . For s = 2 , 1.5 , 1.05 , Figure 1a,b show the evolution of F ( t ) for the reference solution with Δ t = 2 8 and the relative l 2 -errors of ϕ ( x , y , 2 ) for Δ t = 2 6 , 2 5 , , 1 , respectively. Here, the errors are computed by comparison with the reference solution. Figure 1c shows the evolution of Ω ( ϕ ( x , y , t ) ϕ ( x , y , 0 ) ) d x d y for s = 2 and various time steps. It is observed that the method is second-order accurate in time and conserves the mass.

4.2. Energy Stability Test

To demonstrate the energy stability of the proposed method, we take the initial conditions (18) on Ω = [ 0 , 32 ] × [ 0 , 32 ] , and use ϵ = 0.25 , β = 1 , Δ x = Δ y = 1 3 , and t o l = 10 8 . For s = 2 , 1.5 , 1.05 , Figure 2a–c show the evolution of F ˜ ( t ) with different time steps, respectively. The energy curves are all nonincreasing in time, which verifies that the proposed method is unconditionally energy stable (Theorem 3). Moreover, the solution with s = 2 and a sufficiently small time step Δ t = 2 4 is stable at around t = 224 (see the magenta line with dots in Figure 2a); however, the case of s = 1.05 is not stable even at t = 512 (see the magenta line with dots in Figure 2c). These results show that the lifetime of the unstable interface is prolonged as the fractional order s decreases, i.e., a longer time computation is needed to reach a steady state. Figure 3 shows the evolution of ϕ ( x , y , t ) with s = 2 and Δ t = 2 4 .

4.3. Phase Diagram in 2D

In 2D, the phase diagram contains constant ( ϕ c ), striped ( ϕ s ), and triangular ( ϕ t ) states depending on the values of ϕ ¯ (the average density) and ϵ , and these states can be approximated by ϕ c = ϕ ¯ , ϕ s = ϕ ¯ + A s sin ( x ) , and ϕ t = ϕ ¯ + A t cos 3 x 2 cos y 2 1 2 cos ( y ) , where A s and A t are the characteristic amplitudes [23]. To verify that the proposed method does lead to the expected states, we take the initial conditions as ϕ ( x , y , 0 ) = ϕ ¯ + rand , ψ ( x , y , 0 ) = 0 on Ω = [ 0 , 32 ] × [ 0 , 32 ] , where rand is a random number between 0.1 and 0.1 at the grid points. We use s = 2 , Δ x = Δ y = 1 3 , Δ t = 1 4 β , and t o l = 10 9 β , and we compute ϕ ( x , y , t ) for 0 < t t f = 3200 β . To estimate the phase diagram numerically, we calculate the following indicator function:
Λ ( t ) = Ω | ϕ ( x , t ) ϕ ¯ | d x if A = 0 , 1 A Ω | ϕ ( x , t ) ϕ ¯ | d x if A 0 ,
where A = 1 2 max x Ω ϕ ( x , t f ) min x Ω ϕ ( x , t f ) . In the numerical simulations, we set A = 0 if A is less than 10 8 . Figure 4a–c show Λ ( t ) at t = t f with various ϕ ¯ = 0.05 , 0.1 , , 0.5 and ϵ = 0.05 , 0.1 , , 0.5 for β = 0.1 , 1 , 10 , respectively. The results in Figure 4 are consistent with the phase diagram in [1,23]. Sample time evolutions of Λ ( t ) with ϕ ¯ = 0.05 , 0.1 , , 0.5 and ϵ = 0.3 for β = 0.1 , 1 , 10 are shown in Figure 5a–c, respectively. Figure 6 shows ϕ ( x , y , t f ) with various ϕ ¯ and ϵ for β = 1 .

5. Conclusions

In this paper, we introduced the new L 2 -fMPFC equation that conserves the mass precisely without influencing the energy dissipation property. In addition, we proposed the second-order convex splitting method for the L 2 -fMPFC equation. Our mathematical analyses and numerical examples showed that the proposed method inherited the mass conservation and energy stability. The results for pattern formation in 2D that is a standard test in the classical MPFC equation demonstrated the feasibility of the proposed framework.
In this paper, the convergence rate of the proposed method was only demonstrated numerically. We have a plan to carry out a rigorous error analysis for the proposed method in future work.

Funding

The present research was funded by the Research Grant of Kwangwoon University in 2022 and by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2019R1C1C1011112, 2022R1A2C1011708).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks the reviewers for the constructive and helpful comments on the revision of this article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Stefanovic, P.; Haataja, M.; Provatas, N. Phase-field crystals with elastic interactions. Phys. Rev. Lett. 2006, 96, 225504. [Google Scholar] [CrossRef]
  2. Swift, J.; Hohenberg, P.C. Hydrodynamic fluctuations at the convective instability. Phys. Rev. A 1977, 15, 319–328. [Google Scholar] [CrossRef]
  3. Wang, C.; Wise, S.M. An energy stable and convergent finite-difference scheme for the modified phase field crystal equation. SIAM J. Numer. Anal. 2011, 49, 945–969. [Google Scholar] [CrossRef]
  4. Akagi, G.; Schimperna, G.; Segatti, A. Fractional Cahn–Hilliard, Allen–Cahn and porous medium equations. J. Differ. Equ. 2016, 261, 2935–2985. [Google Scholar] [CrossRef]
  5. Ainsworth, M.; Mao, Z. Analysis and approximation of a fractional Cahn–Hilliard equation. SIAM J. Numer. Anal. 2017, 55, 1689–1718. [Google Scholar] [CrossRef]
  6. Weng, Z.; Zhai, S.; Feng, X. A Fourier spectral method for fractional-in-space Cahn–Hilliard equation. Appl. Math. Model. 2017, 42, 462–477. [Google Scholar] [CrossRef]
  7. Bu, L.; Mei, L.; Hou, Y. Stable second-order schemes for the space-fractional Cahn–Hilliard and Allen–Cahn equations. Comput. Math. Appl. 2019, 78, 3485–3500. [Google Scholar] [CrossRef]
  8. Wang, F.; Chen, H.; Wang, H. Finite element simulation and efficient algorithm for fractional Cahn–Hilliard equation. J. Comput. Appl. Math. 2019, 356, 248–266. [Google Scholar] [CrossRef]
  9. Baskaran, A.; Hu, Z.; Lowengrub, J.S.; Wang, C.; Wise, S.M.; Zhou, P. Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation. J. Comput. Phys. 2013, 250, 270–292. [Google Scholar] [CrossRef]
  10. Galenko, P.K.; Gomez, H.; Kropotin, N.V.; Elder, K.R. Unconditionally stable method and numerical solution of the hyperbolic phase-field crystal equation. Phys. Rev. E 2013, 88, 013310. [Google Scholar] [CrossRef] [Green Version]
  11. Dehghan, M.; Mohammadi, V. The numerical simulation of the phase field crystal (PFC) and modified phase field crystal (MPFC) models via global and local meshless methods. Comput. Methods Appl. Mech. Engrg. 2016, 298, 453–484. [Google Scholar] [CrossRef]
  12. Lee, H.G.; Shin, J.; Lee, J.-Y. First- and second-order energy stable methods for the modified phase field crystal equation. Comput. Methods Appl. Mech. Engrg. 2017, 321, 1–17. [Google Scholar] [CrossRef]
  13. Li, Q.; Mei, L.; Yang, X.; Li, Y. Efficient numerical schemes with unconditional energy stabilities for the modified phase field crystal equation. Adv. Comput. Math. 2019, 45, 1551–1580. [Google Scholar] [CrossRef]
  14. Li, X.; Shen, J. Efficient linear and unconditionally energy stable schemes for the modified phase field crystal equation. Sci. China Math. 2021, 1–18. [Google Scholar] [CrossRef]
  15. Shin, J.; Lee, H.G.; Lee, J.-Y. Energy quadratization Runge–Kutta method for the modified phase field crystal equation. Model. Simul. Mater. Sci. Eng. 2022, 30, 024004. [Google Scholar] [CrossRef]
  16. Wise, S.M.; Wang, C.; Lowengrub, J.S. An energy-stable and convergent finite-difference scheme for the phase field crystal equation. SIAM J. Numer. Anal. 2009, 47, 2269–2288. [Google Scholar] [CrossRef]
  17. Hu, Z.; Wise, S.M.; Wang, C.; Lowengrub, J.S. Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation. J. Comput. Phys. 2009, 228, 5323–5339. [Google Scholar] [CrossRef]
  18. Kim, J.; Lee, H.G. Unconditionally energy stable second-order numerical scheme for the Allen–Cahn equation with a high-order polynomial free energy. Adv. Differ. Equ. 2021, 2021, 416. [Google Scholar] [CrossRef]
  19. Shin, J.; Lee, H.G. A linear, high-order, and unconditionally energy stable scheme for the epitaxial thin film growth model without slope selection. Appl. Numer. Math. 2021, 163, 30–42. [Google Scholar] [CrossRef]
  20. Lee, H.G. A non-iterative and unconditionally energy stable method the Swift–Hohenberg equation with quadratic–cubic nonlinearity. Appl. Math. Lett. 2022, 123, 107579. [Google Scholar] [CrossRef]
  21. Lee, H.G.; Shin, J.; Lee, J.-Y. A high-order and unconditionally energy stable scheme for the conservative Allen-Cahn equation with a nonlocal Lagrange multiplier. J. Sci. Comput. 2022, 90, 51. [Google Scholar] [CrossRef]
  22. Li, X.; Han, C.; Wang, Y. Novel patterns in fractional-in-space nonlinear coupled FitzHugh–Nagumo models with Riesz fractional derivative. Fractal Fract. 2022, 6, 136. [Google Scholar] [CrossRef]
  23. Elder, K.R.; Grant, M. Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. Phys. Rev. E 2004, 70, 051605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) Evolution of F ( t ) for the reference solution with ϵ = 0.025 , β = 1 , Δ x = Δ y = 1 3 , and Δ t = 2 8 . (b) Relative l 2 -errors of ϕ ( x , y , 2 ) for Δ t = 2 6 , 2 5 , , 1 . (c) Evolution of Ω ( ϕ ( x , y , t ) ϕ ( x , y , 0 ) ) d x d y for various time steps.
Figure 1. (a) Evolution of F ( t ) for the reference solution with ϵ = 0.025 , β = 1 , Δ x = Δ y = 1 3 , and Δ t = 2 8 . (b) Relative l 2 -errors of ϕ ( x , y , 2 ) for Δ t = 2 6 , 2 5 , , 1 . (c) Evolution of Ω ( ϕ ( x , y , t ) ϕ ( x , y , 0 ) ) d x d y for various time steps.
Fractalfract 06 00472 g001
Figure 2. Evolution of F ˜ ( t ) with different time steps.
Figure 2. Evolution of F ˜ ( t ) with different time steps.
Fractalfract 06 00472 g002
Figure 3. Evolution of ϕ ( x , y , t ) with ϵ = 0.25 , β = 1 , s = 2 , Δ x = Δ y = 1 3 , and Δ t = 2 4 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.6283 , 0.0688 , and 0.4906 , respectively.
Figure 3. Evolution of ϕ ( x , y , t ) with ϵ = 0.25 , β = 1 , s = 2 , Δ x = Δ y = 1 3 , and Δ t = 2 4 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.6283 , 0.0688 , and 0.4906 , respectively.
Fractalfract 06 00472 g003
Figure 4. Values of Λ ( t ) at t = t f with various ϕ ¯ and ϵ for different β .
Figure 4. Values of Λ ( t ) at t = t f with various ϕ ¯ and ϵ for different β .
Fractalfract 06 00472 g004
Figure 5. Evolution of Λ ( t ) with ϕ ¯ = 0.05 , 0.1 , , 0.5 and ϵ = 0.3 for different β .
Figure 5. Evolution of Λ ( t ) with ϕ ¯ = 0.05 , 0.1 , , 0.5 and ϵ = 0.3 for different β .
Fractalfract 06 00472 g005
Figure 6. ϕ ( x , y , t f ) for various ϕ ¯ and ϵ for β = 1 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.5 , 0, and 0.5 , respectively.
Figure 6. ϕ ( x , y , t f ) for various ϕ ¯ and ϵ for β = 1 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.5 , 0, and 0.5 , respectively.
Fractalfract 06 00472 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, H.G. A New L2-Gradient Flow-Based Fractional-in-Space Modified Phase-Field Crystal Equation and Its Mass Conservative and Energy Stable Method. Fractal Fract. 2022, 6, 472. https://doi.org/10.3390/fractalfract6090472

AMA Style

Lee HG. A New L2-Gradient Flow-Based Fractional-in-Space Modified Phase-Field Crystal Equation and Its Mass Conservative and Energy Stable Method. Fractal and Fractional. 2022; 6(9):472. https://doi.org/10.3390/fractalfract6090472

Chicago/Turabian Style

Lee, Hyun Geun. 2022. "A New L2-Gradient Flow-Based Fractional-in-Space Modified Phase-Field Crystal Equation and Its Mass Conservative and Energy Stable Method" Fractal and Fractional 6, no. 9: 472. https://doi.org/10.3390/fractalfract6090472

APA Style

Lee, H. G. (2022). A New L2-Gradient Flow-Based Fractional-in-Space Modified Phase-Field Crystal Equation and Its Mass Conservative and Energy Stable Method. Fractal and Fractional, 6(9), 472. https://doi.org/10.3390/fractalfract6090472

Article Metrics

Back to TopTop