Special Issue

## New Trends in Numerical Simulation and Data Analysis

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 9057387 | https://doi.org/10.1155/2020/9057387

Pengfei Guo, Ariunkhishig Boldbaatar, Dutao Yi, Pengxiang Dai, "Numerical Solution of Sine-Gordon Equation with the Local Kriging Meshless Method", Mathematical Problems in Engineering, vol. 2020, Article ID 9057387, 11 pages, 2020. https://doi.org/10.1155/2020/9057387

# Numerical Solution of Sine-Gordon Equation with the Local Kriging Meshless Method

Accepted13 Aug 2020
Published08 Sep 2020

#### Abstract

This paper develops a local Kriging meshless solution to the nonlinear 2 + 1-dimensional sine-Gordon equation. The meshless shape function is constructed by Kriging interpolation method to have Kronecker delta function property for the two-dimensional field function, which leads to convenient implementation of imposing essential boundary conditions. Based on the local Petrov–Galerkin formulation and the center difference method for time discretization, a system of nonlinear discrete equations is obtained. The numerical examples are presented and the numerical solutions are found to be in good agreement with the results in the literature to validate the ability of the present meshless method to handle the 2 + 1-dimensional sine-Gordon equation related problems.

#### 1. Introduction

The nonlinear sine-Gordon equation (SGE), a type of hyperbolic partial differential equation, is often used to describe and simulate the physical phenomena in a variety of fields of engineering and science, such as nonlinear waves, propagation of fluxons and dislocation of metals . Because sine-Gordon equation leads to different types of soliton solutions, SGE has been receiving an enormous amount of attention. Soliton solution travels without experiencing any deformation through the medium, even when it collides with another soliton. The solitons, identified in many wave and particle systems, are of importance in the theory of nonlinear differential equations.

As one of the crucial equations in nonlinear science, the sine-Gordon equation has been constantly investigated and solved numerically and analytically in recent years . For the one-dimensional sine-Gordon equation, Wazwaz  adopted the tanh method to find the traveling wave solutions. Johnson et al.  provided three exact solutions to deal with the two-dimensional SGE. Zhong and Belić analytically established special two-soliton solution to solve the generalized SGE based on the Hirota bilinear method and self-similar method . For more related research about exactly solving the sine-Gordon equations, the interested readers may refer to . Despite the exact solutions or analytical analyses of the SGEs that can offer deeper understanding and insights of the associated physical phenomena, numerical methods are still indispensable in pursuit those unknown soliton solutions to the SGEs. Dehghan et al. used the collocation method based on the radial basis function to obtain the numerical solutions to the one-dimensional nonlinear sine-Gordon equation . Besides, the difference approach , the boundary integral equation scheme , the homotopy-perturbation method , the modified Adomain decomposition method , and the multilevel augmentation method  have also been applied to solve the one-dimensional SGEs. For dealing with the higher dimensional problems, Argyris et al.  proposed a finite element algorithm to obtain soliton-type solutions. Djidjeli et al.  developed an explicit numerical approach to solve a damped sine-Gordon equation in two space variables. Sheng et al. presented a split cosine scheme , and Bratsos used the method of lines  for seeking the solitary solutions of the two-dimensional sine-Gordon equation.

During recent decades, meshless methods have also been developed for solving the nonlinear sine-Gordon equations, which include the mesh-free kp-Ritz method , the radius basis function (RBF) method [7, 16], the radial point interpolation method (RPIM) , the meshless local Petrov-Galerkin (MLPG) method , and the well-posed moving least-squares (WP-MLS) approximation . Meshless methods are a type of numerical schemes by which the problem domain is represented by a set of scattered nodes, and the mesh or discretization can thus be avoided. They have been usually used to effectively handle those physical and engineering problems or to solve the partial differential equations that might pose challenges of the traditional numerical approaches, such as the boundary element method  and finite element method . A variety of meshless methods have been developed and then applied in scientific and engineering modelling based on different combinations of certain formulation schemes, such as the element-free Galerkin formulation , and the approximation/interpolation techniques, such as moving least-squares (MLS) [31, 32]. Based on MLS, the element-free Galerkin method has been viewed as one of the predominant meshless methods. Cheng et al. proposed the improved moving least-squares (IMLS) approximation to develop the improved element-free Galerkin (IEFG) method  and the interpolating element-free Galerkin method . By introducing the dimension splitting method into the reproducing kernel particle method, a new meshless method was developed to solve three-dimensional potential problems .

In this study, we present the local Kriging meshless method  to seek the numerical solutions of the 2 + 1-dimensional nonlinear sine-Gordon equation. The present meshless approach is based on the Kriging interpolation technique and the local Petrov–Galerkin formulation, which leads to the shape functions having the Kronecker delta property and thus to a convenient implementation required for imposing the essential boundary conditions, like in the traditional finite element analysis. Combining with the Petrov–Galerkin formulation over the local subdomains, the meshless method is featured by avoiding background cells for the numerical integration of the global weak form. Based on the center difference scheme for making time discretization, we have a global system of nonlinear algebraic equations. Numerical tests are performed to validate the convergence and accuracy of the present meshless method for solving the 2 + 1-dimensional sine-Gordon equation. The numerical results are found to agree well with the existing results reported in the published literature.

#### 2. The Local Kriging Meshless Formulation for the 2D Sine-Gordon Equation

##### 2.1. Governing Equation

In the present study, we seek to acquire a meshless approximation to the damped 2 + 1-dimensional sine-Gordon equation ,within the problem domain of for . In equation (1), the parameter is the damping factor for the dissipative term. When , this equation will reduce to an undamped sine-Gordon equation. The functions represent Josephson current density.

The boundary conditions linked to the sine-Gordon equation (equation (1)) enforce a zero normal gradient on the boundary of as follows:

At time , we assume the initial conditions to be in the form ofwhere the functions, and , are specified wave modes and velocity.

##### 2.2. Kriging Interpolation

The problem domain with the boundary can be represented by a set of scattered nodes, the of which is assumed to be located at , where is the total number of nodes. The field function, such as , is approximated/interpolated via the nodal values . Through a weighted linear combination, the field value at point can be evaluated approximately aswhere stands for the nodal value at and represents the support domain associated with the point at , which contains nodes that exert influence on the point at . The weight coefficient will be determined according to the Kriging interpolation method and the Kriging-based shape function can thus be written by

The matrices and are expressed, respectively, aswhere is the unit matrix of order , and is a column vector that contains monomial basis functions, which are given, as an example, in case of two-dimensional linear and quadratic basis functions, as follows:

In terms of the basis function at nodes in the support domain, the matrix is then written as

The matrix and vector are, respectively, expressed aseach entry of which, , is evaluated based on the semivariogram model that determines the relation of any pair of nodes at and as follows:where stands for an expected value of a random function. The semivariogram model in equation (11) is independent of the locations of the two nodes but only depends on their lag vector .

In this numerical investigation, we adopt the Gauss model to construct the shape functions via the Kriging interpolation scheme , which is given as

In equation (12), is the Euclidean distance of any point at and the node at , or the distance of any pair of nodes at and . Thus, for equation (10) and for equation (9). The two factors and are the sill and range that are empirically determined. In the present study, we set to 1 because its value has little influence on the shape functions, whereas is evaluated according to the influence domain size linked to nodesin which is a specified factor . In the local Kriging meshless method, the interpolation domain for any point at is constructed in terms of the influence domain size of that is written aswhere is a dimensional coefficient and represents the average spacing among nodes, which is determined by considering the balance of computational costs and accuracy [49, 50].

In the problem domain , we have the derivatives of the constructed shape functions with respect to and based on equation (5), and their first- and second-order partial derivatives are given as follows:in which and represent and , respectively.

##### 2.3. Local Weak Form of the Governing Equation

In the present study, we apply the local Petrov–Galerkin formulation to construct the weak form of the governing equation over the pre-established local subdomains associated with the nodes in the global problem domain (see Figure 1). The positive definite test functions that are assigned to each node need to be predefined within a test function region . Over one local subdomain centered at node , the local weighted residual weak form of equation (1) is then expressed aswhere we specify a cubic spine function  as the test function in this study, which is defined asin which .

By decomposing the subdomain boundary into three sub-boundaries and applying the Gauss divergence theorem, equation (16) is re-expressed asin which and stand for the outward normal direction cosines of the boundary of the local subdomain and the boundary contains three disjoint parts .

According to the Kriging interpolation method and the constructed shape function (equation (5)) presented in the previous section, the field function to be determined can be approximated in terms of nodal values at the nodes that are included in the support domain associated with the point at as follows:

The time and space derivatives of the filed function are expressed in the following form ofwhere

Substituting the approximate field function in equation (19) for in the local weak form of the governing equation (equation (18)), we obtain the nodal system equation for the node in the matrix form:in which

In order to construct the global system equation, the assembling process similar to the finite difference method has been used for all nodes scattered in the problem domain and it leads to

Apparently, equation (26) is a system of nonlinear algebraic equations, and an iterative technique is then applied to cope with the nonlinearity. In equation (25), the first term on the right-hand side needs to be evaluated according to the latest available approximation of the field function .

##### 2.4. The Time Discretization and Iterative Procedure

In the present study, the time derivatives of the global system equations (equation (25)) is handled by the center difference technique for making time discretization, and the global system equation (equation (25)) is transformed into

The above equation can be rewritten as

When , is determined according to the specified initial condition (equation (3)), that is, , and for the next time level, can be obtained as . For the purpose of handling the nonlinearity of equation (28) using a predictor-corrector technique [4, 26], at a time level of , we need the field function value at the last time level to calculate the term in the equation, and the system of equations is then solved linearly with . Afterwards, we recalculate . Equation (28) can be solved with the latest available to obtain the field function values at the time level of . Performing repeatedly this iterative procedure and setting , we then achieve the nonlinear solutions to equation (28) when if the solutions can converge to a prescribed number . In the present study, the condition for terminating the iteration is given as follows:

We can set and proceed to the next time level, in case the above condition is met.

By conducting the above iterative procedure to solve the equation (equation (28)) till the desired time level, we achieve the numerical solutions to the 2 + 1-dimensional nonlinear sine-Gordon equation.

#### 3. Numerical Examples

In this section, we apply the local Kriging meshless method to several numerical examples regarding two-dimensional line and ring solitons.

##### 3.1. Example 1: Superposition of Two Orthogonal Line Solitons (Undamped and Damped SG Equation) [4, 26, 51]

In the case of an undamped 2D sine-Gordon equation with , we obtain the superposition of two orthogonal line solitons under the given initial conditions,and over a square domain of . The numerical results of are shown in Figure 2 for , and , respectively. For better depicting the breakup of two orthogonal line solitons, we illustrate the contours of the superposition of the two orthogonal line solitons in Figure 3 that move separately from each other along the direction over time. All numerical solutions are in good agreement with those available results in [4, 26, 51].

In order to investigate the effects of the dissipative term on the behavior, the same numerical example is used with and , respectively. Figures 3(d), 4(a), and 4(b) illustrate the contours of the line solitons with varying damping factors when , in which exist apparent propagation delay of the line solitons due to the presence of the dissipative term in the SG equation.

In order to quantitatively compare the numerical results with the literature, we adopt the energy for an undamped sine-Gordon equation , which is conserved and defined as the following form :

The integration of the above equation over the problem domain is calculated according to the trapezoidal rule. Table 1 exhibits the conservation of the energy given by the present method, which agrees well with the literature [4, 26].

  – 175.3027 175.3107 175.3289 175.3743 175.3921 175.4131 175.7865  175.5745 175.2927 175.2761 175.3188 175.3637 175.3899 175.4045 175.6831 (present) 175.7503 175.7515 175.8155 175.9262 176.0511 176.2110 176.3888 176.5849

The convergence analysis of the present method is conducted based on the energy, , shown in Figure 5.

##### 3.2. Example 2: Circular Ring Soliton (Damped) [4, 22]

For the numerical case, in which and , we numerically obtain the circular ring solitons with the specified initial conditions,and in the problem domain of . Figures 6 and 7 show the numerical solutions and their contours in term of for and , respectively. Comparing the surfaces in Figure 6, a complete agreement can be observed with those published in  for the damped SG equation.

##### 3.3. Perturbation of a Static Line Soliton [4, 52]

The perturbation of a line soliton is studied in the numerical example with and and the specified initial conditions,over the problem domain . In Figure 8 which illustrates the symmetric perturbation of single line soliton in terms of for , two symmetric dents firstly move towards each other, then collide at certain time, and proceed to move away from each other.

#### 4. Concluding Remarks

In this study, the local Kriging meshless method has been extended to the 2 + 1-dimensional nonlinear sine-Gordon equation. The Kriging interpolation technique is used to approximate the two-dimensional field function, which leads to convenient implementation of imposing essential boundary conditions. A system of nonlinear discrete equations can be established based on the adoption of the local Petrov–Galerkin formulation and the center difference method for time discretization. The nonlinear algebraic equations are solved by applying the iterative technique and a predictor-corrector scheme. The numerical results are in good agreement with those available in the literature. The present meshless method is thus a potential alternative to other numerical methods, such as the finite element method and finite difference method, for dealing with a 2 + 1-dimensional nonlinear sine-Gordon equation.

#### Data Availability

All results have been obtained by conducting the numerical procedure and the ideas can be shared for the researchers.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was supported by Hainan Provincial Natural Science Foundation of China (Grant no. 119MS039).

1. R. K. Dodd, J. C. Eilbeck, J. D. Gibbon, and H. C. Morris, Solitons and Nonlinear Wave Equations, Academic, London, UK, 1982.
2. A. Barone, F. Esposito, C. J. Magee, and A. C. Scott, “Theory and applications of the sine-Gordon equation,” La Rivista del Nuovo Cimento, vol. 1, no. 2, pp. 227–267, 1971. View at: Publisher Site | Google Scholar
3. A. M. Wazwaz, Partial Differential Equations and Solitary Waves Theory, Springer, Berlin, Germany, 2009.
4. R. J. Cheng and K. M. Liew, “Analyzing two-dimensional sine-Gordon equation with the mesh-free reproducing kernel particle Ritz method,” Computer Methods in Applied Mechanics and Engineering, vol. 245-246, pp. 132–143, 2012. View at: Publisher Site | Google Scholar
5. S. Maitama and Y. F. Hamza, “An analytical method for solving nonlinear sine-Gordon equation,” Sohag Journal of Mathematics, vol. 7, no. 1, pp. 5–10, 2020. View at: Google Scholar
6. L. Su, “Numerical solution of two-dimensional nonlinear sine-Gordon equation using localized method of approximate particular solutions,” Engineering Analysis with Boundary Elements, vol. 108, pp. 95–107, 2019. View at: Publisher Site | Google Scholar
7. M. Dehghan, M. Abbaszadeh, and A. Mohebbi, “An implicit RBF meshless approach for solving the time fractional nonlinear sine-Gordon and Klein-Gordon equations,” Engineering Analysis with Boundary Elements, vol. 50, pp. 412–434, 2015. View at: Publisher Site | Google Scholar
8. M. Dehghan and A. Ghesmati, “Numerical simulation of two-dimensional sine-Gordon solitons via a local weak meshless technique based on the radial point interpolation method (RPIM),” Computer Physics Communications, vol. 181, no. 4, pp. 772–786, 2010. View at: Publisher Site | Google Scholar
9. A.-M. Wazwaz, “The tanh method: exact solutions of the sine-Gordon and the sinh-Gordon equations,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1196–1210, 2005. View at: Publisher Site | Google Scholar
10. S. Johnson, P. Suarez, and A. Biswas, “New exact solutions for the sine-Gordon equation in 2+1 dimensions,” Computational Mathematics and Mathematical Physics, vol. 52, no. 1, pp. 98–104, 2012. View at: Publisher Site | Google Scholar
11. W.-P. Zhong and M. Belić, “Special two-soliton solution of the generalized Sine-Gordon equation with a variable coefficient,” Applied Mathematics Letters, vol. 38, pp. 122–128, 2014. View at: Publisher Site | Google Scholar
12. M. R. Gao and H. T. Chen, “Hybrid solutions of three functions to the (2+1)-dimensional sine-Gordon equation,” Acta Physica Sinica, vol. 61, no. 22, 2012. View at: Google Scholar
13. A. H. Salas, “Exact solutions of coupled sine-Gordon equations,” Nonlinear Analysis: Real World Applications, vol. 11, no. 5, pp. 3930–3935, 2010. View at: Publisher Site | Google Scholar
14. T. Aktosun, F. Demontis, and C. van der Mee, “Exact solutions to the sine-Gordon equation,” Journal of Mathematical Physics, vol. 51, no. 12, p. 123521, 2010. View at: Publisher Site | Google Scholar
15. W.-X. Chen and J. Lin, “Some new exact solutions of (1+2)-dimensional sine-Gordon equation,” Abstract and Applied Analysis, vol. 2014, Article ID 645456, 8 pages, 2014. View at: Publisher Site | Google Scholar
16. M. Dehghan and A. Shokri, “A numerical method for one-dimensional nonlinear sine-Gordon equation using collocation and radial basis functions,” Numerical Methods for Partial Differential Equations, vol. 24, no. 2, pp. 687–698, 2008. View at: Publisher Site | Google Scholar
17. A. G. Bratsos and E. H. Twizell, “A family of parametric finite-difference methods for the solution of the sine-Gordon equation,” Applied Mathematics and Computation, vol. 93, no. 2-3, pp. 117–137, 1998. View at: Publisher Site | Google Scholar
18. M. Dehghan and D. Mirzaei, “The boundary integral equation approach for numerical solution of the one-dimensional sine-Gordon equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 6, pp. 1405–1415, 2008. View at: Publisher Site | Google Scholar
19. M. S. H. Chowdhury and I. Hashim, “Application of homotopy-perturbation method to Klein-Gordon and sine-Gordon equations,” Chaos, Solitons & Fractals, vol. 39, no. 4, pp. 1928–1935, 2009. View at: Publisher Site | Google Scholar
20. S. S. Ray, “A numerical solution of the coupled sine-Gordon equation using the modified decomposition method,” Applied Mathematics and Computation, vol. 175, no. 2, pp. 1046–1054, 2006. View at: Publisher Site | Google Scholar
21. J. Chen, Z. Chen, and S. Cheng, “Multilevel augmentation methods for solving the sine-Gordon equation,” Journal of Mathematical Analysis and Applications, vol. 375, no. 2, pp. 706–724, 2011. View at: Publisher Site | Google Scholar
22. J. Argyris, M. Haase, and J. C. Heinrich, “Finite element approximation to two-dimensional sine-Gordon solitons,” Computer Methods in Applied Mechanics and Engineering, vol. 86, no. 1, pp. 1–26, 1991. View at: Publisher Site | Google Scholar
23. K. Djidjeli, W. G. Price, and E. H. Twizell, “Numerical solutions of a damped Sine-Gordon equation in two space variables,” Journal of Engineering Mathematics, vol. 29, no. 4, pp. 347–369, 1995. View at: Publisher Site | Google Scholar
24. Q. Sheng, A. Q. M. Khaliq, and D. A. Voss, “Numerical simulation of two-dimensional sine-Gordon solitons via a split cosine scheme,” Mathematics and Computers in Simulation, vol. 68, no. 4, pp. 355–373, 2005. View at: Publisher Site | Google Scholar
25. A. G. Bratsos, “The solution of the two-dimensional sine-Gordon equation using the method of lines,” Journal of Computational and Applied Mathematics, vol. 206, no. 1, pp. 251–277, 2007. View at: Publisher Site | Google Scholar
26. D. Mirzaei and M. Dehghan, “Meshless local Petrov-Galerkin (MLPG) approximation to the two dimensional sine-Gordon equation,” Journal of Computational and Applied Mathematics, vol. 233, no. 10, pp. 2737–2754, 2010. View at: Publisher Site | Google Scholar
27. X. Li, “Meshless numerical analysis of a class of nonlinear generalized Klein-Gordon equations with a well-posed moving least squares approximation,” Applied Mathematical Modelling, vol. 48, pp. 153–182, 2017. View at: Publisher Site | Google Scholar
28. D. E. Beskos, “Boundary element methods in dynamic analysis,” Applied Mechanics Reviews, vol. 40, no. 1, pp. 1–23, 1987. View at: Publisher Site | Google Scholar
29. R. Huang, S. Zheng, Z. Liu, and T. Y. Ng, “Recent advances of the constitutive models of smart materials-hydrogels and shape memory polymers,” International Journal of Applied Mechanics, vol. 12, no. 02, p. 2050014, 2020. View at: Publisher Site | Google Scholar
30. T. Belytschko, Y. Y. Lu, and L. Gu, “Element-free Galerkin methods,” International Journal for Numerical Methods in Engineering, vol. 37, no. 2, pp. 229–256, 1994. View at: Publisher Site | Google Scholar
31. J. Cheng, “Data analysis of the factors influencing the industrial land leasing in Shanghai based on mathematical models,” Mathematical Problems in Engineering, vol. 2020, Article ID 9346863, 2020. View at: Publisher Site | Google Scholar
32. J. Cheng, “Analyzing the factors influencing the choice of the government on leasing different types of land uses: evidence from Shanghai of China,” Land Use Policy, vol. 90, p. 104303, 2020. View at: Publisher Site | Google Scholar
33. Y. M. Cheng and M. J. Chen, “A boundary element-free method for linear elasticity,” Acta Mechanica Sinica, vol. 35, no. 2, pp. 181–186, 2003. View at: Google Scholar
34. B. Wang, Y. Ma, and Y. Cheng, “The improved complex variable element-free Galerkin method for bending problem of thin plate on elastic foundations,” International Journal of Applied Mechanics, vol. 11, no. 10, p. 1950105, 2019. View at: Publisher Site | Google Scholar
35. F. Liu and Y. Cheng, “The improved element-free Galerkin method based on the nonsingular weight functions for inhomogeneous swelling of polymer gels,” International Journal of Applied Mechanics, vol. 10, no. 4, p. 1850047, 2018. View at: Publisher Site | Google Scholar
36. Q. Wu, P. P. Peng, and Y. M. Cheng, “The interpolating element-free Galerkin method for elastic large deformation problems,” Science China Technological Sciences, vol. 63, 2020. View at: Google Scholar
37. Q. Wu, F. B. Liu, and Y. M. Cheng, “The interpolating element-free Galerkin method for three-dimensional elastoplasticity problems,” Engineering Analysis with Boundary Elements, vol. 115, pp. 156–167, 2020. View at: Publisher Site | Google Scholar
38. D. Liu and Y. M. Cheng, “The interpolating element-free Galerkin (IEFG) method for three-dimensional potential problems,” Engineering Analysis with Boundary Elements, vol. 108, pp. 115–123, 2019. View at: Publisher Site | Google Scholar
39. Z. J. Meng, H. Cheng, L. D. Ma, and Y. M. Cheng, “The dimension splitting element-free Galerkin method for 3D transient heat conduction problems,” Science China Physics, Mechanics & Astronomy, vol. 62, no. 4, Article ID 40711, 2019. View at: Publisher Site | Google Scholar
40. L. Ma, Z. Meng, J. Chai, and Y. Cheng, “Analyzing 3D advection-diffusion problems by using the dimension splitting element-free Galerkin method,” Engineering Analysis with Boundary Elements, vol. 111, pp. 167–177, 2020. View at: Publisher Site | Google Scholar
41. Z. J. Meng, H. Cheng, L. D. Ma, and Y. M. Cheng, “The dimension split element-free Galerkin method for three-dimensional potential problems,” Acta Mechanica Sinica, vol. 34, no. 3, pp. 462–474, 2018. View at: Publisher Site | Google Scholar
42. H. Cheng, M. Peng, Y. Cheng, and Z. Meng, “The hybrid complex variable element-free Galerkin method for 3D elasticity problems,” Engineering Structures, vol. 219, p. 110835, 2020. View at: Publisher Site | Google Scholar
43. P. P. Peng, Q. Wu, and Y. M. Cheng, “The dimension splitting reproducing kernel particle method for three‐dimensional potential problems,” International Journal for Numerical Methods in Engineering, vol. 121, no. 1, pp. 146–164, 2020. View at: Publisher Site | Google Scholar
44. P. Zhu, L. W. Zhang, and K. M. Liew, “Geometrically nonlinear thermomechanical analysis of moderately thick functionally graded plates using a local Petrov-Galerkin approach with moving Kriging interpolation,” Composite Structures, vol. 107, pp. 298–314, 2014. View at: Publisher Site | Google Scholar
45. P. Zhu and K. M. Liew, “A local Kriging meshless method for free vibration analysis of functionally graded circular plates in thermal environments,” Procedia Engineering, vol. 31, pp. 1089–1094, 2012. View at: Publisher Site | Google Scholar
46. P. Zhu and K. M. Liew, “Free vibration analysis of moderately thick functionally graded plates by local Kriging meshless method,” Composite Structures, vol. 93, no. 11, pp. 2925–2944, 2011. View at: Publisher Site | Google Scholar
47. L. W. Zhang, P. Zhu, and K. M. Liew, “Thermal buckling of functionally graded plates using a local Kriging meshless method,” Composite Structures, vol. 108, pp. 472–492, 2014. View at: Publisher Site | Google Scholar
48. R. Olea, Geostatistics for Engineers and Earth Scientists, Kluwer Academic, Boston, MA, USA, 1999.
49. Z. Liu, G. Wei, and Z. Wang, “The radial basis reproducing kernel particle method for geometrically nonlinear problem of functionally graded materials,” Applied Mathematical Modelling, vol. 85, pp. 244–272, 2020. View at: Publisher Site | Google Scholar
50. Z. Liu, G. Wei, Z. Wang, and J. Qiao, “The meshfree analysis of geometrically nonlinear problem based on radial basis reproducing kernel particle method,” International Journal of Applied Mechanics, vol. 12, no. 4, p. 2050044, 2020. View at: Publisher Site | Google Scholar
51. P. L. Christiansen and P. S. Lomdahl, “Numerical study of 2+1 dimensional Sine-Gordon solitons,” Physica D: Nonlinear Phenomena, vol. 2, no. 3, pp. 482–494, 1981. View at: Publisher Site | Google Scholar
52. A. G. Bratsos, “An explicit numerical scheme for the Sine-Gordon equation in 2+1 dimensions,” Applied Numerical Analysis & Computational Mathematics, vol. 2, no. 2, pp. 189–211, 2005. View at: Publisher Site | Google Scholar

#### More related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.