Contouring the nozzles producing a uniform supersonic flow or a thrust maximum in the presence of a curvilinear sonic line

Within the framework of the ideal, i.e., inviscid and non-heat conducting, gas model we consider the problem of designing the supersonic section of a two-dimensional or axisymmetric nozzle realizing a uniform supersonic flow limitingly similar with a sonic flow when the choked flow involves a curvilinear sonic line. Emphasis is placed on nozzles with abruptly or steeply converging subsonic sections and a strongly curved sonic line formed by the C ⁻-characteristics of the expansion fan with the focus at the lower bend point of the vertical section of the subsonic contour. In the two-dimensional case, the least possible greater-than-unity Mach number M _em at the nozzle exit corresponds to the flow in which the first intersection of the C ⁺-characteristics originated at the closing C ⁻-characteristic of the expansion fan falls on the unknown contour of its supersonic part. For a uniform flow with M _e < M _em the intersection of C ⁺-characteristics beneath the unknown contour make impossible its construction. A part of the contour realizing a uniform flow with M _em > 1 ensures a limitingly rapid flow acceleration and forms the initial region of the supersonic generator of a maximum-thrust nozzle. For this reason, in the case of a curvilinear sonic line the supersonic generators of these nozzles have two, rather than one, bends, which, however, is interesting only for the theory. At least, in the calculated examples the thrusts of the nozzles with one and two bends differ only by a hundredth or even thousandth fractions of per cent.     

On profiling supersonic nozzles to achieve uniform flow in the annular exit section

The problem of profiling supersonic nozzles with a central body in order to achieve uniform supersonic flow with Mach number M _e >1 in the annular exit section is studied. Equal lengths of the profiled inner (central body) and outer nozzle walls can be achieved by displacing downstream the initial bend in the inner contour, i.e., in the central body wall. A displacement of this kind was proposed and tested in [1, 2] for annular nozzles with circular exit sections.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 204–206, March–April, 1995.     

Transonic Shock Problem for the Euler System in a Nozzle

In this paper, we study the well-posedness problem on transonic shocks for steady ideal compressible flows through a two-dimensional slowly varying nozzle with an appropriately given pressure at the exit of the nozzle. This is motivated by the following transonic phenomena in a de Laval nozzle. Given an appropriately large receiver pressure P _r , if the upstream flow remains supersonic behind the throat of the nozzle, then at a certain place in the diverging part of the nozzle, a shock front intervenes and the flow is compressed and slowed down to subsonic speed, and the position and the strength of the shock front are automatically adjusted so that the end pressure at exit becomes P _r , as clearly stated by Courant and Friedrichs [Supersonic flow and shock waves, Interscience Publishers, New York, 1948 (see section 143 and 147)]. The transonic shock front is a free boundary dividing two regions of C ^2,α flow in the nozzle. The full Euler system is hyperbolic upstream where the flow is supersonic, and coupled hyperbolic-elliptic in the downstream region Ω₊ of the nozzle where the flow is subsonic. Based on Bernoulli’s law, we can reformulate the problem by decomposing the 3 × 3 Euler system into a weakly coupled second order elliptic equation for the density ρ with mixed boundary conditions, a 2 × 2 first order system on u ₂ with a value given at a point, and an algebraic equation on (ρ, u ₁, u ₂) along a streamline. In terms of this reformulation, we can show the uniqueness of such a transonic shock solution if it exists and the shock front goes through a fixed point. Furthermore, we prove that there is no such transonic shock solution for a class of nozzles with some large pressure given at the exit. This research was supported in part by the Zheng Ge Ru Foundation when Yin Huicheng was visiting The Institute of Mathematical Sciences, The Chinese University of Hong Kong. Xin is supported in part by Hong Kong RGC Earmarked Research Grants CUHK-4028/04P, CUHK-4040/06P, and Central Allocation Grant CA05-06.SC01. Yin is supported in part by NNSF of China and Doctoral Program of NEM of China.     

Transonic Shocks for the Full Compressible Euler System in a General Two-Dimensional De Laval Nozzle

In this paper, we study the transonic shock problem for the full compressible Euler system in a general two-dimensional de Laval nozzle as proposed in Courant and Friedrichs (Supersonic flow and shock waves, Interscience, New York, 1948): given the appropriately large exit pressure p _e(x), if the upstream flow is still supersonic behind the throat of the nozzle, then at a certain place in the diverging part of the nozzle, a shock front intervenes and the gas is compressed and slowed down to subsonic speed so that the position and the strength of the shock front are automatically adjusted such that the end pressure at the exit becomes p _e(x). We solve this problem completely for a general class of de Laval nozzles whose divergent parts are small and arbitrary perturbations of divergent angular domains for the full steady compressible Euler system. The problem can be reduced to solve a nonlinear free boundary value problem for a mixed hyperbolic–elliptic system. One of the key ingredients in the analysis is to solve a nonlinear free boundary value problem in a weighted Hölder space with low regularities for a second order quasilinear elliptic equation with a free parameter (the position of the shock curve at one wall of the nozzle) and non-local terms involving the trace on the shock of the first order derivatives of the unknown function.     

Jet flow field during screech

Measurements were made of the flow field structure and the near field parameters of a jet exhausting from a sonic nozzle with a 1.27 cm exit diameter. Compressed air was used for obtaining stagnation pressures up to 5 atmospheres. The jet exhausted vertically from a settling chamber into an acoustically insulated room and through an insulated duct out through the roof. Measurements were made with several different reflecting surfaces at the nozzle exit as well as an insulating surface. Schlieren pictures at 500,000 frames/s were taken. Overall sound pressure level, impact pressure level downstream, and sound frequency analyzer measurements were made.It was found that with a reflecting surface there was a radial oscillation of the jet which had the same frequency as the dominant sound (screech) frequency emitted by the jet. No axial motion of the inviscid part of the flow structure was detected. The insulated surface at the nozzle exit appeared to shift the dominant frequencies of the sound generated into the region above the audible (>16 KHz). A reflecting surface yielded pure tones (screech) with one or two harmonics. The dominant (screech) frequency decreased as the stagnation pressure increased. The screech frequency was found to be approximately inversely proportional to the length of the first shock cell.Nomenclature C ₀ speed of sound in ambient gas - D diameter of nozzle exit - f frequency of pure tone (screech frequency) - L ₁ length of first cell, distance between nozzle exit plane and intersection of shock with shear layer - M Mach number based on isentropic expansion to ambient pressure - P ₀ stagnation chamber pressure - P _a ambient pressure - P _i impact pressure - R _LB distance from nozzle centerline to left boundary of jet - R _RB distance from nozzle centerline to right boundary of jet - t time - period of screech, 1/f - X _E axial distance of eddy from nozzle exit plane - X _I axial distance of third cell shock intersection from nozzle exit plane - Y _I transverse distance of third cell shock intersection from nozzle centerline     

UV-mediated coalescence and mixing of inkjet printed drops

In the present study, the characteristics of supersonic rectangular microjets are investigated experimentally using molecular tagging velocimetry. The jets are discharged from a convergent–divergent rectangular nozzle whose exit height is 500 μm. The jet Mach number is set to 2.0 for all tested jets, and the Reynolds number Re is altered from 154 to 5,560 by changing the stagnation pressure. The experimental results reveal that jet velocity decays principally due to abrupt jet spreading caused by jet instability for relatively high Reynolds numbers (Re > ~450). The results also reveal that the jet rapidly decelerates to a subsonic speed near the nozzle exit for a low Reynolds number (Re = 154), although the jet does not spread abruptly; i.e., a transition in velocity decay processes occurs as the Reynolds number decreases. A supersonic core length is estimated from the streamwise distribution of the centerline velocity, and the length is then normalized by the nozzle exit height and plotted against the Reynolds number. As a result, it is found that the normalized supersonic core length attains a maximum value at a certain Reynolds number near which the transition in the velocity decay process occurs.     

Experimental Investigation of Nozzle Exit Reflector Effect on Supersonic Jet

The present study describes an experimental work to investigate the effect of a nozzle exit reflector on a supersonic jet that is discharged from a convergent–divergent nozzle with a design Mach number of 2.0. An annular reflector is installed at the nozzle exit and its diameter is varied. A high-quality spark schlieren optical system is used to visualize detailed jet structures with and without the reflector. Impact pressure measurement using a pitot probe is also carried out to quantify the reflector’s effect on the supersonic jet which is in the range from an over-expanded to a moderately under-expanded state. The results obtained show that for over-expanded jets, the reflector substantially increases the jet spreading rate and reduces the supersonic length of the jet, compared with moderately under-expanded jets. The reflector’s effect appears more significant in imperfectly expanded jets that have strong shock cell structures, but is negligible in correctly expanded jet.     

Aerodynamic experimentation with ducted models as applied to hypersonic air-breathing vehicles

 A methodology of experimentation in high supersonic wind tunnels for studying aerodynamic characteristics of hypersonic flying vehicles powered by air-breathing engines is discussed. Investigations of such total aerodynamic forces as drag, lift and pitching moment at testing the models are implicit when the air flow through the model ducts is accomplished so that to provide the simulation of the external flow around the airplane and flow over the inlets, but the operating engines and, hence, the exhaust jets are not modeled. The methods used for testing such models are based on the measurement of duct stream parameters alongside with the balance measurement of aerodynamic forces acting on the models. In the tests, aerometric tools are used such as narrow metering nozzles (plugs), pitot and static pressure probes, stagnation temperature probes and pressure orifices in walls of the model duct. The aerometric data serve to determine the flow rate and momentum of the stream at the duct exit. The internal non-simulated forces of the model ducts are also determined using the conservation equations for energy, mass flow and momentum, and these forces are eliminated from the aerodynamic test results. The techniques of the said model testing have been well developed as applied to supersonic aircraft, however their application for hypersonic vehicles whose models are tested at high supersonic speeds, Mach number M _∞>4, implies some specific features. In the present paper, the results of experimental and theoretical study of these features are discussed. Some experimental data on aerodynamics of hypersonic aircraft models received in methodological tests are also presented. The tunnel experiments have been carried out in the Mach number range M _∞=2–6. Received: 25 July 1996/ Accepted: 14 December 1998     

Existence of supersonic zones in three-dimensional separation flows

Up till now the region of three-dimensional separation flows which occur with supersonic flow past obstacles has received insufficient study. Supersonic flow with a Mach number of 2.5 past a cylinder mounted on a plate was studied in [1]. A local zone with supersonic velocities was found in the reverse subsonic flow region ahead of the cylinder. Its presence is explained by the three-dimensional nature of the flow. Similar supersonic zones are not observed in the case of supersonic flow over plane and axisymmetric steps.The present paper presents the results of experimental studies whose objective was refinement of the flow pattern ahead of a cylinder on a plate and the study of the local supersonic zones.The experiments were performed in a supersonic wind tunnel with a freestream Mach number M₁=3.11. The 24-mm-diameter cylinder with pressure taps along the generating line was mounted perpendicular to the surface of a sharpened plate. The distance from the plate leading edge to the cylinder axis wasl ₀=140 mm. The plate was pressure tapped along the flow symmetry axis. The Reynolds number was Rl ₀=u₀ l ₀/v ₁, Rl ₀=1.87.10⁷, where u₁ andv ₁ are the freestream velocity and the kinematic viscosity, respectively. The pressures were measured using a Pilot probe with internal and external diameters of 0.15 and 0.9 mm, respectively.The probe was displaced in the flow symmetry plane at a distance of 1.6 mm from the plate surface and at a distance of 1.1 mm along the leading generator of the cylinder. The flow on the surface of the plate and cylinder was studied with the aid of a visualization composition and the flow past the model was photographed with a schlieren instrument. Typical patterns of the visualization composition distribution and the pressure distribution curves over the plate surface, and also photographs of the flow past the model, are shown in [1].     

Experimental investigation of a three-dimensional viscous underexpanded jet

The three-dimensional interaction of jet issuing from two- and four-nozzle systems into ambient space or an outer flow has been investigated experimentally. The range of the important parameters include the following: pressure imbalance n=P_a=/P=10–1.5·10², Mach number at the nozzle exit M_a=3.15, Mach number of the outer flow M=0, 3.1, and 6, the flow is turbulent in the mixing layer (P_a and P are the static pressures at the nozzle exit and in the outer flow). It is shown that the interaction of the jets broadens a multinozzle jet considerably in the plane of interaction, which is a plane of symmetry and which passes through the axis of the system between neighboring nozzles. The cross-sectional shape of a four-nozzle jet is cross-like over the entire length of the initial segment of the jet. The width of the mixing layer in the plane of interaction is considerably larger than in the central plane, which passes through the axis of opposed nozzles.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 21–26, September–October, 1974.     

Experimental investigations of steady and oscillating annular liquid curtains with different surface tensions

Experiments were performed on laminar, vertical, annular, liquid curtains to study the dynamics of steady curtains, and the onset and frequency of oscillating curtains. The experiments were conducted to observe the effects of inertia and pressure on liquid curtains with different surface tensions. For steady curtains, convergence lengths were measured as functions of Froude number and pressure differential for three different surface tensions. The factors causing the onset of oscillations in a pressurized curtain were observed and the frequency of the internal pressure fluctuations were measured for various Froude numbers and two surface tensions.List of symbols b local thickness of curtain sheet - b ₀ initial thickness of curtain or nozzle gap thickness (0.5 mm) - C _P pressure coefficient - Fr Froude number (V ₀ ² /g R ₀) - g gravitational acceleration - g gravitational acceleration - L convergence length of curtain - L ^* dimensionless convergence length (L/R ₀) - N _c convergence number (g ² R ₀ ² b ₀ /2v ₀ ² ) - P _e pressure outside the curtain (ambient) - P _i pressure inside the curtain - P pressure differential (P _i– P _e) - P _cr pressure differential at which curtain begins to oscillate - R local radius of curvature in the horizontal plane - R ₀ initial curtain radius or radius of nozzle exit (50 mm) - r _v local radius of curvature in the vertical plane - V local liquid velocity - V ₀ initial liquid velocity - V ^* dimensionless local liquid velocity (V/V ₀) - z axial distance from the nozzle - z ^* dimensionless axial distance from the nozzel (z/R ₀) - s differential length of curtain - differential angle in the horizontal plane - angle between the direction of the surface tension force in the vertical plane and the direction of r _v - deangle between the direction of the surface tension force in the horizontal plane and the direction of R - angle between r _vand R in the vertical plane - ₀ nozzle exit angle (zero degrees) - surface tension of liquid - liquid density (1.0 gm/cm³)     

Heat transfer of gas-particle flow in a supersonic convergent-divergent nozzle

Summary Heat flux, wall heat transfer coefficients, and wall pressures are determined for high velocity flow of gas-solid mixtures in a converging-diverging nozzle. Flow separation accompanied with oblique shock formation occurs in the diverging section of the nozzle. The shock strength is reduced upon the addition of solid particles. The wall pressure in the convergent section of the nozzle appears unaffected by the presence of solid particles. In the divergent section, however, the wall pressure is slightly lowered. At the maximum ratio of solid to air flow used in the experiments (3.7) increases in the heat transfer rate of up to 20 and 50 percent are obtained in the convergent and separated (divergent) regions of the nozzle, respectively. Slightly larger increases in the wall heat transfer coefficients are also obtained. It is concluded that the wall heat flux and heat transfer coefficients are influenced strongly by the presence of disturbances upstream of the nozzle inlet.Nomenclature W _a air flow rate - W _s solids flow rate - x axial distance from nozzle entrance - L axial length of nozzle - specific heat ratio of fluid - A _e exit cross section of flow - A _* throat cross section of flow - P ₀ inlet pressure - P _s wall separation pressure - P _a ambient exhaust pressure - shock wave angle - shock wave deflection angle - M ₁ Mach number upstream of shock wave - Mach number normal to shock wave - q heat flux - k _f thermal conductivity of fluid - T _wi inside wall temperature - T _wo outside wall temperature - T _ad adiabatic wall temperature - h wall heat transfer coefficient - C nozzle constant - A local cross section of flow - c _p specific heat of fluid - Pr Prandtl number - viscosity of fluid - r _c throat radius of curvature - factor accounting for variation of and Units absolute temperature °R(ankine) °F+459.7 - conductivity 1 BTU (hr ft °F)^–1 4.137×10^–3 cal (s cm °C)^–1 - specific heat 1 BTU (1b °F)^–1 1 cal (g °C)^–1 - absolute pressure 1 psia 0.0680 atm Supported in part by aid provided by the UCLA Space Science Center (Grant NsG 236-62 Libby).Listed for readers not familiar with the units adopted in this paper (editor).     

Effect of primary jet geometry on ejector performance: A cold-flow investigation

The following cold-flow study examines the interaction of the diffracted shock wave pattern and the resulting vortex loop emitted from a shock tube of various geometries, with an ejector having a round bell-shaped inlet. The focus of the study is to examine the performance of the ejector when using different jet geometries (primary flow) to entrain secondary flow through the ejector. These include two circular nozzles with internal diameters of 15 mm and 30 mm, two elliptical nozzles with minor to major axis ratios of a/b = 0.4 and 0.6 with b = 30 mm, a square nozzle with side lengths of 30 mm, and two exotic nozzles resembling a pair of lips with axis ratios of a/b = 0.2 and 0.5 with b = 30 mm. Shock tube driver pressures of P₄ = 4, 8, and 12 bar were studied, with the pressure of the shock tube driven section P₁ being atmospheric. High-speed schlieren photography using the Shimadzu Hypervision camera along with detailed pressure measurements along the ejector and the impulse created by the ejector were conducted.     

Simplified Navier-Stokes Equations for Internal Mixed Flows and a Numerical Method of Solution

Simplified Navier-Stokes equations, of the elliptic and hyperbolic type in the subsonic and supersonic flow regions, respectively, are derived for viscous flows in channels and nozzles with curved walls whose local radii of longitudinal curvature are comparable with the transverse channel dimensions. A new numerical method is developed for the system of equations obtained. This method is of the evolution type along the longitudinal coordinate and includes global iterations of the streamline direction field and the longitudinal pressure gradient field. The effectiveness of the method is illustrated with reference to the solution of the direct Laval nozzle problem for an air flow at Reynolds numbers Re10⁴ and 10⁶ in conical nozzles with throat curvatures K _w=1.0 and 1.6 (K _w is the curvature divided by the inverse radius of the nozzle throat). Two iterations are sufficient to calculate the nozzle flow rate and power correct to 0.01%.     

Gas dynamic ozone generator

The formation of ozone when partially dissociated oxygen flows out of a supersonic nozzle has been investigated experimentally and theoretically. The supersonic flow of a chemically reacting gas mixture containing excess O atoms is calculated in the one-dimensional approximation for a class of plane wedge-shaped nozzles. It is shown that for initial gas pressures ahead of the nozzle inlet of about 10 atm and a temperatureT ₀=1000 K in nozzles with a total vertex angle of 30°C and a throat dimensionh.=1 mm it is possible to obtain an ozone concentration of about 1%, which is comparable with ordinary ozonizers, while the output of the device is two to three orders greater. Experiments on a shock tube fitted with a nozzle to measure the absorption of UV radiation by oxygen recombining in the nozzle under highly nonoptimal conditions revealed the presence in the flow of ozone molecules formed as a result of O+O₂ recombination.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 139–148, November–December, 1994.     

A hybrid pressure–density‐based algorithm for the Euler equations at all Mach number regimes

In the present work, we propose a reformulation of the fluxes and interpolation calculations in the PISO method, a well‐known pressure‐correction solver. This new reformulation introduces the AUSM⁺ ? up flux definition as a replacement for the standard Rhie and Chow method of obtaining fluxes and central interpolation of pressure at the control volume faces. This algorithm tries to compatibilize the good efficiency of a pressure based method for low Mach number applications with the advantages of AUSM⁺ ? up at high Mach number flows. The algorithm is carefully validated using exact solutions. Results for subsonic, transonic and supersonic axisymmetric flows in a nozzle are presented and compared with exact analytical solutions. Further, we also present and discuss subsonic, transonic and supersonic results for the well known bump test‐case. The code is also benchmarked against a very tough test‐case for the supersonic and hypersonic flow over a cylinder. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis of the nozzle exit flow parameter effect on the turbulence characteristics and the noise level in jets issuing from nozzles of different types

The flow, turbulence, and noise parameters of cold and hot jets flowing out of nozzles of different types at subsonic and supersonic velocities are calculated using the high-resolution RANS/ILES method. The effect of the Mach number and the temperature at the nozzle exit on the flow features, the turbulent fluctuations of the velocity, the static pressure, and the temperature, together with the overall noise level is analyzed for all the jets considered. The accuracy of the calculations is confirmed by means of comparing with the available experimental data concerning certain parameters.     

Experimental studies of a steam jet refrigeration cycle: Effect of the primary nozzle geometries to system performance

This paper describes an experimental investigation of a steam jet refrigeration. A 1 kW cooling capacity experimental refrigerator was constructed and tested. The system was tested with various operating temperatures and various primary nozzles. The boiler saturation temperature ranked from 110 to 150 °C. The evaporator temperature was fixed at 7.5 °C. Eight primary nozzles with difference geometries were used. Six nozzles have throat diameters ranked from 1.4 to 2.6 mm with exit Mach number of 4.0. Two remained nozzles have equal throat diameter of 1.4 mm but difference exit Mach number, 3.0 and 5.5. The experimental results show that the geometry of the primary nozzle has strong effects to the ejector performance and therefore the system COP.     

Laminar-turbulent transition in the boundary layer on cones in a hypersonic flow at high Reynolds numbers per meter

The laminar-turbulent transition is experimentally studied in boundary-layer flows on cones with a rectangular axisymmetric step in the base part of the cone and without the step. The experiments are performed in an A-1 two-step piston-driven gas-dynamic facility with adiabatic compression of the working gas with Mach numbers at the nozzle exit M _∞ = 12–14 and pressures in the settling chamber P₀ = 60–600 MPa. These values of parameters allow obtaining Reynolds numbers per meter near the cone surface equal to Re _1e = (53–200) · 10⁶ m ⁻¹. The transition occurs at Reynolds numbers Re _tr = (2.3–5.7) · 10⁶. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 76–83, May–June, 2007.