• Tunnel and Underground Space
• /
• v.12 no.3
• /
• pp.142-151
• /
• 2002

Claesson(2001)'s analytical solution, and two numerical models with Dirichlet and Neuman interior boundary condition respectively were investigated to estimate the transient temperature distribution with distances from the Taejon underground food cold storage pilot cavern. Claesson's solution, which is based on constant temperature boundary condition at the rock wall during a temperature decline step, showed relatively good agreement with temperature measurements in the rock mass in order of average error difference, 0.89$\^{C}$ without any adjustments on laboratory thermal properties to represent the rock mass. For the numerical model with heat flux through the rock wall, a boundary condition setting technique was newly proposed to overcome the difficulty of prescribing variable convective heat tranfer coefficient and far-field air temperature inside the cavern as they may be certainly changed according to the cooling-down time. The results showed also good agreement with measurements in order of average error difference, 1.58$\^{C}$, and were compared to those of the numerical model with fixed temperature at the rock wall. Finally, the most proper procedure to precisely predict the temperature profile around a cavern was proposed as a series of analysis steps including an analytical exact solution and numerical models.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.725-733
• /
• 1990

An analysis of laminar mixed convection flow adjacent to the inclined flat surface which is subjected to a uniform temperature in a uniform free stream is performed. Nonsimilar boundary layed equation are derived by using the mixed convection parameters such that smooth transition from the purely forced convection limit to the purely free convection limit is possible. The governing equations are solved by a finite difference method using the coupled box scheme of sixth order. Numerical results are presented for prandtl numbers of 0.7 and 7 with the angle of inclination ranging from 0 to 90 degree from the vertical. The velocity distributions for the buoyancy assisting flow exhibit a significant overshoot above the free stream value in the region of intense mixed convection and the velocity field is found to be more sensitive to the buoyancy effect than the temperature field. The separation point near the wall was obtained for the buoyancy opposing flow. The local Nusselt number increases for buoyancy assisting flow and decreases for opposing flow with increasing value of the local Grashoff number in the mixed convection parameter. For large Prandtl number, the Nusselt number and the friction factor decrease significantly near the separation point. Present numerical predictions are in good agreement with recent experimental results by Ramachandran.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.8 no.2
• /
• pp.171-177
• /
• 1984

균일하게 분포된 내부발열을 갖는 유체가 든 경사진(수평에서 45.deg.까지)정사각형 단면의 밀폐 공간 내에서의 2차원 자연대류 유동 및 열전달에 관한 수치적인 연구가 수행되었다. 4개 벽면온 도가 동일한 경계조건에서 내부발열로 인한 자연대류 유동이 Rayleigh수 1.5*$10^{5}$까지는 층 류유동 영역의 가정하에서 수치적으로 수렴되었다. 경사진 밀폐 공간에서의 유동형태 및 온도분 포는 수평인 경우에 비하여 윗쪽 벽면 근처에서 그 상이점이 많이 나타났다. 경사각도가 증가함 에 따라 평균 열전달율이 아랫쪽 벽과 오른쪽 벽에서 증가하였고 왼쪽 벽에서 감소하였으며 윗 쪽 벽에서는 거의 일정하였다.다.

• Fire Science and Engineering
• /
• v.19 no.4 s.60
• /
• pp.18-25
• /
• 2005

In this study, numerical analysis of two-dimensional unsteady natural convection of air in a square enclosure heated from below, was performed as a basic research of fire science. SIMPLE algorithm was used to the pressure term of momentum equations in the numerical analysis. The numerical analysis were studied for the two model cases and two heat conditions, respectively, which are different with insulation of enclosures and position of heat applied. Also, the ceiling temperatures of enclosure were measured to compare the accuracy of numerical analysis, and it is found that the temperature predicted by numerical analysis were agreed well with the measurements. Streamline and isotherm of the each model case were acquired for each time step.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.1-11
• /
• 1999

The three-dimensional laminar mixed convection flow between the conductive printed circuit boards. on which the heat generating rectangular blocks are uniformly distributed, has been examined in the present study. The flow and heat-transfer characteristics are assumed to be pseudo periodic in the streamwise direction and symmetric in the cross-stream direction. Using an algorithm of SIMPLER, the continuity equation. the Navier-Stokes equations and the energy equation are solved numerically in the three-dimensional domain Inside the channel. The convective derivative terms are discretized by the QUICK scheme to accurately capture the flow field. The flow and the heat transfer characteristics are thoroughly examined for various Re and Gr.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.2
• /
• pp.39-45
• /
• 2001

With uniform heat generation from the inner surface of the cylindrical heater placed in a cross flow boundary condition, heat flow that is conducted along the wall of the heater creates a non-isothermal surface temperature and non-uniform heat flux distribution. In the present investigation, the effects of circumferential wall heat conduction on convection heat transfer is investigated for the case of forced convection around horizontal circular tube in cross flow of air. The wall conduction number which can be deduced from the governing energy equation should be used to express the effect of circumferential wall heat conduction. It is demonstrated that the circumferential wall heat conduction influences local Nusselt numbers of one-dimensional and two-dimensional solutions.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.4 no.3
• /
• pp.217-232
• /
• 1992

A Numerical study on two-dimensional laminar natural convection with and without surface radiation in fully or partially open square cavity was performed. The cavity has one vertical heated wall facing a vertical opening and two horizontal insulated walls. The pressure boundary condition was applied to the opening instead of the velocity boundary condition. The results of this study showed that the increase of partition length decreased the convective and the radiative Nusselt numbers. It was also found that the increase of wall emissivity decreased the convective Nusselt numbers but increased the radiative Nusselt numbers.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.11 no.5
• /
• pp.14-22
• /
• 2007

A reversed trapezoidal fin with variable lateral surface slope is optimized using a two-dimensional analytic method. For a fin base boundary condition, convection from the inside fluid to the inside wall and conduction from the inside wall to the fin base are considered. Heat loss from the fin tip surface is not ignored. The maximum heat loss at the practical fin length, the corresponding optimum fin efficiency, fin length and fin base height are presented as a function of the fin inside and outside convection characteristic numbers. One of the results shows that the optimum fin shape becomes 'fatter and shorter' as the ratio of fin tip height to base height increases.

• Proceedings of the KIEE Conference
• /
• 2011.07a
• /
• pp.1584-1585
• /
• 2011

분자의 전리(ionization)작용과 이온의 해리(dissociation)작용에 의한 유전체 내 방전해석을 위해서 전자, 양이온, 음이온의 생성과 소멸, 전자부착, 재결합 과정을 포함한 전하연속방정식, 전계에 의한 푸아송 방정식, 유전체 온도에 관한 열확산 방정식을 결합하여 해석하였다. 전계 방출 조건과 열전자 방출 조건이 경계조건으로 부여되었고 에너지 최소화 정류조건을 따르는 유한요소법(finite element method)을 이용하여 해석하였다. 비교적 작은 값을 가지는 확산성을 무시하고 대류성에 역점을 두어 발생하는 수치적 불안정은 인공확산항(artificial diffusion technique)을 도입하여 안정화하였다. 본 논문에서는 IEC standard 60897의 표준규격에 따른 2차원 축대칭 침-구 전극에 적용하여 제안된 다중물리해석기법의 타당성을 입증하였다.