Atomization of liquids by disintegrating thin liquid films using gas jets

Liquid atomization is useful in many applications, such as engineering, science, pharmaceutics, medicine, forensics and others. In the present research, an innovative methodology and a new device for atomization of liquids into mists of micron and submicron droplets have been developed. The new liquid-atomization method exploits the physical phenomenon of fragmentation of thin liquid films into fine micron and submicron droplets by gas jets. For several tested prototypes, the direct observations using a high-speed visualization technique have demonstrated that bubbles were generated within a liquid and their shells have been subsequently destroyed by applying a mechanical impulse (pressure of a compressed air) once the bubbles came over the liquid surface. The main characteristics of the generated tap water mists have been experimentally measured by means of the laser diffraction technique under various conditions for each prototype. One of the prototype devices allowed obtaining mists containing 90–99% of droplets smaller than 1 µm, with the minimum arithmetic and Sauter mean droplet diameters of 1.48 µm and 2.66 µm, and the 2.64 ml/min of droplet flow rate for 3.5 bar manometer pressure of atomizing air. The gas to liquid mass ratios (GLR) in the new device are depending on the atomizing tube length and the number of perforated orifices in the tube: more the tube length, hence more the number of perforated orifices, and therefore more liquid droplets will form for the same gas flow rate. The measured GLR values related to 1 m length of the utilized atomizing tube were in the range of 0.65–1.06, and for the specifically utilized atomizing tube of 72 mm length were among 9.07–14.67. The results of this study demonstrate that the developed method of generation of very fine droplet mists has many advantages over the existing techniques and can be perspective for many practical applications.     

Flow patterns and pressure drop of ionic liquid–water two-phase flows in microchannels

The two-phase flow of a hydrophobic ionic liquid and water was studied in capillaries made of three different materials (two types of Teflon, FEP and Tefzel, and glass) with sizes between 200 μm and 270 μm. The ionic liquid was 1-butyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}amide, with density and viscosity of 1420 kg m⁻³ and 0.041 kg m⁻¹ s⁻¹, respectively. Flow patterns and pressure drop were measured for two inlet configurations (T- and Y-junction), for total flow rates of 0.065–214.9 cm³ h⁻¹ and ionic liquid volume fractions from 0.05 to 0.8. The continuous phase in the glass capillary depended on the fluid that initially filled the channel. When water was introduced first, it became the continuous phase with the ionic liquid forming plugs or a mixture of plugs and drops within it. In the Teflon microchannels, the order that fluids were introduced did not affect the results and the ionic liquid was always the continuous phase. The main patterns observed were annular, plug, and drop flow. Pressure drop in the Teflon microchannels at a constant ionic liquid flow rate, was found to increase as the ionic liquid volume fraction decreased, and was always higher than the single phase ionic liquid value at the same flow rate as in the two-phase mixture. However, in the glass microchannel during plug flow with water as the continuous phase, pressure drop for a constant ionic liquid flow rate was always lower than the single phase ionic liquid value. A modified plug flow pressure drop model using a correlation for film thickness derived for the current fluids pair showed very good agreement with the experimental data.     

Gas temperature in the trace of water droplets streamlined by hot air flow

In order to obtain the knowledge necessary for developing new effective fire extinguishing technologies, we determined experimentally the gas temperature in the trace of water droplets streamlined by hot air flow. It was important to establish how much the temperature in the droplet trace decreases and how fast it recovery to the initial temperature field after the droplet evaporation. The following parameters were varied: droplet size from 1.3 mm to 1.7 mm, velocity from 1 m/s to 5 m/s, initial airflow temperature from 473 K to 773 K, number of droplets (one or two), and the arrangement of droplets relative to the hot inflow (serial or parallel). The study proves the theoretical hypothesis about a significant influence of evaporation on the temperature in the water droplet trace. When a temperature trace of water droplets is formed, irrespective of their arrangement, the role of the evaporation process strengthens with the gas flow temperature rising. Furthermore, the study specifies typical longitudinal dimensions of the aerodynamic and temperature traces of water droplets. It has been established that when droplets are located in series and in parallel, their combined impact on the temperature and velocity of the gas flow in the medium differs rather considerably.     

Drag reduction for liquid flow through micro-apertures

Pressure drops in the flow through micro-orifices and capillaries were measured for silicone oils, aqueous solutions of polyethylene glycol (PEG), and surfactant aqueous solutions. The diameter of micro-orifices ranged from 5 μm to 400 μm. The corresponding length/diameter ratio was from 4 to 0.05 and capillary diameters were 105 μm and 450 μm. The following results were obtained: silicone oils of 10^?6 m²/s and 10^?5 m²/s in kinematic viscosity generated a reduction of pressure drop (RPD), that is, drag reduction, similar to the RPD of water and a glycerol/water mixture reported in the previous paper by the present authors. When RPD occurred, the pressure drop (PD) of silicone oils of 10^?6 m²/s and 10^?5 m²/s had nearly the same magnitude. Namely, the difference in viscosity did not influence RPD. A 10³ ppm aqueous solution of PEG20000 provided almost the same PD as that of PEG8000 for the 400 μm to 15 μm orifices, but a greater PD than that of PEG8000 for the 10 μm to 5 μm orifices. A non-ionic surfactant and a cationic surfactant were highly effective in RPD compared with anionic surfactants: the non-ionic and cationic surfactant solutions had PD one order of magnitude lower than that of water under some flow conditions in the concentration range from 1 ppm to 10⁴ ppm, but the anionic surfactant solutions did not generate RPD except in the case of the smallest orifice of 5 μm in diameter. The PD of the non-ionic surfactant solution showed a steep rise at a Reynolds number (Re_t) for 400 μm to 15 μm orifices. The Re_t provides the relationship Re_t = K/D, where D is the orifice diameter, and K is a constant of 2 × 10^?2 m for the 100–20 μm orifices irrespective of liquid concentration. Capillary flow experiment revealed that the PEG, non-ionic and cationic surfactant solutions generated RPD also in a laminar flow through the capillary of 105 μm in diameter, but not in the flow through the capillary of 450 μm in diameter. In order to clarify the cause of RPD, an additional experiment was carried out by changing the orifice material from metal to acrylic resin. The result gave a different appearance of RPD, suggesting that RPD is related to an interfacial phenomenon between the liquid and wall. The large RPDs found in the present experiment are very interesting from both academic and practical viewpoints.     

Droplet size and velocity characteristics of water-air impinging jet atomizer

A water-air impinging jets atomizer is investigated in this study, which consists of flow visualization using high speed photography and mean droplet size and velocity distribution measurements of the spray using Phase Doppler Anemometry (PDA). Topological structures and break up details of the generated spray in the far and near fields are presented with and without air jet and for an impinging angle of 90°. Spray angle increases with the water jet velocity, air flow rate and impinging angle. PDA results indicate that droplet size is smallest in the spray center, with minimum value of Sauter mean diameter (SMD) of 50 µm at the air flow rate of Q_m = 13.50 g/min. SMD of droplets increases towards the spray outer region gradually to about 120 µm. The mean droplet velocity component W along the air-jet axis is highest in the spray center and decreases gradually with increasing distance from the spray center. SMD normalized by the air nozzle diameter is found firstly to decrease with gas-to-liquid mass ratio (GLR) and air-to-liquid momentum ratio (ALMR) and then remain almost constant. Its increasing with aerodynamic Weber number indicates an exponential variation. The study sheds light on the performance of water-air impinging jets atomizers providing useful information for future CFD simulation works.     

The atomization of water–oil emulsions

The paper presents the results of experimental studies on atomization of the emulsions flowing through twin-fluid atomizers obtained by the use of the digital microphotography method. The main elements of the test installation were: nozzle, reservoir, pump and measurement units of liquid flow. The photographs were taken by a digital camera with automatic flash at exposure time of 1/8000 s and subsequently analyzed using Image Pro-Plus. The oils used were mineral oils 20–90, 20–70, 20–50 and 20–30. The studies were performed at flow rates of liquid phase changed from 0.0014 to 0.011 (dm³/s) and gas phase changed from 0.28 to 1.4 (dm³/s), respectively. The analysis of photos shows that the droplets being formed during the liquid atomization have very different sizes. The smallest droplets have diameters of the order of 10 μm. The experimental results showed that the changes in physical properties of a liquid phase lead to the significant changes in the spray characteristics. The analysis of the photos of water and emulsions atomization process showed that the droplet sizes are dependent on gas and liquid flow rates, construction of nozzle and properties of liquid. The differences between characteristics of atomization for water and emulsions have been observed. Analysis of photos on forming the droplets in air–water and air-emulsions systems showed that droplets are bigger in air-emulsion system (at the same value of gas to liquid mass ratio). The values of Sauter mean diameter (SMD) increased with increase of volume fraction of oil in emulsion. The droplet size increased with emulsion viscosity.     

Atomization of glycerin with a twin-fluid swirl nozzle

Atomization of liquids with high viscosity is always a challenge, especially when small diameter droplets and high liquid flow rates are simultaneously required. In the present research, the performance of a Venturi–vortex twin-fluid swirl nozzle is examined, attending to its capabilities to generate droplets with diameters below 20 µm when atomizing pure glycerin at room temperature. In this nozzle, air is injected tangentially in a central convergent section, and discharges suctioning the liquid fed to a coaxial chamber, here using a gear pump. The resulting spray is visualized and analyzed. Droplet size distributions are measured with a laser diffractometer. As expected, droplet diameter increases with liquid flow rate, and quickly diminishes when air flow rate is increased. Sauter mean diameters (SMD) below 15 µm can be obtained even when atomizing pure glycerin. However, these values are obtained for relatively low glycerin flow rates (∼5 l/h), and with rather wide distributions. For 10 l/h and an air-to-liquid mass flow rate ratio (ALR) of 13.7 more than 26% of the glycerin volume is atomized in droplets smaller than 20 µm. Liquid ligaments are observed near the nozzle exit, but they tend to break up while moving downstream.     

Coalescence of sessile droplets of varying viscosities for line printing

Coalescence of sessile droplets is studied experimentally with water–glycerin mixtures of different viscosities. Effects of viscosity on the dimensionless spreading length (Ψ) and the center-to-center distance (L) are investigated for two droplets; the first droplet (D_s) is stationary on a substrate and the second droplet (D₀) landing at a center-to-center distance L from the first droplet. For a low viscosity fluid, Ψ is maximum when L approaches zero (or λ → 1, where λ = 1 − L/D_s), which represents a head-on collision. For a high viscosity fluid, Ψ is minimum when λ → 0.6. The effect of λ on line printing for various viscosities is also examined by printing multiple droplets. We found that the larger the viscosity, the less the breakup between droplets; viscosities smaller than 60 wt% glycerin yielded line breakup. The overlap ratio of λ > 0.3 produced not a line, but a bigger droplet or puddle because of coalescence. Data obtained in this work can provide insights for the fabrication of conductive microtracks or microinterconnects in printed-electronics applications where a line breakup between droplets would lead to an electrical circuit short.     

Origin of regime transition to turbulent flow in bubble column: Orifice- and column-induced transitions

The possible events during bubble formation on an orifice were investigated using a rectangular bubble column (30 cm × 30 cm × 100 cm). The gas flow rate through a single orifice was adjusted from 0.1 dm³/min to 5.0 dm³/min covering a high flow rate regime. At the high gas flow rate, the bubble formation process was complicated by diverse events, such as wake effect, channeling, and orifice-induced turbulent flow. The detachment period could be used to discern the bubble formation steps because it was strongly affected by the above events. The bubble size distribution around the orifice was also analyzed to gain a clearer understanding of the bubble formation process. Above the rate of 3.0 dm³/min through a single orifice, the detachment period converged to a value of 25 ms irrespective of the orifice diameter. The bubble size distribution also showed little difference in this range of gas flow rate. This could be explained by the development of turbulent flow around the orifice. A 0.15 m in-diameter bubble column was tested to investigate the effect of orifice-induced turbulent flow on the regime transition in which the homogeneous flow regime is converted into the heterogeneous flow regime in the column. Obvious distinction between the orifice- and column-induced transitions was observed.     

Hydrodynamics and energy consumption studies in a three-phase liquid circulating three-phase fluid bed contactor

The hydrodynamics and energy consumption have been studied in a cold flow, bubbling and turbulent, pressurized gas–liquid–solid three-phase fluidized bed (0.15 m ID × 1 m height) with concurrent gas–liquid up flow is proposed with the intention of increasing the gas hold up. The hydrodynamic behaviour is described and characterised by some specific gas and liquid velocities. Particles are easily fluidized and can be uniformly distributed over the whole height of the column. The effect of parameters like liquid flow rate, gas flow rate, particle loading, particle size, and solid density on gas hold up and effect of gas flow rate, solid density and particle size on solid hold up, energy consumption and minimum fluidization velocity has been studied. At the elevated pressures a superior method for better prediction of minimum fluidization velocity and terminal settling velocities has been adopted. The results have been interpreted with Bernoulli’s theorem and Richardson–Zaki equation. Based on the assumption of the gas and liquid as a pretend fluid, a simplification has been made to predict the particle terminal settling velocities. The Richardson–Zaki parameter n′ was compared with Renzo’s results. A correlation has been proposed with the experimental results for the three-phase fluidization.     

An experimental and numerical study of the localization behavior of tantalum and stainless steel

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic–plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material’s local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length is 870–890 μm with deformation imposed through a split-Hopkinson pressure bar system at maximum top-to-bottom velocity in the range of 8–25 m/s. We present experimental results of the deformation response of tantalum and 316L stainless steel samples. The tantalum samples did not form shear bands but the stainless steel sample formed a late stage shear band. We have also modeled these experiments using both conductive and adiabatic continuum models. An anisotropic elasto-viscoplastic constitutive model with damage evolution was used within the finite element code EPIC. A Mie-Gruneisen equation of state and the rate and temperature sensitive MTS flow stress model together with a Gurson flow surface were employed. The models performed well in predicting the experimental data. The numerical results for tantalum suggested a maximum equivalent strain rate on the order of 7 × 10⁴ s⁻¹ in the gage section for an imposed top surface displacement rate of 17.5 m/s. The models also suggested that for an initial temperature of 298 K a temperature in the neighborhood of 900 K was reached within the shear section. The numerical results for stainless steel suggest that melting temperature was reached throughout the shear band shortly after peak load. Due to sample geometry, the stress state in the shear zone was not pure shear; a significant normal stress relative to the shear zone basis line was developed.     

Experimental investigation into transient pressure pulses during pneumatic conveying of fine powders using Shannon entropy

This paper presents the results of an ongoing investigation into transient pressure pulses using Shannon entropy. Pressure fluctuations (produced by gas–solid two-phase flow during fluidized dense-phase conveying) are recorded by pressure transducers installed at strategic locations along a pipeline. This work validates previous work on identifying the flow mode from pressure signals (Mittal, Mallick, & Wypych, 2014). Two different powders, namely fly ash (median particle diameter 45 μm, particle density 1950 kg/m³, loosely poured bulk density 950 kg/m³) and cement (median particle diameter 15 μm, particle density 3060 kg/m³, loosely poured bulk density 1070 kg/m³), are conveyed through different pipelines (51 mm I.D. × 70 m length and 63 mm I.D. × 24 m length). The transient nature of pressure fluctuations (instead of steady-state behavior) is considered in investigating flow characteristics. Shannon entropy is found to increase along straight pipe sections for both solids and both pipelines. However, Shannon entropy decreases after a bend. A comparison of Shannon entropy among different ranges of superficial air velocity reveals that high Shannon entropy corresponds to very low velocities (i.e. 3–5 m/s) and very high velocities (i.e. 11–14 m/s) while low Shannon entropy corresponds to mid-range velocities (i.e. 6–8 m/s).     

CFD modelling and wind tunnel validation of airflow through plant canopies using 3D canopy architecture

The efficiency of pesticide application to agricultural fields and the resulting environmental contamination highly depend on atmospheric airflow. A computational fluid dynamics (CFD) modelling of airflow within plant canopies using 3D canopy architecture was developed to understand the effect of the canopy to airflow. The model average air velocity was validated using experimental results in a wind tunnel with two artificial model trees of 24 cm height. Mean air velocities and their root mean square (RMS) values were measured on a vertical plane upstream and downstream sides of the trees in the tunnel using 2D hotwire anemometer after imposing a uniform air velocity of 10 m s^?1 at the inlet. 3D virtual canopy geometries of the artificial trees were modelled and introduced into a computational fluid domain whereby airflow through the trees was simulated using Reynolds-Averaged Navier–Stokes (RANS) equations and k-ε turbulence model. There was good agreement of the average longitudinal velocity, U between the measurements and the simulation results with relative errors less than 2% for upstream and 8% for downstream sides of the trees. The accuracy of the model prediction for turbulence kinetic energy k and turbulence intensity I was acceptable within the tree height when using a roughness length (y₀ = 0.02 mm) for the surface roughness of the tree branches and by applying a source model in a porous sub-domain created around the trees. The approach was applied for full scale orchard trees in the atmospheric boundary layer (ABL) and was compared with previous approaches and works. The simulation in the ABL was made using two groups of full scale orchard trees; short (h = 3 m) with wider branching and long (h = 4 m) with narrow branching. This comparison showed good qualitative agreements on the vertical profiles of U with small local differences as expected due to the spatial disparities in tree architecture. This work was able to show airflow within and above the canopy in 3D in more details.     

Irreversible thermodynamics of triple junctions during the intergranular void motion under the electromigration forces

A rigorous reformulation of internal entropy production and the rate of entropy flow is developed for multi-component systems consisting of heterophases, interfaces and/or surfaces. The result is a well-posed moving boundary value problem describing the dynamics of curved interfaces and surfaces associated with voids and/or cracks that are intersected by grain boundaries. Extensive computer simulations are performed for void configuration evolution during intergranular motion. In particular we simulate evolution resulting from the action of capillary and electromigration forces in thin film metallic interconnects having a “bamboo” structure, characterized by grain boundaries aligned perpendicular to the free surface of the metallic film interconnects. Analysis of experimental data utilizing previously derived mean time to failure formulas gives consistent values for interface diffusion coefficients and enthalpies of voids. 3.0 × 10⁻⁶ exp(−0.62 eV/kT) m² s⁻¹ is the value obtained for voids that form in the interior of the aluminum interconnects without surface contamination. 6.5 × 10⁻⁶ exp(−0.84 eV/kT) m² s⁻¹ is obtained for those voids that nucleate either at triple junctions or at the grain boundary-technical surface intersections, where the chemical impurities may act as trap centers for hopping vacancies.     

Flow boiling of halogenated refrigerants at high saturation temperature in a horizontal smooth tube

Several correlations are available in the open literature for computing the heat transfer coefficient during flow boiling inside plain channels. With respect to halogenated refrigerants, these correlations are usually compared to data taken in a limited range of evaporation temperature and reduced pressure. More recently, the adoption of new refrigerants, such as high pressure HFCs and carbon dioxide, requires to largely extend the pressure range of application of such correlations. Besides, the design of evaporators for some heat pumping applications, where temperatures are set at higher values as compared to usual evaporating temperatures in air-conditioning equipment, requires proper validation of the computing methods.The present paper aims at comparing four well-known predicting models to a new database collected during flow boiling of HCFC (R22) and HFC refrigerants (R134a, R125 and R410A) in a horizontal 8 mm internal diameter tube. This database is characterized by saturation temperature ranging between 25 °C and 45 °C, reduced pressure spanning between 0.19 and 0.53. Mass velocity ranges between 200 and 600 kg m^?2 s^?1 and heat flux between 9 and 53 kW m^?2.Evaporating heat transfer coefficients of halogenated refrigerants at such high temperatures have not been reported in the open literature so far. The discussion of the results will enlighten some similarities with experimental trends presented in the literature for evaporating carbon dioxide.Two models tested here show good prediction capabilities of the present experimental data, but not for all the data sets in the same way. For the purpose of practical use, a simple modification of the correlation by Gungor and Winterton [1] is proposed, showing that this is able to catch the experimental trends of the present database with good agreement.     

The contribution of small leaks in a baghouse filter to dust emission in the PM2.5 range—A system approach

The contribution of leakage in a baghouse filter (defined as a short circuit between the upstream and downstream sides of the filter) to the emission of fine particles is quantified in comparison to other dust emission sources, and the influence of key operating variables on overall system response is analyzed. The study was conducted on a well-maintained pilot-scale filter unit (9 bags of 500 g/m² calendered polyester needle felt; total surface area 4.2 m²) operated in Δp-controlled mode over a range of pulsing intensities, with two types of test dust (one free-flowing and the other cohesive) at inlet concentrations of 10 and 30 g/m³. Leaks included single holes between 0.5 and 4 mm diameter, intentionally placed in either the plenum plate or one of the filter bags, as well as seamlines from bag confectioning. Emissions were separated by source into a transient contribution due to dust penetration through the filter bags after each cleaning pulse, and a continuous contribution from leaks. This separation was based on a novel method of data processing that relies on time-resolved concentration measurements with a specially calibrated optical particle counter. Tiny leaks on the order of 1 mm generated the same emission level as all the bags combined, and dominated continuous emissions. The equivalent leak cross section (leakage = media emission) was about 1 ppm of the total installed filter surface, independent of upstream dust concentration. Leakage through open seamlines amounted to 75% of media emissions in case of free-flowing test dust. Leakage was restricted to aerodynamic diameters less than ∼5 μm (roughly the PM_2.5 mass fraction). For comparison, time-averaged mass penetration through conventional needle-felt media ranged from about 10⁻⁵ to 10⁻⁶, depending on cohesiveness of the particle material and pulse cleaning intensity, giving emission levels between about 0.02 and 0.2 mg/m³ at the reference concentration of 10 g/m².     

Experimental study of the blockage boundary for dense-phase pneumatic conveying of powders through a horizontal slit

The estimation of the blockage boundary for pneumatic conveying through a slit is of significant importance. In this paper, we investigate the characteristics for blockage of powder (48 μm average diameter) through a horizontal slit (1.6 m × 0.05 m × 0.002 m). The results show that the required critical solid mass flow rate increases as the superficial air velocity increases superficial air velocity. The solid loading ratio and superficial air velocity displayed a decreasing power law relationship. This finding agrees with existing theory and experimental results. However, a minimum inlet solid loading ratio exists. When the air velocity is greater than the corresponding air velocity of the minimum solid loading ratio, the solid loading ratio exhibits an increasing trend in power law. We also found that when the inlet conveying pressure increased, the critical solid mass flow rate required for blockage, the inlet solid loading ratio, and the minimum inlet solid loading ratio increased.     

An investigation of a model of the flow pattern transition mechanism in relation to the identification of annular flow of R134a in a vertical tube using various void fraction models and flow regime maps

In the present study, new experimental data are presented for literature on the prediction of film thickness and identification of flow regime during the co-current downward condensation in a vertical smooth copper tube having an inner diameter of 8.1 mm and a length of 500 mm. R134a and water are used as working fluids in the tube side and annular side of a double tube heat exchanger, respectively. Condensation experiments are done at mass fluxes of 300 and 515 kg m^?2 s^?1. The condensing temperatures are between 40 and 50 °C; heat fluxes are between 12.65 and 66.61 kW m^?2. The average experimental heat transfer coefficient of the refrigerant HFC-134a is calculated by applying an energy balance based on the energy transferred from the test section. A mathematical model by Barnea et al. based on the momentum balance of liquid and vapor phases is used to determine the condensation film thickness of R134a. The comparative film thickness values are determined indirectly using relevant measured data together with various void fraction models and correlations reported in the open literature. The effects of heat flux, mass flux, and condensation temperature on the film thickness and condensation heat transfer coefficient are also discussed for the laminar and turbulent flow conditions. There is a good agreement between the film thickness results obtained from the theoretical model and those obtained from six of 35 void fraction models in the high mass flux region of R134a. In spite of their different valid conditions, six well-known flow regime maps from the literature are found to be predictive for the annular flow conditions in the test tube in spite of their different operating conditions.     

Numerical investigation of particle deposition inside aero-shielded solar cyclone reactor: A promising solution for reactor clogging

Solar cracking of methane is considered to be an attractive option due to its CO₂ free hydrogen production process. Carbon particle deposition on the reactor window, walls and exit is a major obstacle to achieve continuous operation of methane cracking solar reactors. As a solution to this problem a novel “aero-shielded solar cyclone reactor” was created. In this present study the prediction of particle deposition at various locations for the aero-shielded reactor is numerically investigated by a Lagrangian particle dispersion model. A detailed three dimensional computational fluid dynamic (CFD) analysis for carbon deposition at the reactor window, walls and exit is presented using a Discrete Phase Model (DPM). The flow field is based on a RNG k–ε model and species transport with methane as the main flow and argon/ hydrogen as window and wall screening fluid. Flow behavior and particle deposition have been observed with the variation of main flow rates from 10–20 L/min and with carbon particle mass flow rate of 7 × 10⁻⁶ and 1.75 × 10⁻⁵ kg/s. In this study the window and wall screening flow rates have been considered to be 1 L/min and 10 L/min by employing either argon or hydrogen. Also, to study the effect of particle size simulations have also been carried out (i) with a variation of particle diameter with a size distribution of 0.5–234 μm and (ii) by taking 40 μm mono sized particles which is the mean value for the considered size distribution. Results show that by appropriately selecting the above parameters, the concept of the aero-shielded reactor can be an attractive option to resolve the problem of carbon deposition at the window, walls and exit of the reactor.     

CFD simulation of a gas–solid fluidized bed with two vertical jets

A computational fluid dynamics (CFD) model is used to investigate the hydrodynamics of a gas–solid fluidized bed with two vertical jets. Sand particles with a density of 2660 kg/m³ and a diameter of 5.0 × 10^?4 m are employed as the solid phase. Numerical computation is carried out in a 0.57 m × 1.00 m two-dimensional bed using a commercial CFD code, CFX 4.4, together with user-defined Fortran subroutines. The applicability of the CFD model is validated by predicting the bed pressure drop in a bubbling fluidized bed, and the jet detachment time and equivalent bubble diameter in a fluidized bed with a single jet. Subsequently, the model is used to explore the hydrodynamics of two vertical jets in a fluidized bed. The computational results reveal three flow patterns, isolated, merged and transitional jets, depending on the nozzle separation distance and jet gas velocity and influencing significantly the solid circulation pattern. The jet penetration depth is found to increase with increasing jet gas velocity, and can be predicted reasonably well by the correlations of Hong et al. (2003) for isolated jets and of Yang and Keairns (1979) for interacting jets.