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 viscous and heat-conductive fluids governed by the compressible NavierStokes equations with an external potential force, there exist non-trivial stationary solutions with zero velocity. By combining the L^p - L^q estimates for the linearized equations and an elaborate energy method, the convergence rates are obtained in various norms for the solution to the stationary profile in the whole space when the initial perturbation of the stationary solution and the potential force are small in some Sobolev norms. More precisely, the optimal convergence rates of the solution and its first order derivatives in L²-norm are obtained when the L¹-norm of the perturbation is bounded.

The propagation and interaction in more than one space dimension of localized pulse solutions (so-called light bullets) to the sine-Gordon [SG] equation is studied both asymptotically and numerically. Similar solutions and their resemblance to solitons in integrable systems were observed numerically before in vector Maxwell systems. The simplicity of SG allows us to perform an asymptotic analysis of counterpropagating pulses, as well as a fully resolved computation over rectangular domains. Numerical experiments are carried out on single pulse propagation and on two pulse collision under different orientations. The particle nature, as known for solitons, persists in these two space dimensional solutions as long as the amplitudes of initial data range in a finite interval, similar to the conditions on the vector Maxwell systems.

This article is concerned with the time periodic solution to the isentropic compressible Navier-Stokes equations in a periodic domain. Using an approach of parabolic regularization, we first obtain the existence of the time periodic solution to a regularized problem under some smallness and symmetry assumptions on the external force. The result for the original compressible Navier-Stokes equations is then obtained by a limiting process. The uniqueness of the periodic solution is also given.