We study the two-point correlation function of a freely decaying scalar in Kraichnan's model of advection by a Gaussian random velocity field that is stationary and white noise in time, but fractional Brownian in space with roughness exponent 0&lt;<i></i>&lt;2, appropriate to the inertial-convective range of the scalar. We find all self-similar solutions by transforming the scaling equation to Kummer's equation. It is shown that only those scaling solutions with scalar energy decay exponent <i>a</i>(<i>d</i>/<i></i>)+1 are statistically realizable, where <i>d</i> is space dimension and <i></i>=2<i></i>. An infinite sequence of invariants <i>J</i> _p, <i>p</i>=0,1,2,..., is pointed out, where <i>J</i> ₀ is Corrsin's integral invariant but the higher invariants appear to be new. We show that at least one of the invariants <i>J</i> ₀ or <i>J</i> ₁ must be nonzero (possibly infinite) for realizable initial data. Initial datum with a finite, nonzero invariantthe first being <i>J</i> _pconverges at long times to a scaling

We develop an efficient numerical method to compute single slit or double slit diffraction patterns from high frequency wave in inhomogeneous media. We approximate the high frequency asymptotic solution to the Helmholtz equation using the Eulerian Gaussian beam summation proposed in [20, 21]. The emitted rays from a slit are embedded in the phase space using an open segment. The evolution of this open curve is accurately computed using the recently developed Grid Based Particle Method [24] which results in a very efficient computational algorithm. Following the grid based particle method we proposed in [23, 24], we represent the open curve or the open surface by meshless Lagrangian particles sampled according to an underlying fixed Eulerian mesh. The end-points of the open curve are tracked explicitly and consistently with interior particles. To construct the overall wavefield, each of these sampling particles also carry necessary quantities that are obtained by solving advection-reaction equations. Numerical experiments show that the resulting method can model diffraction patterns in inhomogeneous media accurately, even in the occurrence of caustics.

In this paper, our aim is to present (1) an embedded fracture model (EFM) for coupled flow and mechanics problem based on the dual continuum approach on the fine grid and (2) an upscaled model for the resulting fine grid equations. The mathematical model is described by the coupled system of equation for displacement, fracture and matrix pressures. For a fine grid approximation, we use the finite volume method for flow problem and finite element method for mechanics. Due to the complexity of fractures, solutions have a variety of scales, and fine grid approximation results in a large discrete system. Our second focus in the construction of the upscaled coarse grid poroelasticity model for fractured media. Our upscaled approach is based on the nonlocal multicontinuum (NLMC) upscaling for coupled flow and mechanics problem, which involves computations of local basis functions via an energy minimization