• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.6
• /
• pp.2160-2170
• /
• 1991

Three-dimensional turbulent structures in the near field of elliptic jet were experimentally investigated by using a three-color, three-component Laser Doppler Velocimeter. The Reynolds number based on the nozzle exit velocity and nozzle equivalent diameter(De) was about 4*10$^{4}$. The turbulent characteristics of a sharp-edged elliptic nozzle with aspect ratio of 2 were analyzed along major and minor axis at X/De=2,3,5,7 and along the centerline up to X/De=14. Quantities measured at each point with the 3-D LDV system were three orthogonal velocity components, turbulent intensity, skewness, flatness, and Reynolds shear stress. The nondimensional mean velocities coincided well with the Schlichting's empirical curve with going downstream. Elliptic jet of AR=2 had two switching points at about X/De=2 and 16. The turbulent intensity along the minor axis was distributed widely than that along the major axis. In the near field, X/De<5, the Reynolds shear stresses of the inner part of the elliptic jet had negative value, which indicated the enhancement of entrainment toward the inner part.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.10
• /
• pp.1378-1386
• /
• 2002

In the present study, velocity profile and wall shear stress distributions of developing turbulent oscillatory flows in an oscillator connected to straight duct located in exit region of a curved duct was investigated experimentally. The experimental study for air flows was conducted to measure axial velocity profiles, shear stress distributions by using the Laser Doppler Velocimetry(LDV) system with the data acquisition and processing system of Rotating Machinery Resolver(R.M.R) and PHASE software. The results obtained from experimental studies are summarized as follows. The critical Reynolds number for a change from transitional oscillatory flow to turbulent flow was about 7500, in the 60region of dimensionless axial position which was considered as a fully developed flow region. The turbulent oscillatory flow, velocity profiles of the inflow period in the entrance region were gradually developed, but those of the outflow period were not changed nearly. Velocity profiles of inflow and outflow were shown as a symmetric form in a fully developed flow region. The wall shear stress distributions of turbulent oscillatory flow increase rapidly as the flow proceeds to downstream and flow was in good agreement with the theoretically.

• 한국연소학회:학술대회논문집
• /
• 2000.05a
• /
• pp.122-132
• /
• 2000

Computational fluid dynamic simulation is performed for the floating turbulent premixed flames stabilized in stagnation flows of Cho et al. [1] and Cheng and Shepherd [2]. They are both in the wrinkled flamelet regime far from the extinction limit with $u'/S^{0}_{L}$ less than unity. The turbulent flux is given in the first moment closure as a sum of the classical gradient flux due to turbulent motions and the countergradient flux due to thermal expansion. The parameter $N_{B}'s$ are greater than unity with the countergradient flux dominant over the gradient flux. The countergradient flux is assumed to be zero in $\bar{c}<0.05$. The flame surface density is modeled as a symmetric parabolic function with respect to $\bar{c}$. The product of the maximum flame surface density and the mean stretch factor is considered as a tuning constant to match the flame location. Good agreement is achieved with the measured $\tilde{w}$ and $\bar{c}$ profiles along the axis in both flames.

• Journal of the Korean Society for Railway
• /
• v.10 no.2 s.39
• /
• pp.224-230
• /
• 2007

The performance and applicability of FDS code is analyzed for flow simulation in railway tunnel. FDS has been built in NIST(USA) for simulation of fire-driven flow. RANS and DNS's results are compared with FDS's. AJL non-linear ${\kappa}-{\epsilon}$[7,8] model is employed to calculate the turbulent flow for RANS. DNS data by Moser et al.[9] are used to prove the FDS's applicability in the near wall region. Parallel plate is used for simplified model of railway tunnel. Geometrical variables are non-dimensionalized by the height (H) of parallel plate. The length of streamwise direction is 50H and the length of spanwise direction is 5H. Selected Re numbers are 10,667 for turbulent flow and 133 for laminar low. The characteristics of turbulent boundary layer are introduced. AJL model's predictions of turbulent boundary layer are well agreed with DNS data. However, the near wall turbulent boundary layer is not well resolved by FDS code. Slip conditions are imposed on the wall but wall functions based on log-law are not employed by FDS. The heavily dense grid distribution in the near wall region is necessary to get correct flow behavior in this region for FDS.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.2
• /
• pp.218-225
• /
• 1999

Flow visualization study has been conducted to simulate the turbulent dispersion behavior of a crossflow jet physically under the conditions of various thermal stratification in a wind tunnel. A smoke jet with the constant ratio of the jet to freestream velocity is injected normally to the cross flow of the thermally stratified wind tunnel(TSWT) for flow visualization. The typical natures of the smoke dispersion under different thermal stratifications such as neutral, weakly stable, strongly stable, weakly unstable, strongly unstable and inversion layer are successfully reproduced in the TSWT. The Instantaneous velocity and temperature fluctuations are measured by using a cold and hot-wire combination probe. The time averaged dispersion behaviors, the centerline trajectories, the spreading angles and the virtual origins of the cross jet are deduced from the edge detected images with respect to the stability parameter. All the general characteristics of the turbulent dispersion behavior reveal that the definitely different dispersion mechanisms are inherent in both stable and unstable conditions. It is conjectured that the turbulent statistics obtained in the various stability conditions quantitatively demonstrate the vertical scalar flux plays a key role in the turbulent dispersion behavior.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.111-117
• /
• 2008

This paper reports briefly on the computational results of a turbulent Rayleigh-Benard convection with the elliptic-blending second-moment closure (EBM). The primary emphasis of the study is placed on an investigation of accuracy and numerical stability of the elliptic-blending second-moment closure for the turbulent Rayleigh-Benard convection. The turbulent heat fluxes in this study are treated by the algebraic flux model with the temperature variance and molecular dissipation rate of turbulent heat flux. The model is applied to the prediction of the turbulent Rayleigh-Benard convection for Rayleigh numbers ranging from $Ra=2{\times}10^6$ to $Ra=10^9$, and the computed results are compared with the previous experimental correlations, T-RANS and LES results. The predicted cell-averaged Nusselt number follows the correlation by Peng et al.(2006) ($Nu=0.162Ra^{0.286}$) in the 'soft' convective turbulence region ($2{\times}10^6{\leq}Ra{\leq}4{\times}10^7$) and it follows the experimental correlation by Niemela et al. (2000) ($Nu=0.124Ra^{0.309}$) in the 'hard' convective tubulence region ($10^8{\leq}Ra{\leq}10^9$) within 5% accuracy. This results show that the elliptic-blending second-moment closure with an algebraic flux model predicts very accurately the Rayleigh Benard convection.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.359-366
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.7
• /
• pp.905-914
• /
• 2001

A large eddy simulation(LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model is employed and the calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of the inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing them with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structure behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent mixing process behind the flame holder is analyzed by visualizing the sectional views of vortical structure.

• Journal of Mechanical Science and Technology
• /
• v.15 no.7
• /
• pp.895-905
• /
• 2001

This paper presents vector fields, three dimensional mean velocities, turbulent intensities, turbulent kinetic energy and Reynolds shear stresses measured in the X-Y plane of the gas swirl burner with a cone type baffle plate by using an X-type hot-wire probe. This experiment is carried out at the flow rates of 350 and 450ℓ/min which are equivalent to the combustion air flow rate necessary to release 15,000 kcal/hr in a gas furnace. The results show that the maximum axial mean velocity component exists around the narrow slits situated radially on the edge of a burner. Therefore, there is some entrainment of ambient air in the outer region of a burner. The maximum values of turbulent intensities occur around the narrow slits and in front of a burner up to X/R=1.5. Moreover, the turbulent intensity components show a relatively large value in the inner region due to the flow diffusion and mixing processes between the inclined baffle plate and the swirl vane. Consequently, the combustion reaction is expected to occur actively near these regions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.4
• /
• pp.417-427
• /
• 1998

The effect of secondary flow on the heat transfer of a turbulent wake generated by a flat plate was experimentally investigated. The secondary flow was induced in a curved duct in which the flat plate wake generator was installed. All three components of turbulent heat flux were measured in the plane containing the mean radius of curvature of the curved duct. The results showed that mean temperature profiles deviate from the similarity of the straight wake because of the cold fluid transported from the free-stream. The half-width of the mean temperature profile increased rapidly by upwash motion of the secondary flow. The changes to turbulence structure caused by the secondary flow show more pronounced effect on heat transport than on momentum transport. This is because the response to the variation of flow conditions is delayed in temperature field. Negative production of the turbulent heat flux is observed in the inner wake region. From the conditional averaging, it has been found that the negative production of the turbulent heat flux is generated due to a mixing process between the hot and low momentum eddies occupied in the inner wake region and the cold and high momentum eddies in the potential region.