Ray-centred coordinate systems in anisotropic media

Whereas the ray-centred coordinates for isotropic media by Popov and Pšenčík are uniquely defined by the selection of the basis vectors at one point along the ray, there is considerable freedom in selecting the ray-centred coordinates for anisotropic media. We describe the properties common to all ray-centred coordinate systems for anisotropic media and general conditions, which may be imposed on the basis vectors. We then discuss six different particular choices of ray-centred coordinates in an anisotropic medium. This overview may be useful in choosing the ray-centred coordinates best suited for a particular application. The equations are derived for a general homogeneous Hamiltonian of an arbitrary degree and are thus applicable both to the anisotropic-ray-theory rays and anisotropic common S-wave rays.     

Physical splitting of nonlinear effects in high-velocity stable flow through porous media

For a high-velocity stable flow through a periodic corrugated channel representing an element of porous medium, we suggest splitting the overall nonlinear macroscopic effects into two kinds of different physical origin: a pure inertia effect produced by the convective term of Navier–Stokes equations and an inertia–viscous cross effect representing a variation of the viscous dissipation due to a streamline deformation by inertia forces. We will show that the inertia–viscous cross effects may be revealed by simulating a periodic flow, whilst the pure inertia effects are produced by the microscale flow nonperiodicity. We will reveal the individual flow law for each nonlinear component and analyze the relative role of both components numerically by using the finite element method applied to the Navier–Stokes equations. Both the pure inertia and the inertia–viscous cross effects are revealed to be exponential prior to quadratic or cubic ones. The influence of the dead volume is analyzed. The inertia–viscous cross phenomena are shown to be negligible when the flow structure is clearly nonperiodic.     

Numerical simulation of three-dimensional saltwater–freshwater fingering instabilities observed in a porous medium

We consider saltwater–freshwater fingering instabilities in a saturated porous medium. In the first part, we present three-dimensional results obtained from a laboratory experiment using non-invasive imaging. In the second part, we define a set of model problems in which the performed laboratory experiments can be ranged in. Due to its highly non-linear behavior and inevitable modeling errors, a detailed numerical reproduction of the physical concentration measurements cannot be expected. Nevertheless, four criteria have been identified, two quantitative and two qualitative, which facilitate a substantiated comparison of the physical experiment and the numerical simulation. With respect to these criteria a high degree of similarity could be observed. The use of these features allows a deeper understanding of the physical processes and the influence of the initial conditions.     

Effects of mechanical layering on volcano deformation

The migration and accumulation of magma beneath volcanoes often causes surface displacements that can be measured by geodetic techniques. Usually, deformation signals are explained using models with uniform mechanical properties. In this paper, we study surface displacements due to magma chamber inflation, using heterogeneous finite element models. We first present a systematic analysis of the influence of mechanical layering, showing that the stiffness contrast significantly affects the entity and the pattern of vertical and radial displacements. Second, as an example we apply the models to interpret ground displacements at Darwin volcano (Galápagos Islands) as revealed by InSAR data in the period 1992–1998. The considered models suggest that geodetic data interpreted using homogeneous models leads to underestimation of the source depth and volume change. Thus, we propose correction factors for the source parameters estimated by homogeneous models, in order to consider a range of variation due to mechanical layering as analysed in this study. The effect of the mechanical heterogeneities affects the correct understanding of geodetic data and also influences the evaluation of a volcanic hazard potential.     

First-order P-wave ray synthetic seismograms in inhomogeneous weakly anisotropic media

We propose approximate equations for P -wave ray theory Green's function for smooth inhomogeneous weakly anisotropic media. Equations are based on perturbation theory, in which deviations of anisotropy from isotropy are considered to be the first-order quantities. For evaluation of the approximate Green's function, earlier derived first-order ray tracing equations and in this paper derived first-order dynamic ray tracing equations are used.
The first-order ray theory P -wave Green's function for inhomogeneous, weakly anisotropic media of arbitrary symmetry depends, at most, on 15 weak-anisotropy parameters. For anisotropic media of higher-symmetry than monoclinic, all equations involved differ only slightly from the corresponding equations for isotropic media. For vanishing anisotropy, the equations reduce to equations for computation of standard ray theory Green's function for isotropic media. These properties make the proposed approximate Green's function an easy and natural substitute of traditional Green's function for isotropic media.
Numerical tests for configuration and models used in seismic prospecting indicate negligible dependence of accuracy of the approximate Green's function on inhomogeneity of the medium. Accuracy depends more strongly on strength of anisotropy in general and on angular variation of phase velocity due to anisotropy in particular. For example, for anisotropy of about 8 per cent, considered in the examples presented, the relative errors of the geometrical spreading are usually under 1 per cent; for anisotropy of about 20 per cent, however, they may locally reach as much as 20 per cent.     

Time-averaged and time-dependent energy-related quantities of harmonic waves in inhomogeneous viscoelastic anisotropic media

The energy–flux vector and other energy-related quantities play an important role in various wave propagation problems. In acoustics and seismology, the main attention has been devoted to the time-averaged energy flux of time-harmonic wavefields propagating in non-dissipative, isotropic and anisotropic media. In this paper, we investigate the energy–flux vector and other energy-related quantities of wavefields propagating in inhomogeneous anisotropic viscoelastic media. These quantities satisfy energy-balance equations, which have, as we show, formally different forms for real-valued wavefields with arbitrary time dependence and for time-harmonic wavefields. In case of time-harmonic wavefields, we study both time-averaged and time-dependent constituents of the energy-related quantities. We show that the energy-balance equations for time-harmonic wavefields can be obtained in two different ways. First, using real-valued wavefields satisfying the real-valued equation of motion and stress–strain relation. Second, using complex-valued wavefields satisfying the complex-valued equation of motion and stress–strain relation. The former approach yields simple results only for particularly simple viscoelastic models, such as the Kelvin–Voigt model. The latter approach is considerably more general and can be applied to viscoelastic models of unrestricted anisotropy and viscoelasticity. Both approaches, when applied to the Kelvin–Voigt viscoelastic model, yield the same expressions for the time-averaged and time-dependent constituents of all energy-related quantities and the same energy-balance equations. This indicates that the approach based on complex-valued representation of the wavefield may be used for time harmonic waves quite universally. This study also shows importance of joint consideration of time-averaged and time-dependent constituents of the energy-related quantities in some applications.     

Quasi-static deformation of a poroelastic half-space with anisotropic permeability by two-dimensional surface loads

An analytical solution is obtained of the fully coupled diffusion–deformation system of equations governing the quasi-static plane strain deformation of a poroelastic half-space with anisotropic permeability and compressible constituents. The stresses and the pore pressure are taken as the basic state variables. Displacements are obtained by integrating the coupled constitutive relations. The problem of surface loads is discussed in detail. Explicit analytical solutions are derived for normal line loading, shear line loading and normal strip loading. The permeability anisotropy is found to have a significant effect on the quasi-static deformation of the half-space. However, in the drained and undrained limits, the anisotropy has no effect. The stresses in the drained and undrained states are independent of the poroelastic parameters. Numerical computation of the pore pressure indicates that ignoring permeability anisotropy may lead to an overestimation of the pore pressure at points vertically below the point of normal loading. Further, anisotropy in permeability may lead to a dilution in the theoretical prediction of the Mandel–Cryer Effect.