• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.4
• /
• pp.780-788
• /
• 1989

An experimental study has been performed to investigate the process from nonequilibrium state to equilibrium state in redeveloping turbulent boundary layer beyond separation-reattachment using pitot tube and hot-wire anemometer. The model sued in the experiment has the form of a backward facing step which is assembled by a two-dimensional 4:1 half elipse and a plate. Measurements are carried out up to a distance of about 50 step height downstream of the step, where the reattachment observed at about x/h=6.5. The profiles of the shape factor H the Clauser parameter G and the coefficient of friction $C^{f}$ exhibited the characteristics similar to those of the equilibrium turbulent boundary layer from x/h=25, and the profiles of the trubulent quantities did from x/h=35. However, the wake region of the boundary layer does not seem to recover the equilibrium turbulent boundary layer even at x/h=50. By considering the distributions of the intermittency factor it has been noted that the turbulence structure changes gradually from a mixing layer to a turbulent boundary layer along downstream direction after reattachment. This becomes clearer as we analyse the one-dimensional energy spectra and the dissipation energy spectra which are measured and caculated at various downstream positions after the backward facing step.p.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.1
• /
• pp.70-77
• /
• 1986

Turbulent mixing field with recirculating flow which is formed by injecting gaseous fuel on the main air stream is solved numerically by a finite difference method. The turbulence model for obtaining transport properties was k-.epsilon. model, which was obtained from turbulent kinetic energy and its dissipation rate. Considering the effects of streamline curvature, modified k-.epsilon model was used. Generally, Modified k-.epsilon. model makes better predictions than standard model, and from this result, it is recognized that standard model has deficiency when applied to turbulent recirculating flows, and that modified k-.epsilon. model takes into account of streamline curvature effects properly. Meanwhile, A more study will be necessary to find the reason why large differences between predicted and experimental turbulent kinetic energy exist.

• 한국연소학회:학술대회논문집
• /
• 2012.04a
• /
• pp.289-291
• /
• 2012

The present study numerically investigate detailed flame structure of the Syngas diffusion flames. In order to realistically represent the turbulence-chemistry interaction and the spatial inhomogeneity of scalar dissipation rate, the Eulerian Particle Flamelet Model(EPFM) with multiple flamelets has been applied to simulate the combustion processes and NOx formation in the syngas turbulent nonpremixed flames. And level-set approach is also utilized to account for the partially premixing effect at fuel and oxidizer injector in KEPRI nonpremixed combustor. Based on numerical results, the detailed discussion has been made for the precise structure and NOx formation characteristics of the turbulent syngas nonpremixed flames.

• Journal of Advanced Marine Engineering and Technology
• /
• v.23 no.2
• /
• pp.216-225
• /
• 1999

Four turbulent $k-{varepsilon}$models(i.e standard model modified models with streamline curvature modification and/or preferential dissipation modification) are applied in order to analyze the tur-bulent flow of wall-attaching offset jet. The upwind numerical scheme was adopted in the present analyses. The streamline curvature modification results in slightly better prediction while the preferential dissipation modification does not. The obtained analytic results will be used as refer-ences for further study regarding Reynolds stress model. In addition this paper introduced a method of increasing nozzle outlet velocity gradually for numercal convergence. Even though the method was simple it was efficient in view of convergent speed CPU running time computer memory storage programming etc.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.1
• /
• pp.12-19
• /
• 2001

Study of turbulent mixing layers has been a popular subject from the point of view of both practical application and phenomenological importance in engineering field. Turbulent mixing layers can be applied in many fields where rapid transition to turbulence is desirable in order to prevent boundary layer separation or to enhance mixing. The ability to control mixing, structure and growth of the shear flow would obviously have a considerable impact on many engineering applications. In addition to practical applications, free shear flows are one of the simplest flows to understand the fundamental mechanism in the transition process to turbulence. After the discovery of large-scale vortical structure in free shear flows many researchers have investigated the physical mechanism of generation and dissipation processes of the vortical structure. This study investigated the role of the large-scale vortical structures in the turbulent mixing layer using LES(Large-Eddy Simulation). The result shows that the pairing interaction of the vortical structure plays an important role in the growth rate of a mixing layer. It is found that the turbulence quantities depend strongly on the velocity ratio. It is also found that the vorticity in the high-velocity-side can extract energy from the mean flow, while the vorticity in the low-velocity-side lose energy by the viscous dissipation. Finally the results suggest the guideline to obtain the desired flow by control of the velocity ratio.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.1
• /
• pp.137-151
• /
• 1988

Thurbulent thermal convection between two plates, bottom plate is at higher temperature $T_{h}$ and the upper plate is at lower temperature $T_{i}$ is numerically investigated. Model equations are abridged Reynolds stress equations; full Reynolds stress equations are simplified to yield algebraic relations in case of mean square velocity fluctuations in vertical and horizontal directions. Boundary conditions for turbulent kinetic energy k and mean square temperature variance .thera.$^{2}$oner bar at the plate surfaces are set to be zero and those of dissipation rate of turbulent kinetic energy .epsilon. and dissipation rate of mean square temperature variance .epsilon.$_{\theta}$ are assumed at first grid point nearest to the boundary surfaces, whose values are approximated by inviscid estimates. Results show that temperature profiles are in good agreement with experimental data except transition region, in which temperature is over-predicted. Such discrepancy becomes larger as the Rayleigh number becomes smaller. Nusselt numbers, which are calculated from the temperature gradients at the boundary surfaces, are also in good agreement with experimental data.a.a.

• International Journal of Ocean System Engineering
• /
• v.2 no.1
• /
• pp.50-56
• /
• 2012

This study investigates flow fields and energy dissipation due to regular wave interaction with a perforated vertical breakwater, through velocity data measurement in a two-dimensional wave tank. As the waves propagate through the perforated breakwater, the incoming wave energy is reflected back to the ocean, dissipated due to very turbulent flows near the perforations and inside the chamber, and transmitted through the perforations of the breakwater. This transmitted energy is further reduced due to the presence of the perforated back wall. Hence most of the energy is either reflected or dissipated in the vicinity of the structure, and only a small amount of the incoming wave energy is transmitted through the structure. In this study, particle image velocimetry (PIV) technique was employed to measure two-dimensional instantaneous velocity fields in the vicinity of the structure. Measured velocity data was treated statistically, and used to calculate mean flow fields, turbulence intensity and turbulent kinetic energy. For investigation of the flow pattern, time-averaged mean velocity fields were examined, and discussed using the cross-sections through slot and wall for comparison. Flow fields were obtained and compared for various cases with different regular wave conditions. In addition, turbulent kinetic energy was estimated as an approach to understand energy dissipation near the perforated breakwater. The turbulent kinetic energy was distributed against wave height and wave period to see the dependence on wave conditions.

• Journal of Mechanical Science and Technology
• /
• v.15 no.3
• /
• pp.392-402
• /
• 2001

The paper reports a multiple source modeling of low-Reynolds-number dissipation rate equation with aids of DNS data. The key features of the model are to satisfy the wall limiting conditions of the individual source terms in the exact dissipation rate equation using the wall damping functions. The wall damping functions are formulated in term of dimensionless dissipation length scale ι(sup)+(sub)D(≡ι(sub)D($\upsilon$$\xi$)(sup)1/4/$\upsilon$) and the invariants of small and large scale turbulence anisotropy tensors. $\alpha$(sub)ij(=$\mu$(sub)i$\mu$(sub)j/$\kappa$-2$\delta$(sub)ij/3) and e(sub)ij(=$\xi$(sub)ij/$\xi$-2$\delta$(sub)ij/3). The model constants are optimized with aids of DNS data in a plane channel flow. Adopting the dissipation length scale as a parameter of damping function, the applicabilities of $\kappa$-$\xi$ model are extended to the turbulent flow calculation of complex flow passages.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1995.10a
• /
• pp.66-70
• /
• 1995

The oscillatory boundary layer that develops when surface waves propagate over the sea bottom affects many flow-pendent phenomena in the coastal zone. Examples of such phenomena are wave energy dissipation due to bottom friction and the initiation and transport of sediment (Grant and Madsen 1986). In nature the boundary layer under waves will almost always be turbulent (Nielsen 1992). (omitted)