An exact solutin for incompressible flow through a two-dimensional Laval nozzle

A careful examination of the variation of the velocity along the centerline and the contour of a Laval nozzle in the physical plane shows that either the upper or the lower half of the Laval nozzle assumes the same form of a slitted thick airfoil with tandem trailing edges. These two airfoils lie on different Riemann sheets in the hodograph plane. The interior of the airfoil is then mapped onto an infinite strip in the complex potential plane. Making use of these results, we obtained an exact solution for the incompressible potential flow through a two-dimensional Laval nozzle. The solution is applicable for nozzles with any given contraction ratio mexpansion rations, and throat wall radius R^*. As examples of the method, various nozzle contours, the velocity distribution of the flow, and the locations of the fluid particles at different time intervals are presented.     

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%.     

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.     

Symbolic computation of flow in a rotating pipe

A perturbation solution of the fully developed flow through a pipe of circular cross-section, which rotates uniformly around an axis oriented perpendicularly to its own, is considered. The perturbation parameter is given by R = 2Ωa²/ν in terms of the angular velocity Ω, the pipe radius a and the kinematic viscosity ν of the fluid. The two coupled non-linear equations for the axial velocity ω and the streamfunction ? of the transverse (secondary) flow lead to an infinite system of linear equations. This system allows first the computation of a given order ?_n, n ? 1, of the perturbation expansion ? = ∑ Rⁿ?_n in terms of ω_n-1, the (n-1)-th order of the expansion ω = ∑ Rⁿω_n, and of the lower orders ?₁,…,?_{n ? 1}. Then it permits the computation of ω_n from ω₀,…,ω_{n ? 1} and ?₁,…,?;_n. The computation starts from the Hagen–Poiseuille flow ω₀, i.e. the perturbation is around this flow. The computations are performed analytically by computer, with the REDUCE and MAPLE systems. The essential elements for this are the appropriate co-ordinates: in the complex co-ordinates chosen the two-dimensional harmonic (Laplace, Δ) and biharmonic (Δ²) operators are ideally suited for (symbolic) quadratures. Symmetry considerations as well as analysis of the equations for ω_n, ?_n and of the boundary conditions lead to general (polynomial) formulae for these functions, with coeffcients to be determined. Their determination, order by order, implies, in complex co-ordinates, only (symbolic) differentiation and quadratures. The coefficients themselves are polynomials in the Reynolds number c of the (unperturbed) Hagen–Poiseuille flow. They are tabulated in the paper for the orders n ? 6 of the perturbation expansion.     

Flow of a dispersed phase in the Laval nozzle and in the test section of a two-phase hypersonic shock tunnel

An unsteady gas-particle flow in a hypersonic shock tunnel is studied numerically. The study is performed in the period from the instant when the diaphragm between the high-pressure and low-pressure chambers is opened until the end of the transition to a quasi-steady flow in the test section. The dispersed phase concentration is extremely low, and the collisions between the particles and their effect on the carrier gas flow are ignored. The particle size is varied. The time evolution of the particle concentration in the test section is obtained. Patterns of the quasi-steady flow of the dispersed phase in the throat of the Laval nozzle and the flow around a model (sphere) are presented. Particle concentration and particle velocity lag profiles at the test-section entrance are obtained. The particle-phase flow structure and the time needed for it to reach a quasi-steady regime are found to depend substantially on the particle size. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 102–113, September–October, 2008.     

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.     

Oscillation modes of laval nozzle flow

Experimental investigations of Laval nozzle flow show for relatively low supply to exit pressure ratios, which correspond to shock wave positions close to the nozzle throat, three different, oscillatory instabilities.

1. Shock pattern oscillations where the root of a λ-like shock front remains nearly in constant position, but where the proportion between the normal part and the oblique part of the shock changes periodically.
2. Shock wave and separation bubble oscillations where the motion of the shock wave is accompanied by displacements of the separation bubble.
3. Flow rate oscillations where the shock waves leave periodically through the nozzle throat in upstream direction.

     

Numerical study of the spontaneous nucleation of self-rotational moist gas in a converging–diverging nozzle

Spontaneous nucleation is the primary way of droplet formation in the supersonic gas separation technology, and the converging–diverging nozzle is the condensation and separation unit of supersonic gas separation devices. A three-dimensional geometrical model for the generation of self-rotational transonic gas flow is set up, based on which, the spontaneous nucleation of self-rotational transonic moist gas in the converging–diverging nozzle is carried out using an Eulerian multi-fluid model. The simulated results of the main flow and nucleation parameters indicate that the spontaneous nucleation can occur in the diverging part of the nozzle. However, different from the nucleation flow without self-rotation, the distributions of these parameters are unsymmetrical about the nozzle axis due to the irregular flow form caused by the self-rotation of gas flow. The nucleation region is located on the position where gas flows with intense rotation and the self-rotation impacts much on the nucleation process. Stronger rotation delays the onset of spontaneous nucleation and yields lower nucleation rate and narrow nucleation region. In addition, influences of other factors such as inlet total pressure p ₀, inlet total temperature T ₀, the nozzle-expanding ratio ? and the inlet relative humidity ф ₀ on the nucleation of self-rotational moist gas flow in the nozzle are also discussed.     

Shock and wave dynamics in cavitating compressible liquid flows in injection nozzles

Due to the exceptional high inlet pressures up to 2,000 bar flow dynamics and efficiency of modern injection systems are controlled by high frequency wave dynamics of the compressible liquid flow. Corresponding to alternating shock and expansion waves the liquid fluid evaporates and recondenses instantaneously. Here we present CFD simulations of the time accurate evolution of cavitating flows in 2-D plane and in six-hole injection nozzles with focus on the wave dynamics just after initialisation of the flow and within the time scale Δt ≤ 10^?4 s of pilot and multi-point injection. Due to shock reflections at the bottom of the sack hole the instantaneous maximum pressure increases more than three times higher as compared with the prescribed pressure at the nozzle inlet. For instance, in case of an inlet pressure of 600 bar the maximum pressure in the sack and therefore ahead of the nozzle bore holes reaches about 2,100 bar. It is quite reasonable that this amplification of the pressure affects the evolution of the convective flow and therefore the mass flow through the nozzle bore holes.     

Investigation of the effect of nozzle shape on supersonic/hypersonic impactors designed for size discrimination of nanoparticles

In this study the flow field and the nanoparticle collection efficiency of supersonic/hypersonic impactors with different nozzle shapes were studied using a computational modeling approach. The aim of this study was to develop a nozzle design for supersonic/hypersonic impactors with the smallest possible cut-off size d₅₀ and rather sharp collection efficiency curves. The simulation results show that the changes in the angle and width of a converging nozzle do not alter the cut-off size of the impactor; however, using a conical Laval nozzle with an L/D_n ratio less than or equal to 2 reduced d₅₀. The effect of using a cap as a focuser in the nozzle of a supersonic/hypersonic impactor was also investigated. The results show that adding a cap in front of the nozzle had a noticeable effect on decreasing the cut-off size of the impactor. Both flat disks and conical caps were examined, and it was observed that the nozzle with the conical cap had a lower cut-off size.     

Numerical solution of the direct-problem of mixed nozzle flow

A large part of the known results of Laval nozzle theory relates to the inverse problem, in which the velocity distribution on some line (usually the axis of symmetry) is given rather than the nozzle contour. Many important properties of transonic flows have been disclosed as a result of numerous studies, whose basic results were presented together with an extensive bibliography in Ryzhov's monograph [1]. The solution of the inverse problem has recently been used not only to analyze the qualitative characteristics but also to construct nozzles with rather marked variation of the slope of the generator, which are of practical interest. In this connection we note the work of Pirumov [2] and also the studies of Hopkins and Hill [3, 4]. The latter authors, in addition to the classical Laval nozzle, studied several nozzle schemes with a centerbody. Pirumov used a specially developed numerical method for the solution of the inverse problem (we note that in the subsonic part of the nozzle the corresponding Cauchy problem is incorrect), while Hopkins and Hill used a series expansion which was preceded by a change of variables.There are considerably fewer studies devoted to the solution of the direct problem of mixed nozzle flow. Numerical methods have been used by Alikhashkin, Favorskii, and Chushkin [5], Favorskii [6], and Danilov [7], with the method of integral relations being used in the first two studies. Finally, there has recently been extensive development of the method of expansion in powers of ^1/2, where is the ratio of the radius (or half-width of the nozzle to the radius of curvature of the wall, calculated at the throat section. Such expansions have been used by Hall [8] and Kliegel and Quan [9] to study flow in classical Laval nozzles, and by Moore [10] and Moore and Hall [11] to study flow in nozzles with a centerbody. We note that the ^1/2-expansion method is suitable only in those cases in which the wall radii of curvature are large.In the following the asymptotic method is used to solve the direct problem of mixed flow in nozzles. This reduces the very complex boundary value problem for an elliptic-hyperbolic system of equations with two unknown variables to the Cauchy problem (more precisely, to a mixed problem with initial conditions in a bounded two-dimensional region and boundary conditions which are independent of the third variable) for a hyperbolic system with three unknown variables. The integration of the equations describing the two-dimensional (plane of axisymmetric) nonsteady flow was accomplished with the aid of the Godunov-Zabrodin-Prokopov difference scheme [12]. Several types of nozzles with centerbody are calculated as well as the classical Laval nozzle. The contours of the subsonic parts of the nozzles were either closed (finite combustion chamber) or open (nozzle joins an infinite cylindrical tube). In the first case the flow is provided by three-dimensional mass and energy sources which are introduced at some fixed part of the combustion chamber. In the second case there are no mass and energy sources, but a boundary condition is established at a plane perpendicular to the nozzle axis and located at a finite distance from the throat section, and this condition becomes the flow uniformity condition as this plane moves away to infinity.The authors wish to thank I. Yu. Brailovskii for valuable advice in the selection of the difference scheme, U. G. Pirumov for the kind offer of the results of his calculations, and A. M. Konkina and L. P. Frolova for assistance in the calculations.     

The effects of rate of expansion and injection of water droplets on the entropy generation of nucleating steam flow in a Laval nozzle

The entropy generation due to irreversible heat transfer between vapor and liquid phases in a nucleating steam flow in a Laval nozzle is studied. To calculate the entropy generation due to self-condensation in transonic steam flow, a thermodynamic model is presented. The calculations of nucleating steam flow and the predictions of entropy generation rely on one-dimensional two-phase model. This model shows that the most of the thermodynamic losses take place during the nucleation phenomena. The effect of rate of expansion on the exergy losses is considered by decreasing the divergent angle of nozzle. Also micro-sized pure water droplets is injected theoretically to supercooled steam right after the nozzle throat at the onset of divergent section and the effects of injected droplets on thermodynamic losses and nucleation phenomena are investigated. The results indicate that decreasing the divergent angle and also injection of droplets diminishes the pressure rise in transonic steam flow and decreases the thermal entropy generation due to nucleation.     

The effect of magnetic field on two-phase liquid metal flow

In this paper, the basic equations of two-phase liquid metal flow in a magnetic field are derived, and specifically, two-phase liquid metal MHD flow in a rectangular channel is studied, and the expressions of velocity distribution of liquid and gas phases and the ratioK ₀ of the pressure drop in two-phase MHD flow to that in single-phase are derived. Results of calculation show that the ratioK ₀ is smaller than unity and decreases with increasing void fraction and Hartmann number because the effective electrical conductivity in the two-phase case decreases. The Project is supported by the National Natural Science Foundation of China.     

Coefficient of discharge and spray cone angle of a pressure nozzle with combined axial and tangential entry of power-law fluids

Flows of incompressible, time-independent purely viscous power-law fluids through pressure nozzle with combined axial and tangential entry are analysed. Theoretical predictions of coefficient of discharge and spray cone angle are made through an approximate analytical solution of hydrodynamics of flow inside the nozzle. In the converging section of the nozzle, the boundary layer equations have been derived with modified order approximation [O(δ/R)≈1, O(δ ²/R ²)≪1] of Navier-Stokes equations for a better accuracy. Smoother attainment of the free-stream condition at the edge of the boundary layer is ensured by requiring the appropriate shear rate terms, compatible with the above order analysis, to be zero. The pertinent independent input parameters which govern the flow field are the generalized Reynolds number at inlet to the nozzle based on the tangential velocity of injection , the ratio of the axial-to-tangential velocity at the inlet to the nozzle V _R , the flow behaviour index of the fluid n, the length-to-diameter ratio of the swirl chamber L ₁/D ₁, the spin chamber angle 2α and the orifice-to-swirl-chamber-diameter ratio D ₂/D ₁. Experiments reported in the paper corroborate the qualitative trends of analytical results.     

Swirling flows in circular-to-rectangular transition ducts

Miniature axisymmetric supersonic nozzles were produced with exit Mach numbers ranging from 1.0 to 2.8 by forming Pyrex^® capillary tubing of 0.6 and 1.2 mm inside diameter into converging-diverging channels. The nozzle contours were measured and were found to compare favorably to ideal solutions given by the axisymmetric method of characteristics. In addition, the surfaces of these nozzles were quite smooth, providing featureless flows at perfect expansion. Schlieren visualization and pitot pressure measurements of the resulting microjets were compared to the literature available for jets produced by larger-scale nozzles. A postponed transition to turbulence is noted in these microjets due to their low Reynolds number. The pitot pressure on centerline is nearly uniform at perfect expansion over core lengths up to 12 nozzle exit diameters. Supersonic microjet nozzles thus provide a more effective small-scale high-pressure gas delivery device than do sonic nozzles of comparable scale at equivalent mass flow rates. Supersonic microjets may therefore have several industrial applications.List of symbols ^* boundary layer displacement thickness, mm - d diameter of nozzle exit, mm - L length of nozzle diverging section, mm - L _c inviscid core length, mm - L _s supersonic region length, mm - M _c convective Mach number - M _e exit Mach number - P _b backpressure at nozzle exit, (equal to ambient pressure in this experiment) - P _e exit pressure of the supersonic jet - P _be exit pressure ratio (P _b /P _e ) - P _p impingement pressure (pitot pressure) - P ₀ stagnation pressure supplied to nozzle - P _n overall pressure ratio (P ₀/P _b ,) - r radial dimension (cylindrical coordinate system), mm - r ₀ radius of throat, mm - Re _d Reynolds number, based on nozzle exit diameter - V _e exit velocity, m/s - x axial dimension (cylindrical coordinate system), mm This research was sponsored by National Science Foundation Grant DMI 9400119, as part of a study of the assist-gas dynamics of laser cutting.     

Calculating the Thrust Characteristics of Nozzles from Their Outlet Cross-Section Parameters

Within the framework of boundary-layer theory, simple formulas are presented for finding the viscous loss from the gas parameters in the nozzle outlet cross-section. The analysis is performed for ordinary Laval nozzles and spike nozzles. It is found that for nozzles with large expansion ratios the viscous loss is almost independent of the outlet cross-section parameters and is determined only by the parameter values on the nozzle contour. The effect of the longitudinal nozzle curvature on this loss is investigated. It is shown that the viscous losses calculated from the nozzle outlet parameters and by integrating along the nozzle contour with account for the longitudinal curvature fully coincide.     

Experimental and numerical analysis of the structure of pseudo-shock systems in laval nozzles with parallel side walls

Detailed numerical and experimental investigations of pseudo-shock systems in a Laval nozzle with parallel side walls are carried out. The location of the pseudo-shock system is defined in this system of two choked Laval nozzles by the ratio of the critical cross sections A₂^*/A₁^*{{A}_{2}^*/{A}_{1}^*} , the stagnation pressure loss across the shock system and viscous losses. The wall pressure distributions and high-speed schlieren videos recorded in the experiments are compared to the results of a steady and an unsteady numerical simulation. For the steady case, good agreement is found between the calculated and measured shock structure and pressure distribution along the primary nozzle wall, except for a remaining slight deviation in the shock position. For the unsteady case, in which asymmetric shock configurations are observed, deviations of the results with respect to the stochastic wall attachment of the shock system are given which indicate the necessity of further investigations on that topic.     

Calculation of the swirling flow of an ideal gas in a laval nozzle

The numerical solution of the problem of the motion of a swirling flow of an ideal gas in a Laval nozzle in axisymmetric formulation is obtained by the method of stabilization. As a result, a number of effects appear that are essentially not one-dimensional, in particular, the drawing-in of the sonic line into the nozzle, an effect that leads to a decrease in the nozzle's expansion coefficient. The dependence of this coefficient on the intensity of the swirling is obtained. A number of problems connected with the control of the expansion of a gas through a Laval nozzle and with variation of the thrust of a nozzle can be solved successfully in cases where a rotary motion is imparted to the flow of gas exhausted from the nozzle. Investigation of such a swirling flow in [1, 2] and a number of other papers are based on a one-dimensional model of gas flow, which makes it possible in principle to obtain integrated characteristics of the flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 72–76, September–October, 1971.     

3D numerical simulation of compressible swirling flow induced by means of tangential inlets

The objective of the present work is to predict compressible swirl flow in the nozzle of air‐jet spinning using the realizable k–ε turbulence model and discuss the effect of the nozzle pressure. The periodic change of flow patterns can be observed. The recirculation zone near the wall of the injectors upstream increases in size and moves gradually upstream, whereas the vortex breakdown in the injector downstream shifts slowly towards the nozzle outlet during the whole period. A low axial velocity in the core region moves gradually away from the centerline, and the magnitude of the center reverse flow and the area occupied by it increase with axial distance due to the vortex breakdown. From the tangential velocity profile, there is a very small free‐vortex zone. With increasing nozzle pressure, the velocity increases and the location of vortex breakdown is moved slightly downward. However, the increase in the velocity tends to decline at nozzle pressure up to a high level. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical and experimental investigation of a laval nozzle with a cylindrical throat

The results are reported of experimental and numerical investigation of mixed flow and of the parameters of heat transfer in the transonic region of an axisymmetric Laval nozzle whose throat is formed by a cylindrical surface, i.e., the nozzle contour near the minimum cross section contains two bends.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 189–192, September–October, 1984.