No Arabic abstract
In this paper, phase correction and amplitude compensation are introduced to a previously developed mixed domain method (MDM), which is only accurate for modeling wave propagation in weakly heterogeneous media. Multiple reflections are also incorporated with the one-way model to improve the accuracy. The resulting model is denoted as the modified mixed-domain method (MMDM) and is numerically evaluated for its accuracy and efficiency using two distinct cases: a layered medium and a human skull. It is found that the MMDM is significantly more accurate than the MDM for strongly heterogeneous media, especially when the phase aberrating layer is roughly perpendicular to the acoustic beam. Additionally, convergence study suggests that the second-order reflection is sufficient for wave modeling in lossy biological media. The method developed in this work could be used to facilitate therapeutic ultrasound for treating brain-related diseases and disorders.
Biots theory predicts the wave velocities of a saturated poroelastic granular medium from the elastic properties, density and geometry of its dry solid matrix and the pore fluid, neglecting the interaction between constituent particles and local flow. However, when the frequencies become high and the wavelengths comparable with particle size, the details of the microstructure start to play an important role. Here, a novel hydro-micromechanical numerical model is proposed by coupling the lattice Boltzmann method (LBM) with the discrete element method (DEM. The model allows to investigate the details of the particle-fluid interaction during propagation of elastic waves While the DEM is tracking the translational and rotational motion of each solid particle, the LBM can resolve the pore-scale hydrodynamics. Solid and fluid phases are two-way coupled through momentum exchange. The coupling scheme is benchmarked with the terminal velocity of a single sphere settling in a fluid. To mimic a pressure wave entering a saturated granular medium, an oscillating pressure boundary condition on the fluid is implemented and benchmarked with one-dimensional wave equations. Using a face centered cubic structure, the effects of input waveforms and frequencies on the dispersion relations are investigated. Finally, the wave velocities at various effective confining pressures predicted by the numerical model are compared with with Biots analytical solution, and a very good agreement is found. In addition to the pressure and shear waves, slow compressional waves are observed in the simulations, as predicted by Biots theory.
We are interested in the modeling of wave propagation in an infinite bilayered acoustic/poroelastic media. We consider the biphasic Biots model in the poroelastic layer. The first part is devoted to the calculation of analytical solution in two dimensions, thanks to Cagniard de Hoop method. In this second part we consider the 3D case.
In this work, we propose a local multiscale model reduction approach for the time-domain scalar wave equation in a heterogenous media. A fine mesh is used to capture the heterogeneities of the coefficient field, and the equation is solved globally on a coarse mesh in the discontinuous Galerkin discretization setting. The main idea of the model reduction approach is to extract dominant modes in local spectral problems for representation of important features, construct multiscale basis functions in coarse oversampled regions by constraint energy minimization problems, and perform a Petrov-Galerkin projection and a symmetrization onto the coarse grid. The method is expicit and energy conserving, and exhibits both coarse-mesh and spectral convergence, provided that the oversampling size is appropriately chosen. We study the stability and convergence of our method. We also present numerical results on the Marmousi model in order to test the performance of the method and verify the theoretical results.
It has been recently shown that quadruply lensed gravitational-wave (GW) events due to coalescing binaries can be localized to one or just a few galaxies, even in the absence of an electromagnetic counterpart. We discuss how this can be used to extract information on modified GW propagation, which is a crucial signature of modifications of gravity at cosmological scales. We show that, using quadruply lensed systems, it is possible to constrain the parameter $Xi_0$ that characterizes modified GW propagation, without the need of imposing a prior on $H_0$. A LIGO/Virgo/Kagra network at target sensitivity might already get a significant measurement of $Xi_0$, while a third generation GW detector such as the Einstein Telescope could reach a very interesting accuracy.
Context. The frequencies, lifetimes, and eigenfunctions of solar acoustic waves are affected by turbulent convection, which is random in space and in time. Since the correlation time of solar granulation and the periods of acoustic waves ($sim$5 min) are similar, the medium in which the waves propagate cannot a priori be assumed to be time independent. Aims. We compare various effective-medium solutions with numerical solutions in order to identify the approximations that can be used in helioseismology. For the sake of simplicity, the medium is one dimensional. Methods. We consider the Keller approximation, the second-order Born approximation, and spatial homogenization to obtain theoretical values for the effective wave speed and attenuation (averaged over the realizations of the medium). Numerically, we computed the first and second statistical moments of the wave field over many thousands of realizations of the medium (finite-amplitude sound-speed perturbations are limited to a 30 Mm band and have a zero mean). Results. The effective wave speed is reduced for both the theories and the simulations. The attenuation of the coherent wave field and the wave speed are best described by the Keller theory. The numerical simulations reveal the presence of coda waves, trailing the coherent wave packet. These late arrival waves are due to multiple scattering and are easily seen in the second moment of the wave field. Conclusions. We find that the effective wave speed can be calculated, numerically and theoretically, using a single snapshot of the random medium (frozen medium); however, the attenuation is underestimated in the frozen medium compared to the time-dependent medium. Multiple scattering cannot be ignored when modeling acoustic wave propagation through solar granulation.