In the long-time scale, we consider the fluid dynamical limits for the kinetic equations when the fluctuation is decomposed into even and odd parts with respect to the microscopic velocity with different scalings. It is shown that when the background state is an absolute Maxwellian, the limit fluid dynamical equations are the incompressible Navier-Stokes equations with viscous heating. This is different from the case when the even and odd parts of the fluctuation have the same scaling where the standard incompressible Navier-Stokes equations without viscous heating are obtained. On the other hand, when the background is a local Maxwellian, it is shown that the above even-odd decomposition leads to a non-classical fluid dynamical system without viscous heating which has been used to describe the ghost effect in the kinetic theory. In addition, the above even-odd decomposition is justified rigorously for the Boltzmann equation for the former case when the background is an absolute Maxwellian.

In this survey paper, we will present the recent work on the study of the compressible fluids with vacuum states by illustrating its interesting and singular behavior through some systems of fluid dynamics, that is, Euler equations, EulerPoisson equations and NavierStokes equations. The main concern is the well-posedness of the problem when vacuum presents and the singular behavior of the solution near the interface separating the vacuum and the gas. Furthermore, the relation of the solutions for the gas dynamics with vacuum to those of the Boltzmann equation will also be discussed. In fact, the results obtained so far for vacuum states are far from being complete and satisfactory. Therefore, this paper can only be served as an introduction to this interesting field which has many open and challenging mathematical problems. Moreover, the problems considered here are limited to the author's interest and

METHODS AND APPLICATIONS OF ANALYSIS. 2000 International Press Vol. 7, No. 3, pp. 495-510, September 2000 007

We study homogenization of the G-equation with a flow straining term (or the strain G-equation) in two dimensional periodic cellular flow. The strain G-equation is a highly non-coercive and non-convex level set HamiltonJacobi equation. The main objective is to investigate how the flow induced straining (the nonconvex term) influences front propagation as the flow intensity <i>A</i> increases. Three distinct regimes are identified. When <i>A</i> is below the critical level, homogenization holds and the turbulent flame speed <i>s</i> _T (effective Hamiltonian) is well-defined for any periodic flow with small divergence and is enhanced by the cellular flow as <i>s</i> _T <i>O</i>(<i>A</i>/log <i>A</i>). In the second regime where <i>A</i> is slightly above the critical value, homogenization breaks down, and <i>s</i> _T is not well-defined along any direction. Solutions become a mixture of a fast moving part and a stagnant part. When <i>A</i> is

We study wetting front (traveling wave) solutions to the Richards equation that describe the vertical infiltration of water through one-dimensional periodically layered unsaturated soils. We prove the existence, uniqueness, and large time asymptotic stability of the traveling wave solutions under prescribed flux boundary conditions and certain constitutive conditions. The traveling waves are connections between two steady state solutions that form near the ground surface and towards the underground water table. We found a closed form expression of the wave speed. The speed of a traveling wave is equal to the ratio of the flux difference and the difference of the spatial averages of the two steady states. We give both analytical and numerical examples showing that the wave speeds in the periodic soils can be larger or smaller than those in the homogeneous soils which have the same mean diffusivity and conductivity

The minimal speeds (c) of the KolmogorovPetrovskyPiskunov (KPP) fronts at small diffusion ( 1) in a class of time-periodic cellular flows with chaotic streamlines is investigated in this paper. The variational principle of c reduces the computation to that of a principle eigenvalue problem on a periodic domain of a linear advectiondiffusion operator with spacetime periodic coefficients and small diffusion. To solve the advection dominated time-dependent eigenvalue problem efficiently over large time, a combination of spectral methods and finite element, as well as the associated fast solvers, are utilized to accelerate computation. In contrast to the scaling c= O ( 1/4) in steady cellular flows, a new relation c= O (1) as 1 is revealed in the time-periodic cellular flows due to the presence of chaotic streamlines. Residual propagation speed emerges from the Lagrangian chaos which is quantified as a sub

The shape stability of the reaction interface for reactive flow in a layered porous medium is studied. This is done using a complete linearized stability analysis in the setting of a free boundary model of this phenomenon. The spectrum of the linearized problem is obtained in the general case, and it is compared with that obtained for homogeneous media in two typical cases.

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