A novel high-order symmetric and energy-preserving continuous-stage Runge-Kutta-Nyström Fourier pseudo-spectral scheme for solving the two-dimensional nonlinear wave equation
The primary objective of this research is to develop a novel high-order symmetric and energy-preserving method for solving two-dimensional nonlinear wave equations. Initially, the nonlinear wave equation is reformulated as an abstract Hamiltonian ordinary differential equation (ODE) system with separable energy in an appropriate infinite-dimensional function space. Subsequently, an energy-preserving and symmetric continuous-stage Runge-Kutta-Nyström time-stepping scheme is derived. After approximating the spatial differential operator using the two-dimensional Fourier pseudo-spectral method, we derive an energy-preserving fully discrete scheme. A rigorous error analysis demonstrates that the proposed method can achieve at least fourth-order accuracy in time. Finally, numerical examples are provided to validate the accuracy, efficiency, and long-term energy conservation of the method.
Citation: Dongjie Gao, Peiguo Zhang, Longqin Wang, Zhenlong Dai, Yonglei Fang. A novel high-order symmetric and energy-preserving continuous-stage Runge-Kutta-Nyström Fourier pseudo-spectral scheme for solving the two-dimensional nonlinear wave equation[J]. AIMS Mathematics, 2025, 10(3): 6764-6787. doi: 10.3934/math.2025310
Related Papers:
[1]
Patarawadee Prasertsang, Thongchai Botmart .
Improvement of finite-time stability for delayed neural networks via a new Lyapunov-Krasovskii functional. AIMS Mathematics, 2021, 6(1): 998-1023.
doi: 10.3934/math.2021060
[2]
R. Sriraman, P. Vignesh, V. C. Amritha, G. Rachakit, Prasanalakshmi Balaji .
Direct quaternion method-based stability criteria for quaternion-valued Takagi-Sugeno fuzzy BAM delayed neural networks using quaternion-valued Wirtinger-based integral inequality. AIMS Mathematics, 2023, 8(5): 10486-10512.
doi: 10.3934/math.2023532
[3]
Sunisa Luemsai, Thongchai Botmart, Wajaree Weera, Suphachai Charoensin .
Improved results on mixed passive and H∞ performance for uncertain neural networks with mixed interval time-varying delays via feedback control. AIMS Mathematics, 2021, 6(3): 2653-2679.
doi: 10.3934/math.2021161
[4]
Jenjira Thipcha, Presarin Tangsiridamrong, Thongchai Botmart, Boonyachat Meesuptong, M. Syed Ali, Pantiwa Srisilp, Kanit Mukdasai .
Robust stability and passivity analysis for discrete-time neural networks with mixed time-varying delays via a new summation inequality. AIMS Mathematics, 2023, 8(2): 4973-5006.
doi: 10.3934/math.2023249
[5]
Abdulaziz M. Alanazi, R. Sriraman, R. Gurusamy, S. Athithan, P. Vignesh, Zaid Bassfar, Adel R. Alharbi, Amer Aljaedi .
System decomposition method-based global stability criteria for T-S fuzzy Clifford-valued delayed neural networks with impulses and leakage term. AIMS Mathematics, 2023, 8(7): 15166-15188.
doi: 10.3934/math.2023774
[6]
Huahai Qiu, Li Wan, Zhigang Zhou, Qunjiao Zhang, Qinghua Zhou .
Global exponential periodicity of nonlinear neural networks with multiple time-varying delays. AIMS Mathematics, 2023, 8(5): 12472-12485.
doi: 10.3934/math.2023626
[7]
Zirui Zhao, Wenjuan Lin .
Extended dissipative analysis for memristive neural networks with two-delay components via a generalized delay-product-type Lyapunov-Krasovskii functional. AIMS Mathematics, 2023, 8(12): 30777-30789.
doi: 10.3934/math.20231573
[8]
Qinghua Zhou, Li Wan, Hongbo Fu, Qunjiao Zhang .
Exponential stability of stochastic Hopfield neural network with mixed multiple delays. AIMS Mathematics, 2021, 6(4): 4142-4155.
doi: 10.3934/math.2021245
[9]
Nayika Samorn, Kanit Mukdasai, Issaraporn Khonchaiyaphum .
Analysis of finite-time stability in genetic regulatory networks with interval time-varying delays and leakage delay effects. AIMS Mathematics, 2024, 9(9): 25028-25048.
doi: 10.3934/math.20241220
[10]
Yude Ji, Xitong Ma, Luyao Wang, Yanqing Xing .
Novel stability criterion for linear system with two additive time-varying delays using general integral inequalities. AIMS Mathematics, 2021, 6(8): 8667-8680.
doi: 10.3934/math.2021504
Abstract
The primary objective of this research is to develop a novel high-order symmetric and energy-preserving method for solving two-dimensional nonlinear wave equations. Initially, the nonlinear wave equation is reformulated as an abstract Hamiltonian ordinary differential equation (ODE) system with separable energy in an appropriate infinite-dimensional function space. Subsequently, an energy-preserving and symmetric continuous-stage Runge-Kutta-Nyström time-stepping scheme is derived. After approximating the spatial differential operator using the two-dimensional Fourier pseudo-spectral method, we derive an energy-preserving fully discrete scheme. A rigorous error analysis demonstrates that the proposed method can achieve at least fourth-order accuracy in time. Finally, numerical examples are provided to validate the accuracy, efficiency, and long-term energy conservation of the method.
1.
Introduction
Over recent decades, neural networks (NNs) have garnered significant attention and demonstrated success across various engineering and research domains, thereby encompassing image processing, pattern recognition, optimization problems, and associative memory [1,2,3]. Stability properties, which are crucial for effective neural network deployment, include asymptotic and exponential stability. Time delays which are prevalent in numerous control systems [4], poses challenges by potentially destabilizing systems. Consequently, stability analysis, particularly regarding NNs with time delays, is imperative due to the substantial impact of equilibrium point dynamics on practical applications [5,6].
The Lyapunov-Krasovskii functional (LKF) method is a powerful tool for examining the stability of a system. Its effectiveness lies in its ability to identify a positive definite function whose derivative is negative definite along the trajectory of the system [7,8,9]. The choice of an appropriate LKF is crucial to establish the stability criteria. In a notable study [10], a potent methodology known as the delay-product-type functional (DPTF) method was introduced. This method is distinguished by its inclusion of variables that are dependent on the amplitude of the delay. For example, a DPTF V(t) is formulated as V(t)=d(t)νT1(t)P1ν1(t)+(h−d(t))νT2(t)P2ν2(t), where d(t)P1>0, (h−d(t))P2>0 are delay amplitude-dependent matrices, and ν1(t), ν2(t) denote augmented terms related to the state. Importantly, the derivative of V(t) unveils the interconnection among terms related to the state, delay amplitude, and delay derivative, which is subsequently integrated into the final linear matrix inequalities (LMIs). Afterwards, in [11], by using Wirtinger-based integral inequality, the nonintegral terms are connected to the integral terms. As shown in [12], it is efficient for reduction of the conservatism if some double integral terms are introduced in Lyapunov functionals. Nevertheless, there exists an intrinsic conservatism in the LKF due to incomplete state vectors ν1(t) and ν2(t). Various bounding methods have been developed for the stability analysis, including Jensen-based and Wirtinger-based integral inequalities [13,14], as well as slack-matrix-based integral inequalities [15,16,17,18,19], such as the Bessel-Legendre inequality (BLI) [20,21] and the Jacobi-Bessel inequality (JLI) [22]. While BLI offers analytical solutions for constant delay systems, its applicability to time-varying systems is limited due to a reliance on estimated boundaries [23,24,25]. In contrast, affine bessel-Legendre inequality (ABLI) addresses time-varying delay amplitude but suffers from conservatism due to incomplete vectors [26]. To overcome this, generalized free matrix-based integral inequality (GFMBII) was introduced to complement ABLI; however, it still lacked the full incorporation of delay amplitude-dependent slack variables [27]. Although the delay derivative-dependent integral inequality was first introduced in [28], further investigation is warranted as the decision matrices are fixed, thus limiting their flexibility and utilization. This highlights the need for continued research in this area to fully leverage the potential benefits of such integral inequalities.
With the above analysis, this paper focuses on investigating the stability of time-varying delayed NNs. Two kinds of slack matrices with two tunable parameters, which are dependent on both a delay amplitude and a delay derivative, are proposed. These advancements culminate in the formulation of a compound-matrix-based integral inequality (CPBII). By utilizing CPBII, a stability condition tailored for time-varying delayed NNs is developed. Compared to existing literature such as [6,17,19,27], the most significant contribution of this paper is the successful incorporation of both delay amplitude and derivative information into the inequality with the help of a couple of convex parameters. This innovative approach enhances the robustness and accuracy of the analysis. The feasibility of the proposed criterion is demonstrated through a numerical example.
Notation: In this paper, Rn represents the n-dimensional Euclidean space; N represents the nature number; He[X]represents X+XT; Co{⋯} represents a set of points; col[X,Y] represents [XT,YT]T; diag{...} represents a block diagonal matrix; and XT represents the transposition of X.
2.
Problem formation and preliminaries
Let's take the NNs characterized by a time-dependent delays, as depicted by the following equation:
˙u(t)=−Au(t)+F0g(F2u(t))+F1g(F2u(t−τ(t)),
(2.1)
where u(t)=col[u1(t),u2(t),u3(t),u4(t),u5(t),...,un(t)]∈Rn represents the state vector, and g(F2u(t))=col[g1(F21u(t)),g2(F22u(t)),...,gn(F2nu(t))] is the activation function. A=diag{a1,a2,..an} is a positive definite diagonal matrix, and F0, F1, F2 are the appropriately dimensional constant matrices. The delay amplitude τ(t) and the delay derivative ˙τ(t) are bounded by constants d and ρ, respectively, satisfying the following
0≤τ(t)≤d,−ρ≤˙τ(t)≤ρ.
(2.2)
The activation functions gi satisfy gi(0)=0 for all i. Define H− and H+ as diagonal matrices with known constants h−i and h+i, respectively, which can be either positive, negative, or zero. Similarly, let G− and G+ be positive definite diagonal matrices. Given these definitions and m,m1,m2∈R, we can deduce the following inequalities from the previously mentioned activation function properties:
h−i≤fi(a)−fi(b)a−b≤h+i,a≠b.
(2.3)
We can directly acheive the following inequalities from (2.3) with m, m1, m2∈R:
In the existing body of work, such as [18,26,27], the final LMIs often incorporate information about the delay derivative, which is typically derived from the derivatives of the LKFs. Despite this, there has been a noticeable absence of integral inequalities that directly pertain to the delay derivative in the context of time-varying delayed NNs. To address this deficiency, we introduce a novel approach in the form of a CPBII, which is outlined below.
Lemma 1. For any continuously differentiable function u:[−d,0]→Rn, the subsequent inequality is valid for any given parameters γ1 and γ2, R>0, any vector η, and slack variables M and N:
By combining this inequality with the one from the previous lemma, (2.5) is derived. This completes the proof.
Remark 1. In Lemma 1, we present a unique integral inequality, termed as CPBII, which amalgamates slack matrices that are dependent on both the delay amplitude and the delay derivative. This is a pioneering approach in the literature where the delay derivative is factored in [28]. The advantages of CPBII are manifold:
● Through the integration of slack matrices that are influenced by both the delay amplitude and the derivative, the successfully forms a link between vectors related to the system states, the delay amplitude, and the delay derivative. This approach facilitates the retrieval of more interconnected data compared to DPTF, ABLI, and GFMBII, all without the need for extra decision variables.
● The inclusion of parameters γ1 and γ2 aid in circumventing certain incomplete terms. For example, when γ1=0 and γ2=0, the last term dτϑT1ΓdT˜RΓdϑ1+τϑT2ΓhT˜RΓhϑ2 is eliminated. Similarly, when γ1=0, γ2=1, and ˙d=ρ, the first term τMT˜RM+dτNT˜RN disappears. Furthermore, this parameter enhances the systems adaptability.
Remark 2. As noted from [10], there exist two strategies to mitigate the conservatism. The first involves striving to get as close to the left side of the inequality as possible, while the second entails introducing an adequate number of cross terms to ensure sufficient system information within the final conditions. Therefore, this paper opts for the second strategy, albeit at the expense of the first to a certain degree. The most significant challenge resides in the assignment of values to ϵ1 and ϵ2. If these values are not assigned appropriately, it becomes evidently impossible to counterbalance the discrepancy caused by the first strategy.
3.
Main results
In the study [27], it is noted that the delay derivative ˙τ is present in Γ0(τ,˙τ), Γ1(τ,˙τ), and Γ3(˙τ), where it is exclusively coupled with positive definite matrices such as P1h, P2h, Q1, and Q2. Interestingly, the slack matrices that are dependent on the delay derivative are not considered. To exploit the information offered by the delay derivative to its fullest extent, a stability condition for system (2.1) is formulated based on CPBII.
Theorem 1. Provided that there exist positive-definite symmetric matrices P, Q1, Q2, and R, and matrices H1, H2, M, and N, in conjunction with specific scalars γ1, γ2, d, and ρ, that meet the subsequent inequalities, it can be concluded that system (2.1) exhibits asymptotic stability:
where the terms Ψ1(τ,˙τ), Ψ2(τ,˙τ), and Ψ4(˙τ) remain consistent with those outlined in previous part. Upon differentiating v1i(t) and v2i(t) in η1h(t) and η2h(t), the following results are obtained:
Here a2(˙τ) has been defined in Theorem 1, and a1, a0 are the appropriate dimensional matrices. By the Schur complement lemma, the inequalities are equivalent to ˜Ψ(0,˙τ)<0, ˜Ψ(d,˙τ)<0, and −d2a2(˙τ)+˜Ψ(0,˙τ)<0. These correspond to the three conditions f(0)<0, f(d)<0, and −d2a2+f(0)<0 in Lemma 2 of Ref. [4]. Therefore, ˉΨ(τ,˙τ)<0 is ensured for any τ∈[0,d].
Additionally, ˉΨ(τ,˙τ) is affine with respect to ˙τ. Therefore, ˉΨ(τ,˙τ)<0 is ensured for any ˙τ∈[ρ1,ρ2] by ˉΨ(τ,ρ1)<0 and ˉΨ(τ,ρ2)<0.
In conclusion, based on Theorem 1, for a small positive scalar ϵ, it follows that ˙V(t)≤−ϵ||u(t)||<0 for u(t)≠0. This implies the asymptotical stability of NNs (2.1).
Remark 3. In earlier research such as [6,27,28], the inclusion of the delay derivative ˙τ has predominantly been dependent on the derivative of the LKFs. Yet, the full integration of the delay derivative ˙τ remains unaccomplished. In contrast to the LKFs, CPBII effectively incorporates the delay derivative ˙τ. This methodology facilitates the concurrent introduction of ˙τ, the delay amplitude ρ, slack matrices M, N, and the augmented vector ηh(t) along with positive definite matrices. As a result, the system information can be efficiently interconnected. Additionally, the inclusion of parameters γ1 and γ2 assists in circumventing zero terms, thereby enabling the extraction of more coupling information. For example, consider ρd(ρ−γ1−γ2˙τ) in (3.1). If γ1=0, γ2=1, and ˙d=ρ are set, the interrelation between ˙τ and ηh(t) is connected via N would vanish. Furthermore, permitting any parameters where γ1≤0 and γ2≤1 (since ρ−γ1−γ2˙d≥0) enhances the adaptability of Theorem 1, in which demonstrates reduced conservatism.
4.
Numerical example
In this section, we will demonstrate the effectiveness of the proposed stability condition.
We scrutinize NN that is characterized by the structure (2.1), where F2=I. The parameters employed in this analysis are derived from the study [27]:
As depicted in Table 1, with the same LKF in [27], the maximum allowable upper bound of the delay amplitude is presented for a range of ρ values. A clear observation from the table is that the conservatism in the results derived in this study is less pronounced compared to previous studies [6,17,19,27]. A comparative analysis between Theorem 1 of this paper and the results of [27] underscores the efficacy of CPBII in mitigating conservatism. Moreover, it is discerned that the parameters γ1 and γ2 augment the adaptability of the stability condition. It should be highlighted that the values of γ1 and γ2 in Table 1 are random under γ1≤0 and γ2≤1. Thus, the maximum allowable upper bounds of d may be larger by choosing more suitable values, which deserves a further study in the future.
Table 1.
The maximum allowable upper bounds of d for different ρ.
On the other hand, by setting u(t)=[0.50.3−0.3−0.5]T, τ(t)=4.7601+0.1sin(t), g(t)=0.1tanh(u), it can be seen from Figure 1 that the state response is stable, which shows the effectiveness of proposed method.
Figure 1.
The state responses of system (2.1) under τ(t)=4.7601+0.1sin(t).
This study addresseed stability analysis of neural networks with time-varying delays. We introduced to CPBII to incorporate delay derivatives into integral inequalities. Then, a novel stability criterion for such neural networks was derived using CPBII. Notably, CPBII encompassed all augmented vectors and their derivatives from the LKF, thus facilitating comprehensive coupling with the delay amplitude and the delayderivative via slack matrices and tunable parameters. This integration led to less conservative outcomes. The effectiveness of our approach was demonstrated through numerical examples.
Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The author declares that there is no conflict of interest.
References
[1]
T. Aktosun, F. Demontis, C. Der Mee, Exact solutions to the sine-Gordon equation, J. Math. Phys., 51 (2010), 123521. https://doi.org/10.1063/1.3520596 doi: 10.1063/1.3520596
[2]
J. Argyris, M. Haase, J. Heinrich, Finite element approximation to two-dimensional sine-Gordon solitons, Comput. Method. Appl. M., 86 (1991), 1–26. https://doi.org/10.1016/0045-7825(91)90136-T doi: 10.1016/0045-7825(91)90136-T
[3]
Z. Asgari, S. M. Hosseini, Numerical solution of two-dimensional sine-Gordon and MBE models using Fourier spectral and high order explicit time stepping methods, Comput. Phys. Commun., 184 (2013), 565–572. https://doi.org/10.1016/j.cpc.2012.10.009 doi: 10.1016/j.cpc.2012.10.009
[4]
R. Sassaman, A. Biswas, Soliton perturbation theory for phi-four model and nonlinear Klein-Gordon equations, Commun. Nonlinear Sci., 14 (2009), 3239–3249. https://doi.org/10.1016/j.cnsns.2008.12.020 doi: 10.1016/j.cnsns.2008.12.020
[5]
L. Brugnano, F. Iavernaro, Line integral methods for conservative problems, 1 Ed., New York: Chapman and Hall/CRC, 2016. https://doi.org/10.1201/b19319
[6]
L. Brugnano, F. Iavernaro, D. Trigiante, Hamiltonian boundary value methods (energy preserving discrete line integral methods), JNAIAM, 5 (2010), 17–37.
[7]
L. Brugnano, G. Caccia, F. Iavernaro, Energy conservation issues in the numerical solution of the semilinear wave equation, Appl. Math. Comput., 270 (2015), 842–870. https://doi.org/10.1016/j.amc.2015.08.078 doi: 10.1016/j.amc.2015.08.078
[8]
L. Brugnano, F. Iavernaro, D. Trigiante, Energy and quadratic invariants preserving integrators based upon gauss collocation formulae, SIAM J. Numer. Anal., 50 (2012), 2897–2916. https://doi.org/10.1137/110856617 doi: 10.1137/110856617
[9]
H. Y. Cao, Z. Z. Sun, G. H. Gao, A three‐level linearized finite difference scheme for the Camassa-Holm equation, Numer. Meth. Part. D. E., 30 (2014), 451–471. https://doi.org/10.1002/num.21819 doi: 10.1002/num.21819
[10]
J. Chabassier, P. Joly, Energy preserving schemes for nonlinear Hamiltonian systems of wave equations: application to the vibrating piano string, Comput. Method. Appl. M., 199 (2010), 2779–2795. https://doi.org/10.1016/j.cma.2010.04.013 doi: 10.1016/j.cma.2010.04.013
[11]
M. Dehghan, A. Shokri, A numerical method for one-dimensional nonlinear Sine-Gordon equation using collocation and radial basis functions, Numer. Meth. Part. D. E., 24 (2008), 687–698. https://doi.org/10.1002/num.20289 doi: 10.1002/num.20289
[12]
D. Duncan, Sympletic finite difference approximations of the nonlinear Klein-Gordon equation, SIAM J. Numer. Anal., 34 (1997), 1742–1760. https://doi.org/10.1137/S0036142993243106 doi: 10.1137/S0036142993243106
[13]
Y. Gong, Q. Wang, Y. Wang, J. Cai, A conservative Fourier pseudo-spectral method for the nonlinear Schrödinger equation, J. Comput. Phys., 328 (2017), 354–370. https://doi.org/10.1016/j.jcp.2016.10.022 doi: 10.1016/j.jcp.2016.10.022
[14]
P. Guyenne, A. Kairzhan, C. Sulem, Hamiltonian Dysthe equation for three-dimensional deep-water gravity waves, Multiscale Model. Sim., 20 (2022), 349–378. https://doi.org/10.1137/21m1432788 doi: 10.1137/21m1432788
[15]
E. Hairer, Energy-preserving variant of collocation methods, JNAIAM, 5 (2010), 73–84.
[16]
E. Hairer, G. Wanner, C. Lubich, Geometric numerical integration: structure-preserving algorithms for ordinary differential equations, 2 Eds., Berlin: Springer, 2006. https://doi.org/10.1007/3-540-30666-8
[17]
B. Hou, D. Liang, Energy-preserving time high-order AVF compact finite difference schemes for nonlinear wave equations with variable coefficients, J. Comput. Phys., 421 (2020), 109738. https://doi.org/10.1016/j.jcp.2020.109738 doi: 10.1016/j.jcp.2020.109738
[18]
B. Hou, D. Liang, The energy-preserving time high-order AVF compact finite difference scheme for nonlinear wave equations in two dimensions, Appl. Numer. Math., 170 (2021), 298–320. https://doi.org/10.1016/j.apnum.2021.07.026 doi: 10.1016/j.apnum.2021.07.026
[19]
D. Kaya, S. El-Sayed, A numerical solution of the Klein-Gordon equation and convergence of the decomposition method, Appl. Math. Comput., 156 (2004), 341–353. https://doi.org/10.1016/j.amc.2003.07.014 doi: 10.1016/j.amc.2003.07.014
[20]
J. Li, Z. Sun, X. Zhao, A three level linearized compact difference scheme for the Cahn-Hilliard equation, Sci. China Math., 55 (2012), 805–826. https://doi.org/10.1007/s11425-011-4290-x doi: 10.1007/s11425-011-4290-x
[21]
S. Li, L. Vu-Quoc, Finite difference calculus invariant structure of a class of algorithms for the nonlinear Klein-Gordon equation, SIAM J. Numer. Anal., 32 (1995), 1839–1875. https://doi.org/10.1137/0732083 doi: 10.1137/0732083
[22]
C. Liu, J. Li, Z. Yang, Y. Tang, K. Liu, Two high-order energy-preserving and symmetric Gauss collocation integrators for solving the hyperbolic Hamiltonian systems, Math. Comput. Sim., 205 (2023), 19–32. https://doi.org/10.1016/j.matcom.2022.09.016 doi: 10.1016/j.matcom.2022.09.016
[23]
C. Liu, K. Liu, A fourth-order energy-preserving and symmetric average vector field integrator with low regularity assumption, J. Comput. Appl. Math., 439 (2024), 115605. https://doi.org/10.1016/j.cam.2023.115605 doi: 10.1016/j.cam.2023.115605
[24]
C. Liu, A. Iserles, X. Wu, Symmetric and arbitrarily high-order Birkhoff-Hermite time integrators and their long-time behaviour for solving nonlinear Klein-Gordon equations, J. Comput. Phys., 356 (2018), 1–30. https://doi.org/10.1016/j.jcp.2017.10.057 doi: 10.1016/j.jcp.2017.10.057
[25]
W. J. Liu, J. B. Sun, B. Y. Wu, Space-time spectral method for the two-dimensional generalized sine-Gordon equation, J. Math. Anal. Appl., 427 (2015), 787–804. https://doi.org/10.1016/j.jmaa.2015.02.057 doi: 10.1016/j.jmaa.2015.02.057
[26]
C. Liu, Y. Tang, J. Yu, Y. Fang, High-order symmetric and energy-preserving collocation integrators for the second-order Hamiltonian system, J. Math. Chem., 62 (2024), 330–355. https://doi.org/10.1007/s10910-023-01536-x doi: 10.1007/s10910-023-01536-x
[27]
C. Y. Liu, X. Y. Wu, Continuous trigonometric collocation polynomial approximations with geometric and superconvergence analysis for efficiently solving semi-linear highly oscillatory hyperbolic systems, Calcolo, 58 (2021), 6. https://doi.org/10.1007/s10092-020-00394-2 doi: 10.1007/s10092-020-00394-2
[28]
C. Liu, X. Wu, Nonlinear stability and convergence of ERKN integrators for solving nonlinear multi-frequency highly oscillatory second-order ODEs with applications to semi-linear wave equations, Appl. Numer. Math., 153 (2020), 352–380. https://doi.org/10.1016/j.apnum.2020.02.020 doi: 10.1016/j.apnum.2020.02.020
[29]
C. Y. Liu, X. Y. Wu, Arbitrarily high-order time-stepping schemes based on the operator spectrum theory for high-dimensional nonlinear Klein-Gordon equations, J. Comput. Phys., 340 (2017), 243–275. https://doi.org/10.1016/j.jcp.2017.03.038 doi: 10.1016/j.jcp.2017.03.038
[30]
C. Y. Liu, X. Y. Wu, The boundness of the operator-valued functions for multidimensional nonlinear wave equations with applications, Appl. Math. Lett., 74 (2017), 60–67. https://doi.org/10.1016/j.aml.2017.04.026 doi: 10.1016/j.aml.2017.04.026
[31]
C. Y. Liu, X. Y. Wu, W. Shi, New energy-preserving algorithms for nonlinear Hamiltonian wave equation equipped with Neumann boundary conditions, Appl. Math. Comput., 339 (2018), 588–606. https://doi.org/10.1016/j.amc.2018.07.059 doi: 10.1016/j.amc.2018.07.059
[32]
C. Liu, W. Shi, X. Wu, Numerical analysis of an energy-conservation scheme for two-dimensional hamiltonian wave equations with Neumann boundary conditions, Int. J. Numer. Anal. Mod., 16 (2019), 319–339. https://doi.org/2019-IJNAM-12806 doi: 2019-IJNAM-12806
[33]
W. Shi, K. Liu, X. Wu, C. Liu, An energy-preserving algorithm for nonlinear Hamiltonian wave equations with Neumann boundary conditions, Calcolo, 54 (2017), 1379–1402. https://doi.org/10.1007/s10092-017-0232-5 doi: 10.1007/s10092-017-0232-5
[34]
W. Tang, J. Zhang, Symplecticity-preserving continuous-stage Runge-Kutta-Nyström methods, Appl. Math. Comput., 323 (2018), 204–219. https://doi.org/10.1016/j.amc.2017.11.054 doi: 10.1016/j.amc.2017.11.054
[35]
W. Tang, J. Zhang, Symmetric integrators based on continuous-stage Runge-Kutta-Nyström methods for reversible systems, Appl. Math. Comput., 361 (2019), 1–12. https://doi.org/10.1016/j.amc.2019.05.013 doi: 10.1016/j.amc.2019.05.013
[36]
W. Tang, Y. Sun, J. Zhang, High order symplectic integrators based on continuous-stage Runge-Kutta-Nyström methods, Appl. Math. Comput., 361 (2019), 670–679. https://doi.org/10.1016/j.amc.2019.06.031 doi: 10.1016/j.amc.2019.06.031
B. Wang, X. Y. Wu, The formulation and analysis of energy-preserving schemes for solving high-dimensional nonlinear Klein-Gordon equations, IMA J. Numer. Anal., 39 (2019), 2016–2044. https://doi.org/10.1093/imanum/dry047 doi: 10.1093/imanum/dry047
[39]
A. Wazwaz, Exact solutions for the generalized sine-Gordon and the generalized sinh-Gordon equations, Chaos Soliton. Fract., 28 (2006), 127–135. https://doi.org/10.1016/j.chaos.2005.05.017 doi: 10.1016/j.chaos.2005.05.017
[40]
L. Xu, P. Guyenne, Numerical simulation of three-dimensional nonlinear water waves, J. Comput. Phys., 228 (2009), 8446–8466. https://doi.org/10.1016/j.jcp.2009.08.015 doi: 10.1016/j.jcp.2009.08.015
This article has been cited by:
1.
C. Maharajan, C. Sowmiya, Changjin Xu,
Lagrange stability of inertial type neural networks: A Lyapunov-Krasovskii functional approach,
2025,
16,
1868-6478,
10.1007/s12530-025-09681-1
Figure 1. Results for Problem 1 with Mx=My=64. (a) the log-log plot of the global error against different steps h. (b) the log-log plot of the global error against CPU time. (c) the log plot of the relative energy error against integrate time with H(0)=123.9352528130722
Figure 2. Snapshots of the numerical solution of the proposed EP3-FP scheme for solving Problem 2 at different times with the spatial mesh grid scales M=Mx=My=256 and time step size Δt=0.01
Figure 3. Results for Problem 2 with Mx=My=64. (a) the log-log plot of the posterior error against different steps h. (b) the log-log plot of the posterior error against CPU time. (c) the log plot of the relative energy error against integrate time with H(0)=0.377193865036316