We study the convergence rate of Glimm scheme for general systems of hyperbolic conservation laws without the assumption that each characteristic field is either genuinely nonlinear or linearly degenerate. We first give a sharper estimate of the error arising from the wave tracing argument by a careful analysis of the interaction between small waves. With this key estimate, the convergence rate is shown to be o (1) s 1 3| ln s| 1+ , which is sharper compared to o (1) s 1 4| ln s| given in [T. Yang, Convergence rate of Glimm scheme for general systems of hyperbolic conservation laws, Taiwanese J. Math. 7 (2)(2003) 195205]. However, it is still slower than o (1) s 1 2| ln s| given in [A. Bressan, A. Marson, Error bounds for a deterministic version of the Glimm scheme, Arch. Ration. Mech. Anal. 142 (2)(1998) 155176] for systems with each characteristic field being genuinely nonlinear or linearly degenerate. Here s is the

A semi-explicit formula of solution to the boundary layer system for thermal layer derived from the compressible Navier-Stokes equations with the non-slip boundary condition when the viscosity coefficients vanish is given, in particular in three space dimension. In contrast to the inviscid Prandtl system studied by [7] in two space dimension, the main difficulty comes from the coupling of the velocity field and the temperature field through a degenerate parabolic equation. The convergence of these boundary layer equations to the inviscid Prandtl system is justified when the initial temperature goes to a constant. Moreover, the time asymptotic stability of the linearized system around a shear flow is given, and in particular, it shows that in three space dimension, the asymptotic stability depends on whether the direction of tangential velocity field of the shear flow is invariant in the normal direction respective to the boundary.

For the Boltzmann equation with an external potential force depending only on the space variables, there is a family of stationary solutions, which are local Maxwellians with space dependent density, zero velocity and constant temperature. In this paper, we will study the nonlinear stability of these stationary solutions by using the energy method. The analysis combines the analytic techniques used for the conservation laws using the fluid-type system derived from the Boltzmann equation (cf. [14]) and the dissipative effects on the fluid and non-fluid components of the Boltzmann equation through the celebrated H-theorem. To our knowledge, this is the first result on the global classical solutions to the Boltzmann equation with external force and non-trivial large time behavior in the whole space.

The nonlinear stability in the {L^p} -norm, p ≥ 1 , of stationary weak discrete shocks for the Lax-Friedrichs scheme approximating general m × m systems of nonlinear hyperbolic conservation laws is proved, provided that the summations of the initial perturbations equal zero. The result is proved by using both a weighted estimate and characteristic energy method based on the internal structures of the discrete shocks and the essential monotonicity of the Lax-Friedrichs scheme.

We perform a computationl study of front speeds of G-equation models in time dependent cellular flows. The G-equations arise in premixed turbulent combustion, and are Hamilton-Jacobi type level set partial differential equations (PDEs). The curvature-strain G-equations are also non-convex with degenerate diffusion. The computation is based on monotone finite difference discretization and weighted essentially nonoscillatory (WENO) methods. We found that the large time front speeds lock into the frequency of time periodic cellular flows in curvature-strain G-equations similar to what occurs in the basic inviscid G-equation. However, such frequency locking phenomenon disappears in viscous G-equation, and in the inviscid G-equation if time periodic oscillation of the cellular flow is replaced by time stochastic oscillation.