Time reversal super resolution in randomly layered media

In this paper we extend previous work on time reversal in randomly layered media [J.-P. Fouque, J. Garnier, A. Nachbin, K. Sølna, Time reversal refocusing for point source in randomly layered media, Wave Motion 42 (2005) 238–260]. We consider first the case of an active source embedded below the surface in a finely layered random medium. We carry out time reversal with a time reversal mirror placed at the surface and we consider here the case where this mirror is larger than the carrier wavelength. In contrast with the situation addressed in our previous paper, where the size of the mirror was comparable to the wavelength, we show that multi-pathing dramatically enhances the effective aperture of the mirror so that super resolution at the location of the source can be obtained. In other words, the focal spot radius of the refocused field obtained in the case of a multiply scattering medium is much smaller than the spot size obtained in the case of a homogeneous medium. This super resolution effect is obtained by time-reversing the long incoherent waves generated by the multiple scattering due to the thin layers. We also give an application to the problem of focusing on a passive scatterer buried in the random medium and illuminated by a source at the surface.     

Direct numerical modeling of time-reversal acoustic subwavelength focusing

Focusing waves back to their original source position is possible both experimentally and numerically thanks to time reversal mirrors (TRM). For a TRM placed in the far field of the source, the focusing spot of the reversed wavefield is subject to the diffraction limit and cannot be smaller than half the minimum wavelength, even for a very small source. Yet, numerous time reversal experiments in resonating media have shown subwavelength focusing. In this work, we show that it is possible to model these subwavelength focusing observations with simple physics, only the 2-D standard acoustic wave equation, and with specific fine scale heterogeneity. Our work is based on the spectral element method to solve the wave equation and to model time reversal experiments. Such a method makes it possible to propagate very long time series in complex and strongly discontinuous media with high accuracy. The acoustic wave equations are solved at the fine scale in media with one or more split rings of size much smaller than the wavelength. Such split rings produce a Helmholtz resonance effect as well as propagation band-gaps. We show that, in such media, even with a single split ring resonator, subwavelength focusing down to 1/13th of the minimum wavelength can be observed.     

A time-reversal method for an acoustical pulse propagating in randomly layered media

In the recent years a considerable amount of mathematical work has been devoted to the study of reflected signals obtained by the propagation of pulses in randomly layered media. We refer to [M. Asch, W. Kohler, G. Papanicolaou, M. Postel and B. White, “Frequency content of randomly scattered signals”, SIAM Review 33 (4), 519–625 (1991)] for an extensive survey and applications to inverse problems. The analysis is based on separation of scales between the correlation scale of the inhomogeneities present in the medium, the typical wavelengths of the pulse and the macroscopic variations of the medium. On the other hand, in the context of ultrasounds, time-reversal mirrors have been developed and their effects have been studied experimentally by Mathias Fink and his team at the Laboratoire Ondes et Acoustique (ESPCI-Paris). We refer to: [M. Fink, “Time reversal mirrors”, J. Phys. D: Appl. Phys. 26, 1333–1350 (1993)]. Our goal is to present a mathematical analysis of a time-reversal method for analyzing reflected signals in the model described in [M. Asch, W. Kohler, G. Papanicolaou, M. Postel and B. White, “Frequency content of randomly scattered signals”, SIAM Review 33(4), 519–625 (1991)]. We restrict our analysis to the one-dimensional case, the three-dimensional layered case being the content of a forthcoming paper. It is noticeable that we do not introduce new mathematics in the problem but simply put together an already existing mathematical theory and a new device, the time-reversal mirror.     

Statistical theory of radiative transfer in layered media

An equation for the probability density of the wave intensity which takes into account absorption, is obtained with a help of the invariant imbedding method. The limiting case when the medium occupies a half-space, is considered. The field intensity is found for the case of a source inside the medium. The conditions of applicability of the linear theory or radiative transfer are obtained. Numerical solutions of the equations corresponding to the statistical theory of radiative transfer in a layered medium with random inhomogeneities are discussed.     

Invariant imbedding method for wave problems

A new method, which reduces various boundary value problems for a wave equation to initial value problems, is developed. The scalar Helmholtz equation for the one- and three-dimensional cases is considered. The method is extended to the case of nonlinear media. Its applications to the wave equation for different dimensions and various media are described. The source of the wave field may be situated outside or inside the layer occupied by the medium. Governing equations are obtained for cases when one can neglect the backward scattering. The operators that arise are reduced to integral operators. The problem of wave scattering by a weakly rough surface is briefly considered.     

Wave dispersion in randomly layered materials

Wave scattering in materials composed of two kinds of alternating layers with different elastic properties and randomly distributed thicknesses has been modeled. The general form of the dispersion equation is derived for the unbounded layered medium. It defines two basic macroscopic characteristics of the scattered wave: phase velocity and attenuation, which are explicit functions of wave frequency and microscopic parameters of the system: acoustic properties of the layers and stochastic characteristics of their thickness distributions. The analytical expressions are derived for three special cases: for long waves; for a periodic medium composed of layers with constant thicknesses and for random medium with uniform distribution of layer thicknesses. Special attention is paid to the analysis of the frequency dependence of the wave parameters. It was shown that the predictions of the model for long waves and for periodic medium are compatible with the results obtained in the literature.Moreover, comparison of theoretical results for frequency dependent wave parameters with numerical simulations of pulse transmission through the slab of the randomly layered medium shows good qualitative and quantitative agreement in wide frequency range.     

Slip–Brinkman Flow Through Corrugated Microannulus with Stationary Random Roughness

This paper presents a boundary perturbation method of the Brinkman-extended Darcy model to investigate the flow in corrugated microannuli cylindrical tubes with slip surfaces. The stationary random model is used to mimic the surface roughness of the cylindrical walls. The tube is filled with a porous medium. We shall consider the two cases where corrugations are either perpendicular or parallel to the flow, and particular attention is given to the effect of the phase shift. The effects of the corrugations on the flow rate and pressure gradient are investigated as functions of wavelength, the permeability of the medium, the radius ratio and the slip parameter. Particular surface roughnesses are examined as special cases of stationary random surface. It is found that the effect of the partial slip is significant on the corrugation functions. The limiting cases of Stokes and Darcy’s flows and no-slip case are discussed.     

衰变热源作用下饱和多孔介质热固结问题的扩展精细积分法

热源作用下饱和多孔介质热固结效应是土木及能源工程领域的一个重要课题.由于问题的复杂性,已有的研究大多将介质假定为均匀各向同性,且将热源假定为恒定强度.实际工程中,天然饱和多孔介质常表现出明显的分层特性,热源强度也存在衰变性,为此本工作采用扩展精细积分法对衰变热源作用下层状饱和多孔介质的热固结问题进行研究.借助于积分变换,将饱和多孔介质热固结问题的偏微分方程转化为变换域内的常微分方程;然后对饱和多孔介质微层元进行合并消元,并结合边界条件,推导出衰变热源作用下层状饱和多孔介质热固结问题在积分变换域内的扩展精细积分解;对所得解答进行相应的数值积分逆变换,可获得所求温度、超静孔压及竖向位移在物理域内的解答.基于上述求解过程,编制相应的计算程序进行数值计算,通过与已有文献对比,验证本文扩展精细积分法在求解层状饱和多孔介质热固结问题中的适应性和正确性;最后通过几组算例,分析热源衰变周期、热源埋深及介质的成层性对热固结效应的影响.结果表明:热源衰变周期对温度和超静孔压的峰值、以及达到峰值的时间均有明显影响,衰变周期越长,二者峰值均越大,且达到峰值所需时间越长;热源埋深对超静孔压及竖向位移变化影响显著,深埋热源作用时热源两侧竖向位移呈对称分布,而浅埋热源两侧则无此现象;饱和多孔介质的分层特性对热固结效应影响明显.     

Three-dimensional steady-state closed form solution for multilayered fluid-saturated anisotropic finite media due to surface/internal point source

A three-dimensional(3 D)steady-state solution of fluid saturated anisotropic finite media is presented.The eigenequation method and the pseudo-Stroh formalism are used to obtain the exact solution for homogeneous saturated finite media.The propagator matrix method is introduced to deal with the corresponding multilayered poroelastic media.The poroelastic solutions due to surface or internal point fluid source are obtained.The comparison of the results of the saturated isotropic media in a half space and those obtained by the finite element method is given to illustrate the accuracy of the solution in a finite domain.Numerical solutions of a sandwich poroelastic medium are presented to analyze its hydromechanical behaviors.Two ratios of the horizontal permeability to vertical permeability and different source positions are investigated.The results show that the fluid parameters and source positions have great influence on the hydromechanical behaviors of the layered media.     

On the existence of diffusive waves in some non-linear unsteady flows of non-Newtonian fluids through porous media

We investigate the unsteady flow of power law fluids through porous media. We determine the pressure and velocity distributions when fluid is injected into a porous medium of infinite extend. We obtain solutions of progressive-wave type by means of a translation. We determine the necessary conditions for the existence of this type of solution regarding the prescribed pressure of injection and the initial pressure and velocity distributions in the porous medium. Similarity solutions are also obtained for the cases of a prescribed time dependent pressure of injection and a prescribed constant flow rate of injection. In the latter case the resulting ordinary differential equation is solved numerically. Point source solutions are also obtained for the case when an amount of fluid is instantaneously injected into the porous media. In all cases the rheological effects are presented and analyzed.     

基于积分方程方法的弹性波多重散射计算

运用积分方程方法计算了含多个随机分布椭圆柱型孔洞的随机非均匀介质中相干波的速度和衰减系数,分析了这种介质的频散特性。首先,建立了散射位移场满足的积分方程,推导了单个椭圆柱孔洞的散射截面计算公式。接着分析了在含多个随机分布椭圆柱型孔洞的随机非均匀介质中弹性波的多重散射,给出在统计平均意义下的相干波的波速和衰减系数计算公式。然后用Matlab进行了编程,给出了一个数值算例,并将计算结果与波函数展开法进行了比较,分析了随机空隙介质的频散特征及其孔洞椭圆偏心率和材料空隙率的影响。     

The penetrable coated sphere embedded in a point source excitation field

A low frequency acoustic wave field emanates from a given point and fills up the whole space. A penetrable lossy sphere with a coeccentric spherical core, which is also penetrable and lossy but characterized by different physical parameters, disturbs the given point source field. We obtain zeroth- and first-order low frequency solutions of this scattering problem in the interior of the spherical core, within the spherical shell, and in the exterior medium of propagation. We also derive the leading nonvanishing terms of the normalized scattering amplitude, the scattering cross-section as well as the absorption cross-section. The special case of a penetrable sphere is recovered either by equating the physical parameters that characterize the media in the shell and in the exterior, or by reducing the radius of the core sphere to zero. By letting the compressional viscosity of the medium in the interior sphere, or in the shell, go to zero, we obtain corresponding results for the lossless case. The incident point source field is so modified as to be able to obtain the corresponding results for plane wave incidence in the limit as the source point approaches infinity. It is observed that a small scatterer interacts stronger with a point source generated field than with a plane wave. A detailed analysis of the influence that the geometrical and the physical parameters of the problem have on the scattering process is also included. An interesting conclusion is that if the point source is located at a distance more than five radii of the scatterer away from it, then no significant changes with the plane excitation case are observed.     

Precursors for waves in random media

We consider scattering of a pulse propagating through a three-dimensional random media and study the shape of the pulse in the parabolic approximation. We show that, similarly to the one-dimensional O’Doherty–Anstey theory, the pulse undergoes a deterministic broadening. Its amplitude decays only algebraically and not exponentially in time, due to the signal low/midrange frequency component. We also argue that the parabolic approximation captures the front evolution (but not the signal away from the front) correctly even in the fully three-dimensional situation.     

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation‐based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1‐D as well as the 3‐D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1‐D problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Viscoelastic layer under the action of a moving load

Recently, problems concerning the dynamic behavior of imperfect continuous media under various types of actions have been widely investigated. The method of Laplace transformation is very convenient for describing physical processes concerning unsteady phenomena. In viscoelastic media two complications are added: the representation of the properties of a medium depending on time, and the inversion of the obtained solutions containing this additional complication. Certain approximate methods of inversion in the analysis of viscoelastic stresses are discussed in [1]. In [2, 3] a discussion is given for an effective method of constructing the solution of unsteady problems for finite and for infinite imperfect media using auxiliary functions, and a solution is presented for a half-space. Making use of the idea of the inversion of transforms, discussed in [4], in [5] a solution is obtained and a complete picture is presented for the dynamics of the variation of the stress field in a viscoelastic half-space. In the present study we consider the action of a normal moving load that is suddenly applied to the free surface of a viscoelastic layer. By Laplace and Fourier integral transformations we obtain a solution in the form of a uniformly converging series based on longitudinal and transverse waves reflected in the layer. By means of inverting the transforms by the method discussed in [4, 5], we obtain an exact solution for the stress field in the medium under investigation. We consider the special case of a viscoelastic medium of Boltzmann type, for which we obtain a numerical realization of the solution on a digital computer.     

Statistical characteristics of the motion of a fluid particle in a medium with random porosity

In the study of flow of a neutral admixture in a porous medium, it is most often assumed in the stochastic formulation that the porosity is constant and a determinate quantity, and the velocity is a random function [1–4]. The velocity distribution is usually regarded as known. Flow in a porous medium with random porosity has been studied to a far lesser extent. We note [5], which studies the averaged equations obtained within the framework of the correlation approximation. We consider the model problem of one-dimensional motion of a fluid particle (position of the front for flow of a neutral admixture in a porous medium) in a medium with random porosity. For a particular form of random porosity field, expressions are obtained for the one- and two-point densities of the distribution of the position of the particle. A study is made of the dependences of the first four moments and the correlation function of the position of the particle as functions of the time. It is shown that for large values of the time the motion of the particle is asymptotically similar to Brownian motion. It is shown by means of numerical modeling that the results obtained transfer to the case of an arbitrary random porosity field. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 59–65, November–December, 1986.     

Space–time focusing of acoustic waves on unknown scatterers

Consider a propagative medium, possibly inhomogeneous, containing some scatterers whose positions are unknown. Using an array of transmit–receive transducers, how can one generate a wave that would focus in space and time near one of the scatterers, that is, a wave whose energy would confine near the scatterer during a short time? The answer proposed in the present paper is based on the so-called DORT method (French acronym for: decomposition of the time reversal operator) which has led to numerous applications owing to the related space-focusing properties in the frequency domain, i.e., for time-harmonic waves. This method essentially consists in a singular value decomposition (SVD) of the scattering operator, that is, the operator which maps the input signals sent to the transducers to the measure of the scattered wave. By introducing a particular SVD related to the symmetry of the scattering operator, we show how to synchronize the time-harmonic signals derived from the DORT method to achieve space–time focusing. We consider the case of the scalar wave equation and we make use of an asymptotic model for small sound-soft scatterers, usually called the Foldy–Lax model. In this context, several mathematical and numerical arguments that support our idea are explored.     

Propagation of Transient Acoustic Waves in Layered Porous Media: Fractional Equations for the Scattering Operators

Acoustic waves scattering from a rigid air-saturated porous medium is studied in the time domain. The medium is one dimensional and its physical parameters are depth dependent, i.e., the medium is layered. The loss and dispersion properties of the medium are due to the fluid-structure interaction induced by wave propagation. They are modeled by generalized susceptibility functions which express the memory effects in the propagation process. The wave equation is then a fractional telegraphist’s equation. The two relevant quantities are the scattering operators—transmission and reflection operators—which give the scattered fields from the incident wave. They are obtained from Volterra equations which are fractional equations for the scattering operators.     

Overall constitutive description of symmetric layered media based on scattering of oblique SH waves

This papers investigates the scattering of oblique shear horizontal (SH) waves off finite periodic media made of elastic and viscoelastic layers. It further considers whether a Willis-type constitutive matrix (in temporal and spatial Fourier domain) may reproduce the scattering matrix (SM) of such a system. In answering this question the procedure to determine the relevant overall constitutive parameters for such a medium is presented. To do this, first the general form of the dispersion relation and impedances for oblique SH propagation in such coupled Willis-type media are developed. The band structure and scattering of layered media are calculated using the transfer matrix (TM) method. The dispersion relation may be derived based on the eigen-solutions of an infinite periodic domain. The wave impedances associated with the exterior surfaces of a finite thickness slab are extracted from the scattering of such a system. Based on reciprocity and available symmetries of the structure and each constituent layer, the general form of the dispersion and impedances may be simplified. The overall quantities may be extracted by equating the scattering data from TM with those expected from a Willis-type medium. It becomes evident that a Willis-type coupled constitutive tensor with components that are assumed independent of wave vector is unable to reproduce all oblique scattering data. Therefore, non-unique wave vector dependent formulations are introduced, whose SM matches that of the layered media exactly. It is further shown that the dependence of the overall constitutive tensors of such systems on the wave vector is not removable even at very small frequencies and incidence angles and that analytical considerations significantly limit the potential forms of the spatially dispersive constitutive tensors.     

Antiplane elastic wave cloaking using metamaterials,homogenization and hyperelasticity

We consider the problem of how to cloak objects from antiplane elastic waves using two alternative techniques. The first is the use of a layered metamaterial in the spirit of the work of Torrent and Sanchez-Dehesa (2008) who considered acoustic cloaks, motivated by homogenization theories, whilst the second is the use of a hyperelastic cloak in the spirit of the work of Parnell et al. (2012). We extend the hyperelastic cloaking theory to the case of a Mooney–Rivlin material since this is often considered to be a more realistic constitutive model of rubber-like media than the neo-Hookean case studied by Parnell et al. (2012), certainly at the deformations required to produce a significant cloaking effect. Although not perfect, the Mooney–Rivlin material appears to be a reasonable hyperelastic cloak. This is clearly encouraging for applications. We quantify the effectiveness of the various cloaks considered by plotting the scattering cross section as a function of frequency, noting that this would be zero for a perfect cloak.