EXISTENCE AND STABILITY OF VISCOUS SHOCK PROFILES FOR 2-D ISENTROPIC MHD WITH INFINITE ELECTRICAL RESISTIVITY

For the two-dimensional Navier-Stokes equations of isentropic magnetohydrodynamics （MHD） with γ-law gas equation of state, γ≥ 1, and infinite electrical resistivity, we carry out a global analysis categorizing all possible viscous shock profiles. Precisely, we show that the phase portrait of the traveling-wave ODE generically consists of either two rest points connected by a viscous Lax profile, or else four rest points, two saddles and two nodes. In the latter configuration, which rest points are connected by profiles depends on the ratio of viscosities, and can involve Lax, overcompressive, or undercompressive shock profiles. Considered as three-dimensional solutions, undercompressive shocks axe Lax-type （Alfven） waves. For the monatomic and diatomic cases γ= 5/3 and γ=7/5, with standard viscosity ratio for a nonmagnetic gas, we find numerically that the the nodes are connected by a family of overcompressive profiles bounded by Lax profiles connecting saddles to nodes, with no undercompressive shocks occurring. We carry out a systematic numerical Evans function analysis indicating that all of these two-dimensional shock pro- files are linearly and nonlinearly stable, both with respect to two- and three-dimensional perturbations. For the same gas constants, but different viscosity ratios, we investigate also cases for which undercompressive shocks appear; these are seen numerically to be stable as well, both with respect to two-dimensional and （in the neutral sense of convergence to nearby Riemann solutions） three-dimensional perturbations.

Analysis of shock wave reflection from fixed and moving boundaries using a stabilized particle method

In the present paper, the efficiency of an enhanced formulation of the stabilized corrective smoothed particle method （CSPM） for simulation of shock wave propagation and reflection from fixed and moving solid boundaries in compressible fluids is investigated. The Lagrangian nature and its accuracy for imposing the boundary conditions are the two main reasons for adoption of CSPM. The governing equations are further modified for imposition of moving solid boundary conditions. In addition to the traditional artificial viscosity, which can remove numerically induced abnormal jumps in the field values, a velocity field smoothing technique is introduced as an efficient method for stabilizing the solution. The method has been implemented for one- and two-dimensional shock wave propagation and reflection from fixed and moving boundaries and the results have been compared with other available solutions. The method has also been adopted for simulation of shock wave propagation and reflection from infinite and finite solid boundaries.