The role of surface tension in elastic wave scattering in an inhomogeneous poroelastic medium

It is well known that many porous media such as rocks have heterogeneities at nearly all scales. We applied Biot's poroelastic theory to study the propagation of elastic waves in isotropic porous matrix with spherical inclusions. It is assumed that the heterogeneity dimension exceeds significantly the pore size. Modified boundary conditions on poroelastic interface are used to take into account the surface tension effects. The effective wavenumber is calculated using the Waterman and Truell multiple scattering theory, which relates the effective wave number to the amplitude of the wave field scattered by a single inclusion. The calculations were performed for a medium containing fluid-filled cavities or porous inclusions contrasting in saturating fluid elastic properties. The results obtained show that when we consider elastic wave propagation in poroelastic medium containing soft inclusions, it is necessary to take into account the capillary pressure. The influence of the surface tension depends on the diffraction parameter and it is a maximum in the low frequency range.     

The role of surface tension in elastic wave scattering in an inhomogeneous poroelastic medium

It is well known that many porous media such as rocks have heterogeneities at nearly all scales. We applied Biot's poroelastic theory to study the propagation of elastic waves in isotropic porous matrix with spherical inclusions. It is assumed that the heterogeneity dimension exceeds significantly the pore size. Modified boundary conditions on poroelastic interface are used to take into account the surface tension effects. The effective wavenumber is calculated using the Waterman and Truell multiple scattering theory, which relates the effective wave number to the amplitude of the wave field scattered by a single inclusion. The calculations were performed for a medium containing fluid-filled cavities or porous inclusions contrasting in saturating fluid elastic properties. The results obtained show that when we consider elastic wave propagation in poroelastic medium containing soft inclusions, it is necessary to take into account the capillary pressure. The influence of the surface tension depends on the diffraction parameter and it is a maximum in the low frequency range.     

A New Method of Solving Multimode Coupled Equations for Analysis of Uniform and Non-Uniform Fiber Bragg Gratings and Its Application to Acoustically Induced Superstructure Modulation

A new and simple mathematical formulation that is employed to analyze numerically coupled-mode equations modeling uniform and non-uniform gratings in optical fiber is investigated. This method would be straightforward and thus beneficial to solve multimode coupled equations in comparison with a previously used fundamental matrix method, and the Runge--Kutta algorithm. The new formulation proposed in this study is applied to calculate transmission and reflection spectra of core mode and higher-order cladding modes of acoustically induced superstructure modulation caused by microbending through fiber Bragg gratings (FBGs). Co-directional and contra-directional couplings based on acoustically induced modulation in FBGs have been discussed for a variety of induced coupling coefficients.     

A new class of the surface shear magnetoacoustic waves in antiferromagnetic crystals

It is shown that the hybridization of magnetoelastic and dipolar interactions may give rise (even in a zero magnetic field) to a new type of surface shear magnetoacoustic waves near a mechanically free or an acoustically continuous interface between uniaxial ferromagnetic and nonmagnetic media. The effect of the conductivity of the nonmagnetic medium on the localization conditions of this surface mode is studied.     

Poroelastic modeling of seismic boundary conditions across a fracture

Permeability of a fracture can affect how the fracture interacts with seismic waves. To examine this effect, a simple mathematical model that describes the poroelastic nature of wave-fracture interaction is useful. In this paper, a set of boundary conditions is presented which relate wave-induced particle velocity (or displacement) and stress including fluid pressure across a compliant, fluid-bearing fracture. These conditions are derived by modeling a fracture as a thin porous layer with increased compliance and finite permeability. Assuming a small layer thickness, the boundary conditions can be derived by integrating the governing equations of poroelastic wave propagation. A finite jump in the stress and velocity across a fracture is expressed as a function of the stress and velocity at the boundaries. Further simplification for a thin fracture yields a set of characteristic parameters that control the seismic response of single fractures with a wide range of mechanical and hydraulic properties. These boundary conditions have potential applications in simplifying numerical models such as finite-difference and finite-element methods to compute seismic wave scattering off nonplanar (e.g., curved and intersecting) fractures.     

Two-dimensional poroelastic acoustical foam shape design for absorption coefficient maximization by topology optimization method

Optimal shape design of a two-dimensional poroelastic acoustical foam is formulated as a topology optimization problem. For a poroelastic acoustical system consisting of an air region and a poroelastic foam region, two different physical regions are continuously changed in an iterative design process. To automatically account for the moving interfaces between two regions, we propose a new unified model to analyze the whole poroelastic acoustical foam system with one set of governing equations; Biot's equations are modified with a material property interpolation from a topology optimization method. With the unified analysis model, we carry out two-dimensional optimal shape design of a poroelastic acoustical foam by a gradient-based topology optimization setting. The specific objective is the maximization of the absorption coefficient in low and middle ranges of frequencies with different amounts of a poroelastic material. The performances of the obtained shapes are compared with those of well-known wedge shapes, and the improvement of absorption is physically interpreted.     

Dynamic response of a poroelastic stratum to moving oscillating load

The dynamic response of a poroelastic stratum subjected to moving load is studied. The governing dynamic equations for poroelastic medium are solved by using Fourier transform. The general solutions for the stresses and displacements in the transformed domain are established. Based on the general solutions, with the consideration of boundary conditions, the final expressions of stresses and displacements in physical domain are put forward for the three-dimensional single-layer medium. Some numerical solutions for the stresses, displacements and pore fluid pressure are presented and reveal that the response of a poroelastic stratum varies obviously with the moving velocity.     

Ultrasonic measurements on poroelastic slabs: Determination of reflection and transmission coefficients and processing for Biot input parameters

Ultrasonic reflection and transmission measurements on saturated, plane-parallel, poroelastic slabs are presented. A data processing technique is proposed to obtain poroelastic parameters from transmission measurements. A special experimental data acquisition and processing technique is applied to minimize the finite beam effects of the transducers. This technique yields results which can be compared with poroelastic plane-wave theory. Results of normal- and oblique-incidence measurements are presented, and transmission data are processed to yield wave speeds, sample thickness, angle of incidence, tortuosity, and permeability. The results show good agreement with independent measurements, and they are subsequently used as input for a forward modeling of the complete transmitted and reflected waveforms utilizing Biot theory. The agreement between recorded and modeled signals is good, both in time and frequency domain.     

Wave-induced fluid flow in random porous media: attenuation and dispersion of elastic waves

A detailed analysis of the relationship between elastic waves in inhomogeneous, porous media and the effect of wave-induced fluid flow is presented. Based on the results of the poroelastic first-order statistical smoothing approximation applied to Biot's equations of poroelasticity, a model for elastic wave attenuation and dispersion due to wave-induced fluid flow in 3-D randomly inhomogeneous poroelastic media is developed. Attenuation and dispersion depend on linear combinations of the spatial correlations of the fluctuating poroelastic parameters. The observed frequency dependence is typical for a relaxation phenomenon. Further, the analytic properties of attenuation and dispersion are analyzed. It is shown that the low-frequency asymptote of the attenuation coefficient of a plane compressional wave is proportional to the square of frequency. At high frequencies the attenuation coefficient becomes proportional to the square root of frequency. A comparison with the 1-D theory shows that attenuation is of the same order but slightly larger in 3-D random media. Several modeling choices of the approach including the effect of cross correlations between fluid and solid phase properties are demonstrated. The potential application of the results to real porous materials is discussed.     

Symbolic computation on pulse compression without pedestal in an inhomogeneous medium

The coupled nonlinear Schrödinger equations with the harmonic potential and variable coefficients are studied for the pulse propagation in an inhomogeneous medium. With the modified Hirota method and symbolic computation, the bilinear form and analytic one-soliton solutions are obtained. A type of pulse compression technique is proposed, which can have the optical pulses compressed without any external devices. Moreover, the compressed pulses are pedestal free. The influences of the inhomogeneity of the refractive index, Kerr nonlinearity and diffraction are analyzed as well. The proposed technique may provide a different method for the pulse compression.     

Quantum cascade structure for mid-IR four-wave mixing

The design and modeling of a quantum cascade optical amplifier (QCOA) using intra-cavity non-linear interactions to achieve wavelength conversion is proposed. The model is based on the nonlinear equation coupled with Maxwell wave equations for different emission modes. In the proposed structure, four wave mixing (FWM) output exhibits a peak as a function of pump and probe frequency if they are tuned to the energy levels of the QC structure subbands. Results demonstrate that the FWM output signal power significantly depends on how subbands are engineered and interact with optical pulses which propagate in multi layer medium. In addition, we show that by adjusting pump and probe signal frequencies, FWM output power can be tuned.     

A finite difference method for a coupled model of wave propagation in poroelastic materials

A computational method for time-domain multi-physics simulation of wave propagation in a poroelastic medium is presented. The medium is composed of an elastic matrix saturated with a Newtonian fluid, and the method operates on a digital representation of the medium where a distinct material phase and properties are specified at each volume cell. The dynamic response to an acoustic excitation is modeled mathematically with a coupled system of equations: elastic wave equation in the solid matrix and linearized Navier-Stokes equation in the fluid. Implementation of the solution is simplified by introducing a common numerical form for both solid and fluid cells and using a rotated-staggered-grid which allows stable solutions without explicitly handling the fluid-solid boundary conditions. A stability analysis is presented which can be used to select gridding and time step size as a function of material properties. The numerical results are shown to agree with the analytical solution for an idealized porous medium of periodically alternating solid and fluid layers.     

On formulation of a transition matrix for electroporoelastic medium and application to analysis of scattered electroseismic wave

On the basis of Pride's theory (1994) which couples Biot's theory for poroelastic medium (1956) and Maxwell equations via flux/force transport equations, we extend Yeh et al. (2004) approach for poroelastic medium to develop a transition matrix for electroporoelastic medium. The transition matrix, which relates the coefficients of scattered waves to those of incident waves, is then derived through the application of Betti's third identity and the associated orthogonality conditions for the electroporoelastic medium. To illustrate the application, a simple case of the scattering problem of a spherical electroporoelastic inclusion, embedded within the surrounding electroporoelastic medium subjected to an incident plane compressional wave is considered.     

Scaling,propagation, and kinetic roughening of flame fronts in random media

We introduce a model of two coupled reaction-diffusion equations to describe the dynamics and propagation of flame fronts in random media. The model incorporates heat diffusion, its dissipation, and its production through coupling to the background reactant density. We first show analytically and numerically that there is a finite critical value of the background density below which the front associated with the temperature field stops propagating. The critical exponents associated with this transition are shown to be consistent with meanfield theory of percolation. Second, we study the kinetic roughening associated with a moving planar flame front above the critical density. By numerically calculating the time-dependent width and equal-time height correlation function of the front, we demonstrate that the roughening process belongs to the universality class of the Kardar-Parisi-Zhang interface equation. Finally, we show how this interface equation can be analytically derived from our model in the limit of almost uniform background density.     

Scattering of monochromatic longitudinal waves on a planar crack of arbitrary shape in a fluid-saturated poroelastic medium

Scattering of monochromatic longitudinal waves on a planar crack of arbitrary shape in a saturated poroelastic medium is considered. The medium is described by Biot’s constitutive equations, the crack sides are fluid permeable. The problem is reduced to a two-dimensional integral equation for the crack opening vector. Gaussian approximating functions are used for discretization of this equation. For such functions, the elements of the matrix of discretized problem are combinations of four standard one-dimensional integrals that can be tabulated. As a result, numerical integration is not needed. For regular grids of approximating nodes, this matrix has Toeplitz’s structure, and matrix-vector products can be calculated by the fast Fourier transform technique. The latter accelerates substantially the process of iterative solution of the discretized problem. Calculation of crack opening vectors, differential, and total cross-sections of circular and elliptic cracks are performed for longitudinal incident waves orthogonal to the crack surfaces. Dependencies of these characteristics on the medium permeability and wavefrequency are studied. Comparison of a crack in the poroelastic medium and in a dry elastic medium with the same porosity and skeleton elastic properties is presented.     

Features of the phonon localization near the surface of a low-temperature antiferromagnet

Necessary conditions are determined under which elastic surface shear waves can exist in a crystal with magnetic long-range order, even if no account is taken of the magnetic dipole-dipole interaction, in both the case of the mechanically free surface of the crystal and the case of an acoustically continuous interface between the magnetic crystal and a nonmagnetic medium.     

Localization of Excitations near a Thin Structured Spacer between Linear and Nonlinear Crystals

It has been shown that localized and semi-localized stationary states exist near a thin structured defect layer between a linear medium and a Kerr nonlinear medium. Localized states are described by a monotonically decreasing amplitude of the field on the both sides of the interface between the media. Semilocalized states are characterized by the field that has the form of a standing wave in the linear medium and decreases monotonically in the nonlinear medium. Kerr media with self-focusing and defocusing are considered. The proposed model is described by a system of the linear and nonlinear Schrödinger equations with a specific potential simulating a thin structured defect layer. It has been shown that localized and semi-localized states exist in different energy ranges in the case of contact of the linear medium with the self-focusing medium. In the case of contact of the linear medium with the defocusing medium, two types of localized and semi-localized states differing in energy and field profile can exist in different energy ranges. In particular cases, expressions for energies of states of these types have been obtained and conditions of their applicability have been indicated.     

Coupling Biot and Navier–Stokes equations for modelling fluid–poroelastic media interaction

The interaction between a fluid and a poroelastic structure is a complex problem that couples the Navier–Stokes equations with the Biot system. The finite element approximation of this problem is involved due to the fact that both subproblems are indefinite. In this work, we first design residual-based stabilization techniques for the Biot system, motivated by the variational multiscale approach. Then, we state the monolithic Navier–Stokes/Biot system with the appropriate transmission conditions at the interface. For the solution of the coupled system, we adopt both monolithic solvers and heterogeneous domain decomposition strategies. Different domain decomposition methods are considered and their convergence is analyzed for a simplified problem. We compare the efficiency of all the methods on a test problem that exhibits a large added-mass effect, as it happens in hemodynamics applications.     

High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media

This Letter introduces a novel finite-difference time-domain (FDTD) formulation for solving transverse electromagnetic systems in dispersive media. Based on the auxiliary differential equation approach, the Debye dispersion model is coupled with Maxwell's equations to derive a supplementary ordinary differential equation for describing the regularity changes in electromagnetic fields at the dispersive interface. The resulting time-dependent jump conditions are rigorously enforced in the FDTD discretization by means of the matched interface and boundary scheme. High-order convergences are numerically achieved for the first time in the literature in the FDTD simulations of dispersive inhomogeneous media.     

Optical interferometry for measurement of Rayleigh and dilatational waves

Two-beam optical interferometry provides a highly sensitive new technique for the measurement of both Rayleigh waves and dilatational waves. A new optical arrangement overcomes previous objections to interferometry by its differential nature. The acoustically perturbed surface need not be exceptionally flat and the specimen need not be exactly aligned. The technique can be used to measure the phase and amplitude of Rayleigh waves to within a fraction of a wavelength from a discontinuity on any fully or partly reflecting surface. Dilatational waves are measureable in any transparent medium.