Numerical Simulation of Nonlinear Three-Dimensional Waves in Water of Arbitrary Varying Topography

The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.

3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

Reflection of three-dimensional plane waves at the free surface of a rotating triclinic half-space under the context of generalized thermoelasticity

The reflection of three-dimensional(3D) plane waves in a highly anisotropic(triclinic) medium under the context of generalized thermoelasticity is studied. The thermoelastic nature of the 3D plane waves in an anisotropic medium is investigated in the perspective of the three-phase-lag(TPL), dual-phase-lag(DPL), Green-Naghdi-III(GNIII), Lord-Shulman(LS), and classical coupled(CL) theories. The reflection coefficients and energy ratios for all the reflected waves are obtained in a mathematical form. The rotational effects on the reflection characteristics of the 3D waves are discussed under the context of generalized thermoelasticity. Comparative analyses for the reflection coefficients of the waves among these generalized thermoelastic theories are performed. The energy ratios for each of the reflected waves establish the energy conservation law in the reflection phenomena of the plane waves. The highly anisotropic materials along with the rotation may have a significant role in the phenomenon of the reflection behavior of the 3D waves. Numerical computations are performed for the graphical representation of the study.

Laboratory Research on Effective Test Area of Short-Crested Waves Generated by Two-Sided Segmented Wavemakers

The size and shape of the effective test area are crucial to consider when short-crested waves are created by segmented wavemakers. The range of the effective test area of short-crested waves simulated by two-sided segmented wavemakers is analyzed in this paper. The experimental investigation on the wave field distribution of short-crested waves generated by two-sided segmented wavemakers is conducted by using an array of wave gauges. Wave spectra and directional spreading function are analyzed and the results show that when the main direction is at a certain angle with the normal line of wave generators, the wave field of 3D short-crested waves generated by two-sided segmented wavemakers has good spatial uniformity within the model test area. The effective test area can provide good wave environments for seakeeping model tests of various ocean engineering structures in the deep ocean engineering basin.

Seismic analysis of underground tunnels subjected to 3D oblique incident P and SV waves by 2.5D approach

This paper presents,for the first time,the consideration of three-dimensional(3D)oblique incident P and SV waves in calculating the 3D seismic response of a lined tunnel embedded in a half-space by the 2.5D finite/infinite element method(FIEM).Firstly,the applicability of the 2.5D FIEM for 3D seismic analysis is summarized.With the exact solutions obtained for the free field in the Appendix,the equivalent seismic forces are rationally computed for the near-field boundary,considering the horizontal and vertical excitations of the Chi-Chi Earthquake.By performing seismic analysis of the half space embedded with a tunnel using the 2.5D FIEM,the time-domain responses of the tunnel are obtained.The accuracy of the present solutions is verified against those of de Barros and Luco.Conclusions drawn from the parametric study include:(1)Stress concentration for the principal stress under oblique incident seismic waves occurs at the polar angles of 0(vault),90,180(inverted arch),and 270of the lining wall.(2)The vault and inverted arch are the weakest parts of the tunnel during earthquakes.(3)The accelerations of the tunnel during earthquakes can be regarded as of the rigid body type.(4)The responses of the tunnel lining caused by SV waves of an earthquake are much more critical than those by P waves.(5)For arbitrary seismic waves,the maximum longitudinal acceleration azmax is of the same order of magnitude as the maximum horizontal acceleration axmax.

2.5-dimension soil seismic response to oblique incident waves based on exact free-field solution

With the three dimensional(3D)oblique incident waves exactly determined for the free field,the soil seismic responses in both frequency and time domains are studied by the 2.5 dimension(2.5D)finite/infinite element method.First,the free-field responses in frequency domain are solved exactly for 3D arbitrary incident P and SV waves,which requires no coordinate conversion or extra effort for SV waves with super-critical incident angles.Next,the earthquake spectra are incorporated by the concept of equivalent seismic forces on the near-field boundary,based only on the displacements input derived for unit ground accelerations of each frequency using the 2.5D approach.For the asymmetric 2.5D finite/infinite element model adopted,the procedure for soil seismic analysis is presented.The solutions computed by the proposed method are verified against those of Wolf’s and de Barros and Luco’s and for inversely calculated ground motions.Of interest is that abrupt variation in soil response occurs around the critical angle on the wave propagation plane for SV waves.In addition,the horizontal displacements attenuate with increasing horizontal incident angle,while the longitudinal ones increase inversely for 3D incident P and SV waves.

Three-dimensional acoustic wave equation modeling based on the optimal finite-difference scheme

Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.

Hybrid absorbing boundary condition for threedimensional elastic wave modeling

Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.

Tunable three-dimensional nonreciprocal transmission in a layered nonlinear elastic wave metamaterial by initial stresses

In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.

Developing Serpent-Type Wave Generators to Create Solitary Wave Simulations with BEM

Developing serpent-type wave generators to generate solitary waves in a 3D-basin was investigated in this study. Based on the Lagrangian description with time-marching procedures and finite differences of the time derivative, a 3D multiple directional wave basin with multidirectional piston wave generators was developed to simulate ocean waves by using BEM with quadrilateral elements, and to simulate wave-caused problems with fully nonlinear water surface conditions. The simulations of perpendicular solitary waves were conducted in the first instance to verify this scheme. Furthermore, the comparison of the waveform variations confirms that the estimation of 3D solitary waves is a feasible scheme.

Three-dimensional velocity structure and tectonic characteristics of earthquake area in Yibin

In this study,on the basis of absolute first-arrival times of 84756 P-and S-waves from 6085 earthquakes recorded at 56 fixed stations in Yibin and surrounding areas in China from January 2009 to January 2019,focal parameters and three-dimensional(3 D)body-wave high-resolution velocity structures at depths of 0–30 km were retrieved by double-difference tomography.Results show that there is a good correspondence between the spatial distribution of the relocated earthquakes and velocity structures,which were concentrated mainly in the high-velocity-anomaly region or edge of high-velocity region.Velocity structure of P-and S-waves in the Yibin area clearly shows lateral inhomogeneity.The distribution characteristics of the P-and S-waves near the surface are closely related to the geomorphology and geologic structure.The low-velocity anomaly appears at the depth of 15–25 km,which is affected by the lower crust current.The Junlian–Gongxian and Gongxian–Changning earthquake areas,which are the two most earthquake-prone areas in the Yibin region,clearly differ in earthquake distribution and tectonic characteristics.We analyzed the structural characteristics of the Junlian–Gongxian and Gongxian–Changning earthquake areas on the basis of the 3 D bodywave velocity structures in the Yibin region.We found that although most seismicity in the Yibin area is caused by fluid injection,the spatial position of seismicity is controlled by the velocity structures of the middle and upper crust and local geologic structure.Fine-scale 3 D velocity structures in the Yibin area provide important local reference information for further understanding the crustal medium,seismogenic structure,and seismicity.