Investigation of a Gas Flow with Solid Particles in a Supersonic Nozzle

A two-phase flow with high Reynolds numbers in the subsonic, transonic, and supersonic parts of the nozzle is considered within the framework of the Prandtl model, i.e., the flow is divided into an inviscid core and a thin boundary layer. Mutual influence of the gas and solid particles is taken into account. The Euler equations are solved for the gas in the flow core, and the boundary-layer equations are used in the near-wall region. The particle motion in the inviscid region is described by the Lagrangian approach, and trajectories and temperatures of particle packets are tracked. The behavior of particles in the boundary layer is described by the Euler equations for volume-averaged parameters of particles. The computed particle-velocity distributions are compared with experiments in a plane nozzle. It is noted that particles inserted in the subsonic part of the nozzle are focused at the nozzle centerline, which leads to substantial flow deceleration in the supersonic part of the nozzle. The effect of various boundary conditions for the flow of particles in the inviscid region is considered. For an axisymmetric nozzle, the influence of the contour of the subsonic part of the nozzle, the loading ratio, and the particle diameter on the particle-flow parameters in the inviscid region and in the boundary layer is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 65–77, November–December, 2005.     

Effects of gas flow and turbulent diffusion on the centrifugal classification of fine particles

Tthe mass movement of fine particles and changes in their trajectories under the effect of centrifugal, gravitational, and aerodynamic forces is studied by numerical modeling of turbulent swirling flow. The effect of turbulent migration of solid particles due to pulsating gas flow velocities on the aerodynamics of particles and, consequently, on the efficiency of the process of separation of particles by size is analyzed. The studies revealed the behavior laws of probability curve of separation of particles by size and showed the effect of regime and geometric parameters on the process of classification of fine powders.     

Effect of solid mass flux on anisotropic gas–solid flow in risers determined with an LES-SOM model

Within the framework of the two-fluid approach, gas was treated with a large-eddy simulation and a sub-grid-scale (SGS) turbulent kinetic energy model while particles were treated with a second-order-moment method to describe the anisotropy of the fluctuating velocity. A modified Simonin model was derived for the gas–solid interphase fluctuating energy transfer. The anisotropic gas–solid flow in a circulating fluidized bed was investigated. Predictions were in good agreement with experimental data. The distributions of the second-order moment of particles and SGS-turbulent kinetic energy of gas were simulated at different solid mass fluxes. The effects of the solid mass flux on the particle second-order moment, particle anisotropic behavior, gas SGS-turbulent kinetic energy and gas SGS energy dissipation were analyzed for the circulating fluidized bed.     

Investigations of gas and particle dynamics in first generation needle-free drug delivery devices

Abstract. Transdermal powdered drug delivery involves the propulsion of solid drug particles into the skin by means of high-speed gas-particle flow. The fluid dynamics of this technology have been investigated in devices consisting of a convergent-divergent nozzle located downstream of a bursting membrane, which serves both to initiate gas flow (functioning as the diaphragm of a shock tube) and to retain the drug particles before actuation. Pressure surveys of flow in devices with contoured nozzles of relatively low exit-to-throat area ratio and a conical nozzle of higher area ratio have indicated a starting process of approximately 200 s typical duration, followed by a quasi-steady supersonic flow. The velocity of drug particles exiting the contoured nozzles was measured at up to 1050 m/s, indicating that particle acceleration took place primarily in the quasi-steady flow. In the conical nozzle, which had larger exit area ratio, the quasi-steady nozzle flow was found to be overexpanded, resulting in a shock system within the nozzle. Particles were typically delivered by these nozzles at 400 m/s, suggesting that the starting process and the quasi-steady shock processed flow are both responsible for acceleration of the particle payload. The larger exit area of the conical nozzle tested enables drug delivery over a larger target disc, which may be advantageous. Received 12 March 2000 / Accepted 8 June 2000     

The structure of particle-laden,underexpanded free jets

Underexpanded, supersonic gas-particle jets were experimentally studied using the shadowgraph technique in order to examine the influence of the dispersed particles on the shape of the free jet and the structure of the imbedded shock waves. The particle mass loading at the nozzle exit was varied between zero and one, and two sizes of particles (i.e. spherical glass beads) with mean number diameters of 26 and 45 m were used. It was found that the Mach-disc moves upstream towards the orifice with increasing particle loading. The laser light sheet technique was also used to visualize the particle concentration distribution within the particle jet and the spreading rate of the particle jet. Furthermore, the particle velocity along the jet centerline was measured with a modified laser-Doppler anemometer. These measurements revealed that the particles move considerably slower than the gas flow at the nozzle exit. This is mainly the result of the particle inertia, whereby the particles are not accelerated to sonic speed in the converging part of the nozzle.In order to further explore the particle behavior in the free jet, numerical studies were performed by a combined Eulerian/Lagrangian approach for the gas and particle phases, including full coupling between the two phases. The numerical results showed that the application of different particle velocities at the nozzle exit as the inlet conditions, which were below the sonic speed of the gas phase has a significant influence on the free jet shape and the configuration of the shock waves. These results demonstrate that the assumption of equilibrium flow (i.e. zero slip between the phases) at the nozzle exit which has been applied in most of the previous numerical studies is not justified in most cases. Furthermore, the numerical calculations of the free jet shape and the particle velocity along the jet axis were compared with the measurements. Although correlations for rarefaction and compressibility effects in the drag coefficient were taken into account, the particle velocity along the center line was considerably overpredicted.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Experimental study and numerical simulation of the characteristics of a percussive gas–solid separator

The presence of solid particles in the flow of hypersonic wind tunnels damages the appearance of the experiment models in the wind tunnel and influences the accuracy of experimental results. The design of a highly efficient gas–solid separator was therefore undertaken. Particle trajectory imaging methods were used to measure trajectories under different conditions. The flow field and particle movement characteristics for different head angles (HAs) and separation tooth angles (STAs), inlet velocities, and the exhaust gas outlet pressures in the separator, were calculated using simulations based on the discrete phase model. The particle separation efficiency, pressure loss, and flow loss resulting from different structural parameters were also studied. In line with experimental observations, the characteristic angle of particle movements in the separator and the separation efficiency of the separator were found to increase with decreasing HA and with increasing STA. Separation efficiency improves with increasing inlet velocity and with increasing negative pressure of the exhaust gas outlet; however, the corresponding pressure loss and the flow rate of the waste gas also increased.     

Effect of particle size on a two-phase turbulent jet

The effect of particle size on two-phase turbulent jet flow structure is studied in the present experimental investigation. Polystyrene solid particles of 210, 460, and 780 μm were used. The particles' mass loading ratios ranged from 0 to 3.6. The flow Reynolds number was 2 ‘ 10⁴, which was based on the pipe nozzle diameter and the fluid-phase centerline velocity at the nozzle exit. A two-color laser-Doppler anemometer (LDA), combined with the amplitude discrimination method and the velocity filter method, was employed for measurement. The measurement range of the jet flow was from the initial pipe exit to 90D downstream. Results are presented for the mean velocities of particle and fluid phases, the flow's turbulent intensities and the flow's Reynolds stresses. The energy spectra and the correlation functions of the two-phase jet flow were also obtained by using another one-component He-Ne LDA system.     

On the possibility of blowing a gas jet into a supersonic flow without formation of a three-dimensional boundary-layer separation zone

We consider the flow formed by the interaction of a supersonic flow and a transverse sonic or supersonic jet blown at right angles to the direction of the main flow through a nozzle whose exit section is in a flat wall. When a gas jet is blown through a circular opening [1] the pressure rises in front of the jet because of the stagnation of the oncoming flow. This leads to separation of the boundary layer formed on the wall in front of the blowing nozzle. The resulting three-dimensional separation zone leads to a sharp increase in the pressure and the heat fluxes to the wall in front of the blowing nozzle, which is undesirable in many modern applications. The aim of the present investigation was to find a shape of the exit section of the blowing nozzle for which there is no three-dimensional separation zone of the boundary layer in front of the blowing nozzle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 162–165, May–June, 1979.     

Flow of a two-phase medium in a laval nozzle

The back reaction of particles on a gas flow in Laval nozzles was investigated experimentally. Experimental data were obtained that characterize the change produced by the particles of a solid phase in the shape of the sonic line, the pressure distribution on the nozzle profile, and the configuration of the shock waves in the jet. Flow rate coefficients are given for different nozzle profiles and mass fraction and sizes of the particles in the flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 107–111, January–February, 1981.     

MODELLING OF RELEASE OF GAS FROM HIGH-PRESSURE PIPELINES

The problem investigated is the break of a high-pressure pipeline carrying natural single-phase gas which may condensate (retrograde) when the pressure drops. Single-phase non-ideal gas is assumed using a general- ized equation of state. Taking advantage of the choked massflow condition, the break is split into a pipe flow problem and a dispersion flow problem, both solved using a finite difference control volume scheme. The transient flow field from the pipeline break location is expanded analytically, using an approximation of the governing equations, until ambient pressure is reached and matched to the corresponding gas dispersion flow field using as subgrid model a jet box with a time-varying equivalent nozzle area as an internal boundary of the dispersion domain. The turbulence models used for the pipe and dispersion flow fields are an empirical model of Reichard and the k–ϵ model for buoyant flow respectively. The pipe flow simulations indicate that the flow from the pipeline might include dispersed condensate which will affect quantitatively the mass flow rate from the pipeline and qualitatively the gas dispersion if the condensate rains out. The transient dispersion simulation shows that an entrainment flow field develops and mixes supersaturated gas with ambient warmer air to an unsaturated mixture. Because of the inertia of the ambient air, it takes time to develop the entrainment flow field. As a consequence of this and the decay of the mass flow with time, the lower flammability limit of the gas–air mixture reaches its most remote downstream position relatively early in the simulation (about 15 s) and withdraws closer to the break location.     

Numerical simulation of a dense solid particle flow inside a cyclone separator using the hybrid Euler–Lagrange approach

This paper presents a numerical simulation of the flow inside a cyclone separator at high particle loads. The gas and gas–particle flows were analyzed using a commercial computational fluid dynamics code. The turbulence effects inside the separator were modeled using the Reynolds stress model. The two phase gas–solid particles flow was modeled using a hybrid Euler–Lagrange approach, which accounts for the four-way coupling between phases. The simulations were performed for three inlet velocities of the gaseous phase and several cyclone mass particle loadings. Moreover, the influences of several submodel parameters on the calculated results were investigated. The obtained results were compared against experimental data collected at the in-house experimental rig. The cyclone pressure drop evaluated numerically underpredicts the measured values. The possible reason of this discrepancies was disused.     

喷砂冲蚀实验中颗粒轨迹的数值预测

针对冲击磨损实验研究中磨粒群体的运动轨迹难以准确表征的问题,在负压喷射砂粒群冲击Q235钢板的实验中宏观测量了砂粒撞击的速度与位置分布,并使用数值方法模拟了实验砂粒与空气在喷嘴内外的双向耦合过程,以实现负压喷射砂粒群的轨迹预测。计算中提出了非球形粒子在相对马赫数接近1时的曳力模型,以反映空气可压缩引起砂粒表面流动分离的现象,并合理选择Magnus升力模型及壁面反射模型,最终数值预测的砂粒碰撞速度以及撞击位置与实验情况吻合良好。     

Numerical investigation of a perturbed swirling annular two-phase jet

A swirling annular gas–liquid two-phase jet flow system has been investigated by solving the compressible, time-dependent, non-dimensional Navier–Stokes equations using highly accurate numerical methods. The mathematical formulation for the flow system is based on an Eulerian approach with mixed-fluid treatment while an adjusted volume of fluid method is utilised to account for the gas compressibility. Surface tension effects are captured by a continuum surface force model. Swirling motion is applied at the inlet while a small helical perturbation is also applied to initiate the instability. Three-dimensional spatial direct numerical simulation has been performed with parallelisation of the code based on domain decomposition. The results show that the flow is characterised by a geometrical recirculation zone adjacent to the nozzle exit and by a central recirculation zone further downstream. Swirl enhances the flow instability and vorticity and promotes liquid dispersion in the cross-streamwise directions. A dynamic precessing vortex core is developed demonstrating that the growth of such a vortex in annular configurations can be initiated even at low swirl numbers, in agreement with experimental findings. Analysis of the averaged results revealed the existence of a geometrical recirculation zone and a swirl induced central recirculation zone in the flow field.     

Differentially weighted direct simulation Monte Carlo method for particle collision in gas–solid flows

In gas–solid flows, particle–particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles (e.g., agglomeration and fragmentation) in spite of their initial monodisperse particle distribution. The conventional direct simulation Monte Carlo (DSMC) method for particle collision tracks equally weighted simulation particles, which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions. In this study, a new differentially weighted DSMC (DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian–Lagrangian models for hydrodynamics. Three schemes (mass, momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass, momentum or energy of the whole system unchanged respectively. A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency. The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best. Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC, compared with the equally weighted DSMC. Meanwhile, computational cost can be largely reduced in contrast with direct numerical simulation.     

Orientation of cylindrical particles in gas–solid circulating fluidized bed

The orientation of cylindrical particles in a gas–solid circulating fluidized bed was investigated by establishing a three-dimensional Euler–Lagrange model on the basis of rigid kinetics, impact kinetics and gas–solid two-phase flow theory. The resulting simulation indicated that the model could well illustrate the orientation of cylindrical particles in a riser during fluidization. The influences of bed structure and operation parameters on orientation of cylindrical particles were then studied and compared with related experimental results. The simulation results showed that the majority of cylindrical particles move with small nutation angles in the riser, the orientation of cylindrical particles is affected more obviously by their positions than by their slenderness and local gas velocities. The simulation results well agree with experiments, thus validating the proposed model and computation.     

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.     

Numerical solution of an inverse problem of nozzle theory for a two-phase gas-particle mixture

In the present paper gas flows with monodisperse and polydisperse particles in plane and axisymmetric nozzles are calculated by the inverse method [1, 2]. The gas velocity distribution is specified on the axis of symmetry of the nozzle, while the gas and particle parameters are specified in the entrance section. As a result of the numerical integration of a system of equations describing a flow of gas with condensate particles in it we determine the gas and particle parameters, the gas streamlines, and the particle trajectories with allowance for the mutual influence of the gas and particles. One of the gas streamlines is taken as the nozzle contour and the limiting trajectories and pure gas zone are found. A difference method is described which makes it possible to calculate the subsonic, transonic, and supersonic flow regions using a single algorithm, its features are noted, and the results of the calculation for monodisperse mixtures with particle diameters 1 and 5 m and fractions by weight 0.3 are given. A comparison is made with the results of calculations by other methods.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 106–114, July–August, 1986.The authors express their gratitude to N. B. Ponomarev and G. E. Dumnov for their useful discussions and help in carrying out the calculations.     

Oscillating annular liquid membranes

Summary The response of annular liquid membranes to sinusoidal mass flow rate fluctuations at the nozzle exit is analyzed as a function of the amplitude and frequency of the axial velocity fluctuations at the nozzle exit and thermodynamic compression of the gas enclosed by the membrane. It is shown that both the pressure of the gases enclosed by the annular membrane and the axial distance at which the annular membrane merges on the symmetry axis are periodic functions of time which have the same period as that of the mass flow rate fluctuations at the nozzle exit. They are also nearly sinusoidal functions of time for small amplitudes of the mass flow rate fluctuations at the nozzle exit, and exhibit delay and lag times with respect to the sinusoidal axial velocity fluctuations at the nozzle exit. Both the delay and the lag times are functions of the amplitude and frequency of the mass flow rate fluctuations at the nozzle exit and the polytropic exponent. The amplitudes of both the pressure of the gases enclosed by the annular liquid membrane and the convergence length increase and decrease, resp., as the amplitude and frequency of the mass flow rate fluctuations at the nozzle exit, resp., are increased. They also increase as the polytropic exponent is increased.

---

Oszilierende rotationsschalenförmige Flüssigkeitsmembranen

Übersicht Untersucht wird das Verhalten rotationsschalenförmiger Flüssigkeitsmembranen in Abhängigkeit von der Amplitude und Frequenz der axialen Geschwindigkeitsschwankungen an der Düsenmündung und der thermodynamischen Verdichtung des eingeschlossenen Gases, wenn sich der Massestrom an der Düsenmündung sinusförmig ändert. Es wird gezeigt, daß der Druck des eingeschlossenen Gases und der axiale Mündungsabstand des Scheitels der geschlossenen Membran periodische Zeitfunktionen mit der Frequenz der Masseflußschwankung am Düsenaustritt sind. Für kleine Amplituden des Massestroms ist ihr Zeitverhalten ebenfalls fast-sinusförmig, wobei sie bezüglich der sinusförmigen axialen Geschwindigkeitsschwankungen an der Düsenmündung eine Ansprechzeit und Phasenverschiebung aufweisen. Ansprechzeit und Phasenverschiebung sind Funktionen von Amplitude und Frequenz der Massestromschwankung sowie des polytropen Exponenten. Die Amplitude von Gasdruck und Abstand des Membranscheitels von der Düse wächst bzw. fällt mit wachsender Amplitude und Frequenz des Massestroms. Beide nehmen außerdem mit dem polytropen Exponenten zu.