• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.10 no.3
• /
• pp.185-219
• /
• 1981

To analyze two dimensional incompressible laminar natural convection in a rectangular enclosure heated from below and cooled by a horizontal ceiling and two vertical walls, he primitive Navier-Stokes equations and the energy equation were solved numerically in time dependent form by a marker and cell method. A successive over-relaxation method for the elliptic portion of the problem and an explicit method for the parabolic portion were applied for the range of Grashoff number of $5{\times}10^3\;to\;5{\times}10^4$ to get the transient and steady state dimensionless temperature and velocity profiles. For the range of aspect ratio $L/H{\leqq}2.4$ in which only a pair of convection rolls exists mean Nusselt number calculated are as follows : $$N_{NU}0.89\;N_{Gr}^{0.2}(H/L)^{0.45}$$ By path lines drawn by marker particle trajectories roll number of cellular motion were observed for various aspect ratio of the enclosure.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2595-2606
• /
• 1995

Buoyancy-assisted melting of an unconstrained ice in an isothermally heated horizontal enclosure was numerically analyzed in a range of wall temperatures encompassing the density inversion point. The problem as posed here involves two physically distinct domains each of which has its own scales and respective heat transfer mode. These two domains join at the junction where the liquid squeezed out of the film region flushes into the lower melt pool. Both of these domains have been treated separately in the literature by a patching technique which invokes several, otherwise unnecessary, assumptions. The present study eliminates successfully such a superfluous procedure by treating the film and lower melt pool regions as a single domain. As a result of this efficient solution procedure, the interaction of the water stream ejected at the junction and the natural convection in the melt pool could be clarified for different wall temperatures. Though limited by two-dimensionality, the present results conformed indirectly the earlier reported transition of the flow pattern, as the wall temperature was increased over the density inversion point. The transient evolution of the melting surface, the time rate of change in melt volume fraction, the local and temporal variation of the heat transfer coefficients are analyzed and presented.

• Nuclear Engineering and Technology
• /
• v.56 no.2
• /
• pp.611-623
• /
• 2024

The ISP-47 TOSQAN experiment was analyzed with containmentFOAM which is an open-source CFD code based on OpenFOAM. The containment phenomena taking place during the experiment are gas mixing, stratification and wall condensation in a mixture composed of steam and non-condensable gas. The k-ω SST turbulence model was adopted with buoyancy turbulence models. The wall condensation model used is based on the diffusion layer approach. We have simulated the full TOSQAN experiment which had a duration 20000 s. Sensitivity studies were conducted for the buoyancy turbulence models with SGDH and GGDH and there were not significant differences. All the main features of the experiments namely pressure history, temperature, velocity and gas species evolution were well predicted by containemntFOAM. The simulation results confirmed the formation of two large flow stream circulations and a mixing zone resulting by the combined effects of the condensation flow and natural convection flow. It was found that the natural convection in lower region of the vessel devotes to maintain two large circulations and to be varied the height of the mixing zone as result of sensitivity analysis of non-condensing wall temperature. The computational results obtained with the 2D mesh grid approach were comparable to the experimental results.

• Journal of applied mathematics & informatics
• /
• v.27 no.1_2
• /
• pp.301-313
• /
• 2009

This paper presents a Finite Element Method (FEM) application to thermal study of natural three layers of human dermal parts of varying properties. This paper carries out investigation of temperature distributions in these layers namely epidermis, dermis and under lying tissue layer. It is assumed that the outer skin is exposed to atmosphere and the loss of heat due to convection, radiation and evaporation of water have also been taken into account. The computations are carried out for one dimensional unsteady state case and the shape functions in dermal parts have been considered to be quadratic. A Finite Element scheme that uses the Crank-Nicolson method is used to solve the problem and the results computed have been exhibited graphically.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.16 no.2
• /
• pp.204-215
• /
• 1987

Laminar natural convective heat transfer within a two-dimensional rectangular enclosure with horizontal conducting walls and partitions was investigated by numerical analysis and experiment. The enclosure consists of two isothermal vertical walls and two adiabatic horizontal walls. This combined heat transfer problem of conduction and natural convection was solved using finite difference method with SIMPLE algorithm, and temperature distribu-tions in the air filled enclosure was obtained using Mach-Zehnder interferometer. Good agree-ment was obtained between the predicted and measured results. The effect of geometric parameters and thermal properties on heat transfer was studied far Grashof numbers in range, $1\times10^4\;{\leqslant}\;G^r\;{\leqslant}\;6.4\times10^5.$ It was found that both velocity and temperature fields were in-fluenced significantly by thermal conductivity of the conducting walls and the partitions, and by geometry of partitions.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.4
• /
• pp.1511-1517
• /
• 2012

Vacuum chamber or vacuum tube for the exhaust process of vacuum glazing is presently used, where excessive time and expenses are required to make the whole vacuum chamber or tube vacuum. To solve this problem, multifunction equipment for vacuum exhaust and ultrasonic joining at atmospheric pressure has been developed, in which a piezoelectric vibrator experiences excessive temperature rise resulting in optimizing the cooling of the equipment. Therefore, in this study, cooling effects of natural convection and forced convection methods were identified by numerical analysis and experiments, and cooling performance of the multifunction equipment was optimized.

• The Journal of Engineering Research
• /
• v.1 no.1
• /
• pp.23-30
• /
• 1997

Energy transfer by free convection arises in many engineering applications, such as a hot steam radiator for heating a room, refrigeration coils, electric transformers, heating elements and electronic equipments. Generally unsteady natural convection flow in a horizontal channel with arbitrary wall temperatures and the mathematical and physical basis of convection transport has been considered in general. A physically meaningful exact solution of the problem has been obtained in a closed form by the application of the standard finite sine transform technique. Influences of the governing parameters, the Prandtl number and the Rayleigh number, to bring the flow and heat transfer to final steady states have been discussed separately. For constant values of the arbitray wall temperatures and of the function, determining the average axial velocity, the final steady state is approached in different times respectively for the cases when the Prandtl number Pr>1 and Pr<1. It is also seen that the function, representing the axial temperature gradient, is influenced by none of the governing parameters : but the steady state flow is influenced only by the Rayleigh number. There are, of course, many applications. Free convection strongly influences heat transfer from pipes and transmission lines, as well as from various electronic devices. It is also relevant to the environmental sciences, where it is responsible for oceanic and atmospheric motions, as well as related heat transfer processes.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.22-29
• /
• 1998

Parallel implementation is conducted for a SIMPLER finite volume model. The present parallelism is based on domain decomposition and explicit message passing using MPI and SHMEM. Two parallel solvers to tridiagonal matrix equation are employed. The implementation is verified on the Cray T3E system for a benchmark problem of natural convection in a sidewall-heated cavity. The test results illustrate good scalability of the present parallel models. Performance issues are elaborated in view of convergence as well as conventional parallel overheads and single processor performance. The effectiveness of a localized matrix solution algorithm is demonstrated.

• Journal of Electrical Engineering and Technology
• /
• v.4 no.4
• /
• pp.510-514
• /
• 2009

This paper presents the coupled analysis method to calculate the temperature rise in a gas insulated busbar (GIB). Harmonic eddy current analysis is carried out and the power losses are calculated in the conductor and enclosure tank. Two methods are presented to analyze the temperature distribution in the conductor and tank. One is to solve the thermal conduction problem with the equivalent natural convection coefficient and is applied to a single phase GIB. The other is to employ the computational fluid dynamics (CFD) tool which directly solves the thermal-fluid equations and is applied to a three-phase GIB. The accuracy of both methods is verified by the comparison of the measured and calculated temperature in a single phase and three-phase GIB.

• Steel and Composite Structures
• /
• v.24 no.3
• /
• pp.359-367
• /
• 2017

In this work, transient heat transfer analysis of functionally graded (FG) carbon nanotube reinforced nanocomposite (CNTRC) cylinders with various essential and natural boundary conditions is investigated by a mesh-free method. The cylinders are subjected to thermal flux, convection environments and constant temperature faces. The material properties of the nanocomposite are estimated by an extended micro mechanical model in volume fraction form. The distribution of carbon nanotube (CNT) has a linear variation along the radial direction of axisymmetric cylinder. In the mesh-free analysis, moving least squares shape functions are used for approximation of temperature field in the weak form of heat transform equation and the transformation method is used for the imposition of essential boundary conditions. Newmark method is applied for solution time depended problem. The effects of CNT distribution pattern and volume fraction, cylinder thickness and boundary conditions are investigated on the transient temperature field of the nanocomposite cylinders.