As a typical example of hyperbolic conservation laws, the system of Euler equations representing the conservation of mass, momentum and energy has three basic wave patterns. They are two nonlinear waves called shock and rarefaction wave, and one linearly degenerate wave called contact discontinuity. These basic wave patterns have their counterparts both in other physical models for gas motions in equilibrium involving viscosity and heat conductivity, and gas motion in non-equilibrium. The stability of these basic wave patterns in the system of Navier-Stokes equations and the Boltzmann equation has been an active research topic. Even though the stability of the two nonlinear wave patterns has been extensively studied, the stability of the linearly degenerate contact wave was not solved until recently. In this paper, we will briefly present our recent results in [32] on the stability of the contact wave patterns for both the Navier-Stokes equations and the Boltzmann equation.

We investigate the singularity formation of a 3D model that was recently proposed by Hou and Lei (2009) in [15] for axisymmetric 3D incompressible Navier–Stokes equations with swirl. The main difference between the 3D model of Hou and Lei and the reformulated 3D Navier–Stokes equations is that the convection term is neglected in the 3D model. This model shares many properties of the 3D incompressible Navier– Stokes equations. One of the main results of this paper is that we prove rigorously the finite time singularity formation of the 3D inviscid model for a class of initial boundary value problems with smooth initial data of finite energy. We also prove the global regularity of the 3D inviscid model for a class of small smooth initial data.

A new augmented method is proposed for elliptic interface problems with a piecewise variable coefficient that has a finite jump across a smooth interface. The main motivation is to get not only a second order accurate solution but also a second order accurate gradient from each side of the interface. Key to the new method is introducing the jump in the normal derivative of the solution as an augmented variable and rewriting the interface problem as a new PDE that consists of a leading Laplacian operator plus lower order derivative terms near the interface. In this way, the leading second order derivative jump relations are independent of the jump in the coefficient that appears only in the lower order terms after the scaling. An upwind type discretization is used for the finite difference discretization at the irregular grid points on or near the interface so that the resulting coefficient matrix is an M-matrix. A multigrid solver is used to solve the linear system of equations, and the GMRES iterative method is used to solve the augmented variable. Second order convergence for the solution and the gradient from each side of the interface is proved in this paper. Numerical examples for general elliptic interface problems confirm the theoretical analysis and efficiency of the new method.

In this paper, we discuss the Weyl problem in warped product spaces. We apply the method of continuity and prove the openness of the Weyl problem. A counterexample is constructed to show that the isometric embedding of the sphere with canonical metric is not unique up to an isometry if the ambient warped product space is not a space form. Then, we study the rigidity of the standard sphere if we fix its geometric center in the ambient space. Finally, we discuss a Shi-Tam type of inequality for the Schwarzschild manifold as an application of our findings.

Fluid dynamic limit to compressible Euler equations from compressible Navier-Stokes equations and Boltzmann equation has been an active topic with limited success so far. In this paper, we consider the case when the solution of the Euler equations is a Riemann solution consisting two rarefaction waves and a contact discontinuity and prove this limit for both Navier-Stokes equations and the Boltzmann equation when the viscosity, heat conductivity coefficients and the Knudsen number tend to zero respectively. In addition, the uniform convergence rates in terms of the above physical parameters are also obtained. It is noted that this is the first rigorous proof of this limit for a Riemann solution with superposition of three waves even though the fluid dynamic limit for a single wave has been proved.