A comparative study of two stochastic mode reduction methods

We present a comparative study of two methods for the reduction of the dimensionality of a system of ordinary differential equations that exhibits time-scale separation. Both methods lead to a reduced system of stochastic differential equations. The novel feature of these methods is that they allow the use, in the reduced system, of higher order terms in the resolved variables. The first method, proposed by Majda, Timofeyev and Vanden-Eijnden, is based on an asymptotic strategy developed by Kurtz. The second method is a short-memory approximation of the Mori-Zwanzig projection formalism of irreversible statistical mechanics, as proposed by Chorin, Hald and Kupferman. We present conditions under which the reduced models arising from the two methods should have similar predictive ability. We apply the two methods to test cases that satisfy these conditions. The form of the reduced models and the numerical simulations show that the two methods have similar predictive ability as expected.Comment: 35 pages, 6 figures. Under review in Physica

Role of scattering in virtual source array imaging

We consider imaging in a scattering medium where the illumination goes through this medium but there is also an auxiliary, passive receiver array that is near the object to be imaged. Instead of imaging with the source-receiver array on the far side of the object we image with the data of the passive array on the near side of the object. The imaging is done with travel time migration using the cross correlations of the passive array data. We showed in [J. Garnier and G. Papanicolaou, Inverse Problems {28} (2012), 075002] that if (i) the source array is infinite, (ii) the scattering medium is modeled by either an isotropic random medium in the paraxial regime or a randomly layered medium, and (iii) the medium between the auxiliary array and the object to be imaged is homogeneous, then imaging with cross correlations completely eliminates the effects of the random medium. It is as if we imaged with an active array, instead of a passive one, near the object. The purpose of this paper is to analyze the resolution of the image when both the source array and the passive receiver array are finite. We show with a detailed analysis that for isotropic random media in the paraxial regime, imaging not only is not affected by the inhomogeneities but the resolution can in fact be enhanced. This is because the random medium can increase the diversity of the illumination. We also show analytically that this will not happen in a randomly layered medium, and there may be some loss of resolution in this case.Comment: 22 pages, 4 figure