High order modified differential equation of the Beam–Warming method, I. The dispersive features

2020 ◽  
Vol 35 (2) ◽  
pp. 83-94 ◽  
Author(s):  
Yurii Shokin ◽  
Ireneusz Winnicki ◽  
Janusz Jasinski ◽  
Slawomir Pietrek
Keyword(s):  
Differential Equation ◽  
Order Equation ◽  
Advection Equation ◽  
Explicit Difference Scheme ◽  
Eighth Order ◽  
Differential Approximation ◽  
Phase And Group Velocities ◽  
Derivatives Of ◽  
Modified Equation ◽  
Group Velocities

Abstract The analysis of the modified partial differential equation (MDE) of the constant wind speed advection equation explicit difference scheme up to the eighth order with respect to both space and time derivatives is presented. So far, in majority of publications this modified equation has been derived mainly as a fourth-order equation. The MDE is presented in the so-called Π-form of the first differential approximation. This form includes only the space derivatives of higher order p and their coefficients μ(p). Analysis of these coefficients provides indications of the nature of the dissipative and dispersive errors. A fragment of the stencil for determining the modified differential equation up to the eighth-order MDE for the second-order Beam–Warming scheme is included. The derived coefficients μ(p) as well as the analysis of the phase shift errors, the phase and group velocities and dispersive features on the basis of these coefficients have not been published so far. The dissipative features of this method we present in [33].

High order modified differential equation of the Beam–Warming method, II. The dissipative features

2020 ◽  
Vol 35 (3) ◽  
pp. 175-185
Author(s):  
Yurii Shokin ◽  
Ireneusz Winnicki ◽  
Janusz Jasinski ◽  
Slawomir Pietrek
Keyword(s):  
Differential Equation ◽  
Partial Differential Equation ◽  
Wind Speed ◽  
Difference Scheme ◽  
Advection Equation ◽  
Explicit Difference Scheme ◽  
Eighth Order ◽  
Differential Approximation ◽  
Partial Differential ◽  
Linear Advection Equation

AbstractThis paper is a continuation of [38]. The analysis of the modified partial differential equation (MDE) of the constant-wind-speed linear advection equation explicit difference scheme up to the eighth-order derivatives is presented. In this paper the authors focus on the dissipative features of the Beam–Warming scheme. The modified partial differential equation is presented in the so-called Π-form of the first differential approximation. The most important part of this form includes the coefficients μ (p) at the space derivatives. Analysis of these coefficients provides indications of the nature of the dissipative errors. A fragment of the stencil for determining the modified differential equation for the Beam–Warming scheme is included. The derived and presented coefficients μ (p) as well as the analysis of the dissipative features of this scheme on the basis of these coefficients have not been published so far.

Asymptotic derivation of a refined equation for an elastic beam resting on a Winkler foundation

2021 ◽  
pp. 108128652110238
Author(s):  
Barış Erbaş ◽  
Julius Kaplunov ◽  
Isaac Elishakoff
Keyword(s):  
Ad Hoc ◽  
Moving Load ◽  
Low Frequency ◽  
Elastic Beam ◽  
Winkler Foundation ◽  
One Dimensional ◽  
Beam Equations ◽  
Dimensional Equation ◽  
Phase And Group Velocities ◽  
Group Velocities

A two-dimensional mixed problem for a thin elastic strip resting on a Winkler foundation is considered within the framework of plane stress setup. The relative stiffness of the foundation is supposed to be small to ensure low-frequency vibrations. Asymptotic analysis at a higher order results in a one-dimensional equation of bending motion refining numerous ad hoc developments starting from Timoshenko-type beam equations. Two-term expansions through the foundation stiffness are presented for phase and group velocities, as well as for the critical velocity of a moving load. In addition, the formula for the longitudinal displacements of the beam due to its transverse compression is derived.

Compressional‐wave velocities in attenuating media: A laboratory physical model study

Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1162-1167 ◽  
Author(s):  
Joseph B. Molyneux ◽  
Douglas R. Schmitt
Keyword(s):  
Phase Velocity ◽  
Group Velocity ◽  
Wave Velocity ◽  
Propagation Distance ◽  
Compressional Wave ◽  
Arrival Times ◽  
Wave Velocities ◽  
Amplitude Maximum ◽  
Phase And Group Velocities ◽  
Group Velocities

Elastic‐wave velocities are often determined by picking the time of a certain feature of a propagating pulse, such as the first amplitude maximum. However, attenuation and dispersion conspire to change the shape of a propagating wave, making determination of a physically meaningful velocity problematic. As a consequence, the velocities so determined are not necessarily representative of the material’s intrinsic wave phase and group velocities. These phase and group velocities are found experimentally in a highly attenuating medium consisting of glycerol‐saturated, unconsolidated, random packs of glass beads and quartz sand. Our results show that the quality factor Q varies between 2 and 6 over the useful frequency band in these experiments from ∼200 to 600 kHz. The fundamental velocities are compared to more common and simple velocity estimates. In general, the simpler methods estimate the group velocity at the predominant frequency with a 3% discrepancy but are in poor agreement with the corresponding phase velocity. Wave velocities determined from the time at which the pulse is first detected (signal velocity) differ from the predominant group velocity by up to 12%. At best, the onset wave velocity arguably provides a lower bound for the high‐frequency limit of the phase velocity in a material where wave velocity increases with frequency. Each method of time picking, however, is self‐consistent, as indicated by the high quality of linear regressions of observed arrival times versus propagation distance.

scholarly journals Analytical solution to bending of stiffened and continuous antisymmetric laminates

2015 ◽  
Vol 2 (1) ◽  
Author(s):  
Liecheng Sun ◽  
Issam E. Harik
Keyword(s):  
Differential Equation ◽  
Partial Differential Equation ◽  
Analytical Solution ◽  
Displacement Function ◽  
The Other ◽  
Order Partial Differential Equation ◽  
Eighth Order ◽  
Coupling Problem ◽  
Simply Supported ◽  
General Response

AbstractAnalytical Strip Method is presented for the analysis of the bending-extension coupling problem of stiffened and continuous antisymmetric thin laminates. A system of three equations of equilibrium, governing the general response of antisymmetric laminates, is reduced to a single eighth-order partial differential equation (PDE) in terms of a displacement function. The PDE is then solved in a single series form to determine the displacement response of antisymmetric cross-ply and angle-ply laminates. The solution is applicable to rectangular laminates with two opposite edges simply supported and the other edges being free, clamped, simply supported, isotropic beam supports, or point supports.

Errors Introduced in Estimation of Surface Wave Phase and Group Velocities Due to Wrong Assumptions: An Assessment Using a Simple Model For Love Wave

2021 ◽  
Author(s):  
Akash Kharita ◽  
Sagarika Mukhopadhyay
Keyword(s):  
Surface Wave ◽  
Surface Waves ◽  
Love Wave ◽  
Arrival Time ◽  
Epicentral Distance ◽  
Horizontal Distance ◽  
Wave Analysis ◽  
Time Period ◽  
Phase And Group Velocities ◽  
Group Velocities

<p>The surface wave phase and group velocities are estimated by dividing the epicentral distance by phase and group travel times respectively in all the available methods, this is based on the assumptions that (1) surface waves originate at the epicentre and (2) the travel time of the particular group or phase of the surface wave is equal to its arrival time to the station minus the origin time of the causative earthquake; However, both assumptions are wrong since surface waves generate at some horizontal distance away from the epicentre. We calculated the actual horizontal distance from the focus at which they generate and assessed the errors caused in the estimation of group and phase velocities by the aforementioned assumptions in a simple isotropic single layered homogeneous half space crustal model using the example of the fundamental mode Love wave. We took the receiver locations in the epicentral distance range of 100-1000 km, as used in the regional surface wave analysis, varied the source depth from 0 to 35 Km with a step size of 5 km and did the forward modelling to calculate the arrival time of Love wave phases at each receiver location. The phase and group velocities are then estimated using the above assumptions and are compared with the actual values of the velocities given by Love wave dispersion equation. We observed that the velocities are underestimated and the errors are found to be; decreasing linearly with focal depth, decreasing inversely with the epicentral distance and increasing parabolically with the time period. We also derived empirical formulas using MATLAB curve fitting toolbox that will give percentage errors for any realistic combination of epicentral distance, time period and depths of earthquake and thickness of layer in this model. The errors are found to be more than 5% for all epicentral distances lesser than 500 km, for all focal depths and time periods indicating that it is not safe to do regional surface wave analysis for epicentral distances lesser than 500 km without incurring significant errors. To the best of our knowledge, the study is first of its kind in assessing such errors.</p>

Free spheroidal vibrations of the earth at very long periods, Part I— Calculation of periods for several earth models

1963 ◽  
Vol 53 (3) ◽  
pp. 483-501 ◽  
Author(s):  
Leonard E. Alsop
Keyword(s):  
Inner Core ◽  
Computer Programs ◽  
Free Vibrations ◽  
Stoneley Waves ◽  
The Core ◽  
The Earth ◽  
Earth Models ◽  
Shear Modes ◽  
Phase And Group Velocities ◽  
Group Velocities

Abstract Periods of free vibrations of the spheroidal type have been calculated numerically on an IBM 7090 for the fundamental and first two shear modes for periods greater than about two hundred seconds. Calculations were made for four different earth models. Phase and group velocities were also computed and are tabulated herein for the first two shear modes. The behavior of particle motions for different modes is discussed. In particular, particle motions for the two shear modes indicate that they behave in some period ranges like Stoneley waves tied to the core-mantle interface. Calculations have been made also for a model which presumes a solid inner core and will be discussed in Part II. The two computer programs which were made for these calculations are described briefly.

Eulerian partial-differential-equation methods for complex-valued eikonals in attenuating media

Geophysics ◽  
2021 ◽  
pp. 1-53
Author(s):  
Jiangtao Hu ◽  
Jianliang Qian ◽  
Jian Song ◽  
Min Ouyang ◽  
Junxing Cao ◽  
...  
Keyword(s):  
Differential Equation ◽  
Partial Differential Equation ◽  
Ray Tracing ◽  
Seismic Waves ◽  
Real Space ◽  
Advection Equation ◽  
Partial Differential ◽  
The Real ◽  
Eikonal Function ◽  
Complex Valued

Seismic waves in earth media usually undergo attenuation, causing energy losses and phase distortions. In the regime of high-frequency asymptotics, a complex-valued eikonal is an essential ingredient for describing wave propagation in attenuating media, where the real and imaginary parts of the eikonal function capture dispersion effects and amplitude attenuation of seismic waves, respectively. Conventionally, such a complex-valued eikonal is mainly computed either by tracing rays exactly in complex space or by tracing rays approximately in real space so that the resulting eikonal is distributed irregularly in real space. However, seismic data processing methods, such as prestack depth migration and tomography, usually require uniformly distributed complex-valued eikonals. Therefore, we propose a unified framework to Eulerianize several popular approximate real-space ray-tracing methods for complex-valued eikonals so that the real and imaginary parts of the eikonal function satisfy the classical real-space eikonal equation and a novel real-space advection equation, respectively, and we dub the resulting method the Eulerian partial-differential-equation method. We further develop highly efficient high-order methods to solve these two equations by using the factorization idea and the Lax-Friedrichs weighted essentially non-oscillatory (WENO) schemes. Numerical examples demonstrate that the proposed method yields highly accurate complex-valued eikonals, analogous to those from ray-tracing methods. The proposed methods can be useful for migration and tomography in attenuating media.

Sign in / Sign up

Export Citation Format

Share Document