Isotropic Medium Research Articles

Elastic reverse time migration based on nested triangular mesh in 2D isotropic media

Abstract Reverse time migration (RTM) based on two-way wave equation is an important imaging tool, which has been widely used in exploration geophysics. It offers distinct advantages in handling complex geological structures. Moreover, elastic reverse time migration (ERTM) can produce more physically meaningful images for subsurface and has received increasing attention in recent years. However, rugged topography presents challenges in the reconstruction of subsurface structures for ERTM. To overcome this problem, we introduce the grid method for topography ERTM, which is an unstructured mesh based algorithm. However, the conventional gird method for RTM is implementing by discretizing model into small meshes, leading to high computational cost. To address this problem, we introduce the grid method based on nested unstructured meshes to implement ERTM with complex surface topography. Low-order small meshes, due to their fine scale, can accurately depict complex surface topography. However, their generation process is relatively time-consuming. In contrast, high-order large meshes contain many similar small meshes internally, which can be processed in parallel, resulting in higher computational efficiency. Therefore, we propose the following strategy: for the receiver wavefield, several layers on surface discreted by the small meshes are used to place receivers and remaining layers discreted by the large meshes to improve computational efficiency. For the source wavefield, only the large meshes are used because the spacing of shot is relatively large. In this way, we ensure the accurate description of complex surface topography for the implementation of ERTM with high computational efficiency. Finally, the effectiveness of our method is validated through experiments on three models with complex topography.

Negative reflection and total internal reflection at the internal surface of lithium niobate crystal

In this paper we analytically find out the expression for the wave vector and Poynting vector in the reflection domain for negative uniaxial material. We investigate the impact of the optical axis rotation and incident angle on the wave vector and Poynting vector reflected from the internal surface of lithium niobate crystal when the light is propagating from the denser anisotropic medium to the rarer isotropic medium. We derive the Brewster angle condition and discuss the walk-off between TE and TM polarized light. We explore where the wave vector and Poynting vector are parallel after reflection. In further analysis, we look into the phase differences between ordinary and extraordinary reflected light after total internal reflection from the internal surface of lithium niobate crystal. The necessary conditions for linear polarization and circular polarization have been established. The impact of optical axis orientation is also taken into consideration. The refractive index for the isotropic material for the desired polarization response has been found as 2.18, 2.1, and 2.03 at 632 nm, 1550 nm, and 3500 nm, respectively.

A shortest‐path‐aided fast‐sweeping method to improve the accuracy of traveltime calculation in vertically transverse isotropic media

AbstractThe high accuracy and efficiency of traveltime calculation are critical in seismic tomography, migration, static corrections, source locations and anisotropic parameter estimation. The fast‐sweeping method is an efficient upwind finite‐difference approach for solving the eikonal equation. However, the fast‐sweeping method is accurate only along the axis directions. In two‐dimensional or higher dimensional cases, the accuracy is severely decreased in the diagonal directions due to the numerical errors in these directions. These similar numerical errors also arose in higher order fast‐sweeping method and anisotropic fast‐sweeping method. To improve the accuracy of traveltime calculation in two‐dimensional or higher dimensional space, a shortest‐path‐aided fast‐sweeping method is proposed. The shortest‐path‐aided solution is embedded into the sweeping process of the standard fast‐sweeping method to improve the traveltime accuracy in the diagonal directions. Shortest‐path‐aided fast‐sweeping method is very easy to implement nearly without additional computational cost and memory consumption. Furthermore, this method is easy to extend from two‐dimensional to higher dimensional, from low‐order to higher‐order and from isotropic to anisotropic cases.

On uniqueness and Somigliana theorem for a thermoelastic theory

The main concern of the present article is based on unified generalized thermoelasticity in an isotropic and homogeneous medium. The uniqueness theorem and Betti type reciprocal principle are derived. Based on the reciprocal theorem, the generalization of the Somigliana and Green theorem is also established in the present context.

Towards the understanding of convective dissolution in confined porous media: thin bead pack experiments, two-dimensional direct numerical simulations and physical models

We consider the process of convective dissolution in a homogeneous and isotropic porous medium. The flow is unstable due to the presence of a solute that induces a density difference responsible for driving the flow. The mixing dynamics is thus driven by a Rayleigh–Taylor instability at the pore scale. We investigate the flow at the scale of the pores using Hele-Shaw type experiment with bead packs, two-dimensional direct numerical simulations and physical models. Experiments and simulations have been specifically designed to mimic the same flow conditions, namely matching porosities, high Schmidt numbers and linear dependency of fluid density with solute concentration. In addition, the solid obstacles of the medium are impermeable to fluid and solute. We characterise the evolution of the flow via the mixing length, which quantifies the extension of the mixing region and grows linearly in time. The flow structure, analysed via the centreline mean wavelength, is observed to grow in agreement with theoretical predictions. Finally, we analyse the dissolution dynamics of the system, quantified through the mean scalar dissipation, and three mixing regimes are observed. Initially, the evolution is controlled by diffusion, which produces solute mixing across the initial horizontal interface. Then, when the interfacial diffusive layer is sufficiently thick, it becomes unstable, forming finger-like structures and driving the system into a convection-dominated phase. Finally, when the fingers have grown sufficiently to touch the horizontal boundaries of the domain, the mixing reduces dramatically due to the absence of fresh unmixed fluid. With the aid of simple physical models, we explain the physics of the results obtained numerically and experimentally. The solute evolution presents a self-similar behaviour, and it is controlled by different length scales in each stage of the mixing process, namely the length scale of diffusion, the pore size and the domain height.

To the theory of phase transition of a binary solution into an inhomogeneous phase

ABSTRACT In the framework of the theoretical model of the phase transition of binary solutions into spatially inhomogeneous states proposed earlier by the authors [Poluektov Y, Soroka AA. Transition of a binary solution into an inhomogeneous phase. Phase Transition. 2022;95:267–280], which takes into account nonlinear effects, the role of the cubic in concentration term in the expansion of free energy was studied. It is shown that taking into account the cubic term contributes to the stabilization of a homogeneous state. The existence of two inhomogeneous phases in an isotropic medium, considered in [Poluektov Y, Soroka AA. Transition of a binary solution into an inhomogeneous phase. Phase Transition. 2022;95:267–280], proves to be possible only at half the concentration of the solution. The contribution of inhomogeneity effects to thermodynamic quantities is calculated. Phase transitions from a homogeneous state and between inhomogeneous phases are second-order phase transitions.

Deformation and Plasticity of Reinforced Concrete

Statement of the problem.Development of the model of plasticity of reinforced concrete structures is based on the intensity of the "stress-strain" relation to project tensors using transitions, as well as the coefficients of the steps of forces and axes (principal angle and strain-strain, total shear strain, etc.). Results. The dependences of plasticity modulus of concrete, transverse coefficient have been obtained. Complex functions of generalized hypotheses of deplanation of linear and angular deformations and jumps for crack formation breakaway and stiffness are constructed. It includes force flows, i.e., blocks of compressed and tensile concrete for isotropic medium (the first object), "main cracks" from fracture mechanics, functional and two-cone element for the deformation effect of reinforced concrete. The resistance of the tensile concrete is transferred to the working reinforcement using the total average longitudinal and transverse forces for the transversal-isotropic environment, including the average modulus of reinforcement and the reduced tensile concrete coefficient of the second limit state group. The special dowel effect is derived from the second level model of structural mechanics for a rebar with two pinched ends, as well as crack opening (a development of the Thomas –– author hypothesis), crack-track bank shear, and buckling. The main vector of forces in the reinforcement is characterized by two values of longitudinal and transverse displacements and the width of opening and shear of the crack-track (the third object). Conclusions. The energetic theory on the surface of the sphere and the determination of the integral of the mean square value of the tangential stresses of plasticity theory allowed to develop an alternative to the ge-neral model in the form of a paraboloid by summing the volume sectors, radius-levels from the matrix of sliding planes (including octahedral and pure shear) in concrete medium from microcracks to macrocracks.

An Enhanced Vacuum-Assisted Resin Transfer Molding Process and Its Pressure Effect on Resin Infusion Behavior and Composite Material Performance.

In this paper, an enhanced VARTM process is proposed and its pressure effect on resin infusion behavior and composite material performance is studied to reveal the control mechanism of the fiber volume fraction and void content. The molding is vacuumized during the resin injection stage while it is pressurized during the mold filling and curing stages via a VARTM pressure control system designed in this paper. Theoretical calculations and simulation methods are used to reveal the resin's in-plane, transverse, and three-dimensional flow patterns in multi-layer media. For typical thin-walled components, the infiltration behavior of resin in isotropic porous media is studied, elucidating the control mechanisms of fiber volume fraction and void content. The experiments demonstrate that the enhanced VARTM process significantly improves mold filling efficiency and composite's performance. Compared to the regular VARTM process, the panel thickness is reduced by 4% from 1.7 mm, the average tensile strength is increased by 7.3% to 760 MPa, the average flexural strength remains at approximately 720 MPa, porosity is decreased from 1.5% to below 1%, and the fiber volume fraction is increased from 55% to 62%.

Characterization of Rectangular Waveguides Loaded E-Plane Dielectrics Using the Newton-Raphson Method

Homogeneous metallic waveguides have long been used to carry high powers. They are often filled with inhomogeneous, isotropic dielectrics to reduce their size and cut-off frequencies. To characterize these inhomogeneous rectangular waveguides made of homogeneous and isotropic media, the Newton-Raphson method is used in this article. Frequency of cutoff, attenuation, and power flow distribution are all properties of the EM wave that are highly dependent on the physical structure and composition within the guide. This article presents the characterization of an inhomogeneous and isotropic rectangular guide. The analysis of this type of guide is based on the Borgnis potential method for determining the components of the electric field E and the magnetic field H, to obtain the guide's dispersion equations. The modes that were found to exist in these waveguides are hybrid, meaning that they have both axial E- and H-fields. Numerical resolution of these equations using the Newton-Raphson method obtains the guide's propagation constants. A MATLAB program is used to plot these dispersion curves. The propagation constant increases as a function of frequency, and the d/a ratio influences the dispersion curves. Increasing the relative permittivity of the dielectric leads to an increase in the ratio of the propagation constant in the z direction to the wave number.

Radially transverse isotropic inclusions in isotropic elastic media: Local fields, neutral inclusions, effective elastic properties

Radially transverse isotropic inclusions in homogeneous isotropic elastic host media are considered. Mellin transform method is used for solution of the volume integral equation of the problem for an isolated inclusion subjected to a constant external stress (strain) field. The tensor structure of the solution is revealed with precision to three scalar functions of the radial coordinate, and the system of ordinary differential equations for these functions is derived. For multilayered radially transverse isotropic inclusions with constant elastic coefficients inside layers, explicit solution of these equations is obtained. An efficient numerical algorithm of solution for inclusions with an arbitrary number of the layers is proposed. Neutral inclusions that do not disturb homogeneous external fields applied to the medium are considered. It is shown that an inclusion with an isotropic core and radially transverse isotropic external layer can be weak neutral by appropriate choice of the layer elastic constants. The effective field method is used for determination of the effective elastic stiffness tensor of a homogeneous isotropic medium containing a random set of radially transverse isotropic inclusions.

Exact closed forms for the transmittance of electromagnetic waves in one-dimensional anisotropic periodic media

In this work, we obtain closed expressions for the transfer matrix and the transmittance of electromagnetic waves propagating in finite one-dimensional anisotropic periodic stratified media with an arbitrary number of cells. By invoking the Cayley–Hamilton theorem on the transfer matrix for the electromagnetic field in a periodic stratified media formed by N cells, we obtain a fourth-order recursive relation for the matrix coefficients that defines the so-called Tetranacci polynomials (TPs). In the symmetric case, corresponding to a unit-cell transfer matrix with a characteristic polynomial where the coefficients of the linear and cubic terms are equal, closed expressions for the solutions to the recursive relation, known as symmetric TPs, have recently been derived, allowing us to write the transfer matrix and transmittance in a closed form. We show as sufficient conditions that the 4×4 differential propagation matrix of each layer in the binary unit cell, Δ , a) has eigenvalues of the form ±p1 , ±p2 , with p1≠p2 , and b) its off-diagonal 2×2 block matrices possess the same symmetric structure in both layers. Otherwise, the recursive relations are still solvable for any 4×4 matrix and provide an algorithm to compute the Nth power of the transfer matrix without carrying out explicitly the matrix multiplication of N matrices. We obtain analytical expressions for the dispersion relation and transmittance, in closed form, for two finite periodic systems: the first one consists of two birefringent uniaxial media with their optical axis perpendicular to the z-axis, and the second consists of two isotropic media subject to an external magnetic field oriented along the z-axis and exhibiting the Faraday effect. Our formalism applies also to lossy media, magnetic anisotropy or optical activity.

A new highly accurate and efficient pure visco‐acoustic wave equation for tilted transversely isotropic attenuating media

AbstractThe propagation of seismic waves in attenuating anisotropic media exhibits amplitude dissipation and phase dispersion. To describe its effects, the fractional Laplacian pure visco‐acoustic wave equations capable of producing stable and noise‐free wavefields have been derived. However, except for acoustic approximation, previous wave equations utilize the approximations with lower accuracy in simplifying the denominator of the approximate complex‐valued dispersion relation, resulting in reduced accuracy. To address this concern, we use a combination of complex stiffness coefficients to replace the denominator term of the approximate complex‐valued dispersion relation. This approximation effectively reduces the loss of accuracy caused by ignoring the influence of the velocity anisotropy parameter ε and the attenuation anisotropy parameter εQ in the denominator term, leading to a wave equation with high accuracy in media with large anisotropic parameters ε and δ. In addition, the new wave equation only contains two high‐order spatial partial derivatives and has high computational efficiency. Theoretical analysis and numerical examples demonstrate that the proposed pure visco‐acoustic tilted transversely isotropic wave equation outperforms the previous pure visco‐acoustic wave equation in terms of simulation accuracy. The newly developed wave equation is well suited for the application of Q‐compensated reverse time migration and full waveform inversion in attenuating anisotropic media.

Spin-wave self-imaging: Experimental and numerical demonstration of caustic and Talbot-like diffraction patterns

Extending the scope of the self-imaging phenomenon, traditionally associated with linear optics, to the domain of magnonics, this study presents the experimental demonstration and numerical analysis of spin-wave (SW) self-imaging in an in-plane magnetized yttrium iron garnet film. We explore this phenomenon using a setup in which a plane SW passes through a diffraction grating, and the resulting interference pattern is detected using Brillouin light scattering. We have varied the frequencies of the source dynamic magnetic field to discern the influence of the anisotropic dispersion relation and the caustic effect on the analyzed phenomenon. We found that at low frequencies and diffraction fields, the caustics determine the interference pattern. However, at large distances from the grating, when the waves of high diffraction order and number of slits contribute to the interference pattern, the self-imaging phenomenon and Talbot-like patterns are formed. This methodological approach not only sheds light on the behavior of SW interference under different conditions but also enhances our understanding of the SW self-imaging process in both isotropic and anisotropic media.

Linear and Nonlinear Optics of Broad-Band Laser Pulses: Diffraction.

The typical spectrally limited laser pulse in the near-infrared region is narrow-band up to 40-50 fs. Its spectral width Δk is much smaller than the carrying wavenumber k0 (Δk ≪ k0) . For such kinds of pulses, on distances of a few diffraction lengths, the diffraction is of a Fresnel's type and their evolution can be described correctly in the frame of the well-known paraxial evolution equation. The technology established in 1985 of amplification through chirping of laser pulses triggered remarkable progress in laser optics along with the construction of femtosecond (fs) laser facilities producing high intensity fields of the order of 1015-1021 W/cm2. However, the duration of the pulse was quickly shortened from picoseconds down to 5-6 fs, which have a broad-band nature (Δk ∼ k0). The linear and nonlinear propagation dynamics of broad-band pulses is quite different form their narrow-band counterparts. Here, we review the appropriate theoretical approach to study the evolution of the pulse. Moreover, we shed light on the different diffraction regimes inherent to both narrow-band and broad-band laser pulses and compare them to unveil the main differences. Using this very method, in subsequent papers, we will investigate the influence of the dispersion and nonlinearity on the laser pulse propagation in isotropic media.

Truncated Gauss-Newton full-waveform inversion of pure quasi-P waves in vertical transverse isotropic media

Full-waveform inversion (FWI) uses the full information of seismic data to obtain a quantitative estimation of subsurface physical parameters. Anisotropic FWI has the potential to recover high-resolution velocity and anisotropy parameter models, which are critical for imaging the long-offset and wide-azimuth data. We develop an acoustic anisotropic FWI method based on a simplified pure quasi P-wave (qP-wave) equation, which can be solved efficiently and is beneficial for the subsequent inversion. Using the inverse Hessian operator to precondition the functional gradients helps to reduce the parameter tradeoff in the multi-parameter inversion. To balance the accuracy and efficiency, we extend the truncated Gauss-Newton (TGN) method into FWI of pure qP-waves in vertical transverse isotropic (VTI) media. The inversion is performed in a nested way: a linear inner loop and a nonlinear outer loop. We derive the formulation of Hessian-vector products for pure qP-waves in VTI media based on the Lagrange multiplier method and compute the model update by solving a Gauss-Newton linear system via a matrix-free conjugate gradient method. A suitable preconditioner and the Eisenstat and Walker stopping criterion for the inner iterations are used to accelerate the convergence and avoid prohibitive computational cost. We test the proposed FWI method on several synthetic data sets. Inversion results reveal that the pure acoustic VTI FWI exhibits greater accuracy than the conventional pseudoacoustic VTI FWI. Additionally, the TGN method proves effective in mitigating the parameter crosstalk and increasing the accuracy of anisotropy parameters.

Inversion of fracture weakness parameters based on 3C vertical seismic profile data — Part I: Horizontal transversely isotropic media composed of a single fracture set in the isotropic background media

Natural fractures in oil and gas reservoirs are a crucial factor that cannot be ignored because they significantly influence the reservoir’s petrophysical properties and hydrocarbon development. A horizontal transversely isotropic (HTI) medium composed of a single fracture set in an isotropic background is a typical anisotropic medium. Meanwhile, the S-wave splitting is a sensitive response of such anisotropic media, resulting in the generation of fast and slow S waves. The normal and tangential fracture weaknesses are crucial parameters that characterize the anisotropy of fractured media. We develop an inversion method for fracture weakness based on 3C vertical seismic profiling (3CVSP) data. First, assuming weak anisotropy and an HTI medium containing a single fracture set, we derive a first-order linear approximation of the traveltimes of the converted fast and slow S waves (PS1 and PS2 waves) with respect to fracture weakness parameters in the phase velocity domain. By solving for the horizontal projection of the slowness vector, approximate equations of the traveltimes of the PS1 and PS2 waves are converted from the phase velocity domain to the group velocity domain. Furthermore, we devise an inversion workflow consisting of three primary steps: (1) preprocessing the VSP data to derive the traveltimes and azimuth of the HTI medium, (2) constructing a forward model with undetermined fracture weakness parameters, and (3) following the establishment of the objective function, conducting the inversion for the fracture weakness parameters. We demonstrate the reliability of the method through numerical examples and synthetic 3CVSP data. The inversion errors are primarily influenced by the azimuth angle, with minimal influence from the receiver depth. Furthermore, a collective set of inverted results derived from all geophones is more stable and accurate than individual geophones. The application to actual 3CVSP data further confirms the effectiveness of our approach.

Thermal wave reflection and propagation in temperature rate-dependent non-local diffusive solids

The propagation and reflection of thermo-elastic waves through diffusive nonlocal isotropic medium had been studied in this paper. The Green–Lindsay model of thermo-elasticity is incorporated in the context of Temperature Rate Dependent Theory. Using Helmholtz vector decomposition rule, the system of governing equations has been transformed to their respective components. The dispersion relation in frequency indicates the existence of three coupled waves and one independent wave propagating through the medium. The coupled waves are affected by non-locality, temperature field and diffusivity in the medium; anyhow, un-coupled shear vertical wave is only affected by the non-local parameter. The reflection of [Formula: see text]-wave is also studied at the free boundary of the solid and their corresponding amplitude ratios are computed using set of suitable boundary conditions. The obtained results are further discussed graphically for significant physical parameters of interest. The results in the literature are obtained as a special case after ignoring the diffusivity in the solid.

Experimental and theoretical analysis of finely bedded porous rock models with internal cracks

The behavior of ultrasonic waves in synthetic rocks consisting of alternating layers of sand-based (porous and permeable) and acrylic (unporous) materials is investigated. Two synthetic layered media are created, one without cracks and the other with parallel cracks in one of the layers, to increase anisotropy and simulate a heterogeneous transversely isotropic medium with a vertical axis of symmetry. Microcomputed tomography is used to obtain images of the samples, which indicate that there are no bubbles present and that the layers are well connected without any gaps or breaks. The experiment uses ultrasonic measurements to determine the propagation of P and S waves, with transducers at frequencies of 100 kHz, 500 kHz, and 1 MHz for traveltime and waveform measurements. These three sources are chosen to verify the effective behavior of the samples, which depends on the relationship between the dominant wavelength of investigation and the layer thickness ([Formula: see text]). The measured velocities are compared with theoretical models based on Backus, Kennett (traditional models), and Backus (modified for cracked media by Schoenberg and Muir’s theory). The main finding is that the cracked medium exhibits at least twice as large a difference in behavior for P- and S-wave vertical velocities as a function of the number of layers compared with the uncracked medium.

The effect of interstellar medium on LVK’s black holes

ABSTRACT Gravitational radiation alone is not efficient in hardening the orbit of a wide binary black hole (BBH). By employing a toy model for the interstellar medium (ISM) surrounding BBHs, here we discuss the effect of this baryonic medium on BBH dynamics. Depending on the BBH’s mass, we show that a binary surrounded by an isotropic cold neutral medium (i.e. an asymptotic temperature $T_{\infty } \approx 100$ K) with a time-averaged particle density of $\langle n_\mathrm{ H} \rangle = \mathcal {O}(1)$ cm$^{-3}$ can play a significant role in hardening the binary orbit over a $\mathcal {O}(10^9)$ yr time-scale. Additionally, this causes the black hole’s mass to grow at a rate ${\propto} m^2$. We thus discuss the impact of the ISM on the LIGO-Virgo-KAGRA (LVK) observables and quantify the properties of the ISM under which the latter could act as an additional important pathway for driving a subset of LVK’s BBH mergers.

On the singularity point in acoustic orthorhombic media

Abstract The acoustic orthorhombic model is widely used in seismic modeling and processing of P-wave data. However, the anisotropic acoustic models have so called S-wave artifacts (1 artifact in transversely isotropic acoustic medium and two artifacts in orthorhombic acoustic medium). I show that S-wave artifacts can have one singularity point that results in complications in polarization field and the group velocity surface. The conditions of the existence of this point are defined in terms of anellipticity parameters. This singularity point and its group velocity image are the objects of my analysis.