The aim of this paper is to develop a mathematical framework for opto-elastography. In opto-elastography, a mechanical perturbation of the medium produces a decorrelation of optical speckle patterns due to the displacements of optical scatterers. To model this, we consider two optically random media, with the second medium obtained by shifting the first medium in some local region. We derive the radiative transfer equation for the cross-correlation of the wave fields in the media. Then we derive its diffusion approximation. In both the radiative transfer and the diffusion regimes, we relate the correlation of speckle patterns to the solutions of the radiative transfer and the diffusion equations. We present numerical simulations based on our model which are in agreement with recent experimental measurements.

No abstract uploaded!

In this paper, a homotopy continuation method for the computation of nonnegative Z-/H-eigenpairs of a nonnegative tensor is presented. We show that the homotopy continuation method is guaranteed to compute a nonnegative eigenpair. Additionally, using degree analysis, we show that the number of positive Z-eigenpairs of an irreducible nonnegative tensor is odd. A novel homotopy continuation method is proposed to compute an odd number of positive Z-eigenpairs, and some numerical results are presented.

Accurate simulation of seismic waves is of critical importance in a variety of geophysical applications. Based on recent works on staggered discontinuous Galerkin methods, we have developed a new method for the simulations of seismic waves, which has energy conservation and extremely low grid dispersion, so that it naturally provided accurate numerical simulations of wave propagation useful for geophysical applications and was a generalization of classical staggered-grid finite-difference methods. Moreover, it could handle with ease irregular surface topography and discontinuities in the subsurface models. Our new method discretized the velocity and the stress tensor on this staggered grid, with continuity imposed on different parts of the mesh. The symmetry of the stress tensor was enforced by the Lagrange multiplier technique. The resulting method was an explicit scheme, requiring the solutions of a

No abstract uploaded!