Abstract
In the present study, we demonstrate the feasibility to reveal the numerical solution of the multi-dimensional wave equations. A simple semi-analytical meshless method was proposed to obtain the numerical solution of the wave equation with a newly-proposed space-time radial basis function to enhance the numerical stability. The wave equation was discretized into equivalent algebraic equations. By specifying boundary and initial conditions, the wave propagation in a two-dimensional domain can be virtually reconstructed. Our results exhibit that the semi-analytical meshless method is suitable and efficient for solving multi-dimensional wave equations.
1 Introduction
A wide range of physical processes is related to the multi-dimensional wave equation and it has been applied to many practical engineering problems such as underwater sound propagation, motion of vibrating strings and membranes. The ultimately physical model of the wave propagation problem is the time-dependent hyperbolic partial differential equations. For practical engineering problems, only approximate solutions can be obtained by numerical methods. Because of all the problems arose when solving the second-order derivative in time, the numerical solution for this type of equation has been little studied [1].
Several numerical methods have been proposed to get numerical solutions to the multi-dimensional wave equations. The finite-difference-based schemes have gained considerable attention in getting the numerical solutions of different time-dependent partial differential equations [2]. Based on the Houbolt finite difference scheme, the method of the particular solutions and the method of fundamental solutions are combined for the solution of multi-dimensional wave equations [3]. An implicit time difference scheme in conjunction with moving least squares reproducing kernel particle approximation is suggested for time-dependent diffusion-wave equation by Rezvan [4].
The finite-difference-based scheme is another choice. The weak Galerkin finite element method is employed to solve the two-dimensional wave equations [5–7]. A semi-discrete numerical method is introduced for wave equation with the spatial variable discretized by the finite element method [8]. For the other methods, a radial integration boundary element method has been developed for the solution of 2D scalar wave equation. Domain integrals appearing in the integral equations are transformed to the boundary with the help of a modified radial integration method. This technique is accomplished applying two time stepping schemes including Newmark and Houbolt methods [9]. The Laplace transformation is implemented to convert diffusion-wave equation to a series of time-independent nonhomogeneous equations in Laplace domain. A semi-analytical collocation Trefftz scheme is used to obtain the solution of high-order homogeneous equation with boundary-only collocation in Laplace domain [10]. The convolution quadrature method is formulated for the two-dimensional wave equation and the boundary element method is introduced for its spatial discretization [11]. Nevertheless, even when many of these methods obtain satisfactory results, they are based on a two-step solution process. More specifically, the time-dependent problem is treated by the finite difference discretization first which will lead to time-independent equations. Then the other numerical method is employed to solve the time-independent equations.
As is known to all, the radial basis function (RBF) methods perform very well in numerical simulation of mathematical modeling thanks to the features in terms of simple, flexible, and truly meshfree [12–15]. In this paper, we propose a semi-analytical meshfree method with one-step approximation, which is based on newly-proposed RBFs, to analyze the phenomena occurring in the wave propagation.
This paper is briefly organized as follows. Based on the multi-dimensional wave equation, the newly-proposed space-time distance as well as corresponding formulation of RBF is provided in Section 2. Section 3 presents the methodology for multi-dimensional wave equation under initial condition and boundary conditions. Discussions with different wave speed numbers are presented to validate the accuracy and stability of the proposed semi-analytical meshfree method. Section 5 provides some conclusions and future directions.
2 Modeling and methods
The mathematical modeling of wave propagation is one of the earliest well-known multi-dimensional time-dependent wave equationwhere is the space vector with for two-dimensional and for three-dimensional, respectively. is the wave speed. Due to the complexity of practical problems, analytical solutions cannot be obtained for the above wave equation. An alternative is the numerical methods.
Since Eq. 1 is time-dependent, the time-variable is always treated by using the finite difference method, Laplace transformation or the other methods. This will lead to time-nondependent equations. Together with specified boundary and initial conditions, the other numerical methods can be employed to get the approximate solutions of the corresponding mathematical modeling. This procedure is a two-step numerical method. In order to get a one-level numerical method, we propose a semi-analytical meshfree method by using the traditional RBF.
The basic theory of the RBF-based collocation methods lies in that the approximate solution can be written as a linear combination of RBFs. Here, we consider the commonly-used multiquadric RBF as an example [16–18].where is the multiquadric RBF parameter, is the distance between and for two-dimensional cases.
In order to propose a direct meshless method with one-level approximation, we propose a space-time RBF by combination of and as a “space” point . Finally, one can obtain the simple direct radial basis function (DRBF)with for (2 + 1)-dimensional problems andwith for (3 + 1)-dimensional problems. It should be noted that the simple direct radial basis function can be easily extended to higher dimensional cases [19].
3 Methodology for direct meshless method
Based on the definition of space-time radial basis functions, the approximate solution satisfying the Pennes equation Eq. 1 has the formwhere denotes the number of the collocation points. To seek for the unknown coefficients , traditional collocation method can be used, i.e., Eq. 1 is imposed at internal points. The wave equation is reduced to a system of algebraic equations. The system of algebraic equations that corresponds to Eq. 1 is considered when the mixed nonhomogeneous Dirichlet boundary condition and initial conditions are imposed. We can describe the boundary condition as
The initial conditions should be assumed
The following procedure is executed by collocating the boundary conditions Eqs 6–7 at boundary collocation points and initial condition Eq. 8 at initial points , respectively. This procedure yields the following equationswithwhere and are the collocation point numbers on the Dirichlet boundary and Neumann boundary, respectively. is the initial point number which corresponds to time step.
4 Results and discussions
Two examples are considered by implementing the semi-analytical meshless method for the multi-dimensional wave equation. To verify the accuracy and stability of the proposed method in this paper, we consider the mentioned scheme for different values of multiquadric RBF parameter, (the distance between the nodes in space direction and time direction). Numerical solutions obtained from this method are compared with the exact solutions. The root mean square relative error (RMSE) in the following figures is defined aswhere and are the exact and numerical solutions, respectively. NT denotes the total number of testing points.
4.1 Example one
Here, the time-dependent problem in the unit square domain having analytic solutions is considered to validate the capability of the proposed semi-analytical meshless method. The initial boundary conditions can be written aswhile the initial boundary condition is
By using the method of separation of variables, we can obtain the analytical solution for this problem
For the unit wave speed, we consider the effect of RBF parameter to the numerical results. Figure 1 presents the RMSE variation curve for fixed space distance at point . We can see that the RBF parameter performs well in the scope . The optimal choice of RBF parameter in the direct meshless method is similar to the traditional RBF. This is beyond the scope of our investigation, more details related to the optimal choice of RBF parameter can be found in [19–21] and references therein.
FIGURE 1
For fixed RBF parameter, Figure 2 shows the RMSE variation curves for time at three different points. It can be seen that the RMSE are almost smaller than for all the three different points, i.e., these numerical solutions compare well with the analytical solutions. It should be pointed that the method used in [3] requires more fine time step () to achieve the same RMSE while our method only requires .
FIGURE 2
4.2 Example two
In this case, we consider the time-dependent problem in the unit square domain with the analytical solutionwith
We consider the wave speed , the corresponding initial boundary conditions can be written aswhile the initial boundary condition is
Figure 3 presents the RMSE versus RBF parameter for fixed space distance at point . We can see that the RBF parameter performs well in the scope , the quasi-optimal choice is with corresponding . It should be pointed that the method used in [3] requires more fine time step () to achieve RMSE about . Even for the more fine time step , the corresponding RMSE is about . It cannot compete with our method for larger time step . Figure 4 shows the RMSE variation curves for time at three different points.
FIGURE 3
FIGURE 4
5 Conclusion
As presented in the paper, the multiquadric radial basis function is a good base to build an approximate solution of wave propagation problems. The multiquadric radial basis function can be easily generated with the use of any program of Computer Algebra System type executing the symbolic calculations. In the presented paper, the MATLAB program has been used. The semi-analytical meshless method has been presented to find the approximate solution of the multi-dimensional wave equation. Numerical results obtained for the wave equation are similar to/are in accordance with the analytical results.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This work was partially supported by Key Laboratory of Southeast Coast Marine Information Intelligent Perception and Application, Ministry of Natural Resources (KFJJ20220102), the Natural Science Foundation of Xuzhou (No. KC22056) and the Innovation Training Program for College Students (xcx2023119).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
GerardoTHeribertoAJoséAGRicardoRJoséGT. A meshless finite difference scheme applied to the numerical solution of wave equation in highly irregular space regions. Comput Math Appl (2023) 136:25–33. 10.1016/j.camwa.2023.01.035
2.
WuZAlkhalifahTZhangZ. A partial-low-rank method for solving acoustic wave equation. J Comput Phys (2019) 385:1–12. 10.1016/j.jcp.2019.01.054
3.
YoungDLGuMHFanCM. The time-marching method of fundamental solutions for wave equations. Eng Anal Boundary Elem (2009) 33:1411–25. 10.1016/j.enganabound.2009.05.008
4.
RezvanS. A meshless point collocation method for 2-D multi-term time fractional diffusion-wave equation. Numer Algorithms (2017) 74:1145–68. 10.1007/s11075-016-0190-z
5.
NareshK. Supercloseness analysis of a stabilizer-free weak Galerkin finite element method for viscoelastic wave equations with variable coefficients. Adv Comput Math (2023) 49:12. 10.1007/s10444-023-10010-w
6.
HuangYQLiJCLiD. Developing weak Galerkin finite element methods for the wave equation. Numer Methods Partial Differential Equations (2017) 33(3):868–84. 10.1002/num.22127
7.
WangXPGaoFZSunZJ. Weak Galerkin finite element method for viscoelastic wave equations. J Comput Appl Math (2020) 375:112816. 10.1016/j.cam.2020.112816
8.
HanWMSongCHWangFGaoJH. Numerical analysis of the diffusive-viscous wave equation. Comput Math Appl (2021) 102:54–64. 10.1016/j.camwa.2021.10.003
9.
NajarzadehLMovahedianBAzhariM. Numerical solution of scalar wave equation by the modified radial integration boundary element method. Eng Anal Boundary Elem (2019) 105:267–78. 10.1016/j.enganabound.2019.04.027
10.
FuZJYangLWZhuHQXuWZ. A semi-analytical collocation Trefftz scheme for solving multi-term time fractional diffusion-wave equations. Eng Anal Boundary Elem (2019) 98:137–46. 10.1016/j.enganabound.2018.09.017
11.
SilviaFStefanA. The panel-clustering method for the wave equation in two spatial dimensions. J Comput Phys (2016) 305:217–43. 10.1016/j.jcp.2015.10.033
12.
LinJBaiJReutskiySLuJ. A novel RBF-based meshless method for solving time-fractional transport equations in 2D and 3D arbitrary domains. Eng Comput (2023) 39:1905–22. 10.1007/s00366-022-01601-0
13.
ZhengHYangZJZhangCZ. A local radial basis function collocation method for band structure computation of 3D phononic crystals. Appl Math Model (2020) 77:1954–64. 10.1016/j.apm.2019.09.006
14.
GuYFanCMQuWWangFZhangC. Localized method of fundamental solutions for three-dimensional inhomogeneous elliptic problems: Theory and MATLAB code. Comput Mech (2019) 64:1567–88. 10.1007/s00466-019-01735-x
15.
WangFZHouER. A direct meshless method for solving two-dimensional second-order hyperbolic telegraph equations. J Math (2020) 2020:1–9. 10.1155/2020/8832197
16.
WangFZZhangJAhmadIFarooqAAhmadH. A novel meshfree strategy for a viscous wave equation with variable coefficients. Front Phys (2021) 9:701512. 10.3389/fphy.2021.701512
17.
WangFZAhmadIAhmadHAlsulamiMDAlimgeerKSCesaranoCet alMeshless method based on RBFs for solving three-dimensional multi-term time fractional PDEs arising in engineering phenomenons. J King Saud Univ Sci (2021) 33(8):101604. 10.1016/j.jksus.2021.101604
18.
WangFZZhengKHAhmadIAhmadH. Gaussian radial basis functions method for linear and nonlinear convection–diffusion models in physical phenomena. Open Phys (2021) 19(1):69–76. 10.1515/phys-2021-0011
19.
ZhangJWangFZNadeemSSunM. Simulation of linear and nonlinear advection-diffusion problems by the direct radial basis function collocation method. Int Commun Heat Mass Transfer (2021) 130:105775. 10.1016/j.icheatmasstransfer.2021.105775
20.
FuZTangZXiQLiuQGGuYWangFJ. Localized collocation schemes and their applications. Acta Mechanica Sinica (2022) 38:422167. 10.1007/s10409-022-22167-x
21.
ChenWFuZJChenCS. Recent advances in radial basis function collocation methods. In: Springer briefs in applied sciences and Technology. Heidelberg, Germany: Springer (2014).
Summary
Keywords
semi-analytical method, radial basis function, meshless method, wave equation, numerical simulation
Citation
Wang F-Z, Shao M-Y, Li J-L and Zhang Z-L (2023) A space-time domain RBF method for 2D wave equations. Front. Phys. 11:1241196. doi: 10.3389/fphy.2023.1241196
Received
16 June 2023
Accepted
08 August 2023
Published
21 August 2023
Volume
11 - 2023
Edited by
Gang (Gary) Ren, Berkeley Lab (DOE), United States
Reviewed by
Zhuojia Fu, Hohai University, China
Xiao-Wei Gao, Dalian University of Technology, China
Updates
Copyright
© 2023 Wang, Shao, Li and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ming-Yu Shao, 12027@xzit.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.