Search for: All records

We consider random wave coupling along a flat boundary in dimension three, where the coupling is between surface and body modes and is induced by scattering by a randomly heterogeneous medium. In an appropriate scaling regime we obtain a system of radiative transfer equations which are satisfied by the mean Wigner transform of the mode amplitudes. We provide a rigorous probabilistic framework for describing solutions to this system using that it has the form of a Kolmogorov equation for some Markov process. We then prove statistical stability of the smoothed Wigner transform under the Gaussian approximation. We conclude with analyzing the nonlinear inverse problem for the radiative transfer equations and establish the unique recovery of phase and group velocities as well as power spectral information for the medium fluctuations from the observed smoothed Wigner transform. The mentioned statistical stability is essential in monitoring applications where the realization of the random medium may change. 

The shower curtain effect is commonly described as being able to see a person behind a shower curtain better than that person can see us. This asymmetric phenomenon has been observed in numerical simulations in various propagation models and in optics experiments. Here we present an analysis in the paraxial regime to give a novel characterization of the mechanism behind this effect and we discuss applications to imaging. The paraxial regime is for instance appropriate to model the propagation of a laser beam in a turbulent atmosphere. The theory that we present has also applications to tissue imaging. We consider two different measurement and imaging setups (matched field imaging and optical imaging) to clarify the shower curtain mechanism. We give a quantitative description of how the placement of the shower curtain, modeled as a randomly heterogeneous section, affects the optical imaging resolution. We moreover analyze the signal-to-noise ratio of the image. The analysis involves the study of multifrequency fourth-order moments associated with the Itˆo-Schr¨odinger equation and reveals that broadband sources are necessary to ensure statistical stability and high signal-to-noise ratio. 

A theory for the characterization of the fourth-order moment of electromagnetic wave beams is presented in the case when the source is partially coherent. A Gaussian–Schell model is used for the partially coherent random source. The white-noise paraxial regime is considered, which holds when the wavelength is much smaller than the correlation radius of the source, the beam radius of the source, and the correlation length of the medium, which are themselves much smaller than the propagation distance. The complex wave amplitude field can then be described by the Itô-Schrödinger equation. This equation gives closed evolution equations for the wave field moments at all orders and here the fourth-order moment equations are considered. The general fourth-order moment equations are solved explicitly in the scintillation regime (when the correlation radius of the source is of the same order as the correlation radius of the medium, but the beam radius is much larger) and the result gives a characterization of the intensity covariance function. The form of the intensity covariance function derives from the solution of the transport equation for the Wigner distribution associated with the second-order wave moment. The fourth-order moment results for polarized waves are used in an application for imaging of partially coherent sources. 

When waves propagate through a complex medium like the turbulent atmosphere the wave field becomes incoherent and the wave intensity forms a complex speckle pattern. In this paper we study a speckle memory effect in the frequency domain and some of its consequences. This effect means that certain properties of the speckle pattern produced by wave transmission through a randomly scattering medium is preserved when shifting the frequency of the illumination. The speckle memory effect is characterized via a detailed novel analysis of the fourth-order moment of the random paraxial Green's function at four different frequencies. We arrive at a precise characterization of the frequency memory effect and what governs the strength of the memory. As an application we quantify the statistical stability of time-reversal wave refocusing through a randomly scattering medium in the paraxial or beam regime. Time reversal refers to the situation when a transmitted wave field is recorded on a time-reversal mirror then time reversed and sent back into the complex medium. The re-emitted wave field then refocuses at the original source point. We compute the mean of the refocused wave and identify a novel quantitative description of its variance in terms of the radius of the time-reversal mirror, the size of its elements, the source bandwidth, and the statistics of the random medium fluctuations. 

We study the paraxial wave equation with a randomly perturbed index of refraction, which can model the propagation of a wave beam in a turbulent medium. The random perturbation is a stationary and isotropic process with a general form of the covariance that may or may not be integrable. We focus attention mostly on the nonintegrable case, which corresponds to a random perturbation with long-range correlations, that is, relevant for propagation through a cloudy turbulent atmosphere. The analysis is carried out in a high-frequency regime where the forward scattering approximation holds. It reveals that the randomization of the wave field is multiscale: The travel time of the wave front is randomized at short distances of propagation, and it can be described by a fractional Brownian motion. The wave field observed in the random travel time frame is affected by the random perturbations at long distances, and it is described by a Schr\"odinger-type equation driven by a standard Brownian field. We use these results to quantify how scattering leads to decorrelation of the spatial and spectral components of the wave field and to a deformation of the pulse emitted by the source. These are important questions for applications, such as imaging and free space communications with pulsed laser beams through a turbulent atmosphere. We also compare the results with those used in the optics literature, which are based on the Kolmogorov model of turbulence. We show explicitly that the commonly used approximations for the decorrelation of spatial and spectral components are appropriate for the Kolmogorov model but fail for models with long-range correlations. 

We introduce a novel approach to waveform inversion based on a data-driven reduced order model (ROM) of the wave operator. The presentation is for the acoustic wave equation, but the approach can be extended to elastic or electromagnetic waves. The data are time resolved measurements of the pressure wave gathered by an acquisition system that probes the unknown medium with pulses and measures the generated waves. We propose to solve the inverse problem of velocity estimation by minimizing the square misfit between the ROM computed from the recorded data and the ROM computed from the modeled data, at the current guess of the velocity. We give a step by step computation of the ROM, which depends nonlinearly on the data and yet can be obtained from them in a noniterative fashion, using efficient methods from linear algebra. We also explain how to make the ROM robust to data inaccuracy. The ROM computation requires the full array response matrix gathered with colocated sources and receivers. However, we find that the computation can deal with an approximation of this matrix, obtained from towed-streamer data using interpolation and reciprocity on-the-fly. Although the full-waveform inversion approach of nonlinear least-squares data fitting is challenging without low-frequency information, due to multiple minima of the data fit objective function, we find that the ROM misfit objective function has better behavior, even for a poor initial guess. We also find by explicit computation of the objective functions in a simple setting that the ROM misfit objective function has convexity properties, whereas the least-squares data fit objective function displays multiple local minima. 

We present a theory for wave scintillation in the situation of a time-dependent partially coherent source and a time-dependent randomly heterogeneous medium. Our objective is to understand how the scintillation index of the measured intensity depends on the source and medium parameters. We deduce from an asymptotic analysis of the random wave equation a general form of the scintillation index, and we evaluate this in various scaling regimes. The scintillation index is a fundamental quantity that is used to analyze and optimize imaging and communication schemes. Our results are useful to quantify the scintillation index under realistic propagation scenarios and to address such optimization challenges.