scholarly journals One Dimensional Modeling of Jet Diffusion Flame

Author(s):  
Angan Sengupta ◽  
A.K. Gupta ◽  
I.M. Mishra ◽  
S. Suresh

This paper reports on the determination of the flame height of a flare system using theoretical approach based on the laws of conservation of mass, momentum and energy. The set of ordinary differential equations at steady state conditions are solved numerically by fourth order Runge–Kutta method. The extent of reaction between the fuel and the entrained air has been studied by introducing the reaction mixing efficiency parameter, as the reaction rate is fixed by local entrainment rate. The reaction mixing efficiency parameter is a key measure to determine the height of the flame and its variation with the source velocity is limited by the flame width and the maximum vertical flame velocity at the tip of the flaming region. The variation of different parameters as vertical flame velocity, flame geometry and flame temperature with flame height are shown in plots. It is found that the flame geometry undergoes an initial necking up to a certain height, followed by an increase in its spread thereafter. The flame geometry and the flame dynamics depend exclusively on the burner design and the stack exit velocity. The flame height to burner diameter ratio of the jet diffusion flame is found to vary linearly with Froude number with a constant slope in the logarithmic plot.

2020 ◽  
Vol 2020 ◽  
pp. 1-10 ◽  
Author(s):  
Yuchen Ya ◽  
Xiaokang Nie ◽  
Licheng Peng ◽  
Longkai Xiang ◽  
Jialong Hu ◽  
...  

Laminar diffusion flame was used to study the effect of ethanol on n-heptane flame in terms of the morphology and microstructure of soot under atomization combustion. For the same carbon mass flux at the outlet of the burner, the ratio of ethanol doping in n-heptane was changed, and the soot was collected from the axial positions of the flame at different heights using the thermophoresis probe method. The results showed that the flame height increased significantly with the increasing ratio of ethanol doping. When the ratio of ethanol and n-heptane (CE/CN) was 1.5, the flame height increased by 10 mm compared with that of pure n-heptane flame. Besides, the temperature in the center of the flame decreased with the increasing ratio of ethanol doping, but the temperature in the low position was higher than that in the pure n-heptane flame, and the temperature in the high position was lower than that in the pure n-heptane flame. However, the flame temperature was the highest when the proportion of ethanol in the mixture was greater than that of n-heptane. The temperature at the flame center decreased with the increasing ratio of ethanol doping, while the temperature at the flame edge increased with the ratio. The primary particle size of soot (soot size hereafter) in all working conditions increased with the increase of flame height, which was in line with the general growth law of soot. Moreover, the soot size at the same height decreased with the increasing ratio of ethanol doping, and this trend was most obvious at the flame height of 20 mm and 30 mm. Compared with pure n-heptane, when CE/CN was 1.5, the soot size at 20 mm and 30 mm decreased by an average of 34.83%, indicating that ethanol could inhibit the surface growth of soot particle. Furthermore, the density of soot particles collected by a single copper net decreased significantly, indicating that ethanol could reduce the production amount of soot.


Author(s):  
Akira Yamaguchi ◽  
Yuji Tajima

A numerical methodology has been developed to solve a non-premixed diffusion flame under natural convection and applied to the calculation of the liquid sodium pool combustion experiment at various pool temperatures and oxygen concentrations. The computations reproduce the experimental observations concerning burning rate, flame temperature and flame height, consistently. The aerosol release rates are also in good agreement with the measurement. Dominant mechanisms of the mass and heat transfer are identified through the numerical simulation. An intrinsic feature is that the liquid sodium pool combustion is self-limited and negative feedback mechanism is at work Interaction among the thermal-hydraulics, chemical reaction and aerosol behavior plays an important role in the phenomena and it has been successfully analyzed by the numerical simulation. The present method can be used to understand sodium combustion phenomena and applied to the modeling of sodium pool combustion for safety analyses of liquid metal fast reactors. The numerical simulation is a useful tool because it can easily employ various conditions by changing parameters.


Fluids ◽  
2021 ◽  
Vol 6 (1) ◽  
pp. 27
Author(s):  
J. Barry Greenberg ◽  
David Katoshevski

A theoretical investigation of the influence of a standing wave flow-field on the dynamics of a laminar two-dimensional spray diffusion flame is presented for the first time. The mathematical analysis permits mild slip between the droplets and their host surroundings. For the liquid phase, the use of a small Stokes number as the perturbation parameater enables a solution of the governing equations to be developed. Influence of the standing wave flow-field on droplet grouping is described by a specially constructed modification of the vaporization Damkohler number. Instantaneous flame front shapes are found via a solution for the usual Schwab–Zeldovitch parameter. Numerical results obtained from the analytical solution uncover the strong bearing that droplet grouping, induced by the standing wave flow-field, can have on flame height, shape, and type (over- or under-ventilated) and on the existence of multiple flame fronts.


2001 ◽  
Vol 428 ◽  
pp. 349-386 ◽  
Author(s):  
E. J. STRANG ◽  
H. J. S. FERNANDO

The results of a laboratory experiment designed to study turbulent entrainment at sheared density interfaces are described. A stratified shear layer, across which a velocity difference ΔU and buoyancy difference Δb is imposed, separates a lighter upper turbulent layer of depth D from a quiescent, deep lower layer which is either homogeneous (two-layer case) or linearly stratified with a buoyancy frequency N (linearly stratified case). In the parameter ranges investigated the flow is mainly determined by two parameters: the bulk Richardson number RiB = ΔbD/ΔU2 and the frequency ratio fN = ND=ΔU.When RiB > 1.5, there is a growing significance of buoyancy effects upon the entrainment process; it is observed that interfacial instabilities locally mix heavy and light fluid layers, and thus facilitate the less energetic mixed-layer turbulent eddies in scouring the interface and lifting partially mixed fluid. The nature of the instability is dependent on RiB, or a related parameter, the local gradient Richardson number Rig = N2L/ (∂u/∂z)2, where NL is the local buoyancy frequency, u is the local streamwise velocity and z is the vertical coordinate. The transition from the Kelvin–Helmholtz (K-H) instability dominated regime to a second shear instability, namely growing Hölmböe waves, occurs through a transitional regime 3.2 < RiB < 5.8. The K-H activity completely subsided beyond RiB ∼ 5 or Rig ∼ 1. The transition period 3.2 < RiB < 5 was characterized by the presence of both K-H billows and wave-like features, interacting with each other while breaking and causing intense mixing. The flux Richardson number Rif or the mixing efficiency peaked during this transition period, with a maximum of Rif ∼ 0.4 at RiB ∼ 5 or Rig ∼ 1. The interface at 5 < RiB < 5.8 was dominated by ‘asymmetric’ interfacial waves, which gradually transitioned to (symmetric) Hölmböe waves at RiB > 5:8.Laser-induced fluorescence measurements of both the interfacial buoyancy flux and the entrainment rate showed a large disparity (as large as 50%) between the two-layer and the linearly stratified cases in the range 1.5 < RiB < 5. In particular, the buoyancy flux (and the entrainment rate) was higher when internal waves were not permitted to propagate into the deep layer, in which case more energy was available for interfacial mixing. When the lower layer was linearly stratified, the internal waves appeared to be excited by an ‘interfacial swelling’ phenomenon, characterized by the recurrence of groups or packets of K-H billows, their degeneration into turbulence and subsequent mixing, interfacial thickening and scouring of the thickened interface by turbulent eddies.Estimation of the turbulent kinetic energy (TKE) budget in the interfacial zone for the two-layer case based on the parameter α, where α = (−B + ε)/P, indicated an approximate balance (α ∼ 1) between the shear production P, buoyancy flux B and the dissipation rate ε, except in the range RiB < 5 where K-H driven mixing was active.


2019 ◽  
Vol 9 (19) ◽  
pp. 3989 ◽  
Author(s):  
Cheng Wang ◽  
Anthony Chun Yin Yuen ◽  
Qing Nian Chan ◽  
Timothy Bo Yuan Chen ◽  
Qian Chen ◽  
...  

This paper numerically examines the characterisation of fire whirl formulated under various entrainment conditions in an enclosed configuration. The numerical framework, integrating large eddy simulation and detailed chemistry, is constructed to assess the whirling flame behaviours. The proposed model constraints the convoluted coupling effects, e.g., the interrelation between combustion, flow dynamics and radiative feedback, thus focuses on assessing the impact on flame structure and flow behaviour solely attribute to the eddy-generation mechanisms. The baseline model is validated well against the experimental data. The data of the comparison case, with the introduction of additional flow channelling slit, is subsequently generated for comparison. The result suggests that, with the intensified circulation, the generated fire whirl increased by 9.42 % in peak flame temperature, 84.38 % in visible flame height, 6.81 % in axial velocity, and 46.14 % in velocity dominant region. The fire whirl core radius of the comparison case was well constrained within all monitored heights, whereas that of the baseline tended to disperse at 0.5   m height-above-burner. This study demonstrates that amplified eddy generation via the additional flow channelling slit enhances the mixing of all reactant species and intensifies the combustion process, resulting in an elongated and converging whirling core of the reacting flow.


2019 ◽  
Vol 141 (9) ◽  
Author(s):  
Mehdi Bidabadi ◽  
Saman Hosseinzadeh ◽  
Sadegh Sadeghi ◽  
Mostafa Setareh

Due to perspective of biomass usage as a viable source of energy, this paper suggests a potential theoretical approach for studying multiregion nonadiabatic premixed flames with counterflow design crossing through the mixture of air (oxidizer) and lycopodium particles (biofuel). In this research, convective and radiative heat losses are analytically described. Due to the properties of lycopodium, roles of drying and vaporization are included so that the flame structure is created from preheating, drying, vaporization, reaction, and postflame regions. To follow temperature profile and mass fraction of the biofuel in solid and gaseous phases, dimensionalized and nondimensionalized forms of mass and energy balances are expressed. To ensure the continuity and calculate the positions of drying, vaporization, and flame fronts, interface matching conditions are derived employing matlab and mathematica software. For validation purpose, results for temperature profile is compared with those provided in a previous research study and an appropriate is observed under the same conditions. Finally, changes in flame velocity, flame temperature, solid and gaseous fuel mass fractions, and particle size with position measured from the position of stagnation plane, strain rate, and heat transfer coefficient in the presence/absence of losses are evaluated.


2013 ◽  
Vol 70 (1) ◽  
pp. 266-277 ◽  
Author(s):  
Kyongmin Yeo ◽  
David M. Romps

Abstract Lagrangian particle tracking is used in a large-eddy simulation to study an individual cumulus congestus. This allows for the direct measurement of the convective entrainment rate and of the residence times of entrained parcels within the cloud. The entrainment rate obtained by Lagrangian direct measurement is found to be higher than that obtained using the recently introduced method of Eulerian direct measurement. This discrepancy is explained by the fast recirculation of air in and out of cloudy updrafts, which Eulerian direct measurement is unable to resolve. By filtering these fast recirculations, the Lagrangian calculation produces a result in very good agreement with the Eulerian calculation. The Lagrangian method can also quantify some aspects of entrainment that cannot be probed with Eulerian methods. For instance, it is found that more than half of the air that is entrained by the cloud during its lifetime is air that was previously detrained by the cloud. Nevertheless, the cloud is highly diluted by entrained air: for cloudy air above 2 km, its mean height of origin is well above the cloud base. This paints a picture of a cloud that rapidly entrains both environmental air and its own detritus.


2005 ◽  
Author(s):  
Mohsen M. Abou-Ellail ◽  
Karam R. Beshay ◽  
David R. Halka

The present work is a numerical simulation of the, piloted, non-premixed, methane–air flame structure in a new mathematical imaging domain. This imaging space has the mixture fraction of diffusion flame Z1 and mixture fraction of pilot flame Z2 as independent coordinates to replace the usual physical space coordinates. The predications are based on the solution of two–dimensional set of transformed second order partial differential conservation equations describing the mass fractions of O2, CH4, CO2, CO, H2O, H2 and sensible enthalpy of the combustion products which are rigorously derived and solved numerically. A three–step chemical kinetic mechanism is adopted. This was deduced in a systematic way from a detailed chemical kinetic mechanism by Peters (1985). The rates for the three reaction steps are related to the rates of the elementary reactions of the full reaction mechanism. The interaction of the pilot flame with the non-premixed flame and the resulting modifications to the structure and chemical kinetics of the flame are studied numerically for different values of the scalar dissipation rate tensor. The dissipation rate tensor represents the flame stretching along Z1, the main mixture fraction, and in the perpendicular direction, along Z2, the pilot mixture fraction. The computed flame temperature contours are plotted in the Z1-Z2 plane for fixed values of the dissipation rate along Z1 and Z2.These temperature contours show that the flame will become unstable when the dissipate rates along Z1 and Z2 increase, simultaneously, to the limiting value for complete flame extinction of 45 s−1. However, the diffusion flame will extinguish for dissipate rates less than 20 1/s, if unpiloted. It is also noticed that the flame will remain stable if the dissipation rate along Z2 is increased to the limiting value, while the dissipation rate, along Z2, remains constant at a value less than 30 s−1.


Sign in / Sign up

Export Citation Format

Share Document