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

• Journal of the Korean Society of Safety
• /
• v.15 no.1
• /
• pp.19-27
• /
• 2000

This paper numerically deals with combined heat transfer in a enclosure with a block. The block affected by hot wall is located centrally in the enclosure with a radiating gray gas. The discrete ordinate method(DOM) was used for solving the radiative transfer equation. Both laminar and turbulent cases were investigated for various Rayleigh number and standard k-$\varepsilon$ model was adopted to turbulent case. The effects of optical thickness, wall emissivity and fluid-solid thermal conductivity ratio are investigated on the flow and temperature fields. This study shows that as the wall emissivity decreases, the temperature distribution gradually becomes uniform and the heat transfer is reduced in enclosure. It is expected that this study can help to design the energy system related to the combined heat transfer and operate it safely.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.29-37
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• International Journal of Fluid Machinery and Systems
• /
• v.4 no.1
• /
• pp.133-139
• /
• 2011

In this paper, Based on Two-dimensional Flow Theory, adopting quasi-orthogonal method and point-by-point integration method to design the impeller of the low specific speed centrifugal pump by code, and using RANS (Reynolds Averaged N-S) Equation with a standard k-${\varepsilon}$ two-equation turbulence model and log-law wall function to solve 3D turbulent flow field in the impeller of the low specific speed pump. An analysis of the influences of the blade profile on velocity distributions, pressure distributions and pump performance and the investigation of the flow regulation pattern in the impeller of the centrifugal pump are presented. And the result shows that this method can be used as a new way in low speed centrifugal pump impeller design.

• Nuclear Engineering and Technology
• /
• v.30 no.6
• /
• pp.581-591
• /
• 1998

Numerical analysis was peformed for the two-dimensional turbulent natural convection for a long horizontal line with different end temperatures. The turbulent model has been applied a standard k-$\varepsilon$ two equation model of turbulence similar to that the proposed by the Launder and Spalding. The dimensionless governing equations are solved by using SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm which is developed using control volumes and staggered grids. The numerical results are verified by comparison with the operating PWR test data. The analysis focuses on the effects of variation of the heat transfer rates at the pipe surface, the thermal conductivities of the pipe material and the thickness of the pipe wall on the thermal stratification. The results show that the heat transfer rate at the pipe surface is the controlling parameter for mitigating of thermal stratification in the long horizontal pipe. A significant reduction and disappearance of the thermal stratification phenomenon is observed at the Biot number of 4.82$\times$10$^{-1}$ . The results also show that the increment of the thermal conductivity and thickness of the wall weakens a little the thermal stratification and somewhat reduces temperature gradient of y-direction in the pipe wall. These effects are however minor, when compared with those due to the variation of the heat transfer rates at the surface of the pipe wall.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Wind and Structures
• /
• v.23 no.5
• /
• pp.465-483
• /
• 2016

This study aims to enhance the understanding of the surface pressure distribution around rectangular bodies, by considering aspects such as the suction pressure at the leading edge on the top and side faces when the body aspect ratio and wind direction are changed. We carried out wind tunnel measurements and numerical simulations of flow around a series of rectangular bodies (a cube and two rectangular bodies) that were placed in a deep turbulent boundary layer. Based on a modern numerical platform, the Navier-Stokes equations with the typical two-equation model (i.e., the standard $k-{\varepsilon}$ model) were solved, and the results were compared with the wind tunnel measurement data. Regarding the turbulence model, the results of the $k-{\varepsilon}$ model are in overall agreement with the experimental results, including the existing data. However, because of the blockage effects in the computational domain, the pressure recovery region is underpredicted compared to the experimental data. In addition, the $k-{\varepsilon}$ model sometimes will fail to capture the exact flow features. The primary emphasis in this study is on the flow characteristics around rectangular bodies with various aspect ratios and approaching wind directions. The aspect ratio and wind direction influence the type of wake that is generated and ultimately the structural loading and pressure, and in particular, the structural excitation. The results show that the surface pressure variation is highly dependent upon the approaching wind direction, especially on the top and side faces of the cube. In addition, the transverse width has a substantial effect on the variations in surface pressure around the bodies, while the longitudinal length has less influence compared to the transverse width.

• Journal of Energy Engineering
• /
• v.5 no.1
• /
• pp.50-55
• /
• 1996

A numerical simulation on a primary calcinator of porcelain was performed with using Fluent to calculate the heat efficiency by studying velocity vector and temperature profile according to variables such as the location of outlet and porcelain. Control-Volume based Finite Difference Method and Up-wind scheme are used for discretization of differential equation. SIMPLEC Algorithm and standard k-$\varepsilon$ turbulent model are selected to resolve the pressure-velocity coupling and the turbulent. The result of simulation showed that the whole velocity vector field in a calcinator was varied greatly according to the location of outlet. But the whole temperature profile at each zone was still high regardless of the location of outlet because of the radiation. But the temperature of a case with a outlet at sidepart of preheating or cooling zone was little high compared to the case with a outlet on the top of preheating zone. The velocity vector field and temperature profile in a calcinator were almost not affected by the location of porcelain, but the temperature inside a porcelain was much affected according to the place where it was located. The heat efficiency in a calcinator was 44.6% and the gas temperature in the outlet was about 1000 K.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.684-689
• /
• 2001

In this study, three-dimensional viscous flow analysis and optimization are presented for the design of a mixed-flow fan. Steady, imcompressible, three-dimensional Reynolds averaged Navier-Stokes equations are used as governing equations, and standard $k-{\varepsilon}$ turbulence model is chosen as a turbulence model. Governimg equations are discretized using finite volume method. Upwind difference scheme is used for the discretization of the convective term and SIMPLEC algorithm is used as a velocity-pressure correction procedure. The computational results are compared with the results obtained by TASCflow. For the numerical optimization of the design, objective function is defined as a ratio of generation of the turbulent energy to pressure head. Sweep angles are used as design variables.

• Journal of Mechanical Science and Technology
• /
• v.14 no.9
• /
• pp.1013-1019
• /
• 2000

A numerical investigation was performed to determine the effect of the Gurney flap on a NACA 23012 airfoil. A Navier-Stokes code, RAMPANT, was used to calculate the flow field about the airfoil. Fully-turbulent results were obtained using the standard ${\kappa}-{\varepsilon}$ two-equation turbulence model. The numerical solutions showed that the Gurney flap increased both lift and drag. These results suggested that the Gurney flap served to increase the effective camber of the airfoil. The Gurney flap provided a significant increase in the lift-to-drag ratio relatively at low angle of attack and for high lift coefficient. It turned out that 0.6% chord size of flap was the best. The numerical results exhibited detailed flow structures at the trailing edge and provided a possible explanation for the increased aerodynamic performance.