求解可压缩流的二维通量分裂格式

传统的一维通量分裂格式在计算界面数值通量时,只考虑网格界面法向的波系。采用传统的TV格式分别求解对流通量和压力通量。通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的TV通量分裂格式。在计算一维数值算例时,该格式与传统的TV格式具有相同的数值通量计算公式,因此其保留了传统的TV格式精确捕捉接触间断和膨胀激波的优点。在计算二维算例时,该格式比传统的TV格式具有更高的分辨率;在计算二维强激波问题时,消除了传统TV格式的非物理现象,表现出更好的鲁棒性;此外,该格式大大提高了稳定性CFL数,从而具有更高的计算效率。因此,本文方法是一种精确、高效并且具有强鲁棒性的数值方法,在可压缩流的数值模拟中具有广阔的应用前景。     

Identification of concentrated structures in slightly compressible two-dimensional flows

In this paper, different vortex diagnostic methods are compared to obtain a better understanding of boundary layer influence on the transport of vortical structures involving a complete analysis of vorticity, the Vorticity Threshold Criterion (VTC), and the Weiss Criterion (WC). These three techniques are basically confronted to find a suitable understanding of all flow characteristics for a range of laminar to transitional Reynolds numbers. The computations on this dihedral plane are done using a 2D DNS method. The Weiss criterion, coming from the analysis of the incompressible Euler equations is validated and applied to low speed compressible flows (Mach number=0.2).     

A shock tube research facility for high-resolution measurements of compressible turbulence

A new shock tube facility with a 30.5 cm (1 ft) inside diameter is currently in operation that allows for high-spatial-resolution measurements of compressible turbulence. Small scales of turbulence behave very differently from large scales when they interact with shock or expansion waves. Highly resolved measurements can provide new information on the interaction at small scales. Another notable characteristic of the present facility is the ability to control the flow velocity behind the reflected shock through the porosity of the reflecting wall. Tests showed good flow quality with sufficiently long observation times. Measurements of piecewise average skin friction over short segments of the tube indicated strong viscous effects very close to the diaphragm where the shock is developing. The skin friction and the shock propagation speed virtually remained constant inside the working section of the shock tube in all investigated flow cases, even in low Mach number cases where viscous effects are stronger. The experimental results are compared with numerical simulations, including the effects of the reflecting porous wall and viscous effects.     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

Unsteady three-dimensional compressible vortex flows generated during shock wave diffraction

The results of an experimental and numerical investigation into the behaviour of the spiral vortex generated by shock wave diffraction over edges yawed to the incident shock wave are presented. Three-dimensional numerical simulations reveal significant distortion and bending of the free vortex in regions near the boundary of the flow domain, so as to meet it at a right angle. The results of numerical simulations were found to mimic the experimentally obtained photographs very well. The numerical results are used to explain the various features of the resultant flow fields, with particular emphasis placed on the behaviour and properties of the spiral vortex, as it evolves with time. The effects of bending on the structure of the vortex are examined. The rate of circulation production for the three-dimensional shock diffraction cases was calculated, and the trends observed correlated with those for the much published two-dimensional diffraction case.     

Class of two-dimensional nonstationary flows with shock waves

A numerical investigation was made of a two-dimensional plane-parallel nonsteady motion resulting from the decay of an arbitrary discontinuity on the boundaries of several neighboring angular regions filled by gases in different states.     

Planar laser-induced fluorescence imaging in shock tube flows

The planar laser-induced fluorescence (PLIF) imaging method was used to perform flow visualization and quantitative planar thermometry in shock tube flow fields using toluene as a fluorescence tracer in nitrogen. Fluorescence quantum yield values needed to quantify PLIF images were measured in a static cell at low pressures (<1 bar) for various toluene partial pressures in nitrogen bath gas. Images behind incident and reflected shocks were taken in the core flow away from regions affected by boundary layers. Temperature measurements from these images were successfully compared with predicted values using ideal shock equations. Measured temperatures ranged between 296 and 800 K and pressures between 0.15 and 1.5 atm. The average temperature discrepancies between measurements and the predicted values behind the incident and reflected shocks were 1.6 and 3.6%, respectively. Statistical analyses were also conducted to calculate the temperature measurement uncertainty as a function of image resolution. The technique was also applied to the study of more complex supersonic flows, specifically the interaction of a moving shock with a wedge. Measured temperatures agreed well with the results of numerical simulations in all inviscid regions, and all pertinent features of the single Mach reflection were resolved.     

A study on characteristics of shock train inside a shock tube

Shock tubes are devices which are used in the investigation of high speed and high temperature flow of compressible gas. Inside a shock tube, the interaction between the reflected shock wave and boundary layer leads to a complex flow phenomenon. Initially a normal shock wave is formed in the shock tube which migrates toward the closed end of the tube and that in turn leads to the reflection of shock. Due to the boundary layer interaction with the reflected shock, the bifurcation of shock wave takes place. The bifurcated shock wave then approaches the contact surface and shock train is generated. Till date only a few studies have been conducted to investigate this shock train phenomenon inside the shock tube. For the present study a computational fluid dynamics(CFD) analysis has been performed on a two dimensional axi-symmetric model of a shock tube using unsteady, compressible Navier–Stokes equations. In order to investigate the detailed characteristics of shock train, parametric studies have been performed by varying different parameters such as the shock tube length, diameter, pressure ratio used inside the shock tube.     

Analogy between soap film and gas dynamics. II. Experiments on one-dimensional motion of shock waves in soap films

This paper presents an experimental investigation of one-dimensional moving shock waves in vertical soap films. The shock waves were generated by bursting the films with a perforating spark. Images of propagating shock waves and small disturbances were recorded using a fast line scan CCD camera. An aureole and a shock wave preceding the rim of the expanding hole were clearly observed. These images are similar to the x-t diagrams in gas dynamics and give the velocities of shock and sound waves. The moving shock waves cause jumps in thickness. The variations of the induced Mach number, M₂ and the ratio of film thickness across the shock wave, δ ₂/δ ₁, are plotted versus the shock Mach number, M _s. Both results suggest that soap films are analogous to compressible gases with a specific heat ratio of γ≅1.0. Published online: 15 October 2002     

Analysis of the secondary stability of compressible boundary layer flows over a two-dimensional plate

ＩｎｔｒｏｄｕｃｔｉｏｎＲｅｃｅｎｔｌｙ,ｉｎｔｈｅａｎａｌｙｓｉｓｏｆｂｏｕｎｄａｒｙｌａｙｅｒｔｒａｎｓｉｔｉｏｎｆｒｏｍｌａｍｉｎａｒｔｏｔｕｒｂｕｌｅｎｔ,ｔｈｅｉｎｔｅｒａｃｔｉｏｎｏｆｔｗｏ＿ｄｉｍｅｎｓｉｏｎａｌａｎｄｔｈｒｅ＿ｄｉｍｅ...     

A new high-accuracy scheme for compressible turbulent flows

Shock-capturing and broad-bandwidth scale resolutions are two main challenges of compressible turbulent flow simulation. To meet the rigorous requests, a novel fifth-order hybrid scheme based on a uniform hybrid framework is designed. With the help of a continuous weight operator, the new scheme combines an upwind compact scheme for smooth regions and a compact-reconstruction weighted essentially non-oscillatory scheme for discontinuous regions. Numerical analyses and canonical numerical tests confirm that the new scheme has high accuracy, spectral-like resolution property and shock-capturing capability. Besides, the new scheme shows high computational efficiency compared to the related shock-capturing schemes and hybrid ones.     

A weighted least square scheme for compressible flows

The article describes the development of a high order finite volume method for the solution of transonic flow problems. The method is based on a reconstruction procedure similar to the weighted essentially non-oscillatory scheme (WENO). The analysis of accuracy and stability of the method is carried out for the case of smooth data and for simple discontinuity. The computational results demonstrate the performance of the WLSQR method for the solution of several flow problems in 2D and 3D using both structured and unstructured meshes.     

A finite volume scheme for two-phase compressible flows

In gas-particle two-phase flows, when the concentration of the disperesed phase is low, certain assumptions may be made which simplify considerably the equations one has to solve. The gas and particle flows are then linked only via the interaction terms. One may therefore uncouple the full system of equations into two subsystems: one for the gas phase, whose homogeneous part reduces to the Euler equations; and a second system for the particle motion, whose homogeneous part is a degenerate hyperbolic system. The equations governing the gas phase flow may be solved using a high-resolution scheme, while the equations describing the motion of the dispersed phase are treated by a donor-cell method using the solution of a particular Riemann problem. Coupling is then achieved via the right-hand-side terms. To illustrate the capabilities of the proposed method, results are presented for a case specially chosen from among the most difficult to handle, since it involves certain geometrical difficulties, the treatment of regions in which particles are absent and the capturing of particle fronts.     

A segregated method for compressible flow computation Part I: isothermal compressible flows

Traditionally, coupled methods have been employed for the computation of compressible flows, whereas segregated methods have been preferred for the computation of incompressible flows. Compared to coupled methods, segregated solvers present the advantage of reduced computer memory and CPU time requirements, although at the cost of an inferior robustness. Therefore, in a series of papers we present unified computational techniques to compute compressible and incompressible flows with segregated stabilized methods. The proposed algorithms have an increased robustness compared to existing techniques, while possessing additional benefits such as employing standard pressure boundary conditions. In this first part, the thermodynamics of isothermal, thermally perfect compressible flows is set up in the framework of symmetric systems and the corresponding segregated algorithms are introduced. Copyright © 2005 John Wiley & Sons, Ltd.     

Hyperbolic system for viscous fluids and simulation of shock tube flows

A system of field equations for viscous fluids with heat conduction is formulated. Unlike the theory of Navier-Stokes fluids, a hyperbolic system of field equations in conservative form is obtained. Numerical simulation of a viscous shock tube flow is presented as an application.     

A high-order Runge-Kutta discontinuous Galerkin method with a subcell limiter on adaptive unstructured grids for two-dimensional compressible inviscid flows

A robust, adaptive unstructured mesh refinement strategy for high-order Runge-Kutta discontinuous Galerkin method is proposed. The present work mainly focuses on accurate capturing of sharp gradient flow features like strong shocks in the simulations of two-dimensional inviscid compressible flows. A posteriori finite volume subcell limiter is employed in the shock-affected cells to control numerical spurious oscillations. An efficient cell-by-cell adaptive mesh refinement is implemented to increase the resolution of our simulations. This strategy enables to capture strong shocks without much numerical dissipation. A wide range of challenging test cases is considered to demonstrate the efficiency of the present adaptive numerical strategy for solving inviscid compressible flow problems having strong shocks.     

A finite element method for compressible viscous fluid flows

This work presents a mixed three‐dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a ‘stable’ numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf–sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady‐state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright © 2010 John Wiley & Sons, Ltd.     

A lattice Boltzmann model for simulation of compressible flows

This paper presents a new model of lattice Boltzmann method for full compressible flows. On the basis of multi‐speed model, an extra potential energy distribution function is introduced to recover the full compressible Navier–Stokes equations with a flexible specific‐heat ratio and Prandtl number. The Chapman–Enskog expansion of the kinetic equations is performed, and the two‐dimension‐seventeen‐velocity density equilibrium distribution functions are obtained. The governing equations are discretized using the third order monotone upwind scheme for scalar conservation laws finite volume scheme. The van Albada limiter is used to avoid spurious oscillations. In order to verify the accuracy of this double‐distribution‐function model, the Riemann problems, Couette flows, and flows around a NACA0012 airfoil are simulated. It is found that the proposed lattice Boltzmann model is suitable for compressible flows, even for strong shock wave problem, which has an extremely large pressure ratio, 100,000. Copyright © 2014 John Wiley & Sons, Ltd.