Control of Nozzle Flow Using Rectangular Ribs at Sonic and Supersonic Mach Numbers

This study deals with base pressure management in a duct for various values of the Mach number(M),namely,Mach number corresponding to sonic and four supersonic conditions.In addition to the Mach number,the nozzle pressure ratio(NPR),the area ratio,the rib dimension,and the duct length are influential parameters.The following specific values are examined at M=1,1.36,1.64,and 2,and NPRs between 1.5 and 10.The base pressure is determined by positioning ribs of varying heights at predetermined intervals throughout the length of the square duct.When the level of expansion is varied,it is seen that the base pressure initially drops for overexpanded flows and increases for under-expanded flows.When ribs are present,the flow field in the duct and pressure inside the duct fluctuate as the base pressure rises.Under-expanded flows can achieve a base pressure value that is suitably high without experiencing excessive changes in the duct flow in terms of static pressure if a rib height around 10%of the duct height close to the nozzle exit is considered.Rectangular rib passive control does not negatively affect the duct’s flow field.

NUMERICAL SIMULATION OF SUPERSONIC AXISYMMETRIC FLOW OVER MISSILE AFTERBODY WITH JET EXHAUST USING POSITIVE SCHEMES

Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.

PARALLELIZED UPWIND FLUX SPLITTING SCHEME FOR SUPERSONIC REACTING FLOWS ON UNSTRUCTURED HYBRID MESHES

A parallelized upwind flux splitting scheme for supersonic reacting flows on hybrid meshes is presented. The complexity of super/hyper-sonic combustion flows makes it necessary to establish solvers with higher resolution and efficiency for multi-component Euler/N-S equations. Hence, a spatial second-order van Leer type flux vector splitting scheme is established by introducing auxiliary points in interpolation, and a domain decomposition method used on unstructured hybrid meshes for obtaining high calculating efficiency. The numerical scheme with five-stage Runge-Kutta time step method is implemented to the simulation of combustion flows, including the supersonic hydrogen/air combustion and the normal injection of hydrogen into reacting flows. Satisfying results are obtained compared with limited references.

4D-Flow MRI在肥厚型心肌病左室流出道血流评估中的价值探索

目的 探索四维血流(four-dimensional flow,4D-Flow)磁共振成像(magnetic resonance imaging,MRI)技术在左心室腔内应用的可行性。材料与方法 本研究为前瞻性、横断面研究,纳入2022年8月至2023年1月于我院接受心脏MRI检查的21例肥厚型心肌病患者,采用3.0 T MRI扫描仪进行二维血流(tow-dimensional flow,2D-Flow)及4D-Flow成像,收集患者一周内进行的超声心动图检查结果。采用组内相关系数(inter-class correlation coefficient,ICC)、变异系数(coefficients of variation,COV)及Bland-Altman分析比较2D-Flow、4D-Flow评估左室流出道峰值流速的可重复性及一致性,并通过Pearson相关性分析探究二者与超声心动图测量结果的关系。结果 2D-Flow和4D-Flow观察者内/观察者间的ICC分别为0.999/0.999和0.995/0.992,COV分别为0.5%/0.5%和2.4%/2.6%。4D-Flow与超声心动图的测量结果呈中度相关,相关系数r值为0.574(P=0.006),但一致性较差,ICC为0.375(P=0.013)。2D-Flow与4D-Flow和超声心动图间无显著的一致性及相关性。结论 4D-Flow技术能够可视化心腔内血流模式,对左室流出道峰值流速的测量具有高度可重复性,且与超声心动图的测量结果具有显著的一致性。

CHAOTIC MOTIONS AND LIMIT CYCLE FLUTTER OF TWO-DIMENSIONAL WING IN SUPERSONIC FLOW

Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.

Experimental investigations of detonation initiation by hot jets in supersonic premixed flows

A new method to initiate and sustain the detonation in supersonic flow is investigated. The reaction activity of coming flow may influence the result of detonation initiation. When a hot jet initiates a detonation wave successfully, there may exist two types of detonations. If the detonation velocity is greater than the velocity of coming flow, there will be a normal detonation here. Because of the influence of boundary layer separation, the upstream detonation velocity is much greater than the Chapman-Jouguet （C J） detonation velocity. On the other hand, if the detonation velocity is less than the velocity of coming flow, an oblique detonation wave （ODW） will form. The ODW needs a continuous hot jet to sustain itself. If the jet pressure is lower than a certain value, the ODW will decouple. In contrast, the normal detonation wave can sustain itself without the hot jet.

Instantaneous and time-averaged flow structures around a blunt double-cone with or without supersonic film cooling visualized via nano-tracer planar laser scattering

In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scattering （NPLS）, which has a high spatiotemporal resolution. Three experimental cases with different injection mass flux rates were carried out. Many typical flow structures were clearly shown, such as shock waves, expansion fans, shear layers, mixing layers, and turbulent boundary layers. The analysis of two NPLS images with an interval of 5 us revealed the temporal evolution characteristics of flow structures. With matched pressures, the laminar length of the mixing layer was longer than that in the case with a larger mass flux rate, but the full covered region was shorter. Structures like K-H （Kelvin-Helmholtz） vortices were clearly seen in both flows. Without injection, the flow was similar to the supersonic flow over a backward- facing step, and the structures were relatively simpler, and there was a longer laminar region. Large scale structures such as hairpin vortices were visualized. In addition, the results were compared in part with the schlieren images captured by others under similar conditions.

Mode transition and oscillation suppression in supersonic cavity flow

Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation （DNS）. For an uncontrolled cavity, as the thickness of the boundary layer declines, transition of the dominant mode from the steady mode to the Rossiter Ⅱ mode and then to the Rossiter III mode is observed due to the change of vortex-corner interactions. Meanwhile, a low frequency mode appears. However, the wake mode observed in a subsonic cavity flow is absent in the current simulation. The oscillation frequencies obtained from a global dynamic mode decomposition （DMD） approach are consistent with the local power spectral density （PSD） analysis. The dominant mode transition is clearly shown by the dynamic modes obtained from the DMD. A passive control technique of substituting the cavity trailing edge with a quarter-circle is studied. As the effective cavity length increases, the dominant mode transition from the Rossiter Ⅱ mode to the Rossiter Ⅲ mode occurs. With the control, the pressure oscillations are reduced significantly. The interaction of the shear layer and the recirculation zone is greatly weakened, combined with weaker shear layer instability, responsible for the suppression of pressure oscillations. Moreover, active control using steady subsonic mass injection upstream of a cavity leading edge can stabilize the flow.

An efficient unstructured WENO method for supersonic reactive flows

An efficient high-order numerical method for supersonic reactive flows is proposed in this article.The reactive source term and convection term are solved separately by splitting scheme.In the reaction step,an adaptive time-step method is presented,which can improve the efficiency greatly.In the convection step,a third-order accurate weighted essentially non-oscillatory（WENO）method is adopted to reconstruct the solution in the unstructured grids.Numerical results show that our new method can capture the correct propagation speed of the detonation wave exactly even in coarse grids,while high order accuracy can be achieved in the smooth region.In addition,the proposed adaptive splitting method can reduce the computational cost greatly compared with the traditional splitting method.

Nonlinear aeroelastic analysis of a two-dimensional wing with control surface in supersonic flow

Based on the piston theory of supersonic flow and the energy method, a two dimensional wing with a control surface in supersonic flow is theoretically modeled, in which the cubic stiffness in the torsional direction of the control surface is considered. An approximate method of the cha- otic response analysis of the nonlinear aeroelastic system is studied, the main idea of which is that under the condi- tion of stable limit cycle flutter of the aeroelastic system, the vibrations in the plunging and pitching of the wing can approximately be considered to be simple harmonic excita- tion to the control surface. The motion of the control surface can approximately be modeled by a nonlinear oscillation of one-degree-of-freedom. The range of the chaotic response of the aeroelastic system is approximately determined by means of the chaotic response of the nonlinear oscillator. The rich dynamic behaviors of the control surface are represented as bifurcation diagrams, phase-plane portraits and PS diagrams. The theoretical analysis is verified by the numerical results.

Characteristics of supersonic flow in new type externally pressurized spherical air bearings

In order to predict accurately the characteristics of supersonic flow in new type externally pressurized spherical air bearings under large bearing clearance and high air supply pressure, which could decrease their load carrying capacity and stability, a CFD-based analysis was introduced to solve the three-dimensional turbulent complete compressible air flow governing equations. The realizable κ-ε model was used as a turbulent closure illustrate that the interaction exists between shock waves The supersonic flow field near air inlets was analyzed. The flow structures and boundary layer, and the flow separation is formed at the lower comer and the lower wall around the point of a maximum velocity. The numerical results show that the conversion from supersonic flow to subsonic flow in spherical air bearing occurs through a shock region （pseudo-shock）, and the viscous boundary layer results in the flow separation and reverse flow near the shock. The calculation results basically agree with the corresponding experimental data.

Model reduction for supersonic cavity flow using proper orthogonal decomposition(POD)and Galerkin projection

The reduced-order model （ROM） for the two-dimensional supersonic cavity flow based on proper orthogonal decomposition （POD） and Galerkin projection is investigated. Presently, popular ROMs in cavity flows are based on an isentropic assumption, valid only for flows at low or moderate Mach numbers. A new ROM is constructed involving primitive variables of the fully compressible Navier-Stokes （N-S） equations, which is suitable for flows at high Mach numbers. Compared with the direct numerical simulation （DNS） results, the proposed model predicts flow dynamics （e.g., dominant frequency and amplitude） accurately for supersonic cavity flows, and is robust. The comparison between the present transient flow fields and those of the DNS shows that the proposed ROM can capture self-sustained oscillations of a shear layer. In addition, the present model reduction method can be easily extended to other supersonic flows.

Numerical investigation of dynamical behavior of tethered rigid spheres in supersonic flow

The dynamical behavior of two tethered rigid spheres in a supersonic flow is numerically investigated. The tethered lengths and radius ratios of the two spheres are different. The two spheres, which are centroid axially aligned initially, are held stationary first, then released, and subsequently let fly freely in a supersonic flow. The mean qualities of the system and the qualities of the bigger sphere are considered and compared with the situations without the tether. In the separation process, six types of motion caused by the spheres, tether, and fluid interaction are found. The results show that the mean x-velocity of the system changes in a different manner for different radius ratios, and the x-velocity of the bigger sphere is uniformly reduced but through different mechanisms.

NUMERICAL STUDIES ON THE MIXING OF CH_4 AND KEROSENE INJECTED INTO A SUPERSONIC FLOW WITH H_2 PILOT INJECTION

Two-fluid model and divisional computation techniques were used. The multispecies gas fully N-S equations were solved by upwind TVD scheme. Liquid phase equations were solved by NND scheme. The phases-interaction ODE equations were solved by 2nd Runge-Kutta approach. The favorable agreement is obtained between computational results and PLIF experimental results of iodized air injected into a supersonic flow. Then, the numerical studies,were carried out on the mixing of CH, and kerosene injected into a supersonic flow with H-2 pilot injection. The results indicate that the penetration of kerosene approaches maximum when it is injected from the second injector. But the kerosene is less diffused compared with the gas fuels. The free droplet region appears in the flow field. The mixing mechanism of CH4 with H-2 pilot injection is different from that of kerosene. In the staged duct, H-2 can be entrained into both recirculation zones produced by the step mid injectors. But CH, can only be carried into the recirculation between the injectors. Therefore, initiations of H, and CH4 carl occur in those regions. The staged duct is better in enhancing mixing and initiation with H-2 pilot flame.

TIME-PERIODIC ISENTROPIC SUPERSONIC EULER FLOWS IN ONE-DIMENSIONAL DUCTS DRIVING BY PERIODIC BOUNDARY CONDITIONS

We show existence of time-periodic supersonic solutions in a finite interval, after certain start-up time depending on the length of the interval, to the one space-dimensional isentropic compressible Euler equations, subjected to periodic boundary conditions. Both classical solutions and weak entropy solutions, as well as high-frequency limiting behavior are considered. The proofs depend on the theory of Cauchy problems of genuinely nonlinear hyperbolic systems of conservation laws.

Velocity Slip and Interfacial Momentum Transfer in the Transient Section of Supersonic Gas-Droplet Two-Phase Flows

Modelling and simulations are conducted on velocity slip and interfacial momentum transfer for supersonic two-pha.se (gas-droplet) flow in the transient section inside and outside a Laval jet(LJ). The initial velocity slip between gas and droplets causes an interfacial momentum transfer flux as high as (2.0-5.0) x 104 Pa. The relaxation time corresponding to this transient process is in the range of 0.015-0.090ms for the two-phase flow formed inside the LJ and less than 0.5ms outside the LJ. It demonstrates the unique performance of this system for application to fast chemical reactions using electrically active media with a lifetime in the order of 1 ms. Through the simulations of the transient processes with initial Mach number Mg from 2.783 to 4.194 at different axial positions inside the LJ, it is found that Mg has the strongest effect on the process. The momentum flux increases as the Mach number decreases. Due to compression by the shock wave at the end of the LJ, the flow pattern becomes two dimensional and viscous outside the LJ. Laser Doppler velocirneter (LDV) measurements of droplet velocities outside the LJ are in reasonably good agreement with the results of the simulation.

Numerical simulation of gap effect in supersonic flows

The gap effect is a key factor in the design of the heat sealing in super- sonic vehicles subjected to an aerodynamic heat load. Built on S-A turbulence model and Roe discrete format, the aerodynamic environment around a gap on the surface of a supersonic aircraft was simulated by the finite volume method. As the presented results indicate, the gap effect depends not only on the attack angle, but also on the Mach number.

Numerical Analysis of Interaction Between Single-Pulse Laser-Induced Plasma and Bow Shock in a Supersonic Flow

The interaction of laser-induced plasma and bow shock over a blunt body is inves- tigated numerically in an M∞ =6.5 supersonic flow. A ray-tracing method is used for simulating the process of laser focusing. The gas located at the focused zone is ionized and broken down and transformed into plasma. In a supersonic flow the plasma moves downstream and begins to interact with the bow shock when it approaches the surface of the blunt body. The parameters of flowfield and blunt body surface are changed due to the interaction. By analyzing phenomena occurring in the complex unsteady flowfield during the interaction in detail, we can better under- stand the change of pressure on the blunt body surface and the mechanism of drag reduction by laser energy deposition. The results show that the bow shock is changed into an oblique shock due to the interaction of the laser-induced low-density zone with the bow shock, so the wave drag of the blunt body is reduced.

Effects of the jet-to-crossflow momentum ratio on a sonic jet into a supersonic crossflow

Numerical investigation of a transverse sonic jet injected into a supersonic crossflow was carried out using large-eddy simulation for a free-stream Mach number M = 1.6 and a Reynolds number Re = 1.38×10~5 based on the jet diameter.Effects of the jet-to-crossflow momentum ratio on various fundamental mechanisms dictating the intricate flow phenomena,including flow structures, turbulent characters and frequency behaviors,have been studied.The complex flow structures and the relevant flow features are discussed to exhibit the evolution of shock structures,vortical structures and jet shear layers.The strength of the bow shock increases and the sizes of the barrel shock and Mach disk also increase with increasing momentum ratio.Turbulent characters are clarified to be closely related to the flow structures.The jet penetration increases with the increase of the momentum ratio.Moreover,the dominant frequencies of the flow structures are obtained using spectral analysis.The results obtained in this letter provide physical insight in understanding the mechanisms relevant to this complex flow

Thermoelastic stability of closed cylindrical shell in supersonic gas flow

In this paper the problem of linear stability of a closed cylindrical shell under the action of both non-uniform temperature field and supersonic gas flow is considered.The stability conditions for the unperturbed state of the aerothermoelastic system are obtained.It is shown that,by the combined action of the temperature field and the ambient supersonic flow,the process of linear stability can be controlled and the temperature field affects significantly the critical flutter speed.