Next Article in Journal
A Multi-Attribute Decision-Making Algorithm Using Q-Rung Orthopair Power Bonferroni Mean Operator and Its Application
Previous Article in Journal
Cu-Al2O3/Water Hybrid Nanofluid Stagnation Point Flow Past MHD Stretching/Shrinking Sheet in Presence of Homogeneous-Heterogeneous and Convective Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fourth Order Energy Dissipative Scheme for a Traffic Flow Model

1
College of Liberal Arts and Sciences, National University of Defense Technology, Changsha 410073, China
2
State Key Laboratory of High Performance Computing, National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(8), 1238; https://doi.org/10.3390/math8081238
Submission received: 24 June 2020 / Revised: 19 July 2020 / Accepted: 23 July 2020 / Published: 28 July 2020

Abstract

:
We propose, analyze and numerically validate a new energy dissipative scheme for the Ginzburg–Landau equation by using the invariant energy quadratization approach. First, the Ginzburg–Landau equation is transformed into an equivalent formulation which possesses the quadratic energy dissipation law. After the space-discretization of the Fourier pseudo-spectral method, the semi-discrete system is proved to be energy dissipative. Using diagonally implicit Runge–Kutta scheme, the semi-discrete system is integrated in the time direction. Then the presented full-discrete scheme preserves the energy dissipation, which is beneficial to the numerical stability in long-time simulations. Several numerical experiments are provided to illustrate the effectiveness of the proposed scheme and verify the theoretical analysis.

1. Introduction

The Ginzburg–Landau equation (GLE) was first proposed by Ginzburg and Landau in 1950 and applied to model the electrodynamics, quantum mechanics and thermodynamic properties of superconductors. As a state transformation equation, the GLE is closely related to many other state transformation equations, such as Allen–Cahn equation and Chafee–Infante equation. As is well known, the GLE plays an important role in the study of state transformation and unstable wave theory. Moreover, it has been widely used in optical fiber communication. In this paper, we consider the one-dimensional GLE as a model for traffic flow [1]
u t = ( x + ϵ x 2 ) ( 6 u + x 2 u u 3 ) ,
subject to the initial and boundary conditions
u ( x , 0 ) = u 0 ( x ) , x Ω , u ( a , t ) = u ( b , t ) , t 0 ,
where ϵ 0 and Ω = [ a , b ] R is the bounded domain. We denote the L 2 inner product by f , g = Ω f g d x , the L 2 norm f 2 = f , f 1 2 . Under the periodic boundary conditions, the global energy of the GLE is
H ( u ) = Ω ( 3 u 2 1 2 u x 2 1 4 f ( u ) ) d x ,
where f ( u ) = u 4 . Then the GLE in [1] can be rewritten as
u t = ( x + ϵ x x ) m ,
m = δ H ( u ) δ u .
Taking the inner product of the Equation (2) with m and the Equation (3) with u t , we obtain
d d t H ( u ) = x m , m + ϵ x x m , m = ϵ x m 2 0 .
Obviously, the GLE satisfies the energy dissipation law.
In recent years, many researchers have conducted various works about the GLE. Celledoni [1] proved that the average vector field (AVF) method can keep the energy dissipation of the Equation (1). The time-dependent GLE for the car traffic model is derived by Nagatani [2] through the perturbation method. Gunzburger and Hou [3] solved a series of problems for a Ginzburg–Landau model of superconductivity. They proposed finite element approximations for the optimality system and derived corresponding error estimates. The local well-posedness of the initial-boundary value problem for the complex GLE in bounded domains is investigated by Kuroda and Takanori [4]. Meanwhile, there are many effective ways to solve the GLE, such as the discontinuous Galerkin method [5], the semi-implicit method [6] and the trial equation method [7]. In addition, there also exist relevant studies on the fractional GLE [8,9,10,11] and the stochastic GLE [12,13,14,15,16]. However, most of previous researches can only achieve second order accuracy in the time direction. Therefore, the main challenge issue is to construct a high order energy dissipative scheme for the GLE.
The invariant energy quadratization (IEQ) approach [17] is a new method developed in recent years by introducing a Lagrange multiplier and has been successfully applied to a variety of gradient flow models [18,19,20,21,22,23,24]. By introducing a Lagrange multiplier and converting the nonlinear part to a semi-implicit scheme, a new linear system consisting of two coupled equations is calculated at each time step. A time-discretization scheme for the linear systems is constructed by Yang [25], which is easy to implement and achieve unconditionally energy stability. Liu [26] applyed two novel and efficient modified techniques based on IEQ approach to deal with nonlinear terms in gradient flows which succeeded in finding suitable positive preserving functions. Two linear and second order schemes that combine the IEQ approach with the stabilization technique are developed by Xu [27], where several extra linear stabilization terms are added and they are crucial to suppress the non-physical spatial oscillations caused by the strong anisotropy. Gong [28] presented a novel class of arbitrarily high-order and unconditionally energy-stable algorithms for gradient flow models by combining the IEQ approach and a specific class of Runge–Kutta (RK) methods. Since the solution to the Ginzburg–Landau equation shows instability very soon and this high order implicit scheme can unconditionally keep Ginzburg–Landau equation energy dissipation, we adopt this method to solve the Ginzburg–Landau equation. Ref. [28] combines IEQ approach and RK method and apply them to Cahn–Hilliard equation and Allen–Cahn equation to get the arbitrarily high-order schemes, while we apply it to the Ginzburg–Landau equation to construct a fourth order scheme in this paper.
The rest of this paper is organized as follows. In Section 2, the IEQ approach is introduced and applied to the GLE. Then the GLE is reformulated as a new system which admits a quadratic energy dissipation law. In Section 3, the Fourier pseudo-spectral method is used to obtain the semi-discrete the system of the GLE. After that, the implicit fourth order Runge–Kutta (RK4) method is applied for the time discretization to structure a full-discrete scheme. It is proved that the scheme we proposed is energy dissipative. In Section 4, the stability of the proposed scheme is verified by numerical experiments. We make a summary of this work in Section 5.

2. Invariant Energy Quadratization Approach

In this section, we present a new scheme for the GLE to preserve the energy dissipation law. By using the IEQ approach, we transform the Hamlitonian functional into a quadratic functional. We first introduce a Lagrange multiplier r = f ( u ) , then the GLE is transformed into
u t = ( x + ϵ x x ) ( 6 u + x x u u r ) , r t = 2 u u t .
Theorem 1.
The system above preserves the energy dissipation law
d d t H ( u ) 0 ,
where H ( u ) = Ω ( 3 u 2 1 2 | x u 2 | 1 4 r 2 ) d x .
Proof. 
The system (4) can be rewritten as
u t = ( x + ϵ x x ) m ,
m = 6 u + x x u r 4 f ( u ) f ( u ) ,
r t = f ( u ) 2 f ( u ) u t .
By taking the inner products of the Equations (5)–(7) with m, u t and r 2 , respectively, we have the following equalities
u t , m = ( x + ϵ x x ) m , m ,
m , u t = 6 u + x x u r 4 f ( u ) f ( u ) , u t ,
r t , r 2 = f ( u ) 2 f ( u ) u t , r 2 .
By substituting Equations (9) and (10) into Equation (8), we obtain
d d t H ( u ) = d d t ( 1 2 x u 2 + Ω ( 3 u 2 1 4 r 2 ) d x ) = ( x + ϵ x x ) m , m = ϵ x m 2 0 .

3. High Order Energy Dissipative Scheme

3.1. Spatial Discretization by Using Fourier Pseudo-Spectral Method

In this part, we choose the Fourier pseudo-spectral method to discretize the system (4). It is proved that the semi-discrete scheme for the GLE preserves the semi-discrete energy dissipation law.
We primarily divide the domain Ω = [ a , b ] into N equal intervals, where N is an even number. Afterwards, denote h = b a N as the spatial step and x j : = a + j h ( j = 0 , 1 , , N 1 ) as the grid points. The kth-order partial differential optrator x k can be approximated by the Fourier spectral differential matrix D 1 k
( D 1 ) j + 1 , l + 1 = 1 2 ( 1 ) j + l μ c o t ( μ x j a 2 ) , j l , 0 , j = l ,
where μ = 2 π b a , j , l = 0 , 1 , , N 1 . Clearly, D 1 is an N × N skew-symmetric matrix.
By applying Fourier pseudo-spectral method to the system (4), we acquire the semi-discrete scheme
u t = ( D 1 + ϵ D 1 2 ) ( 6 u + ϵ D 1 2 u u · r ) , r t = 2 u · u t ,
where u = [ u 0 , u 1 , , u N 1 ] T , r = [ r 0 , r 1 , , r N 1 ] T and u · r = [ u 0 r 0 , u 1 r 1 , , u N 1 r N 1 ] T .
Theorem 2.
The system (11) preserves the semi-discrete energy dissipation law
d d t H ( u ) 0 ,
where H ( u ) = 1 2 D 1 u 2 + i = 0 N 1 ( 3 u i 2 1 4 r i 2 ) .
Proof. 
The system (11) can be rewritten as follows
u t = ( D 1 + ϵ D 1 2 ) m , m = 6 u 2 + D 1 2 u r 4 f ( u ) f ( u ) , r t = f ( u ) 2 f ( u ) u t .
Taking the inner products of the system (12) with m , u t and r 2 , respectively, we obtain the following semi-discrete energy dissipation law
d d t H ( u , r ) = d d t ( 1 2 D 1 u 2 + i = 0 N 1 3 u i 2 1 4 r i 2 ) = ( D 1 + ϵ D 1 2 ) m , m = ϵ D 1 m 2 0 .
It is known that the Equation (13) is equivalent to
d d t S ω , ω Ω = 2 S F ( ω ) , ω Ω = 2 S ω t , ω Ω 0 , ω Ω
where
S = 3 I + 1 2 D 1 2 0 0 I 4
is a symmetric and constant matrix.

3.2. Full Discretization by Using Implicit High Order Runge–Kutta Method

In previous studies, many methods have been proposed to solve energy dissipation system, such as the second order backward differentiation formula (BDF2), the AVF method and the backward Euler method. However, these methods can only structure first or second order energy dissipative schemes in the time direction meanwhile these schemes always require a lot of iterative computations, which not only affect the efficiency, but also bring cumulative errors in the calculation processes. Therefore, we propose a new scheme to acquire high order energy dissipative schemes for the GLE.
Let b i , a i j , ( i , j = 1 , 2 , , s ) be real numbers and c i = j = 1 s a i j . For given z n = z ( t n ) , when an s-stage RK method is applied to the Equation (1), the following intermediate values are first calculated by
z i = z n + Δ t j = 1 s a i j F j ,
where F j = F ( z j ) . Then the high order energy dissipation scheme of the GLE is updated via
z n + 1 = z n + Δ t i = 1 s b i F i .
Proposition 1.
The RK coefficients are commonly displayed by a Butcher table
cA
b T
=
c 1 a 11 a 1 s
c s a s 1 a s s
       b 1 b s ,
where A R s , s , b R s and c = A I with I = ( 1 , 1 , , 1 ) T R s .
In consequence, we choose RK(4,4) as
b 1 2 b 1 2
b 1 + b 2 2 b 1 b 2 2
b 1 + b 2 + b 3 2 b 1 b 2 b 3 2
b 1 + b 2 + b 3 + b 4 2 b 1 b 2 b 3 b 4 2
       b 1 b 2 b 3 b 4 ,
where b 1 = 2.70309412 , b 2 = 0.53652708 , b 3 = 2.37893931 , b 4 = 1.86068189 .
Definition 1.
Let S be a symmetric, constant matrix. For a dissipative system, the RK method defined by the Butcher array is said to be S -dissipative if applid to the new system (15)
S z n + 1 , z n + 1 S z n , z n
for all n 0 and for all Δ t 0 .
Theorem 3.
([29]) A RK scheme defined by the Butcher array is said to be S -dissipative, if and only if the symmetric coefficient matrix M with elements
m i j = b i a i j + b j a j i b i b j , i , j = 1 , 2 , , s
are nonnegative.
Proof. 
Above all, replace z n + 1 by the new system (15), that can be derived as
Sz n + 1 , z n + 1 Sz n , z n = 2 Δ t i = 1 s b i Sz n , F i + Δ t 2 i , j = 1 s b i b j S F i , F j .
Substituting z n = z i Δ t j = 1 s a i j F j into the first term on the right hand side of the Equation (17), we get
Sz n + 1 , z n + 1 Sz n , z n = 2 Δ t i = 1 s b i Sz i , F i Δ t 2 i , j = 1 s b i a i j S F i , F j Δ t 2 i , j = 1 s b j a j i S F i , F j + Δ t 2 i , j = 1 s b i b j S F i , F j .
Whereafter, it can be arranged as
Sz n + 1 , z n + 1 Sz n , z n = 2 Δ t i = 1 s b i Sz i , F i Δ t 2 i , j = 1 s m i j S F i , F j .
Noticing that S z i , F j 0 and m i j 0 , we have
Sz n + 1 , z n + 1 Sz n , z n .
This completes the proof. □
In this paper, we only consider four stage and fourth order diagonally implicit RK scheme, which belongs to S -dissipative RK scheme can be extended to higher order RK method without any necessary difficulties.

4. Numerical Simulation

In this section, we carry out two numerical experiments by applying the implicit high order energy dissipative scheme. The stability and the convergence of the new scheme (15) are verified by these numerical experiments. Meanwhile, we compare the proposed scheme (15) with the AVF method in energy errors and accuracy. The numerical iteration formula of the AVF method is
u n + 1 u n Δ t = ( 3 E + E 2 D 2 ) ( u n + 1 + u n ) E 4 ( u n ) 3 + ( u n + 1 ) 3 + u n ( u n + 1 ) 2 + ( u n ) 2 u n + 1 ,
where E = D 1 + ϵ D 2 . In the following text, the Newton iteration is mainly used to solve the full-discrete system.
Example 1.
We utilize the initial condition of single solitary wave given by
u ( x , 0 ) = e M 1 ( x 1 2 ) 2 , x [ 10 , 10 ]
to observe how the coefficient ϵ affects the energy dissipation. First, we set M 1 = 0.3 , the number of grid points N = 200 and the time step Δ t = 0.001 . Letting ε = 0 , 0.00025 , 0.0005 , 0.00075 , 0.001 respectively, the energy variations and energy errors of the proposed scheme (15) are presented in Figure 1. Based on Figure 1a, we find the proposed scheme (15) can preserve the discrete energy conservation exactly when ϵ = 0 , and the energy dissipations of the GLE can be observed when ϵ > 0 . Obviously, with bigger ϵ, the energy dissipates faster. Taking the discrete energy obtained from the scheme (15) with Δ t = 0.0001 and N = 200 as the reference discrete energy, we plot the energy errors ( Δ t = 0.001 ) with various ε in Figure 1b. It shows that the proposed new scheme (15) can preserve the energy dissipation of the GLE. Meanwhile, Figure 1b indicates that the energy dissipation is a little bit influenced by the time step size. The waveforms, energy variations and energy errors obtained from the scheme (15) and the scheme (18) with ε = 0.001 , N = 100 , Δ t = 0.001 are pictured in Figure 2. The waveforms obtained from these two schemes at t = 0.1 are almost the same (Figure 2a). Figure 2b shows that they both can keep the energy reducing steadily. According to Figure 2c, comparatively speaking, the energy errors obtained by the proposed scheme (15) is even smaller than that by the scheme (18).
Setting M 1 = 0.5 and ε = 0.001 , we list the convergence rates of the proposed new scheme (15) in Table 1 and Table 2. For the time accuracy test, we select the numerical solution with Δ t = 0.0001 and N = 50 as the reference solution. Based on Table 1, the convergence order of the proposed scheme (15) in time is 4 and the scheme (18) only converges with second-order accuracy. For the space accuracy test, the scheme (15) and the scheme (18) both apply the Fourier pseudo-spectral method in space-discretization, thus the space convergence orders of them decrease in an exponential rate in Table 2 (the numerical solution with N = 560 and Δ t = 0.001 is chosen as the reference solution).
Example 2.
In the second example, we investigate the energy dissipation of double solitary wave with the initial condition,
u ( x , 0 ) = e M 1 ( x 2 ) 2 + M 2 ( x + 2 ) 2 , x [ 10 , 10 ] .
By setting M 1 = M 2 = 1 and ε = 0.001 , we plot the waveform of double solitary wave with N = 200 and Δ t = 0.001 in Figure 3a. The energy dissipation property of the GLE with various time steps can be verified in Figure 3b. It shows the new scheme (15) can keep energy dissipation at a large time step. Taking the discrete energy obtained from the new scheme (15) with Δ t = 0.0001 as the reference discrete energy, the energy error in Figure 3c is only about O ( 10 7 ) . Then using the same grid, we consider the situation of M 1 = 0.2 and M 2 = 2 and plot the energy variation and the energy error in Figure 4. As we can see, the proposed scheme (15) still can maintain the energy dissipation property with various time steps.

5. Conclusions

In this paper, we have developed a stable and high order energy dissipative scheme for solving the GLE. By using the IEQ approach, the GLE is downgraded from fourth order to second order. The full-discrete system of GLE is obtained by using the Fourier pseudo-spectral method in the spatial direction and the diagonally implicit Runge–Kutta method in the time direction. The new scheme (15) we proposed has three main advantages. First of all, it is well-known that the implicit methods enjoy better stability in the simulation of stiff equations. Besides, it has high calculation efficiency due to the application of the diagonally implicit Runge–Kutta method. Last but not least, the high order scheme in time direction can be constructed and the energy dissipation can be preserved. The proposed implicit fourth order scheme can maintain the energy dissipation of the GLE theoretically. Moreover, the numerical results of single solitary wave and double solitary wave have verified the stability and accuracy of the proposed scheme. Compared with the scheme (18), the proposed scheme converges with fourth order in the time direction and can be generalized to construct higher order energy dissipative scheme naturally.

Author Contributions

Conceptualization, X.C.; Project administration, S.S.; Writing—original draft, X.C.; Writing—review editing, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (Grant No. 11901577 11971481) and the National Natural Science Foundation of Hunan (Grant No. S2017JJQNJJ0764), Research Fund of National University of Defense Technology (Grant No. ZK17-03-27), and the fund from Hunan Provincial Key Laboratory of Mathematical Modeling and Analysis in Engineering (Grant No. 2018MMAEZD004).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Celledoni, E.; Grimm, V.; McLachlan, R.I.; McLaren, D.I.; O’Neale, D.; Owren, B.; Quispel, G.R.W. Preserving energy resp. dissipation in numerical PDEs using the “Average Vector Field” method. J. Comput. Phys. 2012, 231, 6770–6789. [Google Scholar] [CrossRef] [Green Version]
  2. Nagatani, T. Time-dependent Ginzburg–Landau equation for the jamming transition in traffic flow. Phys. Stat. Mech. Appl. 1998, 258, 237–242. [Google Scholar] [CrossRef]
  3. Gunzburger, M.D.; Hou, L.S.; Ravindran, S.S. Analysis and approximation of optimal control problems for a simplified Ginzburg–Landau model of superconductivity. Numer. Math. 1997, 77, 243–268. [Google Scholar] [CrossRef]
  4. Kuroda, T.; Ôtani, M. Local well-posedness of the complex Ginzburg–Landau equation in bounded domains. Nonlinear Anal. Real World Appl. 2019, 45, 877–894. [Google Scholar] [CrossRef] [Green Version]
  5. Yildiz, T.A.; Uzunca, M.; Karasozen, B. Structure preserving reduced order modeling for gradient systems. Appl. Math. Comput. 2019, 347, 194–209. [Google Scholar]
  6. Chen, L.Q.; Shen, J. Applications of semi-implicit Fourier-spectral method to phase field equations. Comput. Phys. Commun. 1998, 108, 147–158. [Google Scholar] [CrossRef]
  7. Liu, Y.; Chen, S.Q.; Wei, L.X.; Guan, B. Exact solutions to complex Ginzburg–Landau equation. Pramana J. Phys. 2018, 91. [Google Scholar] [CrossRef]
  8. Zhang, Q.G.; Li, Y.N.; Su, M.L. The local and global existence of solutions for a time fractional complex Ginzburg–Landau equation. J. Math. Anal. Appl. 2019, 469, 16–43. [Google Scholar] [CrossRef]
  9. He, D.D.; Pan, K.J. An unconditionally stable linearized difference scheme for the fractional Ginzburg–Landau equation. Numer. Algorithms 2018, 79, 899–925. [Google Scholar] [CrossRef]
  10. Li, L.; Jin, L.Y.; Fang, S.M. Large time behavior for the fractional Ginzburg–Landau equations near the BCS-BEC crossover regime of Fermi gases. Bound. Value Probl. 2017, 2017, 8. [Google Scholar] [CrossRef] [Green Version]
  11. Lu, H.; Bates, P.W.; Lu, S.J.; Zhang, M.J. Dynamics of the 3D fractional Ginzburg–Landau equation with multiplicative noise on an unbounded domain. Commun. Math. Sci. 2016, 14, 273–295. [Google Scholar] [CrossRef]
  12. Guo, Y.F.; Li, D.L. Random attractor stochastic complex Ginzburg–Landau equation with multiplicative noise ou unbounded domain. Stoch. Anal. Appl. 2017, 35, 409–422. [Google Scholar] [CrossRef]
  13. Shen, T.L.; Huang, J.H. Ergodicity of 2D stochastic Ginzburg–Landau-Newell equations driven by degenerate noise. Math. Method Appl. Sci. 2017, 40, 4812–4831. [Google Scholar]
  14. Shen, T.L.; Xin, J.; Huang, J.H. Time-space fractional stochastic Ginzburg–Landau equation driven by Gaussian white noise. Stoch. Anal. Appl. 2018, 36, 103–113. [Google Scholar] [CrossRef]
  15. Chugreeva, O.; Melcher, C. Vortices in a stochastic parabolic Ginzburg–Landau equation. Stochastics Partial. Differ. Equations Anal. Comput. 2017, 5, 113–143. [Google Scholar] [CrossRef] [Green Version]
  16. Lin, L.; Gao, H.J. A stochastic generalized Ginzburg–Landau equation driven by jump noise. J. Theor. Probab. 2019, 32, 460–483. [Google Scholar] [CrossRef] [Green Version]
  17. Guillen, G.F.; Tierra, G. On linear schemes for a Cahn-Hilliard diffuse interface model. J. Comput. Phys. 2013, 234, 140–171. [Google Scholar] [CrossRef]
  18. Han, D.Z.; Brylev, A.; Yang, X.F.; Tan, Z.J. Numerical analysis of second order, fully discrete energy stable schemes for phase field models of two-phase incompressible flows. J. Sci. Comput. 2017, 70, 965–989. [Google Scholar] [CrossRef]
  19. Zhao, J.; Yang, X.F.; Gong, Y.Z.; Wang, Q. A novel linear second order unconditionally energy stable scheme for a hydrodynamic q-tensor model of liquid crystals. Comput. Method Appl. Mech. Eng. 2017, 318, 803–825. [Google Scholar] [CrossRef] [Green Version]
  20. Yang, X.F. Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends. J. Comput. Phys. 2016, 327, 294–316. [Google Scholar] [CrossRef] [Green Version]
  21. Yang, X.F.; Han, D.Z. Linearly first and second-order, unconditionally energy stable schemes for the phase field crystal model. J. Comput. Phys. 2017, 33, 1116–1134. [Google Scholar] [CrossRef] [Green Version]
  22. Yang, X.F.; Ju, L.L. Efficient linear schemes with unconditional energy stability for the phase field elastic bending energy model. Comput. Method Appl. Mech. Eng. 2017, 315, 691–712. [Google Scholar] [CrossRef] [Green Version]
  23. Yang, X.F.; Ju, L.L. Linear and unconditionally energy stable schemes for the binary fluid-surfactant phase field model. Comput. Method Appl. Mech. Eng. 2017, 318, 1005–1029. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, X.F.; Zhao, J.; Wang, Q. Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method. J. Comput. Phys. 2017, 333, 104–127. [Google Scholar] [CrossRef] [Green Version]
  25. Yang, X.F.; Zhao, J.; Wang, Q.; Shen, J. Numerical approximations for a three-component Cahn-Hilliard phase-field model based on the invariant energy quadratization method. Math. Model. Methods Appl. Sci. 2017, 27, 1993–2030. [Google Scholar] [CrossRef]
  26. Liu, Z.G.; Li, X.L. Efficient modified techniques of invariant energy quadratization approach for gradient flows. Appl. Math. Lett. 2019, 98, 206–214. [Google Scholar] [CrossRef]
  27. Xu, Z.; Yang, X.F.; Zhang, H.; Xie, Z.Q. Efficient and linear schemes for anisotropic Cahn-Hilliard model using the Stabilized-Invariant Energy Quadratization (S-IEQ) approach. Comput. Phys. Commun. 2019, 238, 36–49. [Google Scholar] [CrossRef]
  28. Gong, Y.Z.; Zhao, J. Energy-stable Runge–Kutta schemes for gradient flow models using the energy quadratization approach. Appl. Math. Lett. 2019, 94, 224–231. [Google Scholar] [CrossRef]
  29. Buono, N.D.; Mastroserio, C. Explicit methods based on a class of four stage fourth order Runge–Kutta methods for preserving quadratic laws. J. Comput. Appl. Math. 2002, 140, 231–243. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The energy variations (a) and energy errors (b) of the scheme (15) with various ϵ .
Figure 1. The energy variations (a) and energy errors (b) of the scheme (15) with various ϵ .
Mathematics 08 01238 g001
Figure 2. The waveforms (a), energy variationas (b) and energy errors (c) of the scheme (15) and the scheme (18).
Figure 2. The waveforms (a), energy variationas (b) and energy errors (c) of the scheme (15) and the scheme (18).
Mathematics 08 01238 g002
Figure 3. The waveform (a), energy variation (b) and energy error (c) of the scheme (15) at t = 0.1 with M 1 = 1 and M 2 = 1 .
Figure 3. The waveform (a), energy variation (b) and energy error (c) of the scheme (15) at t = 0.1 with M 1 = 1 and M 2 = 1 .
Mathematics 08 01238 g003
Figure 4. The energy variation (a) and energy error (b) of the scheme (15) at t = 0.1 with M 1 = 0.2 and M 2 = 2 .
Figure 4. The energy variation (a) and energy error (b) of the scheme (15) at t = 0.1 with M 1 = 0.2 and M 2 = 2 .
Mathematics 08 01238 g004
Table 1. Time accuracy test for the scheme (15) and the scheme (18) at t = 0.1 .
Table 1. Time accuracy test for the scheme (15) and the scheme (18) at t = 0.1 .
Δ tNew Scheme (15)Time OrderScheme (18)Time Order
0.002 3.134 × 10 7 - 3.074 × 10 6 -
0.001 1.967 × 10 8 3.994 7.629 × 10 7 2.016
0.0005 1.230 × 10 9 3.999 1.850 × 10 7 2.044
0.00025 7.688 × 10 11 4.000 4.046 × 10 8 2.193
Table 2. Space accuracy test for the scheme (15) and the scheme (18) at t = 0.1 .
Table 2. Space accuracy test for the scheme (15) and the scheme (18) at t = 0.1 .
NNew Scheme (15)Space OrderScheme (18)Space Order
10 6.673 × 10 2 - 6.671 × 10 2 -
20 3.350 × 10 2 0.994 1.858 × 10 2 1.844
40 2.754 × 10 5 10.253 2.092 × 10 5 9.794
80 5.806 × 10 10 15.533 3.242 × 10 10 15.978

Share and Cite

MDPI and ACS Style

Chen, X.; Song, M.; Song, S. A Fourth Order Energy Dissipative Scheme for a Traffic Flow Model. Mathematics 2020, 8, 1238. https://doi.org/10.3390/math8081238

AMA Style

Chen X, Song M, Song S. A Fourth Order Energy Dissipative Scheme for a Traffic Flow Model. Mathematics. 2020; 8(8):1238. https://doi.org/10.3390/math8081238

Chicago/Turabian Style

Chen, Xiaowei, Mingzhan Song, and Songhe Song. 2020. "A Fourth Order Energy Dissipative Scheme for a Traffic Flow Model" Mathematics 8, no. 8: 1238. https://doi.org/10.3390/math8081238

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