• 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 B
• /
• v.20 no.5
• /
• pp.1661-1670
• /
• 1996

Incompressible flow over a backward-facing step is computed by low Reynolds number turbulence models in order to compare with direct simulation results. In this study, selected low Reynolds number 1st and 2nd (Algebraic Stress Model : ASM) moment closure turbulence models are adopted and compared with each other. Each turbulence model predicts different flow characteristics, different re-attachment point, velocity profiles and Reynolds stress distribution etc. Results by .kappa.-.epsilon. turbulence models indicate that predicted re-attachment lengths are shorter than those by standard model. Turbulent intensity and eddy viscosity by low Reynolds number .kappa.-.epsilon. models are still greater than DNS results. The results by algebraic stress model (ASM) are more reasonable than those by .kappa.-.epsilon. models. The convective scheme is QUICK (Quadratic Upstream Interpolation for Convective Kinematics) and SIMPLE algorithm is adopted. Reynolds number based on step height and inlet free stream velocity is 5100.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.1
• /
• pp.212-218
• /
• 1995

Turbulent flow and heat transfer characteristics around staggered tube bundles were studied using a non-orthogonal boundary fitted coordinate system and the low Reynolds .kappa. - .epsilon. turbulence model suggested by Lam and Bremhorst. The predicted flow characteristics for two tube pitches and tube arrangement showed good agreement with the experimental data except the strongly curved region. The predicted Nusselt number was compared with measurements obtained in the staggered rough bundles and it revealed the similar trend to measurements, but the location of the maximum and minimum heat transfer differed somewhat from the measurements.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.272-275
• /
• 2003

A mixed analytical/numerical method is developed here to solve the low Reynolds number kepsilon turbulence model. In this method the advection-diffusion part is solved numerically, while the source terms are split into two parts: one part is solved analytically and the next is solved numerically.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.5
• /
• pp.685-697
• /
• 1998

Series of recent k-.epsilon. model modification have been carried out with the aid of DNS data to include the effect of near wall. Though these methods opened new way of turbulence modelings, newly developed turbulence models of its kind had yet shortcomings in prediction for the turbulent flows with various Reynolds numbers and various geometric conditions. As a remedy for these shortcomings, a new k-.epsilon. model proposed here by improving the dissipation rate equation and the damping function for eddy viscosity model. The new dissipation rate equation was modeled based on the energy spectrum and magnitude analysis. The damping function for eddy viscosity was also formulated on the ground of distribution of dissipation rate length scales near a wall and the DNS data. The new k-.epsilon. model was applied to the fully developed turbulent flows in a channel and a pipe with a wide range of Reynolds numbers. Prediction results showed that the present model represents properly the turbulence properties in all turbulent regions over a wide range of Reynolds numbers.

• Nuclear Engineering and Technology
• /
• v.36 no.3
• /
• pp.248-262
• /
• 2004

Turbulence models are separately assessed for a three dimensional thermal-hydraulic analysis of the integral reactor SMART. Seven models (mixing length, k-l, standard $k-{\epsilon},\;k-{\epsilon}-f{\mu},\;k-{\epsilon}-v2$, RRSM, and ERRSM) are investigated for flat plate channel flow, rotating channel flow, and square sectioned U-bend duct flow. The results of these models are compared to the DNS data and experiment data. The results are assessed in terms of many aspects such as economical efficiency, accuracy, theorization, and applicability. The standard $k-{\epsilon}$ model (high Reynolds model), the $k-{\epsilon}-v2$ model, and the ERRSM (low Reynolds models) are selected from the assessment results. The standard $k-{\epsilon}$ model using small grid numbers predicts the channel flow with higher accuracy in comparison with the other eddy viscosity models in the logarithmic layer. The elliptic-relaxation type models, $k-{\epsilon}-v2$, and ERRSM have the advantage of application to complex geometries and show good prediction for near wall flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.5
• /
• pp.610-620
• /
• 1999

Recently, numerous modifications of low Reynolds number $k-{\epsilon}$ model have boon carried out with the aid of DNS data. However, the previous models made in this way are too intricate to be used practically. To overcome this shortcoming, a new low Reynolds number $k-{\epsilon}$ model has boon developed by considering the distribution of turbulent properties near the wall. This study proposes the revised a turbulence model for prediction of turbulent flow with adverse pressure gradient and separation. Nondimensional distance $y^+$ in damping functions is changed to $y^*$ and some terms modeled for one dimensional flow in $\epsilon$ equations are expanded into two or three dimensional form. Predicted results by the revised model show an acceptable agreement with DNS data and experimental results. However, for a turbulent flow with severe adverse pressure gradient, an additive term reflecting an adverse pressure gradient effect will have to be considered.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.741-750
• /
• 1995

The turbulent buoyancy-driven flow in 2-dimensional enclosed cavities heated from the vertical side is numerically calculated for both cases of a Rayleigh number of 5*10$^{10}$ for air and 2.5*10$^{10}$ for water. Three different turbulence models are considered : standard k-.epsilon. model of Ozoe and low-Reynolds-number model of Lam and Bremhorst, and another low-Reynolds-number model of Davidson. The results indicate that the use of low-Reynolds number models is recommended for the indoor airflow computation, and the results from Davidson model are reasonably close to the reported experimental data. A sensitivity study shows that the amounts of wall-heat transfer and the velocity profiles with the Lam and Bremhorst model largely depend on the choice of the wall function for .epsilon..

• 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.

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

A turbulent flow in a conical diffuser with total divergence angle of 8.deg. was numerically studied. The low Reynolds number k-.epsilon. model(Launder-Sharma model) was adopted to simulate the turbulence. The continuity and time averaged Navier-Stokes equations in a nonorthogonal coordinate system were solved by a finite volume method based on the fully elliptic formulation. The low Reynolds number k-.epsilon. model reasonably simulates the pressure recovery and the mean velocity components. However, there are also considerable discrepancies between predicted and measured shear stress distribution on the wall and turbulent kinetic energy distributions. It is necessary to investigate the flow structure at the entry of the diffuser, numerically as well as experimentally.