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 study the front dynamics of solutions of the initial value problem of the Burgers equation with initial data being the viscous shock front plus the white noise perturbation. In the sense of distribution, the solutions propagate with the same speed as the unperturbed front, however, the front location is random and satisfies a central limit theorem with the variance proportional to the time<i>t</i>, as<i>t</i> goes to infinity. With probability arbitrarily close to one, the front width is<i>O</i>(1) for large time.

Since the work by Guo (Invent Math 153(3):593630, 2003), it has remained an open problem to establish the global existence of perturbative classical solutions around a global Maxwellian to the VlasovMaxwellBoltzmann system with the whole range of soft potentials. This is mainly due to the complex structure of the system, in particular, the degenerate dissipation at large velocity, the velocity-growth of the nonlinear term induced by the Lorentz force, and the regularity-loss of the electromagnetic fields. This paper solves this problem in the whole space provided that initial perturbation has sufficient regularity and velocity-integrability.

We revisit the homogenization problem for the Poisson equation in periodically perforated domains with zero Neumann data at the boundary of the holes and prescribed Dirichlet data at the outer boundary. It is known that, if the periodicity of the holes goes to zero but their volume fraction remains fixed and positive, the limit problem is a Dirichlet boundary value problem posed in the domain without the holes, and the effective diffusion coefficients are non-trivially modified; if that volume fraction goes to zero instead, i.e. the holes are dilute, the effective operator remains the Laplacian (that is, unmodified). Our main results contain the study of a "continuity" in those effective models with respect to the volume fraction of the holes and some new convergence rates for homogenization in the dilute setting. Our method explores the classical two-scale expansion ansatz and relies on asymptotic analysis of the rescaled cell problems using layer potential theory.

We present a domain decomposition method on a semilinear hyperbolic system with multiple relaxation times. In the region where the relaxation time is small, an asymptotic equilibrium equation can be used for computational efficiency. An interface condition based on the sign of the characteristic speed at the interface is provided to couple the two systems in a domain decomposition setting. A rigorous analysis, based on the Laplace Transform, on the $ L^2$ error estimate is presented for the linear case, which shows how the error of the domain decomposition method depends on the smaller relaxation time, and the boundary and interface layer effects. The given convergence rate is optimal. We present a numerical implementation of this domain decomposition method, and give some numerical results in order to study the performance of this method.