A new integral formulation for 3D elastic wave propagation

This paper presents a new integral equation formulation forelastic wave propagation problems which allows for the use ofhigher-order boundary elements with C^0, continuity. This formulation,which is an extension of the one proposed by Panich for acousticproblems, has theoretical proofs of existence and uniquenessof solution for all values of frequency. The computational advantagesof the formulation are verified by numerical applications.     

Efficient simulation of wave propagation in heterogeneous materials

In this contribution we explain the core idea of the recently proposed spectral cell method, which combines a fictitious domain approach with finite elements of high order and noticeably relieves the burden of mesh generation. Moreover, it employs mass lumping techniques and significantly reduces the computational expenditure. Our studies show that the spectral cell method leads to similar results yet with less computational effort as compared to standard techniques. These properties turn the method to a viable tool for the wave propagation analysis of structures that obey a complicated geometry. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of mixed finite element methods on locally refined grids

Summary We consider the mixed finite element method for locally refined triangulations. A local projection operator is defined to satisfy the commutation property that is required in the general theory of mixed methods. Our results can be applied to every known space of arbitrary order over rectangles or triangles. Optimal-order error estimates and superconvergence for the flux along the Gauss lines are established.     

Analysis of a difference scheme for a singular perturbation Cauchy problem on refined grids

A Cauchy problem for a singular perturbation second-order ordinary differential equation is considered. It is proved that the upwind difference scheme on a grid proposed by Shishkin is uniformly convergent. The grid is well known only in application to a boundary value problem. The results of some numerical experiments are discussed.     

Simulation of elastic wave propagation in geological media: Intercomparison of three numerical methods

For wave propagation in heterogeneous media, we compare numerical results produced by grid-characteristic methods on structured rectangular and unstructured triangular meshes and by a discontinuous Galerkin method on unstructured triangular meshes as applied to the linear system of elasticity equations in the context of direct seismic exploration with an anticlinal trap model. It is shown that the resulting synthetic seismograms are in reasonable quantitative agreement. The grid-characteristic method on structured meshes requires more nodes for approximating curved boundaries, but it has a higher computation speed, which makes it preferable for the given class of problems.     

Finite difference on grids with nearly uniform cell area and line spacing for the wave equation on complex domains

A finite difference time-dependent numerical method for the wave equation, supported by recently derived novel elliptic grids, is analyzed. The method is successfully applied to single and multiple two-dimensional acoustic scattering problems including soft and hard obstacles with complexly shaped boundaries. The new grids have nearly uniform cell area (J-grids) and nearly uniform grid line spacing (αγ-grids). Numerical experiments reveal the positive impact of these two grid properties on the scattered field convergence to its harmonic steady state. The restriction imposed by stability conditions on the time step size is relaxed due to the near uniformity cell areas and grid line spacing. As a consequence, moderately large time steps can be used for relatively fine spatial grids resulting in greater accuracy at a lower computational cost. Also, numerical solutions for wave problems inside annular regions of complex shapes are obtained. The use of the new grids results in late time stability in contrast with other classical finite difference time-dependent methods.     

Wave propagation in periodically layered elastic media based on an internal variable model

A phenomenological model of wave propagation in a periodically layered elastic media, with linear elastic constituents, which is based on an internal variable theory, is presented; the main result is a system of coupled second- and fourth-order partial differential equations which describe the motion of each constituent and which, in turn, appear to predict the correct qualitative behavior for both the composite (or main) wave and the precursor wave in the laminate.     

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa-tions with a free surface. The 3D irreg...     

On specific features of the Lebedev scheme in simulating elastic wave propagation in anisotropic media

This paper presents the Lebedev scheme on staggered grids for the numerical simulation of wave propagation in anisotropic elastic media. Primary attention is given to the approximation of the elastic wave equation by the Lebedev scheme. Based on the differential approach, it is shown that the Lebedev scheme approximates a system of equations, which differs from the original equation. It is proved that the approximated system has a set of 24 characteristics, six of them coincide with those of the elastic wave equation and the rest ones are “artifacts.” Requiring the artificial solutions to be equal to zero and the true ones to coincide with those of the elastic wave equation, one comes to the classical definition of the approximation of the initial system on a sufficiently smooth solution. The results obtained and the knowledge of the complete set of characteristics are important for constructing reflectionless boundary conditions during approximation of point sources, etc.     

Energetic boundary element method analysis of wave propagation in 2D multilayered media

The analysis of scalar wave propagation in 2D zonewise homogeneous media with vanishing initial and mixed boundary conditions is carried out. The problem is formulated in terms of time‐dependent boundary integral equations, and then it is set in a weak form, based on a natural energy identity satisfied by the differential problem solution. Several numerical results have been obtained by means of the related energetic Galerkin boundary element method showing accuracy and stability of the method. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite‐difference time‐domain simulation of acoustic propagation in heterogeneous dispersive medium

Accurate modeling of pulse propagation and scattering is a problem in many disciplines (i.e., electromagnetics and acoustics). It is even more tenuous when the medium is dispersive. Blackstock [D. T. Blackstock, J Acoust Soc Am 77 (1985) 2050] first proposed a theory that resulted in adding an additional term (the derivative of the convolution between the causal time‐domain propagation factor and the acoustic pressure) that takes into account the dispersive nature of the medium. Thus deriving a modified wave equation applicable to either linear or nonlinear propagation. For the case of an acoustic wave propagating in a two‐dimensional heterogeneous dispersive medium, a finite‐difference time‐domain representation of the modified linear wave equation can been used to solve for the acoustic pressure. The method is applied to the case of scattering from and propagating through a 2‐D infinitely long cylinder with the properties of fat tissue encapsulating a cyst. It is found that ignoring the heterogeneity in the medium can lead to significant error in the propagated/scattered field. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Effect of the pipe curvature on internal elastic wave propagation

A mathematical model of internal elastic wave propagation in a bent pipe is developed. For slightly bent pipes, a method is devised for reducing the equations to a lower dimension problem. It is shown that the wave dynamics are described by the Korteweg-de Vries and Klein-Gordon-Fock equations. Thus, the class of problems to which these equations apply is expanded.     

A note on the application of the generalized finite difference method to seismic wave propagation in 2D

This paper shows the application of generalized finite difference method (GFDM) to the problem of seismic wave propagation. We investigated stability and star dispersion in 2D.     

3D complex structures imaging based on acoustic wave equation

In this paper, a general expression of the 3D hybrid imaging method based on acoustic wavefield extrapolation is presented. Moreover, the planewave synthesization method is given. The numerical results of 3D shot-profile migration and planewave synthesization migration for the SEG/EAEG 3D benchmark model show the good imaging ability of the hybrid imaging method. This method can be applied in field data processing. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite propagation speed for solutions of the wave equation on metric graphs

We provide a class of self-adjoint Laplace operators ?Δ on metric graphs with the property that the solutions of the associated wave equation satisfy the finite propagation speed property. The proof uses energy methods, which are adaptations of corresponding methods for smooth manifolds.     

Implementation and computational aspects of a 3D elastic wave modelling by PUFEM

This paper presents an enriched finite element model for three dimensional elastic wave problems, in the frequency domain, capable of containing many wavelengths per nodal spacing. This is achieved by applying the plane wave basis decomposition to the three-dimensional (3D) elastic wave equation and expressing the displacement field as a sum of both pressure (P) and shear (S) plane waves. The implementation of this model in 3D presents a number of issues in comparison to its 2D counterpart, especially regarding how S-waves are used in the basis at each node and how to choose the balance between P and S-waves in the approximation space. Various proposed techniques that could be used for the selection of wave directions in 3D are also summarised and used. The developed elements allow us to relax the traditional requirement which consists to consider many nodal points per wavelength, used with low order polynomial based finite elements, and therefore solve elastic wave problems without refining the mesh of the computational domain at each frequency. The effectiveness of the proposed technique is determined by comparing solutions for selected problems with available analytical models or to high resolution numerical results using conventional finite elements, by considering the effect of the mesh size and the number of enriching 3D plane waves. Both balanced and unbalanced choices of plane wave directions in space on structured mesh grids are investigated for assessing the accuracy and conditioning of this 3D PUFEM model for elastic waves.     

Analysis of the influence of Spectral Elements for elastic wave propagation on the CFL-Condition

For the simulation of the interaction of elastic waves in CFRP plates with inhomogeneities and defects the spectral finite element method (SEM) is under investigation. The SEM uses high-order shape functions which are composed of Lagrange polynomials with nodes at the Gauss-Lobatto quadrature (GLq) points. In this way we obtain a diagonal mass matrix which makes an explicit time scheme more efficient. The goal of this work is to investigate the effect of SEM on the CFL-Condition. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of oblique and multidirectional wave propagation and breaking on steep slope based on FEM model of Boussinesq equations

Creating a representative numerical simulation of the propagation and breaking of waves along slopes is an important problem in engineering design. Most studies on wave breaking have focused on the propagation of normal incident waves on gentle slopes. In practice, however, waves on steep slopes are obliquely incident or multidirectional irregular waves. In this paper, the eddy viscosity term is introduced to the momentum equation of the improved Boussinesq equations to model wave dissipation caused by breaking and friction, and a numerical model based on an unstructured finite element method (FEM) is established based on the governing equations. It is applied to simulate wave propagation on a steep slope of 1:5. Parallel physical experiments are conducted for comparative analysis that considered a large number of cases, including those featuring of normal and oblique incident regular and irregular waves, and multidirectional waves. The heights of the incident wave increase for different periods to represent different kinds of waves breaking. Based on examination, the effectiveness and accuracy of the numerical model is verified through a comprehensive comparison between the numerical and the experimental results, including in terms of variation in wave height, wave spectrum, and nonlinear parameters. Satisfactory agreement between the numerical and experimental values shows that the proposed model is effective in representing the breaking of oblique incident regular waves, irregular waves, and multidirectional incident irregular waves. However, the initial threshold of the breaking parameter η_t^(I) takes different values for oblique and multidirectional waves. This needs to be paid attention when the breaking of waves is simulated using the Boussinesq equations.