Euler flow solutions for transonic shock wave–boundary layer interaction

Steady 2D Euler flow computations have been performed for a wind tunnel section, designed for research on transonic shock wave–boundary layer interaction. For the discretization of the steady Euler equations, an upwind finite volume technique has been applied. The solution method used is collective, symmetric point Gauss–Seidel relaxation, accelerated by non-linear multigrid. Initial finest grid solutions have been obtained by nested iteration. Automatic grid adaptation has been applied for obtaining sharp shocks. An indication is given of the mathematical quality of four different boundary conditions for the outlet flow. Two transonic flow solutions with shock are presented: a choked and a non-choked flow. Both flow solutions show good shock capturing. A comparison is made with experimental results.     

Numerical solutions of Euler equations by using a new flux vector splitting scheme

A new flux vector splitting scheme has been suggested in this paper. This scheme uses the velocity component normal to the volume interface as the characteristic speed and yields the vanishing individual mass flux at the stagnation. The numerical dissipation for the mass and momentum equations also vanishes with the Mach number approaching zero. One of the diffusive terms of the energy equation does not vanish. But the low numerical diffusion for viscous flows may be ensured by using higher-order differencing. The scheme is very simple and easy to be implemented. The scheme has been applied to solve the one dimensional (1D) and multidimensional Euler equations. The solutions are monotone and the normal shock wave profiles are crisp. For a 1D shock tube problem with the shock and the contact discontinuities, the present scheme and Roe scheme give very similar results, which are the best compared with those from Van Leer scheme and Liou–Steffen's advection upstream splitting method (AUSM) scheme. For the multidimensional transonic flows, the sharp monotone normal shock wave profiles with mostly one transition zone are obtained. The results are compared with those from Van Leer scheme, AUSM and also with the experiment.     

Stability of Transonic Shock Solutions for One-Dimensional Euler–Poisson Equations

In this paper, both structural and dynamical stabilities of steady transonic shock solutions for one-dimensional Euler–Poisson systems are investigated. First, a steady transonic shock solution with a supersonic background charge is shown to be structurally stable with respect to small perturbations of the background charge, provided that the electric field is positive at the shock location. Second, any steady transonic shock solution with a supersonic background charge is proved to be dynamically and exponentially stable with respect to small perturbations of the initial data, provided the electric field is not too negative at the shock location. The proof of the first stability result relies on a monotonicity argument for the shock position and the downstream density, and on a stability analysis for subsonic and supersonic solutions. The dynamical stability of the steady transonic shock for the Euler–Poisson equations can be transformed to the global well-posedness of a free boundary problem for a quasilinear second order equation with nonlinear boundary conditions. The analysis for the associated linearized problem plays an essential role.     

One-dimensional model for microscale shock tube flow

A one-dimensional model for the numerical simulation of transport effects in small-scale, i.e., low Reynolds number, shock tubes is presented. The conservation equations have been integrated in the lateral directions and three-dimensional effects have been introduced as carefully controlled sources of mass, momentum and energy, into the axial conservation equations. The unsteady flow of gas behind the shock wave is reduced to a quasi-steady flow by choosing a coordinate system attached to the shock. The boundary layer problem is thereby reduced to a laminar solution, similar to the Blasius solution, with the exception that the wall velocity can be nonzero. The resulting one-dimensional equations are then solved numerically using a two-step Lax-Wendroff/ MacCormack scheme with flux correction transport. For validation purposes, comparisons are performed against previously published shock structure and low Reynolds number shock tube experiments; good agreement is observed. The model has been used to predict the performance of a 10μm shock tube and the result of this simulation shows the possibility of shock wave disappearance at lower pressure ratios for a micro-scale shock tube.      

Flux vector splitting solutions for coupling hydraulic transient of gas-liquid-solid three-phase flow in pipelines

The gas-liquid-solid three-phase mixed flow is the most general in multiphase mixed transportation. It is significant to exactly solve the coupling hydraulic transient problems of this type of multiphase mixed flow in pipelines. Presently, the method of characteristics is widely used to solve classical hydraulic transient problems. However, when it is used to solve coupling hydraulic transient problems, excessive interpolation errors may be introduced into the results due to unavoidable multiwave interpolated calculations. To deal with the problem, a finite difference scheme based on the Steger-Warming flux vector splitting is proposed. A flux vector splitting scheme is established for the coupling hydraulic transient model of gas-liquid-solid three-phase mixed flow in the pipelines. The flux subvectors are then discretized by the Lax-Wendroff central difference scheme and the Warming-Beam upwind difference scheme with second-order precision in both time and space. Under the Rankine-Hugoniot conditions and the corresponding boundary conditions, an effective solution to those points located at the boundaries is developed, which can avoid the problem beyond the calculation region directly induced by the second-order discrete technique. Numerical and experimental verifications indicate that the proposed scheme has several desirable advantages including high calculation precision, excellent shock wave capture capability without false numerical oscillation, low sensitivity to the Courant number, and good stability.     

Performance Improvement of a Supercritical Airfoil by a Multi-point Optimized Shock Control Channel

A shock control channel (SCC) is a flow control method introduced here to control the shock wave/boundarylayer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. An SCC transfers an appropriate amount of mass and momentum from downstream of the shock wave location to its upstream to decrease the pressure gradient across the shock wave and as a result the shock-wave strength is reduced. Here, a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of the SCC. This flow control method is implemented on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of the SCC. The solver is validated against experimental works studying effect of cavities in transonic airfoils. It is shown that the application of an SCC improves the average aerodynamic efficiency in a range of off-design conditions by 13.2% in comparison with the original airfoil. The SCC is shown to be an effective tool also for higher angle of attack at transonic flows. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.     

Coupled fluid structure analysis for wing 445.6 flutter using a fast dynamic mesh technology

The flow around wing 445.6 was modelled using Navier–Stokes equations and S-A model. The wing vibration and flow mesh deformation were computed using a fast dynamic mesh technology proposed by our own group. Wing 445.6 flutter was analysed through a strong coupling between the wing vibration and flow. The reduced flutter velocity was predicted and results are in good agreement with the experimental data. It is found that the subsonic flutter is mainly induced by the flow separation and the transonic and supersonic flutter are mainly caused by the oscillating shock wave and its induced flow separation. The positive aerodynamic work increases due to the oscillating shock wave when the subsonic flow becomes transonic reducing the flutter velocity. While the positive aerodynamic work induced by the oscillating shock wave decreases when the transonic flow becomes supersonic increasing the flutter velocity. That is why the transonic dip exists.     

An exact non‐linear Navier–Stokes compressible‐flow solution for CFD code verification

An exact similarity solution of the compressible‐flow Navier–Stokes equations is presented, which embeds supersonic, transonic, and subsonic regions. Describing the viscous and heat‐conducting high‐gradient flow in a shock wave, the solution accommodates non‐linear temperature‐dependent viscosity as well as heat‐conduction coefficients and provides the variation of all the flow variables and their derivatives. Also presented are methods to obtain time‐dependent and/or multi‐dimensional solutions as well as verification benchmarks of increasing severity. Comparisons between the developed analytical solution and CFD solutions of the Navier–Stokes equations, with determination of convergence rates and orders of accuracy of these solutions, illustrate the utility of the developed exact solution for verification purposes. Copyright © 2012 John Wiley & Sons, Ltd.     

由激波诱导产生的二维圆柱上的非定常阻力

激波管是一种具有结构简单、参数可调范围大和使用方便等特点的气体动力学试验装置。因此,它在跨声速和超声速流的实验研究中得以广泛地应用。但由于激波管中气流的持续时间很短(通常为毫秒量级),所以在激波管中直接测量试验模型上的压力分布等气动参数将是十分困难的。本文利用马赫-曾德干涉仪在激波管中拍摄出激波在二维圆柱上的绕射和反射时的干涉图,直接给出了整个流场的密度分布。再由波后气流满足等熵关系式,可以计算得到流场的压力分布和非定常阻力的大小。     

Existence and Stability of Multidimensional Transonic Flows through an Infinite Nozzle of Arbitrary Cross-Sections

We establish the existence and stability of multidimensional steady transonic flows with transonic shocks through an infinite nozzle of arbitrary cross-sections, including a slowly varying de Laval nozzle. The transonic flow is governed by the inviscid potential flow equation with supersonic upstream flow at the entrance, uniform subsonic downstream flow at the exit at infinity, and the slip boundary condition on the nozzle boundary. Our results indicate that, if the supersonic upstream flow at the entrance is sufficiently close to a uniform flow, there exists a solution that consists of a C ^1,α subsonic flow in the unbounded downstream region, converging to a uniform velocity state at infinity, and a C ^1,α multidimensional transonic shock separating the subsonic flow from the supersonic upstream flow; the uniform velocity state at the exit at infinity in the downstream direction is uniquely determined by the supersonic upstream flow; and the shock is orthogonal to the nozzle boundary at every point of their intersection. In order to construct such a transonic flow, we reformulate the multidimensional transonic nozzle problem into a free boundary problem for the subsonic phase, in which the equation is elliptic and the free boundary is a transonic shock. The free boundary conditions are determined by the Rankine–Hugoniot conditions along the shock. We further develop a nonlinear iteration approach and employ its advantages to deal with such a free boundary problem in the unbounded domain. We also prove that the transonic flow with a transonic shock is unique and stable with respect to the nozzle boundary and the smooth supersonic upstream flow at the entrance.     

Focusing of a shock wave in a rarefied gas A numerical study

The process of focusing of a shock wave in a rarefied noble gas is investigated by a numerical solution of the corresponding two dimensional initial–boundary value problem for the Boltzmann equation. The numerical method is based on the splitting algorithm in which the collision integral is computed by a Monte Carlo quadrature, and the free flow equation is solved by a finite volume method. We analyse the development of the shock wave which reflects from a suitably shaped reflector, and we study influence of various factors, involved in the mathematical model of the problem, on the process of focusing. In particular, we investigate the pressure amplification factor and its dependence on the strength of the shock and on the accommodation coefficient appearing in the Maxwell boundary condition modelling the gas-surface interaction. Moreover, we study the dependence of the shock focusing phenomenon on the shape of the reflector, and on the Mach number of the incoming shock. Received 25 May 1998 / Accepted 4 January 2000     

Application of unsteady mixing equations to some aerodynamic problems

Tangential discontinuities [1] are introduced in solving several transient and steady-state problems of gas dynamics. These discontinuities are unstable [2] as a result of the effects of viscosity and thermal conductivity. Therefore it is advisable to replace the tangential discontinuity by a mixing region and account for its interaction with the inviscid flows, establishing on the boundaries of this region the conditions of vanishing friction stress and equality of the velocity and temperature components to the corresponding velocity and temperature components of the inviscid flows. This formulation improves the accuracy of the solution of such problems by posing them as problems with irregular reflection and intersection of shock waves [1].The consideration of the interaction of unsteady turbulent mixing regions with the inviscid flow also permits the formulation of several problems in which the effects of viscosity lead to complete rearrangement of the flow pattern (the lambda-configuration) with the interaction of the reflected shock wave with the boundary layer in the shock tube [3,4], the formation of zones of developed separation ahead of obstacles, etc.).In this connection, §1 presents an analysis of the self-similar solutions of the unsteady turbulent mixing equations (a corresponding analysis of the laminar mixing equations which coincide with the boundary layer equations is presented in [1]). It is shown that these self-similar solutions describe, along with the several problems noted above, the problems of the formation of steady jets and mixing zones in the base wake.As an example, §2 presents, within the framework of the proposed schematization, an approximate solution of the problem of the interaction of a shock wave reflected from a semi-infinite wall with the boundary layer on a horizontal plate behind the incident shock wave. The results obtained are applied to the analysis of reflection in a shock tube. Computational results are presented which are in qualitative agreement with experiment [3, 4].     

A knife-edge load traveling on the surface of an elastic halfspace

A load moving on the surface of an elastic halfspace forms a basic problem that is related to different fields of engineering, such as the subsoil response due to vehicle motion or the ultrasound field due to an angle beam transducer. Many numerical techniques have been developed to solve this problem, but these do not provide the fundamental physical insights that are offered by closed form solutions, which are very rare in comparison. This paper describes the development and analysis of the closed form space-time domain solution for a knife-edge load, i.e. a line segment of normal traction, moving at a constant speed on the surface of an elastic halfspace. The various contributions making up the exact solution, obtained with the Cagniard-De Hoop method, produce all the complicated wave patterns from this distributed type of loading. Examples are the transient wave field at the starting position of the load, angled conical and plane waves propagating into the solid, Rayleigh waves propagating along the surface, and head waves spreading and attenuating in specific directions from the loading path. The influence of the load speed on the wave field is investigated by considering the singularities in the relevant complex domains, for each sonic range relative to the bulk wave velocities. The characteristic wave fronts and wave patterns as exhibited by the particle displacements are evaluated for subsonic, transonic and supersonic load speeds.     

Computation of shock wave diffraction at a ninety degrees convex edge

This paper presents a numerical study of shock wave diffraction at a sharp ninety degrees edge, using an explicit second-order Godunov-type Euler scheme based upon the solution of a generalized Riemann problem (GRP). The Euler computations produce flow separation very close to the diffraction edge, leading to a realistic development of the separated shear layer and subsequent vortex roll-up. The diffracted shock wave, and the secondary shock wave, are both reproduced well. In addition a pair of vortex shocks are shown to form, extending well into the vortex core.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Delta Shock Wave Solution of the Riemann Problem for the Non-homogeneous Modified Chaplygin Gasdynamics

The motivation of the present study is to derive the solution of the Riemann problem for modified Chaplygin gas equations in the presence of constant external force. The analysis leads to the fact that in some special circumstances delta shock appears in the solution of the Riemann problem. Also, the Rankine–Hugoniot relations for delta shock wave which are utilized to determine the strength, position and propagation speed of the delta shocks have been derived. Delta shock wave solution to the Riemann problem for the modified Chaplygin gas equation is obtained. It is found that the external force term, appearing in the governing equations, influences the Riemann solution for the modified Chaplygin gas equation.

     

爆燃波在含内构件管道中传播现象的实验研究

为了评估高温气冷核反应堆热交换器H2泄漏、爆炸的安全性,研究含内构件管道的H2/空气爆燃传播现象,建造了几何相似、尺寸相同的实验管道（真空筒）。分别充入不同初压和当量比H2/空气混合物,在真空筒顶部点火并引发爆燃,利用多通道瞬态压力测量和数据采集系统,记录各测点压力时间曲线。结果表明:对化学计量比H2/空气混合物,在慢化剂室和真空筒顶部空间产生爆燃,邻近测点的压力时间曲线显示了冲击波特征。该冲击波通过慢化剂室和真空筒侧壁的狭缝（2.5 mm）,进入含内构件的扩张管道并形成爆燃。冲击波在真空筒端部反射、向后传播并与火焰相互作用,爆炸流场波系复杂。对富油和低初压化学计量比混合物,在慢化剂室和真空筒顶部空间产生燃烧,高温富油燃气的压力上升速率较慢。当燃气通过上述狭缝时,在真空筒突扩空间内再次点火并形成较强爆燃,压力时间曲线显示了冲击波特征及其在端面的反射。     

Effect of surface gas condensation on the velocity of a reflected shock wave

The time-dependent one-dimensional problem of the normal reflection of a shock wave propagating at constant velocity in a gas (vapor) at rest from the plane surface of its condensed phase under steady-state condensation-evaporation conditions on the interphase plane is considered within the framework of the kinetic equation for a monatomic gas with a model collision operator (S-model). The solution is obtained using a conservative second-order finite-difference method. Attention is concentrated on the steady-state regime of the condensation process. The effect of the condensation (evaporation) coefficient on the velocity of the reflected shock wave is studied.     

An efficient method for solving the mixed direct-inverse problem of the transonic rotational flow in plane casades

The method in [1] has been extended to the case of rotational flow in this paper. A new method for dealing with the shock wave is presented. This method has the advantages of both the shock-fitting and the shock capturing methods. The direct problem and the mixed direct-inverse problem of the rotational flow in a transonic plane cascade at both design and off design conditions are solved, and the results show that the present method has rapid convergence rate and high accuracy even for the flow with moderately strong shocks. The calculations have been carried out on the DPS-8 computer, and for the direct problem, only 50–80 iterations are needed, and 50–80 seconds of CPU time are required.     

An analysis of the stability of the least square finite difference scheme and its shock capturing ability in inviscid flows

A meshless method – The Least Square Finite Difference scheme (LSFD) with diffusion is analyzed and applied to inviscid flows. The scheme is made second-order by using a modified difference in the formulation of LSFD. Several numerical experiments, namely the Sod shock tube and the shallow water problems, are carried out and, in the limelight of the results obtained, the ability of the scheme to resolve shock wave, rarefaction wave, and contact discontinuity is discussed. The conditional stability of the LSFD scheme is established. The LSFD uses weights to diagonalize the least square matrix resulting in the spatial discretization in order to gain computational time. We prove that there exists a unique weight for the resulting optimization problem. The weighted version of LSFD is used to solve the isentropic vortex problem numerically and the results are used to discuss the dissipative nature of the scheme. Five configurations of the two-dimensional Riemann problems are used in our numerical experiments. The capability of the scheme to capture the complex interaction of multiple planar waves is discussed in the limelight of the results on the Riemann problems. The result of the shock reflection problem shows that the scheme is minimally dissipative and leads to sharp and well-resolved shocks.