• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.10a
• /
• pp.480-486
• /
• 1995

Turbulent flow field in a subchannel of bare rod bundles has been numerically simulated using the control volume based finite element method. Launder & Ying model of Reynolds stress and Lam & Bremhorst low-Reynolds number model are implemented in k-$\varepsilon$ equations and momentum equations. Secondary flows are simulated using the stream function and vorticity approach. The control volume based finite element method enable to use the upwind scheme (donor cell scheme). Sensitivity of the constants in the models are studied, and proper values are found to get the close result to the measured flow distributions.

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

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.825-830
• /
• 2004

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc (DCA) compressor cascades. Two types of double-circular-arc cascades were used in this analysis. The appropriate turbulence model for compressor analysis was selected among the conventional turbulence models such as Baldwin-Lomax, k-$\varepsilon$ and k-$\varepsilon$ models. The results of current study were compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE/PWIM algorithm for collocated grid and hybrid scheme for the convective terms were the main features of numerical tools. As commonly known, turbulence modeling is very important for the prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. For selection of turbulence model, 2-D analysis was performed. And then, k-$\varepsilon$ turbulence model with wall function chosen as the reasonable turbulence model for 3-D calculation was used to increase the efficiency of computation times. A reasonable result of 3-D flow pattern passing through the double-circular-arc cascade was obtained.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1569-1580
• /
• 2001

A new nonlinear near-wall turbulence model is developed to predict turbulent flow and heat transfer in strongly nonequilibrium flows. The k-$\varepsilon$-f$\sub$${\mu}$/, model of Park and Sung$\^$(1)/ is extended to a nonlinear formulation. The stress-strain relationship is the thrid-order in the mean velocity gradients. The strain dependent coefficients are obatined from the realizability constraints and the singular behavior at large strains. An improved explicit heat flux model is proposed with the aid of Cayley-Hamilton theorem. This new model includes the quadratic effects of flow deformations. The near-wall asymptotic behavior is incorporated by modifying the f$\sub$λ/ function. The model performance is shown to be satisfactory.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1254-1264
• /
• 1999

A three-dimensional coupled turbulent fluid flow and solidification process were analyzed in a continuous casting process of a steel slab with Electromagnetic Brake(EMBR). A revised low-Reynolds number $k-{\varepsilon}$ turbulence model was used to consider the turbulent effects. The enthalpy-porosity relation was employed to suppress the velocity within a mushy region. The electromagnetic field was described by Maxwell equations. Tile application of EMBR to the mold region results in the decrease of the transfer of superheat to the narrow face, the increase of temperature in free surface region and most liquid of submold region, and the higher temperature gradient near the solidifying shell. The increasing magnetic flux density effects mainly to the surface temperature of the solidifying shell of narrow face, hardly to the one of wide face. It is seen that in the presence of EMBR a thicker solidifying shell is obtained at the narrow face of the slab.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.652-655
• /
• 2008

Despite of the laminar-turbulent transition region co-exist with fully turbulence region around the leading edge of an airfoil, still lots of researchers apply to fully turbulence models to predict aerodynamic characteristics. It is well known that fully turbulent model such as standard k-${\varepsilon}$ model couldn't predict the complex stall and the separation behavior on an airfoil accurately, it usually leads to over prediction of the aerodynamic characteristics such as lift and drag forces. So, we apply correlation based transition model to predict aerodynamic performance of the NREL (National Renewable Energy Laboratory) Phase IV wind turbine. And also, compare the computed results from transition model with experimental measurement and fully turbulence results. Results are presented for a range of wind speed, for a NREL Phase IV wind turbine rotor. Low speed shaft torque, power, root bending moment, aerodynamic coefficients of 2D airfoil and several flow field figures results included in this study. As a result, the low speed shaft torque predicted by transitional turbulence model is very good agree with the experimental measurement in whole operating conditions but fully turbulent model(k-${\varepsilon}$) over predict the shaft torque after 7m/s. Root bending moment is also good agreement between the prediction and experiments for most of the operating conditions, especially with the transition model.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.304-308
• /
• 2000

This paper presents a numerical optimization method to design geometric shape of streamwise periodic ribs mounted on one of the principal walls to enhance turbulent heat transfer in a rectangular channel flow. The golden section method is used for the one dimensional search. The optimization is based on Wavier-Stokes analysis of turbulent forced convection with $k-{\varepsilon}$ turbulence model. The width-to-height ratio of a rib is chosen as a design variable. The object function is defined as an inverse of average Nusselt number. An optimum shape of the rib has been obtained with reasonable computing time.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.6
• /
• pp.824-834
• /
• 1999

Characteristics of two-phase flow and heat transfer were numerically investigated in a submerged gas injection system when temperature of the injected gas was different from that of the liquid. The Eulerian approach was used for both the continuous and dispersed phases. The turbulence in the liquid phase was modeled using the standard $k-\varepsilon$$\varepsilon$ turbulence model. The interphase friction and heat transfer coefficient were calculated from the correlations available in the literature. The turbulent dispersion of the phases was modeled by a "dispersion Prandtl number". In the case with heat transfer where the temperature of the injected gas is higher than the mean liquid temperature, the axial and the radial velocities are lower in comparison with the case of homogeneous temperatures. The results in the present research are of interest in the design and operation of a wide variety of material and chemical processes.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.17 no.4
• /
• pp.386-395
• /
• 2005

With the development of a living standard, the importance of indoor air conditioning system in all kinds of buildings and vehicles has increased. A lot of researches on energy losses in a duct and various kinds of flow pattern in branches or junctions have been carried out over many years, because the primary object of a duct system used in HVAC is to provide equal flow rate in the interior of each room by minimizing pressure drop. In this study, to get equal flow distribution in each branch, a blockage is applied to the rectangular duct system. The flow analysis for flow distribution of a rectangular duct with two branches was performed by CFD. By using SIMPLE algorithm and finite volume method, flow analysis is performed in the case of 3-D, incompressible, turbulent flow. Also, the standard $k-{\varepsilon}$ model and wall function method were used for analysis of turbulent fluid flow. The distribution diagrams of static pressure, velocity vector, turbulent energy and kinetic energy in accordance with variation of Reynolds number and blockages location in a rectangular duct show that flow distribution at duct outlets is improved by a blockage. In this rectangular duct system, mean velocity and flow rate distribution in two branch outlets are nearly constant regardless of variation of Reynolds number, and a flow pattern of the internal duct has a same tendency as well.