Numerical simulation of the flow with a pseudo-shock in an axisymmetric expanding duct with a frontal inlet

The results of the numerical modeling of a flow with a pseudo-shock in an axisymmetric duct are presented. The duct included a frontal inlet with the initial funnel-shaped compression part and the cylindrical throat part as well as the subsequent expanding diffuser. To create a flow with a pseudo-shock, the duct was throttled with the use of the outlet converging insert. Numerical computations of the axisymmetric flow have been conducted on the basis of the solution of the Reynolds-averaged Navier?Stokes equations and with the use of the k-ω SST turbulence model. As a result of computations, such parameters of the flow were determined as the location of the beginning of the pseudo-shock, the length of its supersonic part, the velocity profiles in different cross sections of the pseudo-shock, the pressure distribution on the duct wall, the total pressure recovery factor, and others. The behavior of these parameters at the freestream Mach number М = 6 was analyzed versus the diffuser opening angle and different degree of the inlet duct throttling.     

Numerical and experimental investigations on the mach 2 pseudo-shock wave in a square duct

This paper describes numerical and experimental investigations for the multiple shock wave/turbulent boundary layer interaction in a Mach 2 supersonic square duct. The numerical simulation is carried out with the Harten-Yee second-order accuracy TVD scheme and the Baldwin-Lomax turbulence model. The flow conditions are a free-stream Mach number ofM _∞≈=2.0 and a Reynolds number ofRe _∞;=2.5×10⁷ and the flow confinements are δ_∞/h=0.15 (case A) and δ_∞/h=0.25 (case B), respectively. The computational results for both cases show good agreement with the experimental results. Based on these agreements, the flow quantities, which are very difficult to obtain experimentally, are analyzed by numerical simulation. Moreover, the effect of flow confinement on the pseudo-shock wave characteristics is also presented.     

Numerical simulation of the influence of the control surfaces deflection on the aerodynamics of a slender axisymmetric configuration

A possible influence of the deflection of control surfaces on the aerodynamics of an axisymmetric slender configuration at supersonic flow speeds is considered. A classical configuration consisting from the fuselage in the form of a body of revolution and having cross frontal fins and six-blade trailing stabilizers is considered as the investigation object. The physical flow pattern at the deflection of horizontal fin consoles is investigated and the estimates are obtained for the influence of this deflection on both the characteristics of elements (the body and stabilizers) as well as on the integral aerodynamic characteristics of the entire configuration. Numerical computations of the flow have been done at the freestream Mach number М = 3 in the range of attack angles α = 0?10° and the angles of the control surfaces deflection δ _cs = ±5° on the basis of the averaged Navier?Stokes equations and the SST k-ω turbulence model.     

Investigation of the pseudo-shock wave in a two-dimensional supersonic inlet

Abstract  

This paper describes experimental and numerical investigations into the multiple shock waves/turbulent boundary layer interaction in a supersonic inlet. The test model has a rectangular shape with an asymmetric subsonic diffuser of 5°. Experiments were conducted to obtain the visualization images and static pressure data by using supersonic wind tunnel. Numerical simulation was performed by solving the RANS equations with the Menter’s SST turbulent model. The inflow condition was a free-stream Mach number of 2.5 and a unit Reynolds number of 7.6 × 10⁷/m. Numerical results showed good agreement with the experimental results. Based on this agreement, the flow characteristics which are often very difficult to obtain experimentally alone were analyzed with the aid of numerical simulation. The structures, pressure and velocity distributions, and total pressure loss of the pseudo-shock wave in the supersonic inlet were presented in detail from flow visualization images and static pressures.     

Numerical simulation of the side pylon effect on the aerodynamic characteristics of models at their wind tunnel tests

The possible influence of fastening the models on a side pylon at their tests in wind tunnels on their aerodynamics at supersonic flow speeds has been considered. The physical problem of the pylon and the model interference has been investigated, and the estimates of the pylon influence on integral aerodynamic characteristics have been obtained. The numerical computations of the flow have been done using the averaged Navier-Stokes equations and the SST k-ω turbulence model in the range of freestream Mach numbers M = 2.5-5. As the investigation object the “classical” body of revolution of large aspect ratio is considered, which has a cruciform forward fins and six-blade tail stabilizers.     

Shock wave in the flow field of a source and a vortex source in the energy supply zone at thermal crisis

Thermal crisis of a vortex source outflowing initially in regime I into a rarefied space (into vacuum) with a transition of the supersonic flow into the subsonic flow in the shock wave, and into a stagnant space in regime II with final stagnation is considered in the model of a perfect gas with a constant heat capacity. The shock wave can be located in the energy supply zone or outside the energy release zone depending on the preset total pressure at infinity. In the absence of circulation, a cylindrical source is compared with a spherical source. The dependences of energy parameters and temperature, as well as the total pressure and density, on the coordinates of the shock wave are considered. The dependences of the critical parameters of the flow in the wake behind the zone on the coordinate of the heat supply zone, its length, and gas circulation in the cylindrical vortex source are analyzed.     

Numerical investigation of supersonic flow breakdown at the inlet duct throttling

The work presents the results of investigating the process of supersonic flow deceleration in a duct of the two-dimensional inlet throttled by variation of the outlet cross-sectional area. An inlet with three external compression shock waves designed for the freestream Mach number M_d = 7 was considered as an example for the investigation. A one-dimensional analysis of the conditions for realization of the supersonic flow deceleration regimes in the inlet duct with two throats — in the inlet entrance and at the inlet duct outlet, has been carried out. The parametric numerical computations of two-dimensional inviscid or turbulent flows in the inlet were performed with the use of the Euler and Navier—Stokes codes of the program package FLUENT. The critical conditions for the nonuniform flow in the outlet throat bringing to choking the inlet duct were determined.     

Turbulent structure in supersonic jets with a high Reynolds number

The turbulent properties of a supersonic jet were studied related to a high level of pressure pulsation found in model jets of a reentry flight vehicle approaching the landing ground. This study comprised measurements of total pressure at a small-size target using a dynamic pressure probe placed in a free jet. The most comprehensive data about jet turbulence can be obtained by direct transformation of the pressure reading at the stagnation point near the target into the normalized velocity. The oscillogram of normalized velocity produces the velocity average value, root-mean-square value as well as turbulence intensity and turbulence spectrum. It was demonstrated that a high level of turbulence for a high-head jet retains along the supersonic core length and at the beginning of subsonic interval.     

Design of supersonic three-dimensional inlets using two-dimensional isentropic compression flow

The design of supersonic three-dimensional inlets using the V-shaped body forming a two-dimensional flow including an initial oblique shock wave and a subsequent isentropic compression wave is considered. Such a flow appears attractive for inlets design due to a possibility of obtaining high compression levels of external flow over the inlet ramp with small losses of the total pressure. Numerical computations of the flows around the designed configurations were carried out in design and off-design regimes using Euler code. The flow structure was identified, the aerodynamic characteristics of the inlets were determined. The investigation covers the range of supersonic speeds corresponding to the freestream Mach numbers M_∞= 1.8−2.5.     

Supersonic diffuser for pressure recovery in SCOIL system

Classical supersonic chemical oxygen iodine laser (SCOIL) systems operate under a low total pressure of nearly 18 Torr (2400 Pa) with cavity pressure being in the range 3 Torr (400 Pa) and Mach number of 1.7. These systems handle high flow rates and hence an efficient supersonic diffuser (SD) is a critical first step towards an open-cycle operation, which may be followed by a multi-stage ejector system. The present study discusses the various aspects in the design of a supersonic diffuser for a twin 10 kW COIL module source which employs flow rates of 100 gs⁻¹ in each module. The results of computational studies based on 3-D, viscos compressible flow, k-ε turbulence formulation for the supersonic diffuser geometry have also been discussed. The experimental results from a single-module test of the supersonic diffuser show that a total recovered pressure of nearly 7 Torr is achieved at the diffuser exit.     

On the role of pressure in elasto-inertial turbulence

The dynamics of elasto-inertial turbulence is investigated numerically from the perspective of the coupling between polymer dynamics and flow structures. In particular, direct numerical simulations of channel flow with Reynolds numbers ranging from 1000 to 6000 are used to study the formation and dynamics of elastic instabilities and their effects on the flow. Based on the splitting of the pressure into inertial and polymeric contributions, it is shown that the polymeric pressure is a non-negligible component of the total pressure fluctuations, although the rapid inertial part dominates. Unlike Newtonian flows, the slow inertial part is almost negligible in elasto-inertial turbulence. Statistics on the different terms of the Reynolds stress transport equation also illustrate the energy transfers between polymers and turbulence and the redistributive role of pressure. Finally, the trains of cylindrical structures around sheets of high polymer extension that are characteristics of elasto-inertial turbulence are shown to be correlated with the polymeric pressure fluctuations.     

Influence of the opening angle of a conical supersonic nozzle on the structure of initial interval of non-isobaric jet

The ideal gas exhaustion from an infinite volume into a gas at rest through a supersonic conical Laval nozzle is considered. The problem was solved numerically by steadying in time in a unified formulation for the regions inside the nozzle and in the ambient environment. In such a statement, the nozzle outlet section is no internal boundary of the region under consideration, and there is no need of specifying the boundary conditions here. Local subsonic zones arising in the flow lie inside the region under consideration, which eliminates the possibility of using a marching technique along one of the coordinates. The numerical solution is constructed by a unified algorithm for the entire flow region, which gives a possibility of obtaining a higher accuracy. The computations are carried out in the jet initial interval, where, according to monograph [1], the wave phenomena predominate over the viscous effects. The exhaustion process is described by the system of gas dynamics equations. Their solution is constructed with the aid of a finite difference Harten’s TVD (Total Variation Diminishing) scheme [2], which has the second approximation order in space. The second approximation order in time is achieved with the aid of a five-stage Runge-Kutta method. The solution algorithm has been parallelized in space and implemented on the multi-processor computer systems of the ITAM SB RAS and the MVS-128 of the Siberian Supercomputer Center of SB RAS. The influence of the semi-apex angle of the nozzle supersonic part and the pressure jump between the nozzle outlet section and the ambient environment on the flow in the initial interval of a non-isobaric jet is investigated in the work. A comparison with experimental data is presented. The computations are carried out for the semi-apex angles of the nozzle supersonic part from 0 (parallel flow) to 20 degrees. For all considered nozzles, the Mach number in the nozzle outlet section, which was computed from the one-dimensional theory, equaled three. Computations showed that in the case of flow acceleration in a conical supersonic nozzle, its geometry is one of the main factors determining the formation of the jet initial interval in ambient environment.     

Computation of streamwise vorticity in a compressible flow of a winglet nozzle-based COIL device

Chemical oxygen iodine laser (COIL) is a high-power laser with potential applications in both military as well as in the industry. COIL is the only chemical laser based on electronic transition with a wavelength of 1.315 μm, which falls in the near-infrared (IR) range. Thus, COIL beam can also be transported via optical fibers for remote applications such as dismantling of nuclear reactors. The efficiency of a supersonic COIL is essentially a function of mixing specially in systems employing cross-stream injection of the secondary lasing (I₂) flow in supersonic regime into the primary pumping (O₂¹Δ_g) flow. Streamwise vorticity has been proven to be among the most effective manner of enhancing mixing and has been utilized in jet engines for thrust augmentation, noise reduction, supersonic combustion, etc. Therefore, a computational study of the generation of streamwise vorticity in the supersonic flow field of a COIL device employing a winglet nozzle with various delta wing angles of 5°, 10°, and 22.5° has been carried out. The study predicts a typical Mach number of approximately 1.75 for all the winglet geometries. The analysis also confirms that the winglet geometry doubles up both as a nozzle and as a vortex generator. The region of maximum turbulence and fully developed streamwise vortices is observed to occur close to the exit, at x/λ of 0.5, of the winglets making it the most suitable region for secondary flow injection for achieving efficient mixing. The predicted length scale of the scalloped mixer formed by the winglet nozzle is 4λ. Also, the winglet nozzle with 10° lobe angle is most suitable from the point of view of mixing developing cross-stream velocity of 120 m/s with acceptable pressure drop of 0.7 Torr. The winglet geometry with 5° lobe angle is associated with a low cross-stream velocity of 60 m/s, whereas the one with 22.5° lobe angle is associated with a large static and total pressure drop of 1.87 and 9.37 Torr, respectively, making both the geometries unsuitable for COIL systems. The experimental validation shows a close agreement with the computationally predicted values. The studies for the most suitable 10° lobe angle geometry show an observed Mach number of 1.72 with an improved mixing efficiency of 74% due to the occurrence of predicted streamwise vortices in the flow.     

光在超声速湍流边界层中的传输

采用数值方法分析超声速边界层流场对光传输的影响,并根据光波的变化对流场特征进行分析.用法福尔质量加权平均N-S方程及两方程湍流模型求解三维超声速平板湍流边界层流动;光在流场中的传输采用傍轴近似光波传输方程描述,用相屏法和快速傅里叶变换(FFT)技术求解.利用穿越流场的光波光强和相位的改变,直观分析光波畸变和流场的部分特征信息.     

Flame front analysis of high-pressure turbulent lean premixed methane–air flames

An experimental study on lean turbulent premixed methane–air flames at high pressure is conducted by using a turbulent Bunsen flame configuration. A single equivalence ratio flame at Φ = 0.6 is explored for pressures ranging from atmospheric pressure to 0.9 MPa. LDA measurements of the cold flow indicate that turbulence intensities and the integral length scale are not sensitive to pressure. Due to the decreased kinematic viscosity with increasing pressure, the turbulent Reynolds numbers increase, and isotropic turbulence scaling relations indicate a large decrease of the smallest turbulence scales. Available experimental results and PREMIX code computations indicate a decrease in laminar flame propagation velocities with increasing pressure, essentially between the atmospheric pressure and 0.5 MPa. The u′/S_L ratio increases therefore accordingly. Instantaneous flame images are obtained by Mie scattering tomography. The images and their analysis show that pressure increase generates small scale flame structures. In an attempt to generalize these results, the variance of the flamelet curvatures, the standard deviation of the flamelet orientation angle, and the flamelet crossing lengths have been plotted against which is proportional to the ratio between the integral and Taylor length scales, and which increases with pressure. These three parameters vary linearly with the ratio between large and small turbulence scales and clearly indicate the strong effect of this parameter on premixed turbulent flame dynamics and structure. An obvious consequence is the increase in flame surface density and hence burning rate with pressure, as confirmed by its direct determination from 2D tomographic images.     

����������IJ���������ģ����ֵģ��

采用FLUENT软件分别对外加均匀横向磁场的等截面三维充分发展液态金属管流的层流模型和低雷诺数湍流Lam/Bremhost(LB)模型进行了数值模拟,分析了外加磁场对普通方管LB模型速度分布和压降的影响。比较在相同哈特曼数下,层流和湍流模型方管截面上速度分布和管道中MHD压降。其中,对电流的计算采用磁感应方程来求得。数值模拟结果证明了用低雷诺数LB湍流模型解决方管磁流体流动的可行性。通过层流模型和湍流模型的对比可知,层流模型有较短的入口长度,但管内流体的压降却很大;而湍流模型管内速度更加平均化,管内压降较小,但管内入口长度较长。     

Direct numerical simulation of a supersonic turbulent boundary layer on a flat plate and its analysis

In the past three decades, considerable progress has been made in the investigation of incompressible turbulent boundary layer through experiments, DNS and theoretical works, including: (1) the statistics characteristic and structure of turbulence; (2) the co-herent structures in turbulent flows; (3) turbulence modeling and the large eddy simula-tion (LES). In contrast, the progress was very slow for the compressible, in particular, the super-sonic turbulent boundary layer. Recent works on d…     

Effect of inlet flow turbulence on the combustion instability in a premixed backward-facing step combustor

This paper reports the effect of inlet flow turbulence intensity on the combustion instability characteristics in a backward facing step combustor. The inlet turbulence intensity is varied by a turbulence generator. Unsteady pressure measurements and OH* chemiluminescence images are recorded over a wide range of operating conditions at different inlet turbulence intensities. The study shows an early onset of instability at low turbulence level, i.e., higher turbulence postpones the onset of instability to higher Reynolds number Re and/or higher equivalence ratio Φ. The early onset of instability in the Re and Φ parameter spaces is due to the change in system parameters such as flame speed and size of the recirculation zone downstream of the step at different turbulence levels. Further, the onset is characterized as subcritical bifurcation. At low Re, the hysteresis zone width is small for low turbulence levels and it is large at higher turbulence levels; and at higher Re, the hysteresis width remains constant at all turbulence levels. Investigation of instability characteristics reveals that there are momentary slippages from limit cycle orbit into brief silent regimes in an intermittent manner. The frequency of occurrence of the momentary silent regimes increases with reduction in turbulence, indicating that higher turbulence helps in maintaining the system in a stable limit cycle orbit. High-speed chemiluminescence imaging reveals the necessity of the vortex rollup in the recirculation zone to grow up to the top wall by dilatation from the heat release for the onset of instability. Considerations of the effect of turbulence on both the flame speed and the recirculation zone size together explain all the observed bifurcation trends. These results suggest that inlet flow turbulence should not just be considered as background noise. The turbulence effects on both the flame and flow should be considered in predicting the instability characteristics.     

Effect of cavity location on combustion flow field of integrated hypersonic vehicle in near space

Abstract  

The cavity has been widely employed as the flame holder to prolong the residence time of fuel in supersonic flows since it improves the combustion efficiency in the scramjet combustor, and also imposes additional drag on the engine. In this paper, the two-dimensional coupled implicit Reynolds Average Navier–Stokes equations, the RNG k–ε turbulence model and the finite-rate/eddy-dissipation reaction model have been employed to numerically simulate the combustion flow field of an integrated hypersonic vehicle. The effect of cavity location on the combustion flow field of the vehicle has been investigated, and the fuel, namely hydrogen, was injected upstream of the cavity on the walls of the first stage combustor. The obtained results show that the viscous lift force, drag force and pitching moment of the vehicle are nearly unchanged by varying the cavity location over the location range and designs considered in this article, namely the configurations with single cavity, double cavities in tandem and double cavities in parallel. The variation of the fuel injection strategy affects the separation of the boundary layer, and the viscous effect on the drag force of the vehicle is remarkable, but the viscous effects on the lift force and the pitching moment are both small and they can be neglected in the design process of hypersonic vehicles. In addition to varying the location of the cavities, three fuel injection configurations were considered. It was found that one particular case can restrict the inlet unstart for the scramjet engine.     

Numerical investigation of the physics of rotating-detonation-engines

Rotating-detonation-engines (RDE’s) represent an alternative to the extensively studied pulse-detonation-engines (PDE’s) for obtaining propulsion from the high efficiency detonation cycle. Since it has received considerably less attention, the general flow-field and effect of parameters such as stagnation conditions and back pressure on performance are less well understood than for PDE’s. In this article we describe results from time-accurate calculations of RDE’s using algorithms that have successfully been used for PDE simulations previously. Results are obtained for stoichiometric hydrogen–air RDE’s operating at a range of stagnation pressures and back pressures. Conditions within the chamber are described as well as inlet and outlet conditions and integrated quantities such as total mass flow, force, and specific impulse. Further computations examine the role of inlet stagnation pressure and back pressure on detonation characteristics and engine performance. The pressure ratio is varied between 2.5 and 20 by varying both stagnation and back pressure to isolate controlling factors for the detonation and performance characteristics. It is found that the detonation wave height and mass flow rate are determined primarily by the stagnation pressure, whereas overall performance is closely tied to pressure ratio. Specific impulses are calculated for all cases and range from 2872 to 5511 s, and are lowest for pressure ratios below 4. The reason for performance loss is shown to be associated with the secondary shock wave structure that sets up in the expansion portion of the RDE, which strongly effects the flow at low pressure ratios. Expansion to supersonic flow behind the detonation front in RDE’s with higher pressure ratios isolate the detonation section of the RDE and thus limit the effect of back pressure on the detonation characteristics.