Evaluation of stochastic particle dispersion modeling in turbulent round jets

ODT (one-dimensional turbulence) simulations of particle-carrier gas interactions are performed in the jet flow configuration. Particles with different diameters are injected onto the centerline of a turbulent air jet. The particles are passive and do not impact the fluid phase. Their radial dispersion and axial velocities are obtained as functions of axial position. The time and length scales of the jet are varied through control of the jet exit velocity and nozzle diameter. Dispersion data at long times of flight for the nozzle diameter (7 mm), particle diameters (60 and 90 µm), and Reynolds numbers (10, 000–30, 000) are analyzed to obtain the Lagrangian particle dispersivity. Flow statistics of the ODT particle model are compared to experimental measurements. It is shown that the particle tracking method is capable of yielding Lagrangian prediction of the dispersive transport of particles in a round jet. In this paper, three particle-eddy interaction models (Type-I, -C, and -IC) are presented to examine the details of particle dispersion and particle-eddy interaction in jet flow.     

A new analytical model for thermal stresses in multi-phase materials and lifetime prediction methods

Based on the fundamental equations of the mechanics of solid continuum, the paper employs an analytical model for determination of elastic thermal stresses in isotropic continuum represented by periodically distributed spherical particles with different distributions in an infinite matrix, imaginarily divided into identical cells with dimensions equal to inter-particle distances, containing a central spherical particle with or without a spherical envelope on the particle surface. Consequently, the multi-particle-（envelope）- matrix system, as a model system regarding the analytical modelling, is applicable to four types of multi-phase materials. As functions of the particle volume fraction v, the inter-particle distances dl, d2, d3 along three mutually per- pendicular axes, and the particle and envelope radii, R1 and R2, respectively, the thermal stresses within the cell, are originated during a cooling process as a consequence of the difference in thermal expansion coefficients of phases rep- resented by the matrix, envelope and particle. Analytical-（experimental）-computational lifetime prediction methods for multi-phase materials are proposed, which can be used in engineering with appropriate values of parameters of real multi-phase materials.     

Statistical Analysis of Turbulent Flame-Droplet Interaction: A Direct Numerical Simulation Study

Turbulent combustion of mono-disperse droplet-mist has been analysed based on three-dimensional Direct Numerical Simulations (DNS) in canonical configuration under decaying turbulence for a range of different values of droplet equivalence ratio (?_d), droplet diameter (a_d) and root-mean-square value of turbulent velocity (u^′). The fuel is supplied in liquid phase and the evaporation of droplets gives rise to gaseous fuel for the flame propagation into the droplet-mist. It has been found that initial droplet diameter, turbulence intensity and droplet equivalence ratio can have significant influences on the volume-integrated burning rate, flame surface area and burning rate per unit area. The droplets are found to evaporate predominantly in the preheat zone, but some droplets penetrate the flame front, reaching the burned gas side where they evaporate and some of the resulting fuel vapour diffuses back towards the flame front. The combustion process in gaseous phase takes place predominantly in fuel-lean mode even for ?_d > 1. The probability of finding fuel-lean mixture increases with increasing initial droplet diameter because of slower evaporation of larger droplets and this predominantly fuel-lean mode of combustion exhibits the attributes of low Damköhler number combustion and gives rise to thickening of flame with increasing droplet diameter. The chemical reaction is found to take place under both premixed and non-premixed modes of combustion and the relative contribution of non-premixed combustion to overall heat release increases with increasing droplet size. The statistical behaviours of the flame propagation and mode of combustion have been analysed in detail and detailed physical explanations have been provided for the observed behaviour.     

Investigation of turbulent vapor-air jets in the presence of condensation and the injection of foreign particles

A study is made of flow in turbulent jets when there is condensation of the water vapor contained in them. A necessary condition for condensation in vapor-air jets is formulated. Relations are obtained for the regime of equilibrium condensation. An experimental investigation was made of the local characteristics of an isobaric turbulent vapor jet exhausting into air at rest when condensation develops in the jet and foreign condensation nuclei (smoke particles) and charged particles (ions produced in a corona discharge) are introduced into the flow. Measurements were made of the local characteristics of the condensed disperse phase — the Sauter diameter d₃₂ of the drops and their volume concentration c_s — using the optical method of an integrating diaphragm. It is shown that d₃₂ and₃₂ c_s increase downstream in the main section of the jet. Specific features of temperature measurements using an end-type microthermocouple were established. Quantitative data were obtained about the influence on the condensation of the thermal conditions and the presence of the foreign particles. The conditions under which there is an intensification of the condensation in vapor-air jets in the presence of ions were determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 53–61, May–June, 1984.     

Numerical investigation of nonequilibrium flows of mixtures of fusible metal particles and gases in a laval nozzle

There are numerous papers [1–11] on the determination of the parameters of condensed oxide particles which are formed during combustion of metallized fuels. The ambiguity, and sometimes the contradictoriness, of the test results obtained [3–5, 9–11] indicate the difficulties in conducting correct experimental investigations. In this connection, numerical studies using mixtures of calibrated liquid-metal particles and different gases are of practical interest. Different probes can be calibrated by using calibrated two-phase flows, the two-phase flow around models and probes can be studied, as can the interaction between liquid-metal particles and the front of an aerodynamic compression shock, their intrusion in different entraining media, the interaction between fine particles (particle-projectiles) and large size particles (particle-targets), etc. In many cases, the prehistory of the flow and the parameters of the gas mixture with the particles in the area of the nozzle exit section must be known to investigate the above-mentioned phenomena. The parameters of different nonequilibrium flows of mixtures of gallium particles and gases in a Laval nozzle are investigated numerically in this paper; the maximum diameter (upper boundary of the spectrum) of the particles (d_s = 30 ) which are not destroyed in the nozzle under the effect of the aerodynamic forces and are suitable for use in a calibrated two-phase stream is determined. The computations were carried out in a one-dimensional approximation according to [12–14].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 86–91, March–April, 1976.The authors are grateful to V. K. Starkov and U. G. Pirumov for discussing the results of the research and to N. M. Alekseev for aid in constructing the graphs.     

Flame Structure and Propagation in Turbulent Flame-Droplet Interaction: A Direct Numerical Simulation Analysis

Three-dimensional Direct Numerical Simulations (DNS) in canonical configuration have been employed to study the combustion of mono-disperse droplet-mist under turbulent flow conditions. A parametric study has been performed for a range of values of droplet equivalence ratio ?_d, droplet diameter a_d and root-mean-square value of turbulent velocity u^′. The fuel is supplied entirely in liquid phase such that the evaporation of the droplets gives rise to gaseous fuel which then facilitates flame propagation into the droplet-mist. The combustion process in gaseous phase takes place predominantly in fuel-lean mode even for ?_d>1. The probability of finding fuel-lean mixture increases with increasing initial droplet diameter because of slower evaporation of larger droplets. The chemical reaction is found to take place under both premixed and non-premixed modes of combustion: the premixed mode ocurring mainly under fuel-lean conditions and the non-premixed mode under stoichiometric or fuel-rich conditions. The prevalence of premixed combustion was seen to decrease with increasing droplet size. Furthermore, droplet-fuelled turbulent flames have been found to be thicker than the corresponding turbulent stoichiometric premixed flames and this thickening increases with increasing droplet diameter. The flame thickening in droplet cases has been explained in terms of normal strain rate induced by fluid motion and due to flame normal propagation arising from different components of displacement speed. The statistical behaviours of the effective normal strain rate and flame stretching have been analysed in detail and detailed physical explanations have been provided for the observed behaviour. It has been found that the droplet cases show higher probability of finding positive effective normal strain rate (i.e. combined contribution of fluid motion and flame propagation), and negative values of stretch rate than in the stoichiometric premixed flame under similar flow conditions, which are responsible for higher flame thickness and smaller flame area generation in droplet cases.     

A single-camera technique for simultaneous measurement of large solid particles transported in rapid shallow channel flows

This paper describes a measurement technique that was successfully applied in a study of bed load transport of large spherical solid particles in a shallow and supercritical flow (Fr?=?2.59–3.17) down a steep slope. The experimental condition was characterized by the relatively large solid particle size compared to the flow depth (d _p /h?=?0.23–0.35), and compared to the tracer diameter (d _p /d _t ?≈?130). The technique incorporated particle image velocimetry and particle tracking velocimetry (PTV) to simultaneously measure the characteristics of the two phases. In order to detect true solid particles and to distinguish them from each other and the unwanted objects, a particle characterization (PCR) algorithm based on Hough transform was employed. The output from the PCR process was utilized for PTV, as well as to generate the corresponding tracer images for special needs. Validation tests have confirmed the pixel accuracy and high reliability of the combined technique. Experimental results obtained with the developed technique include flow velocities, particle velocities, and concentration. The analysis has shown that the particle concentration profile followed an exponential relationship of the form similar to that of Rouse’s profiles, despite the large d _p /h ratio. It also revealed the effect of phase interaction, as a low loading rate of light particles on the order of O(10^?3) could yield a noticeable slowdown in the streamwise fluid velocity.     

Study of airborne particles produced by normal impact of millimetric droplets onto a liquid film

The aim of this work is to carry out an experimental investigation into the generation of airborne microparticles when millimetric droplets of aqueous solutions impact onto a liquid film. Impact experiments using 3.9 mm diameter droplets were carried out for Weber numbers between 159 and 808, with a fixed Ohnesorge number of 2 × 10⁻³ and film parameters S _f (the ratio between the thickness of the liquid film h _film and the diameter of the impacting droplet d _i) between 0.3 and 1. Observed results show that the deposition/splashing threshold is independent of the parameter S _f in agreement with the data in the literature. The aerosol measurement results demonstrate the production of solid particles from the evaporation of secondary microdroplets with diameters less than 30 μm formed when splash occurs. The median diameter of these microdroplets is around 20 μm, corresponding to a value of d ₅₀/d _i = 5 × 10⁻³. Taken together, the results show that the mass and the number of particles emitted increase as the Weber number increases. Moreover, at a Weber number of 808, the results show that the mass and number of particles emitted increases as the parameter S _f decreases. In this case, the mean number of microdroplets emitted per impact is equal to 14 for S _f = 1 and equal to 76 for S _f = 0.3.     

Numerical modelling of film cooling with and without mist injection

Two-phase CFD calculations, using a Lagrangian model and commercial code Fluent 6.2.16, were employed to calculate the gas and droplet flows and film cooling effectiveness with and without mist on a flat plate. Two different three dimensional geometries are generated and the effects of the geometrical shape, size of droplets, mist concentration in the coolant flow and temperature of mainstream flow for different blowing ratios are studied. A cylindrical and laterally diffused hole with a streamwise angle of 30° and spanwise angle of 0° are used. The diameter of film cooling (d) hole, and the hole length to diameter ratio (L/d) for both of geometries are 10 mm and 4, respectively. Also the blowing ratio ranges from 1.0 to 2.0, and the mainstream Reynolds number based on the mainstream velocity and hole diameter (Re _d) is 6,219. The results are shown for different droplets diameters (1–10 μm), concentrations (1–5%) and mainstream temperatures (350–500 K). The centreline effectiveness and distribution of effectiveness on the surface of cooling wall are presented.     

Lattice–Boltzmann simulations for analysing the detachment of micron-sized spherical particles from surfaces with large-scale roughness structures

Fully resolved numerical simulations of a micron-sized spherical particle residing on a surface with large-scale roughness are performed by using the Lattice–Boltzmann method. The aim is to investigate the influence of surface roughness on the detachment of fine drug particles from larger carrier particles for transporting fine drug particles in a DPI (dry powder inhaler). Often the carrier surface is modified by mechanical treatments for modifying the surface roughness in order to reduce the adhesion force of drug particles. Therefore, drug particle removal from the carrier surface is equivalent to the detachment of a sphere from a rough plane surface. Here a sphere with a diameter of 5 μm at a particle Reynolds number of 1.0, 3.5 and 10 are considered. The surface roughness is described as regularly spaced semi-cylindrical asperities (with the axes oriented normal to the flow direction) on a smooth surface. The influence of asperity distance and size ratio (i.e. the radius of the semi-cylinder to the particle radius, R_c/R_d) on particle adhesion and detachment are studied. The asperity distance is varied in the range 1.2 < L/R_d < 2 and the semi-cylinder radius between 0.5 < R_c/R_d < 0.75. The required particle resolution and domain size are appropriately selected based on numerical studies, and a parametric analysis is performed to investigate the relationship between the contact distance (i.e. half the distance between the particle contact points on two neighbouring semi-cylinders), the asperity distance, the size ratio, and the height of the particle centroid from the plane wall. The drag, lift and torque acting on the spherical particle are measured for different particle Reynolds numbers, asperity distances and sizes or diameters. The detachment of particles from rough surfaces can occur through lift-off, sliding and rolling, and the corresponding detachment models are constructed for the case of rough surfaces. These studies will be the basis for developing Lagrangian detachment models that eventually should allow the optimisation of dry powder inhaler performance through computational fluid dynamics.     

Particle transport in Poiseuille flow in narrow channels

Particle-tracking experiments were performed to validate a model [Staben, M.E., Zinchenko, A.Z., Davis, R.H., 2003. Motion of a particle between two parallel plane walls in low-Reynolds-number Poiseuille flow. Phys. Fluids 15, 1711–1733] for neutrally buoyant spherical particles convected by a Poiseuille flow in a thin microchannel for particles as large as d_p/H = 0.95, where d_p is the particle diameter and H is the channel width (narrow dimension). The measured and predicted velocities agree within experimental error and show that a particle’s velocity is more retarded when it is larger and/or closer to a channel wall. The particle distribution across the channel for a blunt entrance shows a focusing of small particles away from the walls and towards the center of the channel, whereas the particle distribution for an offset-angled entrance is slightly skewed towards the wall encountered first in the entrance region. As a result, the average particle velocities for the blunt entrance exceed those for the angled entrance. Moreover, due to the depletion of particles from the slow-moving region within one radius of the wall, the average particle velocity exceeds the average fluid velocity unless the particle diameter exceeds about 80% of the channel width.     

Experimental investigation of a confined flow downstream of a circular cylinder centred between two parallel walls

In this work, we present an experimental study of the wall confinement effect on the wake formation behind a circular cylinder of diameter d_c=10 mm and of length L_c=30d_c. The experiments were performed in a water tunnel with the dimensions (length=300d_c, height=3d_c, span L_c=30d_c). The confinement rate was r=1/3 and the Reynolds number was in the range of 30–277. The experiments were done using 2-D PIV measurements. The first instability was delayed by the confinement and the von Kármán vortices characteristics are different from the unconfined case. Proper orthogonal decomposition (POD) of the flow was used for a filtering purpose and to extract the energetic contribution of the different modes. A low-order representation of the flow, constructed from the first pair of modes in the well-defined region of the flow, shows that von Kármán vortices are equivalent to vanishing progressive waves. Measurements done above the cylinder show the presence of 3-D span instabilities showing great similarities with “Mode A” and “Mode B” found in the unconfined case.     

Influence of Internal Structure and Medium Length on Transport and Deposition of Suspended Particles: A Laboratory Study

A laboratory study was undertaken on the transport and the deposition of suspended particles (silt of modal diametre 6 μm) in three columns of different length, filled with glass beads or gravel. Tracer tests were carried out at various flow velocities by short pulses of a mixture of suspended particles/dissolved tracer. The breakthrough curves were competently described with the analytical solution of a convection dispersion equation with a first-order deposition rate and the hydro-dispersive parameters were deduced. For the same experimental conditions, the results showed a difference in the behaviour of the suspended particles transport and deposition rates within the two porous media tested. The internal structure of both media governs the particle-grain collision frequency as well as the particles trapping. The scale effect was highlighted and affects the dispersivity, the size exclusion effect, the recovery rates and the deposition rates. Longitudinal dispersion increases with mean pore velocity and is described with a nonlinear relationship. The dispersivity increases with the column length. The size exclusion effect is more important in the short column. The recovery rate increases with flow velocity and decreases while increasing column length. The deposition rates increases until a critical flow velocity then decreases. This critical velocity is also sensitive to the scale effect, and increases with the column length.     

Thermocapillary migration of liquid droplets in a temperature gradient in a density matched system

An experimental investigation of thermocapillary flow in droplets of a vegetable oil (partially hydrogenated soybean oil) immersed in silicone oil was conducted in a test cell with a heated top wall and a cooled bottom wall. The liquids are nearly immiscible and have equal densities at a temperature below the room temperature, thus providing a simulation of low-gravity conditions by reducing the buoyancy forces. The interfacial tension between the two oils was measured in the temperature range 20–50°C using a capillary tube and d/dT was determined to be negative. Droplets ranging in sizes from 3 mm to 1 cm diameter were injected into the silicone oil. The vertical temperature profile in the bulk liquid (silicone oil) produces temperature variations along the interface which induce variations in the interfacial tension. The flow inside the droplet driven by the resulting interfacial shear stresses was observed using a laser light-sheet flow visualization technique. The flow direction is consistent with the sign of d/dT. The observed maximum surface velocities are compared to the theoretical predictions of Young et al. (1959).For short times after injection, the droplets were driven by this flow towards the hot wall above the matched-density temperature until the droplets reached a point where the forces due to the flow and buoyancy were equal. After longer times, the droplets moved to the cold side due to suspected density changes caused by mass transfer from the droplets to the silicone oil. This was confirmed by tests under isothermal conditions, where it was observed that droplets of all sizes fell to the cold bottom eventually.Thus, even though the thermocapillary flow inside the droplets persisted for long times in spite of the mass transfer, the migration of droplets towards the hot side was eventually affected by uncontrolled buoyancy forces resulting from density changes due to mass transfer. While additional liquids are being tried, it is suggested from the present experience that reduced gravity experiments will probably be necessary to provide unambiguous data for the migration of droplets.     

High-throughput generation of uniform droplets from parallel microchannel droplet generators and the preparation of polystyrene microsphere material

To prepare uniform polystyrene particles with ten microns of diameter, a parallel scaling-up strategy for the capillary-assembled stepwise microchannel was developed, which created uniform droplets with high-throughput and formed a large amount of emulsion templates for the polymerization of styrene and cross-linker. The microchannel droplet generator was robust for the flow rate deviation of the continuous phase in the jetting flow, and droplet generation frequency up to 2.8 × 10⁴ Hz was achieved with only four parallel droplet generators, which were much more efficient than the parallelly scaled microfluidic devices working in dripping flow. 32–52 μm average diameter droplets with 4.5%–8.4% diameter variation coefficients were successfully prepared from the microchannel device fabricated by low-cost 3D-print method, and the droplets were subsequently turned to solid particles via a two-step polymerization in the platform. The polystyrene particles were further reduced to 16.9–23.5 μm with 5.0%–8.6% diameter variation coefficients due to the accompanying emulsion polymerization, and the working capacity of the platform reached hundred milligrams of particles per hour.     

Determination of particle charge and size from the propagation of ultrasound in suspensions

The propagation of small perturbations in raulticomponent disperse media consisting of an uncharged dispersion fluid, positive and negative ions and charged particles or droplets of another fluid is investigated. When weak waves pass through emulsions and suspensions, because of the difference in the velocities of the ions and charged particles a non-uniform distribution of electric potential develops in the medium [1–3]. Expressions relating the amplitude of the electric potential and the amplitude of the fluid velocity in the wave, the particle charge and the parameters characterizing the medium are derived. Relations are obtained for the phase shift between the values of the electric potential and the fluid velocity. It is proposed to use the expressions obtained, which describe the propagation of ultrasound, for the experimental determination of the particle charge and other parameters of the disperse medium, in particular, the particle size.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 122–128, January–February, 1988.     

Simultaneous air/fuel-phase PIV measurements in a dense fuel spray

A new diagnostic has been developed that is capable of obtaining simultaneous two-phase velocity measurements in a gasoline direct-injection fuel spray. This technique utilizes a two-laser (double-pulse) two-camera (double-frame) setup to simultaneously image the injected fuel and entrained air to determine the 2D velocity vector fields of both phases using cross-correlation particle image velocimetry (PIV). The air phase is visualized through fluorescence from seeding particles introduced into the static measurement volume while Mie scattering signals are collected from the fuel droplets. The combination of different laser wavelengths and a spectral signal shift for the air phase allows spectral separation of the signals. Independent timing of the laser pulses permits optimized adaptation of the velocity dynamic range for the two phases to account for the large difference in velocities between air and fuel droplets.     

Transient behavior of heat removal from a cylindrical heat storage vessel packed with spherical porous particles

The transient heat transfer behavior in the case of heat removal from a cylindrical heat storage vessel packed with spherical particles was investigated experimentally for various factors (flow rate, diameter of spherical particles packed, temperature difference between flowing cold air and spherical particles accumulating heat, and physical properties of spherical particles). The experiments were covered in ranges of Reynolds number based on the mean diameter of spherical particles packed Re_d = 10.3–2200, porosity?=0.310 to 0.475, ratio of spherical particle diameter to cylinder diameterd/D = 0.0075–0.177 and ratio of length of the cylinder to cylinder diameterL/D=2.5–10. It was found that especially the flow rate and the dimension of spherical particles played an important role in estimating the transient local heat transfer characteristics near the wall of the cylindrical vessel in the present heat storage system. As flow rate and diameter of spherical particles were increased under a given diameter of the cylinder heat storage vessel, the mean heat transfer coefficient between the flow cold air and the hot spherical particles increased and the time period to finish removing heat from the vessel reduced. In addition, the useful experimental correlation equations of mean heat transfer coefficient between both phases and the time period to finish removing heat from the vessel were derived with the functional relationship of Nusselt numberNu _d=f [modified Prandtl numberPr ^* (d/D), Re_d) and Fourier numberFo = f(d/D, L/D, Pr^*, Re_d).     

Solids behavior in dilute zone of a CFB riser under turbulent conditions

The behavior of the solid phase in the upper zone of a circulating fluidized bed riser was studied using a phase Doppler anemometer. Glass particles of mean diameter 107 μm and superficial gas velocities U_g covering the turbulent and the beginning of the fast fluidization regime were investigated. Three static bed heights were tested. Ascending and descending particles were found co-existing under all operating conditions tested, and at all measurement locations. Superficial gas velocity proved/happened to have a larger effect on descending particles at the wall and on ascending particles in the central region. Transversal particle velocities in both directions (toward the center and toward the wall) behaved relatively equivalently, with only slight difference observed at the wall. However, observation of the number of particles moving in either transversal direction showed a change in bed structure when increasing U_g. Furthermore, a balance was constantly observed between the core zone and the annulus zone where the mutual mass transfer between these two zones occurred continuously. Transition from a slow to a fast particle motion was accompanied by a transition to high levels of velocity fluctuations, and was found corresponding to the appearance of significant solid particle flow rate.     

Simultaneous velocity field measurements in two-phase flows for turbulent mixing of sprays by means of two-phase PIV

This paper describes a novel derivative of the PIV method for measuring the velocity fields of droplets and gas phases simultaneously using fluorescence images rather than Mie scattering images. Two-phase PIV allows the simultaneous and independent velocity field measurement of the liquid phase droplets and ambient gas in the case of two-phase flow mixing. For phase discrimination, each phase is labelled by a different fluorescent dye: the gas phase is seeded with small liquid droplets, tagged by an efficient fluorescent dye while the droplets of the liquid phases are tagged by a different fluorescent dye. For each phase, the wavelength shift of fluorescence is used to separate fluorescence from Mie scattering and to distinguish between the fluorescence of each phase. With the use of two cross-correlation PIV cameras and adequate optical filters, we obtain two double frame images, one for each phase. Thus standard PIV or PTV algorithms are used to obtain the simultaneous and independent velocity fields of the two phases. Because the two-phase PIV technique relies on the ability to produce two simultaneous and independent images of the two phases, the choice of the labelling dyes and of the associated optical filter sets is relevant for the image acquisition. Thus a spectroscopic study has been carried out to choose the optimal fluorescent dyes and the associated optical filters. The method has been evaluated in a simple two-phase flow: droplets of 30–40 μm diameter, produced by an ultrasonic nozzle are injected into a gas coflow seeded with small particles. Some initial results have been obtained which demonstrate the potential of the method.