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.

Let U(n) be the unitary group, and u(n)¤ the dual of its Lie algebra, equipped with the Kirillov Poisson structure. In their 1983 paper, Guillemin-Sternberg introduced a densely defined Hamiltonian action of a torus of dimension (n−1)n/2 on u(n)¤, with moment map given by the Gelfand-Zeitlin coordinates. A few years later, Flaschka-Ratiu described a similar, ‘multiplicative’ GelfandZeitlin system for the Poisson Lie group U(n)¤. By the Ginzburg-Weinstein theorem, U(n)¤ is isomorphic to u(n)¤ as a Poisson manifold. Flaschka-Ratiu conjectured that one can choose the Ginzburg-Weinstein diffeomorphism in such a way that it intertwines the linear and nonlinear Gelfand-Zeitlin systems. Our main result gives a proof of this conjecture, and produces a canonical Ginzburg-Weinstein diffeomorphism.

In this paper, we develop the high order weighted essentially non-oscillatory (WENO) schemes for solving the Degasperis-Procesi (DP) equation, including finite volume (FV) and finite difference (FD) methods. The DP equation contains nonlinear high order derivatives, and possibly discontinuous or sharp transition solutions. The finite volume method is de- signed based on the total variation bounded property of the DP equation. And the finite difference method is constructed based on the L2 stability of the DP equation. Due to the adoption of the WENO reconstruction, both schemes are arbitrary high order accuracy and shock capturing. The numerical simulation results for different types of solutions of the nonlinear Degasperis-Procesi equation are provided to illustrate the accuracy and capability of the methods.

No abstract uploaded!

No abstract uploaded!