• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.783-788
• /
• 2001

Choosing the most suitable mathematical model and relating this to turbulent tangential tensions model are very important in the investigations of turbulent two-phase flow. This paper considers two-fluid scheme. According to it, two phases have their own densities, velocities, and temperatures at any spatial point and at any moment. The equations of motion and heat transfer for each phase are linked with the forces of interaction between two phases. These forces are considered as predominant for the flow. As a closure in the system of motion equations, one modification of $K - {\epsilon}$ turbulent model is worked out. The modification uses two equations for turbulent kinetic energy of the phases and one - for the turbulent energy loss of main phase. This model can be set as a $K_g - K_p -{\epsilon}$ model. The modified model has been tested for both a two-phase non-isothermal flat jet and axially symmetrical jet. The numerical results are compared with the reference data revealing a good agreement between them.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1399-1408
• /
• 2000

A 2-dimensional temperature field measurement technique using PLIF (Planar Laser Induced Fluorescence) was developed and it was applied to an axisymmetric buoyant jet. Rhodamine B was used as a fluorescent dye. Laser light sheet illuminated a two-dimensional cross section of the jet. The intensity variations of LIF signal from Rhodamine B molecules scattered by the laser light were captured with an optical filter and a CCD camera. The spatial variations of temperature field of buoyant jet were derived using the calibration data between the LIF signal and real temperature. The measured results show that the turbulent jet is more efficient in mixing compared to the transition and laminar jet flows. As the initial flow condition varies from laminar to turbulent flow, the entrainment from ambient fluid increases and temperature decay along the jet center axis becomes larger. In addition to the mean temperature field, the spatial distributions of temperature fluctuations were measured by the PLIF technique and the result shows the shear layer development from the jet nozzle exit.

• Journal of Energy Engineering
• /
• v.10 no.3
• /
• pp.272-277
• /
• 2001

An experimental study of jet impinging the non-isothermal heating surface with linear temperature gradient is conducted with the presentation of the turbulent flow characteristics and the heat transfer rate, represented by the Nusselt number. The jet Reynolds number ranges from 15,000 to 30,000, the temperature gradient of the plate is 2~4.2$^{\circ}C$/cm and the dimensionless nozzle to plate distance (H/D) is from 2 to 10. The results show that the peak of heat transfer rate occurs at the stagnation point, and the heat transfer rate decreases as the radial distance from the stagnation point increases. A remarkable feature of the heat transfer rate is the existence of the second peak. This is due to the turbulent development of the wall jet. Maximum heat transfer rate occurs when the axial distance from the nozzle to nozzle diameter (H/D) is 6 or 8. The heat transfer rate can be correlated as a power function of Prandtl number, Reynolds number, the dimensionless nozzle to plate distance (H/D) and temperature gradient (dT/dr). It has been found that the heat transfer rate increases with increasing turbulent intensity. The wall jet is influenced by temperature gradient and the effect becomes more important at higher radii.

• Journal of Power System Engineering
• /
• v.10 no.4
• /
• pp.17-24
• /
• 2006

The turbulent flow field characteristics of a gun-type gas burner with and without a duct were investigated under the isothermal condition of non-combustion. Vectors and mean velocities were measured by hot-wire anemometer system with an X-type hot-wire probe in this paper. The turbulent flow field with a duct seems to cause a counter-clockwise recirculation flow from downstream to upstream due to the unbalance of static pressure between a main jet flow and a duct wall. Moreover, the recirculation flow seems to expand the main jet flow to the radial and to shorten it to the axial. Therefore, the turbulent flow field with a duct increases a radial momentum but decreases a axial momentum. As a result, an axial mean velocity component with a duct above the downstream range of about X/R=1.5 forms a smaller magnitude than that without a duct in the inner part of a burner, but it shows the opposite trend in the outer part.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.729-732
• /
• 2010

The present study has been motivated by the development of a reliable numerical methodology for simulation of kerosene/LOx coaxial swirl injectors. To deal with thermodynamic non-ideality and anomalies of transport properties pronounced at supercritical pressures, a set of subroutine libraries has been constructed based on the cubic equations of state, and applied to an existing flamelet analysis code. For computational efficiency, two-dimensional axisymmetric RANS formulation with swirl was adopted and validated successfully against an isothermal coaxial swirling jet. For the actual problem with high pressure combustion, however, numerical results show that the RANS models yield excessive production of turbulence probably due to high density gradient magnitude in the vicinity of mixing layer of swirling film flow, and imply strongly further improvement of the turbulence models.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2386-2396
• /
• 1996

A numerical simulation has been performed for isothermal and reacting flows in an exisymmetric, bluff-body research combustor. The present formulation is based on the density-weighted averaged Navier-Stokes equations together with a k-epsilon. turbulence model and a modified eddy-breakup combustion model. The PISO algorithm is employed for solution of thel Navier-Stokes system. Comparison between measurements and predictions are made for a centerline axial velocities, location of stagnation points, strength of recirculation zone, and temperature profile. Even though the numerical simulation gives acceptable agreement with experimental data in many respects, the present model is defictient in predicting the recoveryt rate of a central near-wake region, the non-isotropic turbulence effects, and variation of turbulent Schmidt number. Several possible explanations for these discrepancies have been discussed.