• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.821-827
• /
• 2003

Fully developed turbulent flow in a square duct is numerically predicted with two nonlinear low-Reynolds-number ${\kappa}-{\varepsilon}$ models. Typical predicted quantities such as axial and secondary velocities, turbulent kinetic energy and Reynolds stresses are compared in detail with each other. It is found that the nonlinear low-Reynolds-number ${\kappa}-{\varepsilon}$ model adopted in a commercial code is unable to predict accurately duct flows involving turbulence-driven secondary motion with the prediction level of secondary flows one order less than that of the experiment.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2816-2821
• /
• 2008

This paper investigates the flow field and organism concentration in a UV disinfection channel in which vertical ultraviolet lamps are arranged in a staggered configuration. Turbulence is described by low Reynolds number ${\kappa}-{\varepsilon}$ turbulence model and standard ${\kappa}-{\varepsilon}$ turbulence model, respectively. P-1 method has been employed to solve the radiative transfer equation. The obtained incident radiation is used to compute the inactivation term in the species equation. The CFD results are in good agreement with the existing experimental data for the UV channel. For the flow field, the low-Reynolds number ${\kappa}-{\varepsilon}$ model is superior to the standard ${\kappa}-{\varepsilon}$ model. The approach velocity has a significant effect on the disinfection efficiency. The organism concentration at the outlet decreases fast to a low inlet velocity.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.1
• /
• pp.304-320
• /
• 1996

A multiple production .epsilon. equation model was developed in the low Reynolds number $\kappa$-$\varepsilon$ model with the aids of DNS data. We derived the model theoretically and avoided the use of empirical correlations as much as possible in order for the model to have generality in the prediction of complex turbulent flow. Unavoidable model constants were, however, optimized with the aids of DNS data. All the production and dissipation models in the $\varepsilon$ equation were modified with damping functions to satisfy the wall limiting behavior. A new $f_{\mu}$ function, turbulent diffusion and pressure diffusion model for the k and .epsilon. equations were also proposed to satisfy the wall limiting behavior. By, computational investigation on the plane channel flows, we found that the multiple production model for .epsilon. equation could improve the near wall turbulence behavior compared with the standard production model without the complicated empirical modification. Satisfication of the wall limiting conditions for each turbulence model term was found to be most important for the accurate prediction of near wall turbulence behaviors.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.6 no.4
• /
• pp.325-336
• /
• 1994

A numerical simulation of velocity and temperature fields and Nusselt number distributions is performed by using the algebraic stress model (ASM) for the velocity profiles and low Reynolds number ${\kappa}-{\varepsilon}$ model and the algebraic heat flux model(AHFM) for turbulent heat transfer in a $180^{\circ}$ bend with a constant wall heat flux. In the low Reynolds number ${\kappa}-{\varepsilon}$ model, turbulent Prandtl number is modified by considering the streamline curvature effect and the non-equilibrium effect between turbulent kinetic energy production and dissipation rate. Every heat flux term presented in the transport equation of turbulent heat flux is reduced to algebraic expressions in a way similar to algebraic stress model. Also. in the wall region, low Reynods number algebraic heat flux model(AHFM) is applied.

• Journal of Mechanical Science and Technology
• /
• v.17 no.7
• /
• pp.1054-1062
• /
• 2003

The numerical simulations of flow and pollutant particle dispersion are described for two-dimensional bell shaped hills with various aspect ratios. The Reynolds-averaged incompressible Navier-Stokes equations with low Reynolds number $\kappa$-$\varepsilon$ turbulent model are used to simulate the flowfield. The gradient diffusion equation is used to solve the pollutant dispersion field. The code was validated by comparison of velocity, turbulent kinetic energy, Reynolds shear stress, speed-up ratio, and ground level concentration with experimental and numerical data. Good agreement has been achieved and it has been found that the pollutant dispersion pattern and ground level concentration have been strongly influenced by the hill shape and aspect ratio, as well as the location and height of the source.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.14 no.12
• /
• pp.1004-1013
• /
• 2002

Turbulent carbon dioxide flows and cooling heat transfers under supercritical state in a straight duct with a square cross-section are numerically analyzed employing low Reynolds number $\kappa-\varepsilon$ model and algebraic stress model. The flow is assumed to be (quasi-incompressible. Predicted Nusselt number and friction factor are compared with the experimental data, Blasius correlation for friction factor and Dittus-Boelter correlation for Nusselt number. Computational results for the Fanning's friction factor agree well with the all Rohsenow and Choi's correlation, Liou and Hwang's experimental data and Blasius correlation. The results obtained by algebraic stress model agree more with the Liou and Hwang's experimental data, while the results obtained by low Reynolds number $\kappa-\varepsilon$ model agree more with Blasius correlation. In the computation of Nusselt number, Dittus-Boelter correlation can not exactly fit the computational results. Therefore we propose the new correlation$Nu=0.053Re^{0.73}Pr^{0.4}$for the turbulent cooling heat transfer of carbon dioxide under supercritical state.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.321-328
• /
• 2009

Free surface impinging jet on a moving plate, which is applicable to cooling of hot metals in a steel-making process, is investigated numerically by solving the Navier-Stokes equations in the liquid and gas phases. The free surface of liquid-gas interface is tracked by a level-set method which is improved by incorporating the ghost fluid approach based on a sharp-interface representation. The method is further improved by employing a nonequilibrium $\kappa-\varepsilon$ turbulence model including the effect of low Reynolds number. The computations are made to investigate the effects of the nozzle pitch, moving velocity of plate and jet velocity on the interfacial motion and the associated flow and temperature fields.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2294-2303
• /
• 1993

This paper investigated the characteristics of the turbulent incompressible flow past the orifice ring in an axi-symmetric pipe. The flow field was the turbulent pulsatile flow for Reynolds number of $2{\times}10^{5}$ which was defined based on the maximum velocity and the pipe diameter at the inlet, with oscillating frequence $(f_{os})=1/4{\pi}$ which was considered as quasi-steady state frequence. In the present investigation, finite analytic method was used to solve the governing equations in Navier Stokes and turbulent transport formulations. Particularly at high Reynolds number and low oscillation frequency, the effects of orifice ring on the flow were numerically investigated. The separation zone behind the orifice ring during the acceleration phase was found to be decreased. However, during the deceleration phase, the separation behind the orifice ring for pulsatile flow continuously grow to a size even larger than that in steady flow. The pressure drop in steady flow was found to be constant and always positive while for pulsatile flow the pressure drop change with time. And large turbulent kinetic energy, dissipation rate were found to be located in the region where the flow passes through the orifics ring. The maximum turbulent kinetic energy, generally occurs along the shear layer where the velocity gradient is large.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.9-25
• /
• 1994

In the present study, the Reynolds-averaged Navier-Stokes equations, together with the equations of the $k-{\varepsilon}$ model of turbulence, were solved numerically in a general body-fitted coordinate system for three-dimensional turbulent flows around the six basic shapes of the magnetically levitated train(MAGLEV). The numerical computations were conducted on the MAGLEV model configurations to provide information on shapes of this type very near the elevated track at a constant Reynolds number of $1.48{\times}10^{6}$ based on the body length. The coordinate system was generated by numerically solving a set of Poisson equations. The convective transport equations were discretized using the finite-analytic scheme which employed analytic solutions of the locally-linearized equations. A time marching algorithm was employed to enable future extensions to be made to handle unsteady and fully-elliptic problems. The pressure-velocity coupling was treated with the SIMPLER-algorithm. Of particular interests were wall effect by the elevated track on the aerodynamic forces and flow characteristics of the six models calculated. The results indicated that the half-circle configuration with extended sides and with smooth curvature of sides was desirable because of the low aerodynamic forces and pitching moment. And it was found that the separation bubble was occured at wake region in near the elevated track.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2315-2328
• /
• 1993

The turbulence effect of particle deposition on a horizontal free-standing wafer in a vertical flow has been studied numerically by using the low-Reynolds-number k-.epsilon. turbulence model. For both the upper and lower surfaces of the wafer, predictions are made of the averaged particle deposition velocity and its radial distribution. Thus, it is now possible to obtain local information about the particle deposition on a free-standing wafer. The present result indicates that the particle deposition velocity on the lower surface of wafer is comparable to that on the upper one in the diffusion controlled deposition region in which the particle sizes are smaller than $0.1{\mu}m$. And it is found in this region that, compared to the laminar flow case, the averaged deposition velocity under the turbulent flow is about two times higher, and also that the local deposition velocity at the center of wafer is high equivalent to that the wafer edge.