• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.5
• /
• pp.1253-1261
• /
• 1993

Numerical analysis and measurements of the velocity and temperature distributions in buoyancy assisting laminar mixed convection flow over a vertically located, two-dimensional backward-facing step are reported. Laser-Doppler Velocimeter and Constant Temperature Anemometer operated in constant current were used to measure simultaneously the velocity and temperature distributions in the recirculation region downstream of the step. The reattachment length was measured by using flow visualization technique for different inlet velocities, wall temperatures and step heights. While the reattachment length $X_r$ increases as the inlet velocity or step height increase, it decreases as the buoyancy force increases, causing the size of the recirculation region to decrease. For the experimental range of $Gr_s$/$Re_{s}^{2}$$\times$$10^3$<17, a correlation equation for the reattachment length can be given by $X_{r}=1.05(2.13+0.021 Re_{s})exp$ $(-33.7_s^{-0.186}/Gr_{s}/Re_{s}^2).$ The Nusselt number is found to increase and the location of its maximum value moves closer to the step as the buoyancy force increases. The location of the maximum Nusselt number occurs downstream of the reattachment point, and distance between the reattachment point and the location of the maximum Nusselt mumber increases as the buoyancy force increases. Computational prediction agrees favorably well with measured results.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.450-455
• /
• 2001

An analytical study is made of transient adjustment process of an initially stationary, stably-stratified fluid in a square. The boundary walls are highly-conducting. The overall Rayleigh number $R_a$ is large. Considerations are given to both opposing (${\delta}w/{\delta}T>0$) and cooperating (${\delta}w/{\delta}T<0$) configurations. The flow character in opposing configuration can be classified into (a) a forced-convection dominaut mode (${\delta}w/{\delta}T>1/\sqrt{2}$), (b) a buoyancy-convection dominant mode ($0<{\delta}w/{\delta}T<1/\sqrt{2}$), and (c) a static mode (${\delta}w/{\delta}T{\cong}1/\sqrt{2}$). Global evolutionary processes are depicted. and physical rationalizations are provided.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.46-53
• /
• 2016

Direct numerical simulation is conducted to observe the behavior of microbubbles in isotropic turbulence. Navier-Stokes equation and the motion of equation for microbubbles are solved with periodic boundary condition in a cube domain. Vorticity contour, enstrophy ratio, relative reduction of bubble rise velocity, and the closest distance of particles are investigated for various Stokes numbers and gravity factors to understand clustering of microbubbles. Also, clustering due to the effect of the lift force is investigated.

• Journal of Biosystems Engineering
• /
• v.18 no.2
• /
• pp.144-157
• /
• 1993

환기중인 실험축사내에서 가축의 현열과 환기공기의 온도차에 의한 열부력(熱浮力)(thermal buoyancy)이 공기유동 및 온도분포에 미치는 영향을 구명(究明)하기 위하여 TEACH 컴퓨터프로그램($k-{\varepsilon}$ 난류모형 및 SIMPLE계열 Algorithm)을 Curvilinear Coordinates에 맞게 변형하였다. 계산한 축사내 공기유통 및 온도분포의 유의성(有意性) 검증은 Boon(1978)의 실험결과를 이용하였다. 열부력의 크기에 따른 유동의 변화를 관찰하기 위하여 유입공기의 온도를 $17^{\circ}C$와 $10^{\circ}C$ 두 수준으로 입력하였으며, 가축의 현열플릭스(flux)는 실내온도에 따라 변화하므로 유압공기의 온도가 $17^{\circ}C$일 때는 130W/$m^2$, $10^{\circ}C$일 때는 170W/$m^2$을 경계조건으로 입력하였다. 예측한 공기유동의 형태는 실험값(Boon, 1978)과 비교하여 대체로 만족할만한 결과를 얻었다. 그러나 유입공기의 온도가 $10^{\circ}C$인 경우, 예측 공기유동은 실험 유동형태와 차이가 있었다. 즉, 실험에서는 수평슬롯으로 유입된 공기가 바로 아래로 굴절되어 유동(流動)하였으나, 계산의 결과는 일정 거리로 수평방향으로 유동하다가 아래로 굴절하였다. 이런 유동의 차이는 경험적으로 열부력(熱浮力)에 민감하게 반응하지 않는 k-${\varepsilon}$ 난류(亂流)모형의 적용이 원인이 되거나 실험의 부적절한 수행이 원인이 될 수도 있다. 이 유동(流動)의 Reynolds 수(數) (Re)는 약 3,300, 수정Ar수(修正Ar數)(Corrected Archimedes Number : $Ar_c$)64로써, $Ar_c$ <30 이거나 $Ar_c$ >75이면 유입공기의 제트는 수평유동한다는 Randall & Battams(1979)의 연구결과와는 일치하였다. 그러나 공기제트의 굴절은 유동의 특성이 같다하더라도 유체의 성질, 축사의 기하학적 형태에 따라서 매우 민감하게 반응하므로 실제 실험을 통한 재검정과정을 거쳐야 할 것으로 판단된다. Fig. 9와 Fig. 10의 기하학적 형태의 지점별 예측온도와 측정온도(Boon, 1978)와의 편차는 대부분의 지점에서는 $1^{\circ}C$ 미만으로 상당히 정확하였으며, 최대의 온도차는 Fig. 10의 지점 13에서 $1.7^{\circ}C$이었다.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.121-126
• /
• 1998

A numerical optimization procedure is developed to find the inlet velocity profile that yields the most uniform epitaxial layer in a vertical MOCVD reactor. It involves the solution of fully elliptic equations of motion, temperature, and concentration; the finite volume method based on SIMPLE algorithm has been adopted to solve the Navier-Stokes equations. The overall optimization process is highly nonlinear and has been efficiently treated by the sequential linear programming technique that breaks the non-linear problem into a series of linear ones. The optimal profile approximated by a 6th-degree Chebyshev polynomial is very successful in reducing the spatial non-uniformity of the growth rate. The optimization is particularly effective to the high Reynolds number flow. It is also found that a properly constructed inlet velocity profile can suppress the buoyancy driven secondary flow and improve the growth-rate uniformity.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.220-225
• /
• 2007

This paper is a continuation of the recent development on the hermite-based divergence free basis function and deals with a non-isothermal fluid flow thru the buoyancy driven flow in a square cavity with temperature difference across the two sides. The basis functions for the velocities consist of the hermite function and its curl. However, the basis for the temperature are the hermite function and its gradienst. Hence, the number of degrees of freedom at a node becomes 6, which are the stream function, two velocities, the temperature and its x- and y-derivatives. Numerical results for the streamlines, the temperatures, the x-velocities and the y-velocities show good agreements with those of De vahl Davis[7].

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1657-1662
• /
• 2003

In this study, to research characteristics of heat flow of magnetic fluid, it's studied about numerical and experimental method of natural convections change and characteristics of heat transfer in Concentric double pipe annuli as analysis model. In the result, natural convection of magnetic fluid is controlled by direction and strength of the impressed magnetic field. Especially, according to average Nusselt number, heat transfer is the smallest on the balancing point between body force and buoyancy.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.64-71
• /
• 2016

This research is to numerically investigate the convective cooling performance in the Disk brake. Research concentrates on the heat transfer coefficient and cooling performance which are selected with cooling local locations. Cooling performance of the Hole disk has been compared by Ventilated Disk. According to the results of heat transfer on the disk brake, activated velocity distributions more appear in the Hole disk. This is due to the fact that a number of hole units have exactly 120 on the surface of the hole disk. Therefore, velocity distributions of hole disk brake is better activated than Ventilated disk. According to the calculations of Nusselt number between surface and atmosphere in the interested cooling area, average value of cooling effect has been increased 13.5% by the hole disk at driving of speed 65 km/h situation and grown 18% by the hole disk at driving speed of 100 km/h. Due to the flow of air through the hole route, cooling performance of the hole disk was very excellent. In addition, cooling effect on edge of the bottom is better than the vicinity of center.