No abstract uploaded!

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.

No abstract uploaded!

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.