The wave equation in generalized coordinates

Geophysics ◽  
1994 ◽  
Vol 59 (12) ◽  
pp. 1911-1919 ◽  
Author(s):  
José M. Carcione

This work introduces a spectral collocation scheme for the viscoelastic wave equation transformed from Cartesian to generalized coordinates. Both the spatial derivatives of field variables and the metrics of the transformation are calculated by the Chebychev pseudospectral method. The technique requires a special treatment of the boundary conditions, which is based on 1-D characteristics normal to the boundaries. The numerical solution of Lamb’s problem requires two 1-D stretching transformations for each Cartesian direction. The results show excellent agreement between the elastic numerical and analytical solutions, demonstrating the effectiveness of the differential operator and boundary treatment. Another example, requiring 1-D transformations, tests the propagation of a Rayleigh wave around a corner of the numerical mesh. Two‐dimensional transformations adapt the grid to topographic features: a syncline, and an anticlinal structure formed with fine layers.

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. A53-A56 ◽  
Author(s):  
José M. Carcione

The Fourier pseudospectral (PS) method is generalized to the case of derivatives of nonnatural order (fractional derivatives) and irrational powers of the differential operators. The generalization is straightforward because the calculation of the spatial derivatives with the fast Fourier transform is performed in the wavenumber domain, where the operator is an irrational power of the wavenumber. Modeling constant-[Formula: see text] propagation with this approach is highly efficient because it does not require memory variables or additional spatial derivatives. The classical acoustic wave equation is modified by including those with a space fractional Laplacian, which describes wave propagation with attenuation and velocity dispersion. In particular, the example considers three versions of the uniform-density wave equation, based on fractional powers of the Laplacian and fractional spatial derivatives.


Geophysics ◽  
2002 ◽  
Vol 67 (5) ◽  
pp. 1602-1611 ◽  
Author(s):  
Andrew Curtis ◽  
Johan O. A. Robertsson

“Volumetric recording” of the seismic wavefield implies that the local receiver group or array approximately encloses a volume of the earth. We show how volumetric recording can be used to measure several spatial derivatives of the wavefield. By making use of the full elastic wave equation, the free surface condition on elastic wavefields, and derivative centering techniques analagous to Lax‐Wendroff corrections used in synthetic finite‐difference modeling, these derivative estimates can be inverted for P‐ and S‐velocities in the near surface directly beneath the receiver group. The quantities estimated are the effective velocities of the P‐ and S‐components experienced by the wavefield at any point in time. Hence, the velocity estimates may vary with both wave type and wavelength. The estimates may be useful to aid statics estimation and are exactly the effective velocities required for separation of the wavefield into P‐ and S‐, and up‐ and down‐going components.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. T335-T345
Author(s):  
Enjiang Wang ◽  
José M. Carcione ◽  
Jing Ba ◽  
Mamdoh Alajmi ◽  
Ayman N. Qadrouh

We have applied the nearly perfectly matched layer (N-PML) absorber to the viscoelastic wave equation based on the Kelvin-Voigt and Zener constitutive equations. In the first case, the stress-strain relation has the advantage of not requiring additional physical field (memory) variables, whereas the Zener model is more adapted to describe the behavior of rocks subject to wave propagation in the whole frequency range. In both cases, eight N-PML artificial memory variables are required in the absorbing strips. The modeling simulates 2D waves by using two different approaches to compute the spatial derivatives, generating different artifacts from the boundaries, namely, 16th-order finite differences, where reflections from the boundaries are expected, and the staggered Fourier pseudospectral method, where wraparound occurs. The time stepping in both cases is a staggered second-order finite-difference scheme. Numerical experiments demonstrate that the N-PML has a similar performance as in the lossless case. Comparisons with other approaches (S-PML and C-PML) are carried out for several models, which indicate the advantages and drawbacks of the N-PML absorber in the anelastic case.


2018 ◽  
Vol 24 (23) ◽  
pp. 5650-5664 ◽  
Author(s):  
Shang–Teh Wu ◽  
Shan-Qun Tang ◽  
Kuan–Po Huang

This paper investigates the vibration control of a two-link flexible manipulator carried by a translational stage. The first and the second links are each driven by a stage motor and a joint motor. By treating the joint motor as a virtual spring, the two-link manipulator can be regarded as an integral flexible arm driven by the stage motor. A noncollocated controller is devised based on feedback from the deflection of the virtual spring, which can be measured by a shaft encoder. Stability of the closed-loop system is analyzed by examining the spatial derivatives of the modal functions. By including a bandpass filter in the feedback loop, residual vibrations can be attenuated without exciting high-frequency vibrations. The control method is simple to implement; its effectiveness is confirmed by simulation and experimental results.


1983 ◽  
Vol 29 (2) ◽  
pp. 243-253 ◽  
Author(s):  
Tomikazu Namikawa ◽  
Hiromitsu Hamabata

The ponderomotive force generated by random Alfvén waves in a collisionless plasma is evaluated taking into account mean magnetic and velocity shear and is expressed as a series involving spatial derivatives of mean magnetic and velocity fields whose coefficients are associated with the helicity spectrum function of random velocity field. The effect of microscale random Alfvén waves through ponderomotive and mean electromotive forces generated by them on the propagation of large-scale Alfvén waves is also investigated.


2021 ◽  
Vol 5 (4) ◽  
pp. 203
Author(s):  
Suzan Cival Buranay ◽  
Nouman Arshad ◽  
Ahmed Hersi Matan

We give fourth-order accurate implicit methods for the computation of the first-order spatial derivatives and second-order mixed derivatives involving the time derivative of the solution of first type boundary value problem of two dimensional heat equation. The methods are constructed based on two stages: At the first stage of the methods, the solution and its derivative with respect to time variable are approximated by using the implicit scheme in Buranay and Arshad in 2020. Therefore, Oh4+τ of convergence on constructed hexagonal grids is obtained that the step sizes in the space variables x1, x2 and in time variable are indicated by h, 32h and τ, respectively. Special difference boundary value problems on hexagonal grids are constructed at the second stages to approximate the first order spatial derivatives and the second order mixed derivatives of the solution. Further, Oh4+τ order of uniform convergence of these schemes are shown for r=ωτh2≥116,ω>0. Additionally, the methods are applied on two sample problems.


2019 ◽  
Vol 23 (6 Part B) ◽  
pp. 3739-3744
Author(s):  
Feng Gao

In this article we investigate the general fractional-order derivatives of the Riemann-Liouville type via Lorenzo-Hartley kernel, general fractional-order integrals and the new general fractional-order wave equation defined on the definite domain with the analytical soluton.


2018 ◽  
Vol 184 ◽  
pp. 01023
Author(s):  
Gordana V. Jelić ◽  
Vladica Stanojević ◽  
Dragana Radosavljević

One of the basic equations of mathematical physics (for instance function of two independent variables) is the differential equation with partial derivatives of the second order (3). This equation is called the wave equation, and is provided when considering the process of transverse oscillations of wire, longitudinal oscillations of rod, electrical oscillations in a conductor, torsional vibration at waves, etc… The paper shows how to form the equation (3) which is the equation of motion of each point of wire with abscissa x in time t during its oscillation. It is also shown how to determine the equation (3) in the task of electrical oscillations in a conductor. Then equation (3) is determined, and this solution satisfies the boundary and initial conditions.


Sign in / Sign up

Export Citation Format

Share Document