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.     

Efficient Simulation of Elastic Waves Propagation in Anisotropic Media

The paper deals with finite–difference (f-d) approach to simulation of elastic waves' propagation in anisotropic elastic media with general symmetry. Any implementation of this approach claims resolution of two key problems:

• construction of an effective f-d scheme itself; we propose to use the Lebedev's scheme (LS) being a natural generalization of Virieux staggered grid scheme (VS) widely used for isotropy; we prove that LS possesses better dispersion properties in comparison with well known Rotated Staggered Grids Scheme (RSGS).
• stable domain distension. The Perfectly Matched Layers(PML) useful for isotropic problems can be unstable in the case of anisotropy. Lebedev scheme allows one to use Optimal Grids (OG) which gives a possibility to implement efficient and low cost domain distension.

(© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Reflectionless truncation of target area for axially symmetric anisotropic elasticity

This paper presents an approach to numerical simulation of seismic wave propagation in anisotropic elastic media in cylindrical coordinates by means of conservative finite difference scheme on staggered grids. An original low cost domain distension based on Optimal Grids is proposed to restrict computational domain.     

Stability and convergence of staggered Runge-Kutta schemes for semilinear wave equations

A staggered Runge-Kutta (staggered RK) scheme is a Runge-Kutta type scheme using a time staggered grid, as proposed by Ghrist et al. in 2000 [6]. Afterwards, Verwer in two papers investigated the efficiency of a scheme proposed by Ghrist et al. [6] for linear wave equations. We study stability and convergence properties of this scheme for semilinear wave equations. In particular, we prove convergence of a fully discrete scheme obtained by applying the staggered RK scheme to the MOL approximation of the equation.     

A new type of boundary treatment for wave propagation

Numerical approximation of wave propagation can be done very efficiently on uniform grids. The Yee scheme is a good example. A serious problem with uniform grids is the approximation of boundary conditions at a boundary not aligned with the grid. In this paper, boundary conditions are introduced by modifying appropriate material coefficients at a few grid points close to the embedded boundary. This procedure is applied to the Yee scheme and the resulting method is proven to be \(L^2\)-stable in one space dimension. Depending on the boundary approximation technique it is of first or second order accuracy even if the boundary is located at an arbitrary point relative to the grid. This boundary treatment is applied also to a higher order discretization resulting in a third order accurate method. All algorithms have the same staggered grid structure in the interior as well as across the boundaries for efficiency. A numerical example with the extension to two space dimensions is included.     

On energy preserving consistent boundary conditions for the Yee scheme in 2D

The Yee scheme is one of the most popular methods for electromagnetic wave propagation. A main advantage is the structured staggered grid, making it simple and efficient on modern computer architectures. A downside to this is the difficulty in approximating oblique boundaries, having to resort to staircase approximations.     

A sweeping preconditioner for Yee’s finite difference approximation of time-harmonic Maxwell’s equations

This paper is concerned with the fast iterative solution of linear systems arising from finite difference discretizations in electromagnetics. The sweeping preconditioner with moving perfectly matched layers previously developed for the Helmholtz equation is adapted for the popular Yee grid scheme for wave propagation in inhomogeneous, anisotropic media. Preliminary numerical results are presented for typical examples.     

An explicit staggered-grid method for numerical simulation of large-scale natural gas pipeline networks

We present an explicit second order staggered finite difference (FD) discretization scheme for forward simulation of natural gas transport in pipeline networks. By construction, this discretization approach guarantees that the conservation of mass condition is satisfied exactly. The mathematical model is formulated in terms of density, pressure, and mass flux variables, and as a result permits the use of a general equation of state to define the relation between the gas density and pressure for a given temperature. In a single pipe, the model represents the dynamics of the density by propagation of a non-linear wave according to a variable wave speed. We derive compatibility conditions for linking domain boundary values to enable efficient, explicit simulation of gas flows propagating through a network with pressure changes created by gas compressors. We compare our staggered grid method with an explicit operator splitting method and a lumped element scheme, and perform numerical experiments to validate the convergence order of the new discretization approach. In addition, we perform several computations to investigate the influence of non-ideal equation of state models and temperature effects on pipeline simulations with boundary conditions on various time and space scales.     

地形构造中地震波传播的非对称交错网格模拟

提出了一种新的三维空间对称交错网格差分方法,模拟地形构造中弹性波传播过程．通过具有二阶时间精度和四阶空间精度的不规则网格差分算子用来近似一阶弹性波动方程,引入附加差分公式解决非均匀交错网格的不对称问题．该方法无需在精细网格和粗糙网格间进行插值,所有网格点上的计算在同一次空间迭代中完成．使用精细不规则网格处理海底粗糙界面、 断层和空间界面等复杂几何构造, 理论分析和数值算例表明, 该方法不但节省了大量内存和计算时间, 而且具有令人满意的稳定性和精度．在模拟地形构造中地震波传播时,该方法比常规方法效率更高．     

Efficient Simulation of Wave Propagation with Implicit Finite Difference Schemes

Finite difference method is an important methodology in the approximation of waves. In this paper, we will study two implicit finite difference schemes for the simulation of waves. They are the weighted alternating direction implicit (ADI) scheme and the locally one-dimensional (LOD) scheme. The approximation errors, stability conditions, and dispersion relations for both schemes are investigated. Our analysis shows that the LOD implicit scheme has less dispersion error than that of the ADI scheme. Moreover, the unconditional stability for both schemes with arbitrary spatial accuracy is established for the first time. In order to improve computational efficiency, numerical algorithms based on message passing interface (MPI) are implemented. Numerical examples of wave propagation in a three-layer model and a standard complex model are presented. Our analysis and comparisons show that both ADI and LOD schemes are able to efficiently and accurately simulate wave propagation in complex media.     

Simulation of planar ionization wave front propagation on an unstructured adaptive grid

The paper deals with the numerical solution of a basic 2D model of the propagation of an ionization wave. The system of equations describing this propagation consists of a coupled set of reaction–diffusion-convection equations and a Poissons equation. The transport equations are solved by a finite volume method on an unstructured triangular adaptive grid. The upwind scheme and the diamond scheme are used for the discretization of the convection and diffusion fluxes, respectively. The Poisson equation is also discretized by the diamond scheme. Numerical results are presented. We deal in more detail with numerical tests of the grid adaptation technique and its influence on the numerical results. An original behavior is observed. The grid refinement is not sufficient to obtain accurate results for this particular phenomenon. Using a second order scheme for convection is necessary.     

High-order difference methods for waves in discontinuous media

Second-, fourth- and sixth-order one-step methods have been constructed for the solution of wave propagation problems. The method is based on the first-order system formulation of the wave equation, and uses a staggered grid both in space and time. The method is applied with good results to a problem with discontinuous coefficients without using any special procedure across the discontinuity. The behavior of the truncation error has been investigated for one-dimensional problems and stability criteria have been derived for one- and two-dimensional cases.     

An explicit fourth-order staggered finite-difference time-domain method for Maxwell's equations

A new explicit fourth-order accurate staggered finite-difference time-domain (FDTD) scheme is proposed and applied to electromagnetic wave problems. It is fourth-order accurate in both space and time, conditionally stable, and highly efficient (with respect to Yee's scheme) and still retains much of the original simplicity of Yee's scheme. Both extension to perfectly matched layers and modification to deal with dielectric interfaces and perfectly conducting boundaries of the scheme have also been presented. Numerical examples are shown to illustrate the efficiency of the method.     

A short note on the asymptotic stability of an oscillation-free eikonal splitting method

Different beam propagation methods (BPMs) have been fundamental in modern electromagnetical wave simulations. Challenges of the numerical strategy include the computational efficiency and stability, in particular when highly oscillatory optical waves are present. This paper concerns an eikonal splitting BPM scheme for two-dimensional paraxial Helmholtz equations together with transparent boundary conditions in slowly varying envelope approximations of active laser beams. It is shown that the finite difference method investigated is not only oscillation-free, but also asymptotically stable. This ensures the high efficiency and applicability in highly oscillatory wave applications.     

Development and analysis of Crank‐Nicolson scheme for metamaterial Maxwell's equations on nonuniform rectangular grids

In this paper, we develop a fully implicit scheme on staggered grids to solve the Maxwell's equations when Drude metamaterial is involved. Unconditional stability and optimal error estimate of the scheme are proved. Numerical results are provided to support the theoretical analysis, and used to demonstrate the applicability of the scheme to simulate the complicated backward wave propagation phenomenon occurring in metamaterials.     

Microscale Investigations of Highfrequency Wave Propagation Through Highly Porous Media

Wave propagation in highly porous materials has a well established theoretical background. Still there are parameters which require complex laboratory experimentation in order to find numerical values. This paper presents an effective method to calculate the tortuosity of aluminum foam using numerical simulations. The work flow begins with the acquisition of the foam geometry by means of a micro-CT scanner and further image segmentation and analysis. The elastodynamic wave propagation equation is solved using a velocity-stress rotated staggered finite-difference technique. The effective wave velocities are calculated and using the fluid and, aluminum effective properties, the tortuosity is determined. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Sur la magnétohydrodynamique anisotrope relativiste

Summary This paper is concerned with the relativistic phenomenological theory of anisotropic magnetohydrodynamics. An anisotropic fluid scheme is defined and studied. The main system of anisotropic magnetohydroldynamics is deduced. This system may describe a collisionless anisotropic plasma embedded in a strong magnetic field. The main system is shown to yield to three types of waves as in isotropic (perfect) magnetohydrodynamics: the entropic waves, the magnetosonic waves and the Alfven waves. For the rays associated respectively to the magnetosonic and Alfven waves the fundamental property concerning the propagation of infinitesimal discontinuities of variables is established. The conditions under which the velocities of propagation of magnetosonic and Alfven waves are real are derived: these conditions imply as in the classical theory the absence of fire hose and mirror instabilities in the fluid. The study of wave cones allows, on the one hand to point out some particularities of the propagation of waves in anisotropic magnetohydrodynamics, and on the other hand to clear up the hyperbolicity character of differential operators associated to various waves.

---

---

Entrata in Redazione il 23 aprile 1975.     

用CT-VOF法在数值波试验槽中生成线性和非线性波

为渡水槽中波的模拟和传播提出了二维的数值模型．假设流动的流体为粘性、不可压缩的,并将Navier-Stokes方程和连续性方程作为控制方程．用标准的k-ε模型来模拟紊流流动;用交错网格的有限差分法,离散化Navier-Stokes方程;并用简化的标识和单元(SMAC)方法进行求解．使用活塞型波发生器生成并传播波;数值渡水槽的端部采用敞开式的边界条件．为了证明模型的有效性,进行了一些标准的试验,如顶盖驱动的方腔测试试验、单向的常速度场试验以及干燥河床上的溃坝试验．为了论证方法的性能及其精度,将所生成波的结果与已有波理论的结果进行比较．最后,采用群集技术(CT)生成网格,并提出最佳的网格生成条件．     

The Well-Posedness of Dynamical Equations of Magneto-electro-elasticity

A mathematical model of wave propagation in magneto-electro-elastic materials is obtained in the form of a symmetric hyperbolic system of the first-order partial differential equations. This model is a result of the qualitative analysis of the coupled time-dependent Maxwell’s equations and equations of elastodynamics which are considered together with constitutive relations in non-homogeneous anisotropic magneto-electro-elastic materials. Applying the theory and methods of symmetric hyperbolic systems, we have proved that the reported model of wave propagation in magneto-electro-elastic materials satisfies the Hadamards correctness requirements: solvability, uniqueness and stability with respect to perturbation of data.     

Numerical model for the shallow water equations on a curvilinear grid with the preservation of the Bernoulli integral

Methodological aspects concerning the construction of a two-dimensional numerical model for reservoir flows based on the shallow water equations are considered. A numerical scheme is constructed by applying the control volume method on staggered grids in combination with the Bernoulli integral, which is used to interpolate the desired fields inside a grid cell. The implementation of the method yields a monotone numerical scheme. The results of numerical integration are compared with the exact solution.