Mathematical Problems in Engineering
Volume 2020 (2020), Article ID 7180489, 16 pages
https://doi.org/10.1155/2020/7180489

Research Article

Transient Wave Propagation Dynamics with Edge-Based Smoothed Finite Element Method and Bathe Time Integration Technique

Yingbin Chai1,2,3 and Yongou Zhang¹

¹School of Transportation, Wuhan University of Technology, Wuhan 430063, China
²School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
³Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration (CISSE), Shanghai 200240, China

Correspondence should be addressed to Yongou Zhang; zhangyo@whut.edu.cn

Received 1 May 2020; Revised 3 July 2020; Accepted 8 July 2020; Published 31 July 2020

Academic Editor: Marek Lefik

Copyright © 2020 Yingbin Chai and Yongou Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this work, the edge-based smoothed finite element method (ES-FEM) is incorporated with the Bathe time integration scheme to solve the transient wave propagation problems. The edge-based gradient smoothing technique (GST) can properly soften the “overly soft” system matrices from the standard finite element approach; then, the spatial numerical dispersion error of the calculated solutions for wave problems can be significantly reduced. To effectively solve the transient wave propagation problems, the Bathe time integration scheme is employed to perform the involved time integration. Due to the appropriate “numerical dissipation effects” from the Bathe time integration method, the spurious oscillations in the relatively large wave numbers (high frequencies) can be effectively suppressed; then, the temporal numerical dispersion error in the calculated solutions can also be notably controlled. A number of supporting numerical examples are considered to examine the capabilities of the present approach. The numerical results show that ES-FEM works very well with the Bathe time integration technique, and much more numerical solutions can be reached for solving transient wave propagation problems compared to the standard FEM.

1. Introduction

The wave propagation problems are usually encountered in real engineering applications. For several simple wave problems (such as a single wave propagating in a one-dimensional space), the exact solutions can be obtained using the analytical approach. When it comes to the relatively complex wave propagation problems, we have to resort to the numerical methods.

Actually, a lot of different numerical techniques can be employed to deal with the wave problems, such as the finite element method (FEM) [1, 2], finite volume method (FVM) [3], the boundary element [4–6] or boundary-type numerical method [7–11], and meshless techniques [12–22].

Among them, the classical finite element method is one of the most widely used and well-developed numerical approaches for analyzing wave problems. Based on the finite element method, several powerful and versatile commercial software packages have also been successfully developed for engineering computation. However, the finite element solutions for wave problems usually suffer from the numerical dispersion error issue for large wave numbers (namely, high frequencies) [23–27]. As a result, the calculated numerical solutions of wave problems from FEM are only reliable for relatively small wave numbers. In the large wave number range, the spurious oscillations always arise in the numerical results, and very inaccurate results can be frequently obtained. Quite importantly, when the low-order linear elements (such as triangular or quadrilateral element) are used, this troublesome dispersion issue will be even more severe.

In order to address the numerical dispersion error issue. Much research effort has been made and a variety of advanced or modified finite element schemes have been proposed (see, e.g., [28–34]), including the smoothed finite element method (S-FEM) (see, e.g., [35–40]). The S-FEM is developed by combining the classical finite element concepts and the generalized gradient smoothing technique (GGST) which is frequently employed in the meshless methods. Since the GGST can provide proper “softening effects” to the discretized model, then the inherent “overly stiff” property of the FEM model can be relieved to some extent and then the stiffness of obtained system matrices is closer to the real system. In consequence, the numerical dispersion error in solving wave propagation problems can be significantly reduced.

The present work mainly focuses on tackling the transient wave propagation problems. Many direct time integration techniques can be employed to solve the transient wave propagation dynamics, such as the central difference method [41], Wilson’s θ method [42], Houbolt’s method [43], Newmark’s method [44], and so on. However, these different time integration schemes usually show several undesirable properties. For example, the central difference method is only a conditional stable method for transient dynamic analysis; the well-known Newmark method is unconditional stable scheme, but several extra adjustable parameters should be carefully determined for desirable computation accuracy; the Wilson θ method and Houbolt method always lead to a relatively large numerical error in period elongation and amplitude decay [45].

In this work, the Bathe implicit time integration technique is combined with the edge-based FEM (ES-FEM) to solve the transient wave propagation problems. The Bathe time integration scheme is a typical two-substep method. The Newmark trapezoidal rule is used in the first substep and the 3-point Euler backward method is employed in the second substep. In the Bathe method, the appropriate numerical dissipation effects are introduced to suppress the spurious wave modes in high frequency range; then spurious oscillations in the calculated numerical solutions can be effectively eliminated and the solution accuracy can be significantly improved. The numerical examples in this paper show that the edge-based smoothed FEM (ES-FEM) works very well with the Bathe implicit time integration technique in solving transient wave propagation problems. The numerical dispersion error in the calculated numerical results can be significantly suppressed and much more accurate numerical solutions can be obtained compared to the conventional FEM.

The rest of the present work is structured as follows: in the next section the formulation of the ES-FEM for wave problems is briefly retrospected. A comprehensive dispersion analysis of the present ES-FEM with Bathe time integration scheme for transient wave propagation problems is given in Section 3; both the spatial dispersion error and temporal dispersion error are investigated in detail. A number of supporting numerical experiments and the related conclusions are then summarized in the remaining sections.

2. Formulation of the Edge-Based Smoothed FEM for Wave Problems

For a typical scalar wave problem (such as the wave propagation in a prestressed membrane), the governing equation is given by [46]

in which denotes the field variable (such as the displacement), is the wave speed, stands for the Laplace operator, and the overdot represents the derivative with respect to time.

Following the basic finite element discretizations [47, 48], we can obtain the corresponding matrix form of (1):

in which represents the unknown nodal variables; and are the obtained mass matrix and stiffness matrix, respectively:

in which is the used nodal shape function matrix.

In this work, the edge-based smoothed FEM [49], in which the standard FEM is combined with the edge-based gradient smoothing technique, is used to deal with the transient wave propagation problems. In the ES-FEM model, the standard triangular mesh is used, and each edge-based smoothing domain is formed by sequentially connecting the endpoints of each edge and the centers of the neighboring elements (as shown in Figure 1). For the interior edges, two different elements will participate in constructing the smoothing domain, while for the edges on the global boundary, the corresponding smoothing domain only involves one element. In consequence, each edge corresponds to each smoothing domain uniquely and the gradient smoothing technique will be operated over these obtained smoothing domains.

Figure 1: The description of forming the edge-based smoothing domains.

In this work, the gradient smoothing operation is performed by smoothing the particle velocity , and the smoothed particle velocity can be obtained by [49]

in which is the obtained smoothing domain for edge k and is a predefined smoothing function given by

in which stands for the area of the smoothing domain.

For the wave propagation in a prestressed membrane, the relationship between the particle velocity and the displacement is given by

in which stands for the outward unit normal vector.

Using (6) and (4), it can be rewritten by

in which is the boundary of the corresponding smoothing domain.

Following the conventional finite element interpolation steps, we have

in which denotes the number of involved nodes which participate in forming the smoothed gradient matrix .

In this work, the Gauss integration scheme is used to perform the related numerical integration along the boundary of smoothing domain which consists of segments; then smoothed gradient matrix can be obtained by

in which is the number of Gauss points in each segment and it is determined by the order of the used nodal shape functions, is the weighting coefficients, and and are the components of the outward unit normal vector.

Once the smoothed gradient matrix is obtained, then the smoothed element stiffness matrix and smoothed global stiffness matrix can be obtained as in the standard FEM scheme:

in which is the smoothed global stiffness matrix and is the smoothed element stiffness matrix for edge k.

3. Dispersion Analysis

The numerical solutions of the transient wave propagation problems usually suffer from the dispersion error issue. Both the spatial discretization and temporal discretization are able to lead to dispersion error. In this section, both the spatial discretization error and temporal discretization error will be investigated in detail by using the uniform node distributions with average nodal space h (see Figure 2).

(a)

(b)

Figure 2: The employed uniform mesh pattern and the illustration of involved nodes for dispersion analysis. (a) FEM. (b) ES-FEM.

3.1. Dispersion Error from Spatial Discretization

For the time-independent form of the wave equation, the following matrix equation can be obtained without considering the boundary condition [47, 48]:

in which is the wave number which is given by

in which is the angular frequency.

For an interior node (here p and q represent the row number and column number, resp.), Figure 2 presents all the nodes (the blue solid nodes) which have contributions in forming the system matrix equations associated with the central node (the red solid node) for both FEM and ES-FEM models. We can see that in the present ES-FEM model more nodes are used to construct the system matrix equation compared to those in the standard FEM model. The reason for this is that the edge-based gradient smoothing operations are employed in the ES-FEM model.

Since no additional boundary conditions are considered here, so the numerical solutions corresponding to the involved nodes in Figure 2 should have the same amplitude and can be represented by

in which is the amplitude of the numerical solution, , is the numerical wave number which corresponds to the exact wave number shown in (11), is a unit vector denoting the wave propagation direction, and is the position vector of the considered point.

By referring to the uniform mesh pattern shown in Figure 2 and substituting (13) into (11), we have

in which and are the matrices corresponding to the system stiffness matrix and system mass matrix .

For the standard FEM and ES-FEM models, the matrices and are obtained by

To ensure that the nontrivial solutions to (14) exist, the following equation should be satisfied [46]:

From (15)–(18), it is found that the matrices and are the functions of the numerical wave number , so the relationship between the exact wave number and the numerical wave number can be determined by (19). In other words, for any given numerical wave number , we can obtain the corresponding exact wave number by using (19); then we can obtain the spatial discretization error which is defined by .

Figure 3 gives the spatial discretization error in different wave propagation directions as a function of the normalized numerical wave number (here is the numerical wavelength) for both the standard FEM and the present ES-FEM. It is found that both the FEM and ES-FEM can produce very small spatial discretization error for the relatively small normalized wave numbers, and the spatial discretization error will become larger with the increase of the considered normalized wave number. However, it is clearly seen that the spatial discretization error resulting from the present ES-FEM is much smaller than that from the standard FEM in the whole considered normalized wave number range. The reason for this is that the edge-based gradient smoothing operations in the ES-FEM are indeed helpful to suppress the spatial discretization error in solving wave propagation problems. Therefore, it is reasonable to expect that the present ES-FEM can produce more accurate numerical solutions for solving transient wave propagation problems than the standard FEM.

(a)

(b)

Figure 3: Comparisons of the spatial discretization error as a function of the normalized numerical wave number in different wave propagation directions for the standard FEM and the present ES-FEM. (a) FEM. (b) ES-FEM.

In addition, we also can find that the spatial discretization error results from the standard FEM are strongly affected by the wave propagation directions; namely, the standard FEM shows clear “numerical anisotropy” property in solving wave propagation problems even if the uniform mesh pattern is used. From the results shown in Figure 3, it is seen that the present ES-FEM also suffers from the “numerical anisotropy” issue; however, it has been significantly relieved compared to the standard FEM.

3.2. Dispersion Error from Temporal Discretization

The temporal discretization is also an indispensable ingredient for solving the transient wave propagation problems. Many direct time integration schemes can be used to deal with the transient wave propagation dynamics. However, the temporal discretization always can result in additional temporal discretization error which can degrade the quality of the calculated numerical solutions. In this section, we mainly focus on investigating the additional dispersion effects from the temporal discretization. Since the Bathe time integration technique is a very effective approach in dynamic analysis and shows several evident advantages compared to other time integration schemes [45], in this work the Bathe time integration method is used for temporal discretization.

If the additional boundary conditions are not considered, the fundamental solution to the matrix equation shown in (2) has the following form:

in which denotes the numerical angular frequency.

By using (20) and referring to the node distributions shown in Figure 2, the following equation can be obtained from (2):

Using the standard Bathe time integration scheme [45], (21) can be rewritten by

in which and , is the exact angular frequency which can be obtained from (19) for a given numerical wave number , and is the used time step for time integration.

Based on the standard Bathe time integration scheme [45], the following equations can be directly obtained:

in which is the numerical damping ratio, is the critical time step and it satisfies , and AD denotes the percentage amplitude decay (AD) per period .

For any given numerical wave number , we can obtain the corresponding numerical angular frequency by using (23); then the total dispersion error for transient wave propagation problems can be obtained by

in which is the numerical wave speed and is the numerical period.

From (26), we can find the following interesting points:

(1)The total dispersion error of the numerical solutions for transient wave propagation problems can be split into two different parts: the first part, which is defined by , is from the spatial discretization and the second part, which is defined by , is induced by the temporal discretization(2)The spatial discretization error is determined by (19) and it is a function of the average nodal space h; when the used nodal space h trends to zero, the obtained spatial discretization error will also trend to zero (namely, will trend to 1)(3)For the given wave speed c and nodal space h, the temporal discretization error is mainly determined by the used time step ; when the used time step trends to zero, the obtained temporal discretization error will also trend to zero (namely, will trend to 1).

For a fixed time step  = 0.5, Figure 4 gives the calculated total dispersion error results (which is measured by ) in various wave propagation directions as a function of the normalized wave number for both FEM and ES-FEM. It is found that the magnitude of the dispersion error from the standard FEM is clearly much larger than that from the present ES-FEM. Both the standard FEM and the present ES-FEM suffers from the “numerical anisotropy” issue; namely, the dispersion error results are different in different considered wave propagations; however, it is also found that this issue has been successfully suppressed to a certain extent by the present ES-FEM. In addition, the additional dispersion error induced by the temporal discretization (which is measured by ) and the percentage amplitude decay (AD) per period for the above-mentioned two different numerical techniques are displayed in Figures 5 and 6. Likewise, the very similar findings can be obtained from the figures. This means that the present ES-FEM with the Bathe time integration technique possesses better predictive capabilities in tackling the transient wave propagation problems than the standard FEM, and much more accurate numerical solutions can be reached.

(a)

(b)

Figure 4: Comparisons of the total dispersion error results as a function of the normalized wave number in various wave propagation directions for both FEM and ES-FEM. (a) FEM. (b) ES-FEM.

(a)

(b)

Figure 5: Comparisons of the additional dispersion error results induced by the temporal discretization as a function of the normalized wave number in various wave propagation directions for both FEM and ES-FEM. (a) FEM. (b) ES-FEM.

(a)

(b)

Figure 6: Comparisons of the percentage amplitude decay (AD) as a function of the normalized wave number in various wave propagation directions for both FEM and ES-FEM. (a) FEM. (b) ES-FEM.

4. Numerical Examples

In this section, a number of supporting numerical examples are given to examine the effectiveness of the present ES-FEM with Bathe time integration scheme in solving transient wave propagation problems. It should be pointed out that the nonreflecting boundary conditions [50] (such as the Dirichlet-to-Neumann map, the absorbing boundary condition, and the perfect matched layer) are not considered in this work because all the wave components do not reach the boundary of the considered problem domain for the considered simulation time.

4.1. Two-Dimensional Scalar Wave Propagation in a Square Domain

The first considered numerical example is the two-dimensional scalar wave propagation in a square prestressed membrane. This problem is described in Figure 7. The length of the square problem domain is L = 2.4 m and the wave propagation speed is c = 1 m/s. A concentrated force is prescribed at the center of the square prestressed membrane. Due to symmetry, only a quarter of this in the involved problem domain is modelled for computation. The computation domain is discretized into uniform triangular elements with average nodal space h = 0.015 m and the time step  = 0.1 s is used for temporal discretization. For this scalar wave propagation problem, the transverse displacement is governed by the following equation:

(a)

(b)

Figure 7: The problem description and used mesh patterns for the two-dimensional scalar wave propagation. (a) The geometry description of the square prestressed membrane. (b) The used uniform mesh.

The first concentrated force is a Ricker wavelet which is defined by [46]

in which the magnitude A = 0.4 N, the peak frequency  = 5 Hz, and the time shift  = 0.25 s are used here.

Figure 8 gives the calculated displacement results from the standard FEM and the present ES-FEM along the wave propagation directions at the observation time t = 0.9 s. The exact solution shown in the figure is obtained by using the Green function method. From the results shown in the figure, it is very clear to see that the present ES-FEM can produce better numerical solutions than the standard FEM and there are almost no spurious oscillations in the ES-FEM results. However, the FEM results deviate from the exact solution very much, and several spurious peaks can also be found. These findings indicate that the edge-based gradient smoothing technique used in the present ES-FEM is indeed effective to suppress the dispersion error, and then the accuracy of the calculated numerical solutions can be significantly improved.

Figure 8: The calculated displacement results from the standard FEM and the present ES-FEM along the wave propagation directions for the first concentrated force case.

Figures 9 and 10 show the calculated displacement results from the standard FEM and the present ES-FEM along three different wave propagation directions (, , and ) at the observation time t = 0.9 s. It is seen that the numerical results from the standard FEM are strongly dependent on the wave propagation directions; namely, the standard FEM exhibits clear “numerical anisotropy” property in solving transient wave problems, while the present ES-FEM can successfully restrain this property to some degree. From Figure 10, we can find that the calculated solutions from the present ES-FEM are almost the same in different wave propagation directions. This numerical example demonstrates that the present ES-FEM works very well with the Bathe time integration scheme and indeed has obvious advantages in solving transient wave propagation problems compared to the standard FEM.

Figure 9: Comparisons of the calculated displacement results from the standard FEM along various propagation directions for the first concentrated force case.

Figure 10: Comparisons of the calculated displacement results from the present ES-FEM along various propagation directions for the first concentrated force case.

Next, we consider another concentrated force which is defined by [46]

For this case, the displacements from the standard FEM and the present ES-FEM along the wave propagation direction at observation time t = 0.8 s are firstly calculated and shown in Figure 11. Similarly, we also investigated how the numerical solutions from the two mentioned methods are affected by the wave propagation directions, and the related numerical results are plotted in Figures 12 and 13. We again find that the present ES-FEM shows weaker “numerical anisotropy” property than the standard FEM in solving transient wave propagation problems, and better numerical results can be obtained.

Figure 11: The calculated displacement results from the standard FEM and the present ES-FEM along the wave propagation directions for the second concentrated force case.

Figure 12: Comparisons of the calculated displacement results from the standard FEM along various propagation directions for the second concentrated force case.

Figure 13: Comparisons of the calculated displacement results from the present ES-FEM along various propagation directions for the second concentrated force case.

4.2. Two-Dimensional Scalar Wave Propagation Scattering Problem

The second considered numerical example is still the two-dimensional scalar wave propagation in a square prestressed membrane. In this case, four same holes are located uniformly in the problem domain. This problem is described in Figure 14. The external excitation load at the center of the membrane is still a Ricker wavelet with magnitude A = 0.4 N, the peak frequency  = 10 Hz, and the time shift  = 0.1 s. It is clear that both reflected wave and transmitted wave will arise when the original incident wave reach the boundary of the holes. In this numerical example, the abilities of the two different methods (FEM and ES-FEM) in predicting the behaviors of these different wave components will be investigated in detail. Likewise, only a quarter of the problem is modelled due to symmetry and the computation problem domain is discretized using standard triangular mesh with average nodal space h = 0.125 m, and the time step  = 0.1 s is used for temporal discretization. Note that the corresponding exact solutions to this problem are not easy to obtain; here the numerical solutions from the traditional FEM with a very fine mesh are provided as the reference solutions.

(a)

(b)

Figure 14: The problem description and used mesh patterns for the two-dimensional scalar wave scattering problem. (a) The geometry description of the considered problem. (b) The used mesh pattern.

At the observation time t = 0.9 s, the comparisons of the calculated displacement results from the standard FEM and the ES-FEM along the wave propagation direction are plotted in Figure 15. In this wave propagation direction, there only exists the original incident wave, and the ES-FEM results are very close to the reference solution, while the FEM results are much worse than the ES-FEM results and many spurious peaks can be found in the FEM results.

Figure 15: Comparisons of the calculated displacement results from the standard FEM and the ES-FEM along the wave propagation direction for two-dimensional scalar wave propagation scattering problem.

Then, the corresponding displacement results along the other two different wave propagation directions ( and ) at the observation time t = 1 s are presented in Figures 16 and 17. In these wave propagation directions, both the reflected waves and the transmitted waves can be seen, and it is found that the present ES-FEM is indeed superior to the standard FEM in predicting the behaviors of the reflected waves and the transmitted waves because the ES-FEM results are much closer to the reference solutions than the FEM results.

Figure 16: Comparisons of the calculated displacement results from the standard FEM and the ES-FEM along the wave propagation direction for two-dimensional scalar wave propagation scattering problem.

Figure 17: Comparisons of the calculated displacement results from the standard FEM and the ES-FEM along the wave propagation direction for two-dimensional scalar wave propagation scattering problem.

Furthermore, several snapshots of the displacement distributions in the global computation domain from the standard FEM and ES-FEM at several observation times are shown in Figures 18 and 19. From the results shown in the figures, we can see that the dispersion error in the ES-FEM solutions is much smaller than that in the standard FEM solutions; hence, much smoother displacement distribution results can be reached and all the physical behaviors of the different wave components can be accurately predicted.

(a)

(b)

(c)

(d)

Figure 18: Snapshots of the displacement distributions from the FEM for the two-dimensional scalar wave scattering problem at different observation times: (a) t = 0.6 s; (b) t = 0.7 s; (c) t = 0.8 s; (d) t = 0.9 s.

(a)

(b)

(c)

(d)

Figure 19: Snapshots of the displacement distributions from the ES-FEM for the two-dimensional scalar wave scattering problem at different observation times: (a) t = 0.6 s; (b) t = 0.7 s; (c) t = 0.8 s; (d) t = 0.9 s.

5. Conclusions

This paper focuses on presenting an edge-based smoothed FEM (ES-FEM) with the Bathe time integration technique for solving two-dimensional transient wave propagation problems. The dispersion error properties of the present approach for solving transient wave problems are investigated detailedly. It is found that the present ES-FEM works quite well with the Bathe time integration technique, and the dispersion error can be significantly reduced compared to the standard FEM. Quite importantly, the troublesome “numerical anisotropy” property of the standard FEM for transient wave problems can be clearly suppressed by the present ES-FEM; hence, the obtained numerical results are almost the same in the different wave propagation directions. The numerical results show that the present ES-FEM indeed surpasses the standard FEM in solving transient wave problems and can provide much more accurate numerical results in predicting the behaviors of the waves; hence, it can be regarded as a very promising numerical approach for solving practical transient wave dynamic problems.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors wish to express their gratitude to the National Natural Science Foundation of China (Grant no. 51809208), the Fundamental Research Funds for the Central Universities (Grant no. WUT: 2019IVB012), and the China Postdoctoral Science Foundation (Grant nos. 2018M632866 and 2018M642940).

References

I. Harari, “A survey of finite element methods for time-harmonic acoustics,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 13–16, pp. 1594–1607, 2006.

X. Wu, M. Yun, M. Wang et al., “Self-imaging in multi-walled carbon nanotube arrays at visible wavelengths,” Carbon, vol. 108, pp. 47–51, 2016.

T. R. Fogarty and R. J. Leveque, “High-resolution finite-volume methods for acoustic waves in periodic and random media,” The Journal of the Acoustical Society of America, vol. 106, no. 1, pp. 17–28, 1999.

J. Li, W. Chen, and Q. Qin, “A modified dual-level fast multipole boundary element method based on the Burton-Miller formulation for large-scale three-dimensional sound field analysis,” Computer Methods in Applied Mechanics and Engineering, vol. 340, pp. 121–146, 2018.

L. Chen and X. Li, “Boundary element-free methods for exterior acoustic problems with arbitrary and high wavenumbers,” Applied Mathematical Modelling, vol. 72, pp. 85–103, 2019.

L. Chen, X. Liu, and X. Li, “The boundary element-free method for 2D interior and exterior Helmholtz problems,” Computers & Mathematics with Applications, vol. 77, no. 3, pp. 846–864, 2019.

Z.-J. Fu, Q. Xi, W. Chen, and A. H.-D. Cheng, “A boundary-type meshless solver for transient heat conduction analysis of slender functionally graded materials with exponential variations,” Computers & Mathematics with Applications, vol. 76, no. 4, pp. 760–773, 2018.

J. Li, W. Chen, Z. Fu, and Q.-H. Qin, “A regularized approach evaluating the near-boundary and boundary solutions for three-dimensional Helmholtz equation with wideband wavenumbers,” Applied Mathematics Letters, vol. 91, pp. 55–60, 2019.

J. Li and W. Chen, “A modified singular boundary method for three-dimensional high frequency acoustic wave problems,” Applied Mathematical Modelling, vol. 54, pp. 189–201, 2018.

J. Li, W. Chen, Q.-H. Qin, and Z. Fu, “A modified multilevel algorithm for large-scale scientific and engineering computing,” Computers & Mathematics with Applications, vol. 77, no. 8, pp. 2061–2076, 2019.

F. J. Wang, C. M. Fan, Q. S. Hua, and Y. Gu, “Localized MFS for the inverse Cauchy problems of two-dimensional Laplace and biharmonic equations,” Applied Mathematics and Computation, vol. 364, Article ID 124658, 2020.

Z.-J. Fu, J. Zhang, P.-W. Li, and J.-H. Zheng, “A semi-Lagrangian meshless framework for numerical solutions of two-dimensional sloshing phenomenon,” Engineering Analysis with Boundary Elements, vol. 112, pp. 58–67, 2020.

Z. J. Fu, Z. Y. Xie, S. Y. Ji, C. C. Tsai, and A. L. Li, “Meshless generalized finite difference method for water wave interactions with multiple-bottom-seated-cylinder-array structures,” Ocean Engineering, vol. 195, Article ID 106736, 2020.

W. Qu, C.-M. Fan, and X. Li, “Analysis of an augmented moving least squares approximation and the associated localized method of fundamental solutions,” Computers & Mathematics with Applications, vol. 80, no. 1, pp. 13–30, 2020.

X. Y. You, W. Li, and Y. B. Chai, “A truly meshfree method for solving acoustic problems using local weak form and radial basis functions,” Applied Mathematics and Computation, vol. 365, Article ID 124694, 2020.

T. Zhang and X. L. Li, “Variational multiscale interpolating element-free Galerkin method for the nonlinear Darcy–Forchheimer model,” Computers & Mathematics with Applications, vol. 79, no. 2, pp. 363–377, 2019.

J. Lin, Y. Xu, and Y. H. Zhang, “Simulation of linear and nonlinear advection-diffusion-reaction problems by a novel localized scheme,” Applied Mathematics Letters, vol. 99, Article ID 106005, 2020.

X. You, W. Li, and Y. Chai, “Dispersion analysis for acoustic problems using the point interpolation method,” Engineering Analysis with Boundary Elements, vol. 94, pp. 79–93, 2018.

Y. Gu and H. Sun, “A meshless method for solving three-dimensional time fractional diffusion equation with variable-order derivatives,” Applied Mathematical Modelling, vol. 78, pp. 539–549, 2020.

Y. Wang, Y. Gu, and J. Liu, “A domain-decomposition generalized finite difference method for stress analysis in three-dimensional composite materials,” Applied Mathematics Letters, vol. 104, Article ID 106226, 2020.

F. Wang, C. M. Fan, C. Zhang et al., “A localized space-time method of fundamental solutions for diffusion and convection-diffusion problems,” Advances in Applied Mathematics and Mechanics, vol. 12, no. 4, pp. 940–958, 2020.

F. Wang, Y. Gu, W. Qu et al., “Localized boundary knot method and its application to large-scale acoustic problems,” Computer Methods in Applied Mechanics and Engineering, vol. 361, Article ID 112729, 2020.

A. Deraemaeker, I. Babuka, and P. Bouillard, “Dispersion and pollution of the FEM solution for the Helmholtz equation in one, two and three dimensions,” International Journal for Numerical Methods in Engineering, vol. 46, no. 4, pp. 471–499, 1999.

X. Hu, X. Cui, Q. Zhang, G. Wang, and G. Li, “The stable node-based smoothed finite element method for analyzing acoustic radiation problems,” Engineering Analysis with Boundary Elements, vol. 80, pp. 142–151, 2017.

X. Hu, X. Y. Cui, Q. Y. Zhang, G. Wang, and G. Y. Li, “Dispersion error analysis of stable node-based finite element method for the Helmholtz equation,” Communications in Computational Physics, vol. 23, no. 3, pp. 795–821, 2018.

Y. B. Chai, Z. X. Gong, W. Li et al., “Application of smoothed finite element method to two-dimensional exterior problems of acoustic radiation,” International Journal of Computational Methods, vol. 15, no. 5, Article ID 1850029, 2018.

Y. B. Chai, X. Y. You, and W. Li, “Dispersion reduction for the wave propagation problems using a coupled “FE-Meshfree” triangular element,” International Journal of Computational Methods, Article ID 1950071, 2019.

I. Harari and T. J. R. Hughes, “Galerkin/least-squares finite element methods for the reduced wave equation with non-reflecting boundary conditions in unbounded domains,” Computer Methods in Applied Mechanics and Engineering, vol. 98, no. 3, pp. 411–454, 1992.

Z. C. He, E. Li, G. R. Liu, G. Y. Li, and A. G. Cheng, “A mass-redistributed finite element method (MR-FEM) for acoustic problems using triangular mesh,” Journal of Computational Physics, vol. 323, pp. 149–170, 2016.

Y. Xu, G. Zhang, B. Zhou, H. Wang, and Q. Tang, “Analysis of acoustic radiation problems using the cell-based smoothed radial point interpolation method with Dirichlet-to-Neumann boundary condition,” Engineering Analysis with Boundary Elements, vol. 108, pp. 447–458, 2019.

G. Zhang, Z. Chen, Z. Sui et al., “A cell-based smoothed radial point interpolation method with virtual nodes for three-dimensional mid-frequency acoustic problems,” International Journal for Numerical Methods in Engineering, vol. 119, no. 6, pp. 548–566, 2019.

W. Li, Y. Chai, M. Lei, and T. Li, “Numerical investigation of the edge-based gradient smoothing technique for exterior Helmholtz equation in two dimensions,” Computers & Structures, vol. 182, pp. 149–164, 2017.

D. Komatitsch, C. Barnes, and J. Tromp, “Simulation of anisotropic wave propagation based upon a spectral element method,” Geophysics, vol. 65, no. 4, pp. 1251–1260, 2000.

G. Seriani and S. P. Oliveira, “Dispersion analysis of spectral element methods for elastic wave propagation,” Wave Motion, vol. 45, no. 6, pp. 729–744, 2008.

G. R. Liu and T. T. Nguyen-Thoi, “Smoothed Finite Element Method”, CRC Press, New York, NY, USA, 2010.

G. Wang, X. Y. Cui, H. Feng, and G. Y. Li, “A stable node-based smoothed finite element method for acoustic problems,” Computer Methods in Applied Mechanics and Engineering, vol. 297, pp. 348–370, 2015.

Y. Chai, W. Li, Z. Gong, and T. Li, “Hybrid smoothed finite element method for two-dimensional underwater acoustic scattering problems,” Ocean Engineering, vol. 116, pp. 129–141, 2016.

C. Jiang, Z.-Q. Zhang, G.-J. Gao, and G. R. Liu, “A modified immersed smoothed FEM with local field reconstruction for fluid-structure interactions,” Engineering Analysis with Boundary Elements, vol. 107, pp. 218–232, 2019.

Y. Chai, Z. Gong, W. Li, T. Li, and Q. Zhang, “A smoothed finite element method for exterior Helmholtz equation in two dimensions,” Engineering Analysis with Boundary Elements, vol. 84, pp. 237–252, 2017.

W. Li, Z. X. Gong, Y. B. Chai et al., “Hybrid gradient smoothing technique with discrete shear gap method for shell structures,” Computers & Mathematics with Applications, vol. 74, no. 8, pp. 1826–1855, 2017.

B. Wu, H. Bao, J. Ou, and S. Tian, “Stability and accuracy analysis of the central difference method for real-time substructure testing,” Earthquake Engineering & Structural Dynamics, vol. 34, no. 7, pp. 705–718, 2005.

K.-J. Bathe and G. Noh, “Insight into an implicit time integration scheme for structural dynamics,” Computers & Structures, vol. 98-99, pp. 1–6, 2012.

A. Soroushian and J. Farjoodi, “A unified starting procedure for the Houbolt method,” Communications in Numerical Methods in Engineering, vol. 24, no. 1, pp. 1–13, 2008.

D. Roy and M. K. Dash, “A stochastic newmark method for engineering dynamical systems,” Journal of Sound and Vibration, vol. 249, no. 1, pp. 83–100, 2002.

K.-J. Bathe, “Conserving energy and momentum in nonlinear dynamics: a simple implicit time integration scheme,” Computers & Structures, vol. 85, no. 7-8, pp. 437–445, 2007.

K.-T. Kim, L. Zhang, and K.-J. Bathe, “Transient implicit wave propagation dynamics with overlapping finite elements,” Computers & Structures, vol. 199, pp. 18–33, 2018.

Y. B. Chai and K. J. Bathe, “Transient wave propagation in inhomogeneous media with enriched overlapping triangular elements,” Computers & Structures, vol. 237, Article ID 106273, 2020.

Y. Chai, C. Cheng, W. Li, and Y. Huang, “A hybrid finite element-meshfree method based on partition of unity for transient wave propagation problems in homogeneous and inhomogeneous media,” Applied Mathematical Modelling, vol. 85, pp. 192–209, 2020.

Y. Chai, X. You, W. Li, Y. Huang, Z. Yue, and M. Wang, “Application of the edge-based gradient smoothing technique to acoustic radiation and acoustic scattering from rigid and elastic structures in two dimensions,” Computers & Structures, vol. 203, pp. 43–58, 2018.

J. B. Keller and D. Givoli, “Exact non-reflecting boundary conditions,” Journal of Computational Physics, vol. 82, no. 1, pp. 172–192, 1989.