• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1262-1272
• /
• 2003

A two-dimensional direct numerical simulation is performed to investigate the dynamic behaviors of a single vortex in counter reacting and non-reacting flow field. A predictor-corrector-type numerical scheme with a low Mach number approximation is used in this simulation. A 16-step augmented reduced mechanism is adopted to treat the chemical reaction. The budget of the vorticity transport equation is examined to reveal a mechanism leading to the formation, destruction and transport of a single vortex according to the direction of vortex generation in reacting and non-reacting flows. The results show that air-side vortex has more larger strength than that of fuel-side vortex in both non-reacting and reacting flows. In reacting flow, the vortex is more dissipated than that in non-reacting flow as the vortex approach the flame. The total circulation in reacting flow, however, is larger than that in non-reacting flow because the convection transport of vorticity becomes much large by the increased velocity near the flame region. It is also found that the stretching and the convection terms mainly generate vorticity in non-reacting and reacting flows. The baroclinic torque term generates vorticity, while the viscous and the volumetric expansion terms attenuate vorticity in reacting flow. Furthermore, the contribution of volumetric expansion term on total circulation for air-side vortex is much larger than that of fuel-side vortex. It is also estimated that the difference of total circulation near stagnation plane according to the direction of vortex generation mainly attributes to the convection term.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.35-40
• /
• 2002

The dynamic behaviors of the single vortex in a reacting and non-reacting methane-air jet flow were investigated numerically. The numerical method was based on a predictor-corrector for low Mach number flow A two-step global reaction mechanism was adopted as a combustion model. After fuel and air were developed entirely in computational domain, the single vortex was generated by an axisymmetric jet that was impulsed to emit a cold fuel. Through comparisons of single vortex in reacting and non-reacting jet flow, it was found that global dynamic behaviors and the mechanisms leading to the formation, transport processes of vortex ring were influenced significantly by heat release from reaction. In addition, the interaction between a single vortex and flame bulge generated by buoyance effect in a reacting jet flow was found.

• Journal of Mechanical Science and Technology
• /
• v.18 no.3
• /
• pp.499-512
• /
• 2004

The effects of swirl intensity on non-reacting and reacting flow characteristics in a flat flame burner (FFB) with four types of swirlers were investigated. Experiments using the PIV method were conducted for several flow conditions with four swirl numbers of 0, 0.26, 0.6 and 1.24 in non-reacting flow. The results show that the strong swirling flow causes a recirculation, which has the toroidal structures, and spreads above the burner exit plane. Reacting flow characteristics such as temperature and the NO concentrations were also investigated in comparison with non-reacting flow characteristics. The mean flame temperature was measured as the function of radial distance, and the results show that the strong swirl intensity causes the mean temperature distributions to be uniform. However the mean temperature distributions at the swirl number of 0 show the typical distribution of long flames. NO concentration measurements show that the central toroidal recirculation zone caused by the strong swirl intensity results in much greater reduction in NO emissions, compared to the non-swirl condition. For classification into the flame structure interiorly, the turbulence Reynolds number and the Damkohler number have been examined at each condition. The interrelation between reacting and non-reacting flows shows that flame structures with swirl intensity belong to a wrinkled laminar-flame regime.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.449-452
• /
• 2008

Flow difference between reacting and non-reacting case in a swirl stabilized annular combustor is investigated using 3D Large Eddy Simulation with flamelet turbulent combustion model. The combustor of concern is the LM6000, lean premixed dry low-NOx annular combustor, developed by GEAE. Boundary conditions are based on experimental data. Heat release as a result of combustion put the dilatation of density in primary combustion zone highly increased so that the main swirl stream behind of a swirl cup stretched further downstream than that of non-reacting case. The oval shape of core flow in cross-section to flow direction, which clearly observed in non-reacting case, tends to be circle, and small vorticities in wide range in non-reacting case disappears, but the size of iso-vorticity increase in reacting case.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.99-105
• /
• 1997

A finite-difference method based on conservative supra characteristic method type upwind flux difference splitting has been developed to study the nonequilibrium chemically reacting inviscid flow. For nonequilibrium air, NS-1 species equations were strongly coupled with flowfield equations through convection and species production terms. Inviscid nonequilibrium chemically reacting air mixture flows over Blunt-body were solved to demonstrate the capability of the current method. At low altitude flight conditions the nonequilibrium air models predicted almost the same temperature, density and pressure behind the shock as equilibrium flow: however, at high altitudes they showed substantial differences due to nonequilibrium chemistry effect. The new nonequilibrium chemically reacting upwind flux difference splitting mettled can be extended to viscous flow and multi-dimensional flow conditions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.10
• /
• pp.28-37
• /
• 2004

A two-dimensional direct numerical simulation was performed to investigate the evolution and vortical structure of a single vortex in reacting and non-reacting jet flow fields. A predictor-corrector-type numerical scheme with a low Mach number approximation was used, and a two-step global reaction mechanism was adopted as the combustion model. Through the comparisons of single vortex behaviors in reacting and non-reacting jet flow fields, it was found that the evolution characteristics and vortical structure of the single vortex were significantly influenced by a outer vortex that was generated from the buoyance effect as well as the chemical heat release. Furthermore, it was also identified that the differences of the vortical structure in reacting and non-reacting jet flow fields were mainly attributed to the thermal expansion, Baroclinic torque and buoyance effect.

• 한국연소학회:학술대회논문집
• /
• 2002.11a
• /
• pp.49-56
• /
• 2002

A two-dimensional direct numerical simulation is performed to investigate the formation characteristics of a single vortex interacting with $CH_4/N_2$-Air counterflow nonpremixed flame. The numerical method was based on a predictor-corrector scheme for a low Mach number flow. The detailed transport properties and a 16-step augmented reduced mechanism are adopted in this calculation. The budgets of the vorticity transport equation arc examined to reveal the mechanisms leading to the formation, evolution and dissipation of a single vortex interacting with counterflow nonpremixed flame. It is found that the stretching term, which depends on the azimuthal component of vorticity, and radial velocity, mainly generates vortieitv in non-reacting and reacting flows. The viscous and baroclinic torque term destroy the vorticity in non-reacting flow. In addition, the baroclinic torque term due to density and pressure gradient generates vorticity, while viscous and the volumetric expansion terms due to density gradient destroy vorticity in reacting flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.2
• /
• pp.252-259
• /
• 2001

Computational study has been performed to observe the flow characteristics of combustion chamber for geometrical modification in municipal solid waste incinerator. A series of geometrical modification has been carried out as an attempt to reduce the size of recirculation zone, to obtain uniform flow field in the secondary combustion chamber and to improve the mixing of combustion gas. Two dimensional non-reacting turbulent flow has been studied as the first step to get such goals and the result of design optimization is presented. In addition, three dimensional non-reacting and reacting flow analyses were performed to verify the validity of two dimensional approach.

• 한국연소학회:학술대회논문집
• /
• 2013.06a
• /
• pp.123-126
• /
• 2013

Combustor development requires high fidelity simulation capable of predicting recirculation zone (RZ), temperature field, and pollutant emission. Swirling flow is widely used in combustor for its benefits in efficient mixing and flame stabilization by RZ. Large eddy simulation (LES) is used to calculate swirling flow in an expanding pipe [1], and shows higher accuracy than RANS. Reactive flow modeling using LES and flamelet model is validated with experiments by Barlow et al. [4] and Masri et al. [3]. Finally, heat transfer simulation of Samsung Techwin's combustor liner is presented.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.