Next Article in Journal
Strong Tolerance and Strong Universality of Interval Eigenvectors in a Max-Łukasiewicz Algebra
Next Article in Special Issue
Numerical Investigation of Freely Falling Objects Using Direct-Forcing Immersed Boundary Method
Previous Article in Journal
Progressive Iterative Approximation with Preconditioners
Previous Article in Special Issue
Fourier-Spectral Method for the Phase-Field Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Efficient and Accurate Method for the Conservative Swift–Hohenberg Equation and Its Numerical Implementation

Department of Mathematics, Kwangwoon University, Seoul 01897, Korea
Mathematics 2020, 8(9), 1502; https://doi.org/10.3390/math8091502
Submission received: 21 August 2020 / Revised: 30 August 2020 / Accepted: 2 September 2020 / Published: 4 September 2020
(This article belongs to the Special Issue Open Source Codes for Numerical Analysis)

Abstract

:
The conservative Swift–Hohenberg equation was introduced to reformulate the phase-field crystal model. A challenge in solving the conservative Swift–Hohenberg equation numerically is how to treat the nonlinear term to preserve mass conservation without compromising efficiency and accuracy. To resolve this problem, we present a linear, high-order, and mass conservative method by placing the linear and nonlinear terms in the implicit and explicit parts, respectively, and employing the implicit-explicit Runge–Kutta method. We show analytically that the method inherits the mass conservation. Numerical experiments are presented demonstrating the efficiency and accuracy of the proposed method. In particular, long time simulation for pattern formation in 2D is carried out, where the phase diagram can be observed clearly. The MATLAB code for numerical implementation of the proposed method is provided in Appendix.

1. Introduction

The phase-field crystal (PFC) model describes the microstructure of two-phase systems on atomic length and diffusive time scales and has been used to study grain growth, dendritic and eutectic solidification, and epitaxial growth [1,2]. The PFC model is the H 1 -gradient flow for the Swift–Hohenberg (SH) energy functional [3]:
E ( ϕ ) : = Ω Φ ( ϕ ) | ϕ | 2 + 1 2 ( Δ ϕ ) 2 d x ,
where Ω is a domain in R d ( d = 1 , 2 , 3 ) , ϕ is the density field, Φ ( ϕ ) = 1 4 ϕ 4 g 3 ϕ 3 + 1 ϵ 2 ϕ 2 , and g 0 and ϵ > 0 are positive constants with physical significance.
Recently, conservative SH equations were introduced to reformulate the PFC model [4,5]. In [4], Zhang and Yang derived the following equation:
ϕ t = Φ ( ϕ ) + 2 Δ ϕ + Δ 2 ϕ + β ˜ ( t ) ,
where β ˜ ( t ) is a nonlocal Lagrange multiplier and β ˜ ( t ) = 1 | Ω | Ω Φ ( ϕ ( x , t ) ) d x , and developed a second-order energy stable scheme by combining the invariant energy quadratization idea with the stabilization technique. However, the scheme involves solving a linear system with complicated variable coefficients. In [5], Lee introduced the following equation:
ϕ t = Φ ( ϕ ) + 2 Δ ϕ + Δ 2 ϕ + Φ ( ϕ ) r β ( t ) ,
where Φ ( ϕ ) r β ( t ) is a nonlocal and local Lagrange multiplier and β ( t ) = Ω Φ ( ϕ ( x , t ) ) d x Ω Φ ( ϕ ( x , t ) ) r d x , and proposed mass conservative first- and second-order operator splitting methods. However, the methods lead to the necessity nonlinear equations to be solved at each time step thus require an iterative solver for solving the nonlinear equations.
Therefore, the aim of this paper is to present an efficient and accurate method that preserves mass conservation for solving the conservative SH Equation (3). We place the linear and nonlinear terms in the implicit and explicit parts, respectively, where an extra linear stabilizing term is added to improve the stability while preserving the simplicity. And we employ the implicit-explicit Runge–Kutta (RK) method [6]. As a result, our method is linear, high-order accurate in time, and mass conservative. We show analytically that the method inherits the mass conservation. In addition, the Fourier spectral method [5,7,8,9,10] is used for the spatial discretization. The MATLAB code for numerical implementation of the method in 2D is provided in Appendix A.
This paper is organized as follows. In Section 2, we construct the linear, high-order, and mass conservative method and show analytically that the method inherits the mass conservation. Numerical examples showing the efficiency and accuracy of the proposed method are presented in Section 3. Finally, conclusions are drawn in Section 4. In Appendix A, we provide the MATLAB code for numerical implementation of the proposed method in 2D.

2. Linear, High-Order, and Mass Conservative Method

For simplicity and clarity of exposition, we consider Equation (3) in one-dimensional space Ω = [ 0 , L ] with a periodic boundary condition:
ϕ ( x , t ) t = Φ ( ϕ ( x , t ) ) + 2 2 ϕ ( x , t ) x 2 + 4 ϕ ( x , t ) x 4 + Φ ( ϕ ( x , t ) ) r β ( t ) ,
where β ( t ) = Ω Φ ( ϕ ( x , t ) ) d x Ω Φ ( ϕ ( x , t ) ) r d x . Two- and three-dimensional cases are defined analogously. Let M be a positive integer, Δ x = L M be the space step size, and Δ t be the time step size. Let ϕ m n be an approximation of ϕ ( x m , t n ) , where x m = m Δ x for m = 0 , 1 , , M 1 and t n = n Δ t . The discrete Fourier transform and its inverse transform are
ϕ ^ k = m = 0 M 1 ϕ m e i x m ξ k
and
ϕ m = 1 M k = 0 M 1 ϕ ^ k e i x m ξ k ,
where ξ k = 2 π k L .
To develop a linear, high-order (up to third-order), and mass conservative method for solving Equation (4), we treat s ϕ ( x , t ) + 2 2 ϕ ( x , t ) x 2 + 4 ϕ ( x , t ) x 4 implicitly and Φ ( ϕ ( x , t ) ) s ϕ ( x , t ) + Φ ( ϕ ( x , t ) ) r β ( t ) explicitly, where s is a non-negative number, and employ the implicit-explicit RK method. First- (S1), second- (S2), and third- (S3) order methods are as follows:
S 1 : ϕ m n + 1 = ϕ m n + Δ t p ( ϕ m n + 1 ) + q ( ϕ m n ) ,
S 2 : ϕ m ( 1 ) = ϕ m n + Δ t γ p ( ϕ m ( 1 ) ) + γ q ( ϕ m n ) ,
ϕ m n + 1 = ϕ m n + Δ t γ p ( ϕ m n + 1 ) + ( 1 γ ) p ( ϕ m ( 1 ) ) + ( 1 δ ) q ( ϕ m ( 1 ) ) + δ q ( ϕ m n ) ,
S 3 : ϕ m ( 1 ) = ϕ m n + Δ t 1 2 p ( ϕ m ( 1 ) ) + 1 2 q ( ϕ m n ) ,
ϕ m ( 2 ) = ϕ m n + Δ t 1 2 p ( ϕ m ( 2 ) ) + 1 6 p ( ϕ m ( 1 ) ) + 1 18 q ( ϕ m ( 1 ) ) + 11 18 q ( ϕ m n ) ,
ϕ m ( 3 ) = ϕ m n + Δ t ( 1 2 p ( ϕ m ( 3 ) ) + 1 2 p ( ϕ m ( 2 ) ) 1 2 p ( ϕ m ( 1 ) )   + 1 2 q ( ϕ m ( 2 ) ) 5 6 q ( ϕ m ( 1 ) ) + 5 6 q ( ϕ m n ) ) ,
ϕ m n + 1 = ϕ m n + Δ t ( 1 2 p ( ϕ m n + 1 ) + 1 2 p ( ϕ m ( 3 ) ) 3 2 p ( ϕ m ( 2 ) ) + 3 2 p ( ϕ m ( 1 ) )   7 4 q ( ϕ m ( 3 ) ) + 3 4 q ( ϕ m ( 2 ) ) + 7 4 q ( ϕ m ( 1 ) ) + 1 4 q ( ϕ m n ) ) ,
where p ( ϕ m ( · ) ) = s ϕ m ( · ) + 2 2 ϕ m ( · ) x 2 + 4 ϕ m ( · ) x 4 , q ( ϕ m ( · ) ) = Φ ( ϕ m ( · ) ) s ϕ m ( · ) + Φ ( ϕ m ( · ) ) r β ( · ) , β ( · ) = m = 0 M 1 Φ ( ϕ m ( · ) ) m = 0 M 1 Φ ( ϕ m ( · ) ) r , γ = 2 2 2 , and δ = 1 1 2 γ .
For the method S1, Equation (7) can be transformed into the discrete Fourier space using (6):
ϕ ^ k n + 1 = ϕ ^ k n + Δ t q ^ k n 1 + Δ t s 2 ξ k 2 + ξ k 4 ,
where q ^ k ( · ) = F [ q ( ϕ m ( · ) ) ] and F denotes the discrete Fourier transform. After updating ϕ ^ k n + 1 with ϕ ^ k n , we recover ϕ m n + 1 from ϕ ^ k n + 1 using (6). To satisfy the mass conservation property, we should have m = 0 M 1 ϕ m n + 1 = m = 0 M 1 ϕ m n for n = 0 , 1 , . From Equation (14), we get
m = 0 M 1 ϕ m n + 1 = ϕ ^ 0 n + 1 = ϕ ^ 0 n + Δ t q ^ 0 n 1 + s Δ t = ϕ ^ 0 n + s Δ t ϕ ^ 0 n 1 + s Δ t = ϕ ^ 0 n = m = 0 M 1 ϕ m n
since
q ^ 0 ( · ) = m = 0 M 1 q ( ϕ m ( · ) ) = m = 0 M 1 ( Φ ( ϕ m ( · ) ) s ϕ m ( · ) ) + β ( · ) m = 0 M 1 Φ ( ϕ m ( · ) ) r = m = 0 M 1 ( Φ ( ϕ m ( · ) ) s ϕ m ( · ) ) + m = 0 M 1 Φ ( ϕ m ( · ) ) m = 0 M 1 Φ ( ϕ m ( · ) ) r m = 0 M 1 Φ ( ϕ m ( · ) ) r = s m = 0 M 1 ϕ m ( · ) = s ϕ ^ 0 ( · ) .
Thus, the method S1 inherits the mass conservation. Next, for the method S2, we have
ϕ ^ 0 ( 1 ) = ϕ ^ 0 n + γ Δ t q ^ 0 n 1 + γ s Δ t = ϕ ^ 0 n + γ s Δ t ϕ ^ 0 n 1 + γ s Δ t = ϕ ^ 0 n
from Equation (8) and
m = 0 M 1 ϕ m n + 1 = ϕ ^ 0 n + 1 = ϕ ^ 0 n + Δ t ( 1 γ ) s ϕ ^ 0 ( 1 ) + ( 1 δ ) q ^ 0 ( 1 ) + δ q ^ 0 n 1 + γ s Δ t = ϕ ^ 0 n + Δ t ( 1 γ ) s ϕ ^ 0 n + ( 1 δ ) s ϕ ^ 0 n + δ s ϕ ^ 0 n 1 + γ s Δ t = ϕ ^ 0 n = m = 0 M 1 ϕ m n
from Equation (9). For the method S3, we have ϕ ^ 0 ( 3 ) = ϕ ^ 0 ( 2 ) = ϕ ^ 0 ( 1 ) = ϕ ^ 0 n from Equations (10)–(12) and
m = 0 M 1 ϕ m n + 1 = ϕ ^ 0 n + 1 = ϕ ^ 0 n + Δ t s 2 ϕ ^ 0 ( 3 ) + 3 s 2 ϕ ^ 0 ( 2 ) 3 s 2 ϕ ^ 0 ( 1 ) 7 4 q ^ 0 ( 3 ) + 3 4 q ^ 0 ( 2 ) + 7 4 q ^ 0 ( 1 ) + 1 4 q ^ 0 n 1 + s Δ t 2 = ϕ ^ 0 n = m = 0 M 1 ϕ m n
from Equation (13). Thus, the methods S2 and S3 also inherit the mass conservation.

3. Numerical Experiments

3.1. Convergence Test

We demonstrate the convergence of the proposed methods with the initial condition [11,12]
ϕ ( x , 0 ) = 0.07 0.02 cos 2 π ( x 12 ) 32 + 0.02 cos 2 π ( x + 10 ) 32 0.01 sin 2 4 π x 32
on Ω = [ 0 , 32 ] . We set ϵ = 0.25 , g = 0 , r = 0 , and s = 0 , and compute ϕ ( x , t ) for 0 < t 96 . The grid size is fixed to Δ x = 1 3 which provides enough spatial accuracy. To estimate the convergence rate with respect to Δ t , simulations are performed by varying Δ t = 2 10 , 2 9 , , 2 4 . We take the quadruply over-resolved numerical solution using the method S3 as the reference solution. Figure 1a shows the relative l 2 -errors of ϕ ( x , 40 ) for various time steps. Here, the errors are computed by comparison with the reference solution. In addition, Figure 1b–d show the evolution of Ω ( ϕ ( x , t ) ϕ ( x , 0 ) ) d x using the methods S1S3, respectively. Here, Ω ( ϕ ( x , t ) ϕ ( x , 0 ) ) d x is approximated by m = 0 M 1 ( ϕ m n ϕ m 0 ) Δ x . It is observed that the methods give desired order of accuracy in time and conserve the total mass.

3.2. Efficiency of the Proposed Method

To show the efficiency of the proposed method, we take the initial condition (20) and parameter values used to create Figure 1. Figure 2 presents the CPU time (in seconds, averaged over 10 trials performed on Intel Core i5-7500 CPU at 3.40 GHz with 8 GB RAM) consumed using the methods S1S3 for various time steps. The results suggest that the CPU time is almost linear with respect to the number of steps and the methods S2 and S3 are about two and four times more expensive than the method S1, respectively.

3.3. Phase Diagram in 2D

In 2D, the phase diagram contains striped, hexagonal, and constant states depending on the values of ϕ ¯ and ϵ [1] (see Figure 3a). To verify that the proposed method does lead to the expected states, we take an initial condition as ϕ ( x , y , 0 ) = ϕ ¯ + rand on Ω = [ 0 , 32 ] × [ 0 , 32 ] . Here, rand is a random number between 0.1 and 0.1 at the grid points, and we use g = 0 , r = 0 , Δ x = Δ y = 1 3 , s = 2 , and the method S3. For saving computational time, we choose different time steps as the solution evolves from random noisy stage to smooth stage: Δ t = 1 4 for 0 < t 128 and Δ t = 4 for 128 < t 2176 . To estimate the phase diagram numerically, we calculate the indicator function defined similarly in [13]:
Λ ( t ) = Ω | ϕ ( x , t ) ϕ ¯ | d x if Ω | ( ϕ ( x , t ) ϕ ¯ ) | d x 0 Ω | ϕ ( x , t ) ϕ ¯ | d x Ω | ( ϕ ( x , t ) ϕ ¯ ) | d x otherwise .
Here, we set Ω | ( ϕ ( x , t ) ϕ ¯ ) | d x 0 if Ω | ( ϕ ( x , t ) ϕ ¯ ) | d x is less than 10 10 . Figure 3b shows Λ ( t ) at t = 2176 with various ϕ ¯ = 0.02 , 0.04 , , 0.3 and ϵ = 0.02 , 0.04 , , 0.3 . Results in Figure 3b are consistent with the phase diagram in Figure 3a. Sample time evolutions of ϕ ( x , y , t ) with ( ϕ ¯ , ϵ ) = ( 0.02 , 0.1 ) , ( 0.14 , 0.1 ) , and ( 0.26 , 0.1 ) are shown in Figure 4, Figure 5 and Figure 6, respectively. Figure 7 shows evolutions of E ( t ) and Λ ( t ) with ( ϕ ¯ , ϵ ) used in Figure 4, Figure 5 and Figure 6. We remark that the solution ϕ ( x , y , t ) with ( ϕ ¯ , ϵ ) = ( 0.26 , 0.1 ) has a small but measurable perturbation from a constant state until t = 128 (see Figure 6). In this case, Ω | ( ϕ ( x , t ) ϕ ¯ ) | d x 0 and Λ ( t ) is calculated using Ω | ϕ ( x , t ) ϕ ¯ | d x / Ω | ( ϕ ( x , t ) ϕ ¯ ) | d x . Afterward, the solution evolves to a constant state, i.e., Ω | ϕ ( x , t ) ϕ ¯ | d x 0 and Ω | ( ϕ ( x , t ) ϕ ¯ ) | d x 0 . Thus, Λ ( t ) is calculated using Ω | ϕ ( x , t ) ϕ ¯ | d x and approximately zero.

3.4. Comparison with Other Method

To compare the proposed method with other method, we solve the conservative SH Equation (3) using the proposed method S2 and the second-order operator splitting (OS2) method in [5] with the initial condition and parameter values used to create Figure 4 except for Δ t . Figure 8 and Figure 9 show evolutions of ϕ ( x , y , t ) using the method OS2 with Δ t = 1 4 and 2, respectively. The method OS2 with a smaller time step Δ t = 1 4 leads to the expected striped state, whereas a constant state is observed for Δ t = 2 . Figure 10 and Figure 11 show evolutions of ϕ ( x , y , t ) using the method S2 with Δ t = 1 4 and 2, respectively. The method S2 gives the striped state even for a large time step. Evolutions of E ( t ) for Figure 8, Figure 9, Figure 10 and Figure 11 are shown in Figure 12.

4. Conclusions

In this paper, we developed linear, first-, second-, and third-order, and mass conservative methods for the conservative SH equation by placing the linear and nonlinear terms in the implicit and explicit parts, respectively, and employing the implicit-explicit RK method. We confirmed that the proposed methods give desired order of accuracy in time, inherit the mass conservation, and are efficient (the CPU time was almost linear with respect to the number of steps and of stages). And we performed long time simulation for pattern formation in 2D, where the phase diagram can be observed clearly.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2019R1C1C1011112).

Acknowledgments

The corresponding 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.

Appendix A. Matlab Code

MATLAB code for numerical implementation of the methods S1S3 (7)–(13) in 2D. Mathematics 08 01502 i001a Mathematics 08 01502 i001b

References

  1. Elder, K.R.; Katakowski, M.; Haataja, M.; Grant, M. Modeling elasticity in crystal growth. Phys. Rev. Lett. 2002, 88, 245701. [Google Scholar] [CrossRef] [Green Version]
  2. 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]
  3. Swift, J.; Hohenberg, P.C. Hydrodynamic fluctuations at the convective instability. Phys. Rev. A 1977, 15, 319–328. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, J.; Yang, X. Numerical approximations for a new L2-gradient flow based Phase field crystal model with precise nonlocal mass conservation. Comput. Phys. Commun. 2019, 243, 51–67. [Google Scholar] [CrossRef]
  5. Lee, H.G. A new conservative Swift-Hohenberg equation and its mass conservative method. J. Comput. Appl. Math. 2020, 375, 112815. [Google Scholar] [CrossRef]
  6. Ascher, U.M.; Ruuth, S.J.; Spiteri, R.J. Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations. Appl. Numer. Math. 1997, 25, 151–167. [Google Scholar] [CrossRef]
  7. Lee, H.G. A semi-analytical Fourier spectral method for the Swift-Hohenberg equation. Comput. Math. Appl. 2017, 74, 1885–1896. [Google Scholar] [CrossRef]
  8. Lee, H.G. An energy stable method for the Swift-Hohenberg equation with quadratic-cubic nonlinearity. Comput. Methods Appl. Mech. Eng. 2019, 343, 40–51. [Google Scholar] [CrossRef]
  9. Chen, X.; Song, M.; Song, S. A fourth order energy dissipative scheme for a traffic flow model. Mathematics 2020, 8, 1238. [Google Scholar] [CrossRef]
  10. Yoon, S.; Jeong, D.; Lee, C.; Kim, H.; Kim, S.; Lee, H.G.; Kim, J. Fourier-spectral method for the phase-field equations. Mathematics 2020, 8, 1385. [Google Scholar] [CrossRef]
  11. 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] [Green Version]
  12. 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] [Green Version]
  13. Shin, J.; Lee, H.G.; Lee, J.-Y. Long-time simulation of the phase-field crystal equation using high-order energy-stable CSRK methods. Comput. Methods Appl. Mech. Eng. 2020, 364, 112981. [Google Scholar] [CrossRef]
Figure 1. (a) Relative l 2 -errors of ϕ ( x , 40 ) for various time steps with ϵ = 0.25 and Δ x = 1 3 . (b)–(d) Evolution of Ω ( ϕ ( x , t ) ϕ ( x , 0 ) ) d x for various time steps using the methods S1S3.
Figure 1. (a) Relative l 2 -errors of ϕ ( x , 40 ) for various time steps with ϵ = 0.25 and Δ x = 1 3 . (b)–(d) Evolution of Ω ( ϕ ( x , t ) ϕ ( x , 0 ) ) d x for various time steps using the methods S1S3.
Mathematics 08 01502 g001
Figure 2. CPU time versus time step. Each line segment is obtained by least squares fitting of all corresponding points.
Figure 2. CPU time versus time step. Each line segment is obtained by least squares fitting of all corresponding points.
Mathematics 08 01502 g002
Figure 3. (a) Phase diagram (Reprinted with permission from [1]). Here, ψ ¯ denotes the averaged field ( = ϕ ¯ ). (b) Values of Λ ( t ) at t = 2176 with various ϕ ¯ and ϵ .
Figure 3. (a) Phase diagram (Reprinted with permission from [1]). Here, ψ ¯ denotes the averaged field ( = ϕ ¯ ). (b) Values of Λ ( t ) at t = 2176 with various ϕ ¯ and ϵ .
Mathematics 08 01502 g003
Figure 4. Evolution of ϕ ( x , y , t ) using the method S3 with ( ϕ ¯ , ϵ ) = ( 0.02 , 0.1 ) . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Figure 4. Evolution of ϕ ( x , y , t ) using the method S3 with ( ϕ ¯ , ϵ ) = ( 0.02 , 0.1 ) . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Mathematics 08 01502 g004
Figure 5. Evolution of ϕ ( x , y , t ) using the method S3 with ( ϕ ¯ , ϵ ) = ( 0.14 , 0.1 ) . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3926 , 0.0108 , and 0.3710 , respectively.
Figure 5. Evolution of ϕ ( x , y , t ) using the method S3 with ( ϕ ¯ , ϵ ) = ( 0.14 , 0.1 ) . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3926 , 0.0108 , and 0.3710 , respectively.
Mathematics 08 01502 g005
Figure 6. Evolution of ϕ ( x , y , t ) using the method S3 with ( ϕ ¯ , ϵ ) = ( 0.26 , 0.1 ) . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.2610 , 0.2600 , and 0.2590 , respectively.
Figure 6. Evolution of ϕ ( x , y , t ) using the method S3 with ( ϕ ¯ , ϵ ) = ( 0.26 , 0.1 ) . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.2610 , 0.2600 , and 0.2590 , respectively.
Mathematics 08 01502 g006
Figure 7. Evolutions of (a) E ( t ) and (b) Λ ( t ) with ( ϕ ¯ , ϵ ) used in Figure 4, Figure 5 and Figure 6.
Figure 7. Evolutions of (a) E ( t ) and (b) Λ ( t ) with ( ϕ ¯ , ϵ ) used in Figure 4, Figure 5 and Figure 6.
Mathematics 08 01502 g007
Figure 8. Evolution of ϕ ( x , y , t ) using the method OS2 with Δ t = 1 4 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Figure 8. Evolution of ϕ ( x , y , t ) using the method OS2 with Δ t = 1 4 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Mathematics 08 01502 g008
Figure 9. Evolution of ϕ ( x , y , t ) using the method OS2 with Δ t = 2 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Figure 9. Evolution of ϕ ( x , y , t ) using the method OS2 with Δ t = 2 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Mathematics 08 01502 g009
Figure 10. Evolution of ϕ ( x , y , t ) using the method S2 with Δ t = 1 4 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Figure 10. Evolution of ϕ ( x , y , t ) using the method S2 with Δ t = 1 4 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Mathematics 08 01502 g010
Figure 11. Evolution of ϕ ( x , y , t ) using the method S2 with Δ t = 2 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Figure 11. Evolution of ϕ ( x , y , t ) using the method S2 with Δ t = 2 . In each snapshot, the yellow, green, and blue regions indicate ϕ = 0.3810 , 0.0195 , and 0.3420 , respectively.
Mathematics 08 01502 g011
Figure 12. Evolutions of E ( t ) for Figure 8, Figure 9, Figure 10 and Figure 11.
Figure 12. Evolutions of E ( t ) for Figure 8, Figure 9, Figure 10 and Figure 11.
Mathematics 08 01502 g012

Share and Cite

MDPI and ACS Style

Lee, H.G. An Efficient and Accurate Method for the Conservative Swift–Hohenberg Equation and Its Numerical Implementation. Mathematics 2020, 8, 1502. https://doi.org/10.3390/math8091502

AMA Style

Lee HG. An Efficient and Accurate Method for the Conservative Swift–Hohenberg Equation and Its Numerical Implementation. Mathematics. 2020; 8(9):1502. https://doi.org/10.3390/math8091502

Chicago/Turabian Style

Lee, Hyun Geun. 2020. "An Efficient and Accurate Method for the Conservative Swift–Hohenberg Equation and Its Numerical Implementation" Mathematics 8, no. 9: 1502. https://doi.org/10.3390/math8091502

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop