Stability of viscous shock wave for compressible Navier-Stokes equations with free boundary

A free boundary problem for the one-dimensional compressible Navier-Stokes equations is investigated. The asymptotic stability of the viscous shock wave is established under some smallness conditions. The proof is given by an elementary energy estimate.Comment: 20 page

A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting

An explicit moving boundary method for the numerical solution of time-dependent hyperbolic conservation laws on grids produced by the intersection of complex geometries with a regular Cartesian grid is presented. As it employs directional operator splitting, implementation of the scheme is rather straightforward. Extending the method for static walls from Klein et al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme calculates fluxes needed for a conservative update of the near-wall cut-cells as linear combinations of standard fluxes from a one-dimensional extended stencil. Here the standard fluxes are those obtained without regard to the small sub-cell problem, and the linear combination weights involve detailed information regarding the cut-cell geometry. This linear combination of standard fluxes stabilizes the updates such that the time-step yielding marginal stability for arbitrarily small cut-cells is of the same order as that for regular cells. Moreover, it renders the approach compatible with a wide range of existing numerical flux-approximation methods. The scheme is extended here to time dependent rigid boundaries by reformulating the linear combination weights of the stabilizing flux stencil to account for the time dependence of cut-cell volume and interface area fractions. The two-dimensional tests discussed include advection in a channel oriented at an oblique angle to the Cartesian computational mesh, cylinders with circular and triangular cross-section passing through a stationary shock wave, a piston moving through an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil profile.Comment: 30 pages, 27 figures, 3 table

Cavitation Induction by Projectile Impacting on a Water Jet

The present paper focuses on the simulation of the high-velocity impact of a projectile impacting on a water-jet, causing the onset, development and collapse of cavitation. The simulation of the fluid motion is carried out using an explicit, compressible, density-based solver developed by the authors using the OpenFOAM library. It employs a barotropic two-phase flow model that simulates the phase-change due to cavitation and considers the co-existence of non-condensable and immiscible air. The projectile is considered to be rigid while its motion through the computational domain is modelled through a direct-forcing Immersed Boundary Method. Model validation is performed against the experiments of Field et al. [Field, J., Camus, J. J., Tinguely, M., Obreschkow, D., Farhat, M., 2012. Cavitation in impacted drops and jets and the effect on erosion damage thresholds. Wear 290–291, 154–160. doi:10.1016/j.wear.2012.03.006. URL http://www.sciencedirect.com/science/article/pii/S0043164812000968 ], who visualised cavity formation and shock propagation in liquid impacts at high velocities. Simulations unveil the shock structures and capture the high-speed jetting forming at the impact location, in addition to the subsequent cavitation induction and vapour formation due to refraction waves. Moreover, model predictions provide quantitative information and a better insight on the flow physics that has not been identified from the reported experimental data, such as shock-wave propagation, vapour formation quantity and induced pressures. Furthermore, evidence of the Richtmyer-Meshkov instability developing on the liquid-air interface are predicted when sufficient dense grid resolution is utilised

Absolute and convective instabilities of an inviscid compressible mixing layer: Theory and applications

This study aims to examine the effect of compressibility on unbounded and parallel shear flow linear instabilities. This analysis is of interest for industrial, geophysical, and astrophysical flows. We focus on the stability of a wavepacket as opposed to previous single-mode stability studies. We consider the notions of absolute and convective instabilities first used to describe plasma instabilities. The compressible-flow modal theory predicts instability whatever the Mach number. Spatial and temporal growth rates and Reynolds stresses nevertheless become strongly reduced at high Mach numbers. The evolution of disturbances is driven by Kelvin -Helmholtz instability that weakens in supersonic flows. We wish to examine the occurrence of absolute instability, necessary for the appearance of turbulent motions in an inviscid and compressible two-dimensional mixing layer at an arbitrary Mach number subject to a three-dimensional disturbance. The mixing layer is defined by a parametric family of mean-velocity and temperature profiles. The eigenvalue problem is solved with the help of a spectral method. We ascertain the effects of the distribution of temperature and velocity in the mixing layer on the transition between convective and absolute instabilities. It appears that, in most cases, absolute instability is always possible at high Mach numbers provided that the ratio of slow-stream temperature over fast-stream temperature may be less than a critical maximal value but the temporal growth rate present in the absolutely unstable zone remains small and tends to zero at high Mach numbers. The transition toward a supersonic turbulent regime is therefore unlikely to be possible in the linear theory. Absolute instability can be also present among low-Mach-number coflowing mixing layers provided that this same temperature ratio may be small, but nevertheless, higher than a critical minimal value. Temperature distribution within the mixing layer also has an effect on the growth rate, this diminishes when the slow stream is heated. These results are applied to the dynamics of mixing layers in the interstellar medium and to the dynamics of the heliopause, frontier between the interstellar medium, and the solar wind. (C) 2009 American Institute of Physics

A Compact Third-order Gas-kinetic Scheme for Compressible Euler and Navier-Stokes Equations

In this paper, a compact third-order gas-kinetic scheme is proposed for the compressible Euler and Navier-Stokes equations. The main reason for the feasibility to develop such a high-order scheme with compact stencil, which involves only neighboring cells, is due to the use of a high-order gas evolution model. Besides the evaluation of the time-dependent flux function across a cell interface, the high-order gas evolution model also provides an accurate time-dependent solution of the flow variables at a cell interface. Therefore, the current scheme not only updates the cell averaged conservative flow variables inside each control volume, but also tracks the flow variables at the cell interface at the next time level. As a result, with both cell averaged and cell interface values the high-order reconstruction in the current scheme can be done compactly. Different from using a weak formulation for high-order accuracy in the Discontinuous Galerkin (DG) method, the current scheme is based on the strong solution, where the flow evolution starting from a piecewise discontinuous high-order initial data is precisely followed. The cell interface time-dependent flow variables can be used for the initial data reconstruction at the beginning of next time step. Even with compact stencil, the current scheme has third-order accuracy in the smooth flow regions, and has favorable shock capturing property in the discontinuous regions. Many test cases are used to validate the current scheme. In comparison with many other high-order schemes, the current method avoids the use of Gaussian points for the flux evaluation along the cell interface and the multi-stage Runge-Kutta time stepping technique.Comment: 27 pages, 38 figure