• Journal of computational fluids engineering
• /
• v.7 no.4
• /
• pp.8-18
• /
• 2002

The primary objective of the present study is evaluation of the k-ε-vv-f turbulence model for prediction of natural convection in a rectangular cavity. As a comparative study, the two-layer k-ε model is also considered. Both models, with and without algebraic heat flux model, are applied to the analysis of natural convection in a rectangular cavity. The performances of turbulence models are investigated through comparison with available experimental data. The predicted results of vertical velocity component, turbulent heat fluxes, turbulent shear stress, local Nusselt number and wall shear stress are compared with experimental data. It is shown that, among the turbulence models considered in the present study, the k-ε-vv-f model with an algebraic heat flux model predicts best the vertical mean velocity and velocity fluctuation, and the inclusion of algebraic heat flux model slightly improves the accuracy of results.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.1-10
• /
• 2016

The analysis of characteristics of turbulent flow and thermal boundary layer for natural convection caused by fire along vertical wall is performed. The 4m-high vertical copper plate is heated and kept at a uniform surface temperature of $60^{\circ}C$ and the surrounding fluid (air) is kept at $16.5^{\circ}C$. The flow and temperature is solved by large eddy simulation(LES) of FDS code(Ver.6), in which the viscous-sublayer flow is calculated by Werner-Wengle wall function. The whole analyzed domain is assumed as turbulent region to apply wall function even through the laminar flow is transient to the turbulent flow between $10^9$<$Gr_z$<$10^{10}$ in experiments. The various grids from $7{\times}7{\times}128$ to $18{\times}18{\times}128$ are applied to investigate the sensitivity of wall function to $x^+$ value in LES simulation. The mean velocity and temperature profiles in the turbulent boundary layer are compared with experimental data by Tsuji & Nagano and the results from other LES simulation in which the viscous-sublayer flow is directly solved with many grids. The relationship between heat transfer rate($Nu_z$) and $Gr_zPr$ is investigated and calculated heat transfer rates are compared with theoretical equation and experimental data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.1
• /
• pp.112-121
• /
• 1992

A numerical analysis was carried out to study two-dimensional turbulent natural convection in a square enclosure containing fluid of Prandtl number 6.05 within internal energy sources. The square enclosure was bounded by four rigid planes of constant equal temperature. Inclination angles of 0, 15, 30 and 45 deg. from the horizon for Rayleigh numbers from 1 * 10$^{6}$ to 1 * 10$^{9}$ were studied. Local and average Nusselts numbers are obtained on all four walls. If inclination angle exists, the average Nusselt number appears in increasing order at bottom, left, right and top wall.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.7 no.1
• /
• pp.27-40
• /
• 1997

The heat in Czochralski method is transfered by all transport mechanisms such as convection, conduction and radiation and convection is caused by the temperature difference in the molden pool, the rotations of crystal or crucible and the difference of surface tension. This study delvelops the simulation model of Czochralski growth by using the finite difference method with fixed grids combined with new latent heat treatment model. The radiative heat transfer occured in the surfce of the system is treated by calculating the view factors among surface elements. The model shows that the flow is turbulent, therefore, turbulent modeling must be used to simulate the transport phenomena in the real system applied to 8" Si single crystal growth process. The effects of a cusp magnetic field imposed on the Czochralski silicon melt are studied by numerical analysis. The cusp magnetic field reduces the natural and forced convection due to the rotation of crystal and crucible very effectively. It is shown that the oxygen concentration distribution on the melt/crystal interface is sensitively controlled by the change of the magnetic field intensity. This provides an interesting way to tune the desired O concentration in the crystal during the crystal growing.

• 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
• /
• v.17 no.6
• /
• pp.1609-1620
• /
• 1993

A Study on the trombe wall system, a kind of passive solar systems, has been performed numerically. The system is treated as a two-dimensional steady turbulent natural convection including constant heat source per unit area. The numerical code, "PHOENICS, " was employed to analyze this conduction-convection conjugated heat transfer. The general mode of the flow field was examined, and the exchange of mass between two recirculating flows is found to be the major mechanism of the heat transfer. It is shown that the performance is affected by the changes in the geometrical factors-the thickness of the wall, the width between the windowand the wall, and size of the vents. Further analysis has been performed to show the optimal geometry with regard to the last two factors.o factors.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.4
• /
• pp.556-565
• /
• 1995

The physical model of interest is based upon the concentric cylinder, where the outside cylinder is filled with optically thick and high temperature phase change material(PCM). The fluid is flowing through the inside cylinder to transfer the appropriate energy. The fluid is flowing through the inside cylinder to transfer the appropriate energy. The governing equations for the phase change material including internal thermal radiation and for the turbulent transfer fluid have been employed and numerically solved. The optically thick phase change justifies the P-l spherical harmonics approximation, which is believed to be appropriate choice particularly for the much coupled problem like in this study. The solid/liquid interface, temperature distribution within the PCM and the heat flux from the PCM to the transfer fluid have been obtained and compared with those of laminar transfer fluid. The numerical results show that the turbulent transfer fluid accelerates the solid/liquid interface and results in the increase of heat transfer rate from the PCM. The internal thermal radiation within the PCM, however, does not always playa role to increase the heat transfer rate throughout the inside cylinder. It is believed that the combined heat flux has been picked up more in the inflowing area than in the pure conductive phase change material.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.12 no.2
• /
• pp.150-160
• /
• 2000

Turbulent natural convective flow and heat transfer in a square enclosure with horizontal partition are investigated numerically. The enclosure is composed of a lower hot and a upper cold horizontal walls and adiabatic vertical walls. Partitions carried with the upward, downward, and both control plates are attached perpendicularly to the one of the vertical insulated walls, respectively. The low Reynolds number $k-\varepsilon$ model is adopted to calculate the turbulent thermal convection. The governing equations are solved by using the finite element method with Galerkin method. The computations have been carried out by varying the length of partition, the position of control plates, and the Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height for water(Pr=4.95). When the control plates are attached at the edge of partition, the stability of oscillating flow grows wrose with the increase of Rayleigh number and the partition length. The heat transfer rate has been reducer than that of no control plate due to the restraint of control plates with the increase of Rayleigh number.

• Fisheries and Aquatic Sciences
• /
• v.7 no.2
• /
• pp.90-95
• /
• 2004

An analytical model for predicting the convection-diffusion of solute dumped in a homogeneous open sea of constant water depth has been developed in a time-integral form. The model incorporates spatially uniform, uni-directional, mean and oscillatory currents for horizontal convection, the settling velocity for the vertical convection, and the anisotropic turbulent diffusion. Two transformations were introduced to reduce the convection-diffusion equation to the Fickian type diffusion equation, and then the Galerkin method was then applied via the expansion of eigenfunctions over the water column derived from the Sturm-Liouville problem. A series of calculations has been performed to demonstrate the applicability of the model.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.15 no.4
• /
• pp.129-134
• /
• 2005

A numerical analysis was performed on magnetic field effects of silicon melt motion in Czochralski crystal puller. The turbulent modeling was used to simulate the transport phenomena in 18' single crystal growing process. For small crucible angular velocity, the natural convection is dominant. As the crucible angular velocity is increased, the forced convection is increased and the distribution of temperature profiles is broadened. The cusp magnetic field reduces effectively the natural and forced convection near the crucible and the temperature profiles of the silicon fluids is similar in the case of conduction.