• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

• 한국신재생에너지학회:학술대회논문집
• /
• 2008.05a
• /
• pp.642-645
• /
• 2008

A pin fin with variable fin base thickness is analysed by using the one dimensional analytic method. Heat loss and fin thermal resistance are presented as a function of the fin base thickness, pin fin outer radius and convection characteristic numbers ratio. The relationship between the fin outer radius and fin base thickness for the same amount of heat loss is shown. One of the results indicates the fin thermal resistance decreases as the fin outer radius and/or convection characteristic numbers ratio increase whereas the fin thermal resistance is independent on the variation of fin base thickness.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.4
• /
• pp.742-747
• /
• 2009

Many works on the temperature distribution of power apparatus have usually done by coupled magneto-thermal analysis. Such a method can not consider the internal gas or oil flow in the power apparatus such as gas insulated switchgear, GIS bus bar, and power transformer. Moreover it can not show the internal temperature distribution of the power apparatus exactly. This paper proposes a coupled magneto-thermal-flow analysis considering Navier-Stokes equations. The convection heat transfer coefficient is calculated analytically by applying Nusselt number for natural convection and is applied to the boundary condition of proposed method. Temperature distribution of the GIS bus bar model considering thermal flow is obtained by the proposed method and shows good agreement with the experimental data.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.5
• /
• pp.365-375
• /
• 1991

Confined natural convection due to lateral temperature difference in rectangular enclosures was studied numerically and experimentally for the insulated and the constant temperature enclosures. In the case of insulated enclosure, the flow pattern and heat transfer modes are rather simple depending mainly upon Rayleigh number. In the case of isothermal enclosure, however, the phenomena of flow and heat transfer are somewhat complex and interesting due to the stable thermal gradients and various circumstances resulted from four wall temperature conditions. As a dimensionless variable, to describe properly the flow and heat transfer phenomena in the isothermal enclosure, temperature difference ratio ${\Delta}T_v/{\Delta}T_H$ is newly introduced and this parameter seems to be appropriate in the analysis of results on the effect of stabilizing thermal gradient.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.4
• /
• pp.459-479
• /
• 2015

Two-dimensional steady laminar natural convection around a heat conducting and generating solid body inside a square enclosure with different thermal boundaries is performed. The mathematical model is governed by the coupled equation of mass, momentum and energy. These equations are discretized by finite volume method with power-law scheme and solved numerically by SIMPLE algorithm with under-relaxation technique. Effect of Rayleigh number, temperature difference ratio of solid-fluid, aspect ratio of solid-enclosure and the thermal conductivity ratio of solid-fluid are investigated numerically for Pr = 0.7. The flow and heat transfer aspects are demonstrated in the form of streamlines and isotherms respectively.

• Journal of the Korea Concrete Institute
• /
• v.17 no.4 s.88
• /
• pp.551-558
• /
• 2005

The setting and hardening of concrete is accompanied with nonlinear temperature distribution caused by development of hydration heat of cement. Especially at early ages, this nonlinear distribution has a large influence on the crack evolution. As a result, in order to predict the exact temperature history in concrete structures it is required to examine thermal properties of concrete. In this study, the convection heat transfer coefficient which presents thermal transfer between surface of concrete and air, was experimentally investigated with variables such as velocity of wind, curing condition and ambient temperature. At initial stage, the convection heat transfer coefficient is overestimated by the evaporation quantity. So it is essential to modify the thermal equilibrium considered with the boiling effect. From experimental results, the convection heat transfer coefficient was calculated using equations of thermal equilibrium. Finally, the prediction model for equivalent convection heat transfer coefficient including effects of velocity of wind, curing condition, ambient temperature and boiling effects was theoretically proposed. The convection heat transfer coefficient in the proposed model increases with velocity of wind, and its dependance on wind velocity is varied with curing condition. This tendency is due to a combined heat transfer system of conduction through form and convection to air. From comparison with experimental results, the convection heat transfer coefficient by this model was well agreed with those by experimental results.

• Journal of computational fluids engineering
• /
• v.10 no.3 s.30
• /
• pp.36-42
• /
• 2005

The present study numerically investigates the natural convection heat transfer in a 2-D square cavity containing a centered heat conducting body. Special emphasis is given to the influences of the Rayleigh number, the dimensionless conducting body size, and the ratio of the thermal diffusivity of the body to that of the fluid on the natural convection heat transfer in overall concerned region. The analysis reveals that the fluid flow and heat transfer processes are governed by all of them. Results for isotherms, vector plots and wall Nusselt numbers are reported for Pr = 0.71 and relatively wide ranges of the other parameters. Heat transfer across the cavity, in comparison to that in the absence of a body, are enhanced (reduced) in general by a body with a thermal diffusivity ratio less (greater) than unity. It is also found that the heat transfer attains a minimum as the body size is increased with a thermal diffusivity ratio greater than unity.

• Korean Journal of Materials Research
• /
• v.14 no.12
• /
• pp.870-875
• /
• 2004

Heat flow characteristics in the newly developed Wave Heat Sink were analyzed under natural and forced convections by Icepak program using the finite volume method. Temperature distribution and thermal resistance of Wave Heat Sink with/without air vent hole on the top of fin were compared with those of a commercial Al extruded heat sink(Intel Heat Sink). Under the natural convection, the maximum temperature was $45.1^{\circ}C$ in the air vent hole typed Wave Heat Sink, which was superior to that of Intel Heat Sink. The thermal resistance was $2.51^{\circ}C/W$ in the air vent hole typed Wave Heat Sink, and it changed to $2.65^{\circ}C/W\;and\;2.16^{\circ}C/W$ with changes of gravity direction and fin height, respectively. Under the forced convection, the maximum temperature became lower than that under the natural convection. In addition, the thermal resistance lowered in the air vent hole typed Wave Heat Sink with higher fin height and it decreased with increasing the air flux.

• Proceedings of the KIPE Conference
• /
• 2002.11a
• /
• pp.205-208
• /
• 2002

Thermal design is required with considering thermal stability to verify the reliability of electric power device with using IGBT. Numerical analysis is performed to analyzed the change in thermal resistance with respect to the various thermal density of heating element. Correlations between thermal resistance and heat generation density are established. With using these correlations, performance curve is composed with respect to the change in thermal resistance of cooling conditions for natural convection and forced convection. Thermal fatigue is occurred at the Inside and outside of IGBT by repeated heat load. The crack is occurred between base plate and ceramic substrate for the inside. When the crack length is 4mm, the failure is occurred. Therefore, Thermal design method considering thermal density, thermal fatigue resistance is presented on this study and it is expected to thermal design with considering life prediction.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.6
• /
• pp.1330-1338
• /
• 1989

본 연구에서는 직각밀폐용기의 수평경계면이 단열인 경우뿐만 아니라 등온 조건을 갖는 경우에 대해 실험적으로 연구하여 경계조건의 변화가 직각밀폐용기내 흐름 및 열전달에 미치는 영향, 특히 등온조건을 갖는 경우 수직 온도차에 따르는 안정온도 구배효과로 예상되는 흐름의 억제, 지연효과를 작종 물리적 변수들의 영향과 함께 세밀히 조사하였다.