• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.10
• /
• pp.941-951
• /
• 2013

A three dimensional numerical simulation of laminar flow and heat transfer in fully developed curved pipe flow has been performed to study the effects of Dean number and pipe curvature on the flow and temperature fields under the thermal boundary condition of axially uniform wall heat flux. The Reynolds number under consideration ranges from 100 to 4000, and the Prandtl number is 0.71. The curvature ratios are 0.01, 0.025, 0.05 and 0.1. The axial velocity and temperature profiles and the local Nusselt number obtained from the present study are in good agreement with the previous numerical and experimental results currently available. To show the effects of pipe curvature on the flow and heat transfer, the resistance coefficients and heat transfer coefficients are computed and compared with the results of the previous theoretical and experimental studies. The averaged Nusselt number is correlated with Dean and Prandtl numbers. Furthermore, the critical Reynolds number for transition to turbulent flow is observed to depend upon the curvature ratio.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.12
• /
• pp.3329-3343
• /
• 1995

The interaction of mixed convection and surface radiation in a printed circuit board(PCB) is investigated numerically. The electronic equipment is modeled by a two-dimensional channel with three hot blocks. In order to calculate the turbulent flow characteristics, the low Reynolds number k-.epsilon. model which is proposed by Launder and Sharma is applied. The S-4 approximation is used to solve the radiative transfer equation. The effects of the Reynolds number and geometric configuration variation of PCB on the flow and heat transfer characteristics are analyzed. As the results of this study, it is found that the thermal boundary layer occured at adiabatic wall in case with thermal radiation included, and the effect of radiation is also found to be insignificant for high Reynolds numbers. It is found, as well, that the heat transfer increases as the Reynolds number and block space increase and the channel height decreases and the heat transfer of vertical channel is greater than that of horizontal channel.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.92-99
• /
• 2011

In order to protect the turbine blade from working fluid of high temperature, many cooling techniques such as internal convection cooling, film cooling, impinging jet cooling and thermal barrier coating have been developed. With all other things, film-cooling has been widely used as the important alternative. In the present work, numerical analysis has been performed to investigate and to compare the film-cooling performance of various film-cooling hole schemes such as cylindrical, crescent, louver, and dumbbell holes. To analyze the turbulent flow and the film-cooling mechanism, three-dimensional Reynolds-averaged Navier-Stokes analysis has been performed with shear stress transport turbulence model. The validation of numerical results has been assessed in comparison with experimental data. The characteristics of fluid flow and the film-cooling performance for each shaped hole have been investigated and evaluated in terms of centerline, laterally averaged and spatially averaged film-cooling effectivenesses. Among the film cooling holes, the dumbbell shaped hole shows better film-cooling effectiveness than the other shaped holes. And the louver and cylindrical shaped hole show the worst film cooling performance, and concentrated flows on near the centerline only.

• The KSFM Journal of Fluid Machinery
• /
• v.3 no.4 s.9
• /
• pp.22-29
• /
• 2000

Numerical analysis using three dimensional finite volume method for the discretization, adaptive grid method for the numerical accuracy, multiple rotating frame method for the rotating body and the standard $k-{\epsilon}$ model for the turbulent flow was performed to understand the heat transfer phenomena and to improve the efficiency of the scroll compressor. The temperature measurement was carried out under ARI condition. It was found that the fluid temperature in the compressor was predicted accurately while the temperature of the motor coil showed large discrepancy between the calculation and experiment due to the large anisotropy of the conductivity and non homogeneity. We found that the efficiency of the compressor depends on the inlet temperature of the compressing part and the flow pattern around the inlet region of the compressing part influences the inlet temperature due to high surface temperature of the main frame. The efficiency of the compressor using Coanda effect is higher than the previous one because the smooth suction at the inlet region of the compressing part leads to low heat transfer to the refrigerant of the compressor.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.5
• /
• pp.604-611
• /
• 2004

Simulation of three-dimensional turbulent flow with LES and DNS lakes much time and expense with currently available computing resources and requires big computing resources especially for high Reynolds number. The emerging alternative to provide the required computing power and working environment is the Grid computing technology. We developed the CFD code which carries out the parallel computing under the Grid environment. We constructed the Grid environment by connecting different PC-cluster systems located at two different institutes of Pusan National University in Busan and KISTI in Daejeon. The specification of PC-cluster located at two different institutes is not uniform. We run our parallelized computer code under the Grid environment and compared its performance with that obtained using the homogeneous computing environment. When we run our code under the Grid environment, the communication time between different computer nodes takes much larger time than the real computation time. Thus the Grid computing requires the highly fast network speed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.5
• /
• pp.381-388
• /
• 2009

In some portions of nuclear piping systems, stratification phenomena may occur due to the density difference between hot and cold stream. When the temperature difference is large, the stratified flow under diverse operating conditions can produce high thermal stress, which leads to unanticipated piping integrity issues. The objectives of this research are to examine controvertible numerical factors such as model size, grid resolution, turbulent parameters, governing equation, inflow direction and pipe wall. Parametric three-dimensional computational fluid dynamics analyses were carried out to quantify effects of these parameters on the accuracy of temperature profiles in a typical nuclear piping with complex geometries. Then, as a key finding, it was recommended to use optimized mesh of real piping with the conjugated heat transfer condition for accurate thermal stratification analyses.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.19 no.7
• /
• pp.500-511
• /
• 2007

The Optimal solutions of the design variables in internally finned tubes have been obtained for three-dimensional periodically fully developed turbulent flow and heat transfer. For a trapezoidal fin profile, performances of the heat exchanger are determined by considering the heat transfer rate and pressure drop, simultaneously, that are interdependent quantities. Therefore, Pareto frontier sets of a heat exchanger can be acquired by integrating CFD and a multi-objective optimization technique. The optimal values of fin widths $(d_1,\;d_2)$, fin height(h) and helix angle$(\gamma)$ are numerical1y obtained by minimizing the pressure loss and maximizing the heat transfer rate within ranges of $d_1=0.5\sim1.5mm$, $d_2=0.5\sim1.5mm$, $h=0.5\sim1.5mm$, and $\gamma=0\sim20^{\circ}$. For this, a general CFD code and a global genetic algorithm(GA) are used. The Pareto sets of the optimal solutions can be acquired after $30^{th}$ generation.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.19 no.2
• /
• pp.191-201
• /
• 2007

It was noted by the several researchers that the voltage outputs in response to a single yawed hot-wire sensor in a flow perpendicular to the axis deviate from the theoretical voltage output by King's law and Jorgensen's relation. This study noticed that the calibration coefficients of original Grande's method are not constant and fairly sensitive to the radial angle (${\alpha}_{R}$). For more accuracy, this study interpolated the parameters of the Grande relation as a function of radial angle and compared velocity components with ones by Jorgensen and original Grande relation in the calibration jet flow. Finally, as a test case, 3-dimensional turbulent flows of the inlet plane of 180 degree bend are measured and compared the velocity components by above three methods and showed the characteristics of the flows.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.11a
• /
• pp.363-366
• /
• 2006

The effectiveness of two passive control techniques for alleviating the pressure oscillation generated in a supersonic cavity flow is investigated numerically. The passive devices suggested in the present research include a triangular bump and a sub-cavity installed near the upstream edge of a rectangular cavity. The supersonic cavity flow characteristics are examined by using the three-dimensional, unsteady Wavier-Stokes computation based on a finite volume scheme. Large eddy simulation (LES) is carried out to properly predict the turbulent features of cavity flow. The results show that the pressure oscillation near the downstream edge dominates overall time-dependent cavity pressure variations. Such an oscillation is attenuated more considerably using the sub-cavity compared with other methods, and a larger sub-cavity leads to better control performance.