• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.24 no.3
• /
• pp.247-255
• /
• 2012

Numerical calculations are carried out for the natural convection in a square enclosure with two hot cylinders induced by temperature difference between a cold outer rectangular cylinder and two hot circular cylinders. A two-dimensional solution for unsteady natural convection is obtained, using the immersed boundary method (IBM) to model two inner circular cylinders based on finite volume method, for different Rayleigh numbers varying over the range of $10^3$ to $10^5$. The study goes further to investigate the effect of the location of two cylinders on the heat transfer and fluid flow. The location of inner circular cylinders is changed vertically along the center-line of square enclosure. The changes of heat transfer quantities have been presented.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.139-150
• /
• 2005

The present study has numerically investigated two-dimensional laminar flow over a circular cylinder with a uniform planar shear, where the free-stream velocity varies linearly across the cylinder. Numerical simulations using the immersed boundary method are performed for the ranges of $50{\le}Re{\le}160,\;K{\le}0.2$, and B=0.1 and 0.05 where Re, K and B are the Reynolds number, the non-dimensionalized velocity gradient and the blockage ratio, respectively. Results show that the flow depends significantly on B as well as Re and K. It is found, especially, that the blockage effect accounts for some causes of apparent discrepancies among previous studies on the flow. With increasing K, the vortex shedding frequency and the mean drag stay nearly constant or slightly decrease whereas the mean lift, acting from the higher-velocity side to the lower, increases linearly. Flow statistics as well as instantaneous flow fields are presented to identify the characteristics of the flow and then to understand the underlying mechanism.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.22 no.11
• /
• pp.767-772
• /
• 2010

The incompressible viscous flow past a sphere under forced oscillation is numerically investigated at a Reynolds number of 300. The immersed boundary method is used to handle the sphere oscillating vertically to the streamwise direction. There are two important variables to characterize the oscillating state of a sphere. One is an oscillating amplitude normalized by the sphere diameter is set as a fixed number of 0.2. Another is the frequency ratio which is defined by $f_e/f_o$, where fe and fo are the excited frequency and the natural frequency of vortex shedding for the stationary sphere. In this study, three different frequency ratios of 0.8, 1.0 and 1.2 are considered. The results show a periodic flow with hairpin vortices shedding from upper and lower positions as well as vortical legs obliquely extended by oscillating motion of sphere. The enveloping vortical structure experience rupture twice in one period of oscillation. As the frequency of oscillation is increased, the vortical legs are getting shorter and eventually the hairpin vortices are much closer to the adjacent one.

• Journal of Advanced Research in Ocean Engineering
• /
• v.4 no.4
• /
• pp.146-162
• /
• 2018

The unsteady problem of a rigid body impact onto a floating plate is studied. Both the plate and the water are at rest before impact. The plate motion is caused by the impact force transmitted to the plate through an elastic layer with viscous damping on the top of the plate. The hydrodynamic force is calculated by using the second-order model of plate impact by Iafrati and Korobkin (2011). The present study is concerned with the deceleration experienced by a rigid body during its collision with a floating object. The problem is studied also by a fully-nonlinear computational-fluid-dynamics method. The elastic layer is treated with a moving body-fitted grid, the impacting body with an immersed boundary method, and a discrete-element method is used for the contact-force model. The presence of the elastic layer between the impacting bod- ies may lead to multiple bouncing of them, if the bodies are relatively light, before their interaction is settled and they continue to penetrate together into the water. The present study is motivated by ship slamming in icy waters, and by the effect of ice conditions on conventional free-fall lifeboats.

• Computers and Concrete
• /
• v.30 no.5
• /
• pp.323-337
• /
• 2022

This paper proposes a numerical model to simulate the rotational behavior of steel fiber in fresh cement-based materials in the presence of a magnetic field. The results indicate that as the aspect ratio of fiber increases, the required minimum magnetic field intensity to make fiber rotate in viscous fluid increases. The optimal magnetic field intensity is 0.03 T for aligning steel fiber in fresh cement-based materials to ensure that the applying time of the magnetic field can be conducted concurrently with the vibrating process to increase the aligning efficiency. The orientation factor of steel fiber in cement mortar can exceed 0.85 after aligning by 0.03 T of the uniform magnetic field. When the initial angle of the fiber to the magnetic field direction is less than 10°, the magnetic field less than 0.03 T cannot make the fiber overcome the yield stress of fluid to rotate. The coarse aggregate in steel fiber-reinforced concrete is detrimental to the rotation and alignment of the steel fiber. But the orientation factor of ASFRC under the 0.03T of the magnetic field can also exceed 0.8, while the orientation factor of SFRC without magnetic field application is around 0.6.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.12
• /
• pp.1240-1246
• /
• 2010

Enthalpy probe with the inner and outer diameters of 1.5 mm and 4.8 mm, respectively, is designed and used to measure the temperatures and velocities along the centerline of Ar arc plasma flow until the probe was destroyed. For this purpose, Ar arc plasma flow is generated by non-transferred type DC arc heater with the power level of 17 kW. From this experiment, it is shown that the designed enthalpy probe can measure the temperature and velocity of arc plasma flow up to 12,000 K and 600 m/s, respectively, without destroy of probe tip. In this extreme case, the arc plasma flow is calculated to transfer the heat flux of ${\sim}5{\times}10^7\;W/m^2$ to the probe based on the heat and thermal boundary equations near the forward stagnation point of a body immersed in arc plasma flow. Consequently, the designed enthalpy probe can measure the wide ranges of plasma temperatures, velocities and concentrations simultaneously, which are generated by various types of arc heaters within the heat flux ranges of $0{\sim}5{\times}10^7\;W/m^2$ on the probe tip.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.69-69
• /
• 2011

본 논문에서는 비정상 상태의 비압축성 유동장을 해석하기 위하여 물체맞춤격자방법이 아닌 가상경계법을 사용하였다. 가상경계법은 구조격자를 사용하여 구조물 경계면에서 Momentum Forceing을 사용하여 가상의 경계를 만들어 유동장을 해석하는 방법이다. Navier-Stoke 방정식의 수치 이산화 방법으로 Kim et al(1985)이 사용한 Fractional Step Method(FSM)을 사용하였다. 시간에 대하여 semi-implicit FSM를 사용하였고, 확산항에 대해서는 2차 정확도의 Crank-Nicolson Method를 대류항은 3차 정확도의 Runge-Kutta Method를 사용 하였다. 본 연구에서는 가상경계법을 이용한 유동장 해석이 교량 단면에 대하여 수치해석이 가능한지 검토하였다. 가상경계법은 현재 많은 연구가 유선형의 구조물에 대하여 수행되어 오고 있다. 교량 단면과 같은 각 진 구조물에 대한 검토는 아직 미비한 실정이다. 가상경계법에서 다루고 있는 구조물 경계면에서의 Momentum Forcing 방법이 유선형의 구조물에 맞추어 연구가 진행되었기 때문이다. 먼저 본 연구의 프로그램을 검증하기 위하여 원형 실린더에 대하여 가상경계법을 적용한 결과 Re 수 200에서 Strouhal Number, 양력계수, 항력계수를 이전 연구 결과와 비교하였다. Williamson(1988)과 Zhang(1995)의 연구결과와 유사한 결과를 얻을 수 있다. 그리고 교량의 단면과 같은 각진 구조물(Bluff Body)에 대하여 가상경계법 적용하였다. 본 논문의 연구에서 평가 대상으로 하고 있는 2차원 교량 단면에 대하여 유동장 해석을 하였다. 본 논문에서 정량적인 유체력과 유동장에 대한 비교 및 검토가 이루어지지 못했지만 압력장과 유선의 형태가 이론적인 값을 벗어나지 않고 있는 것으로 확인 되었다. Re 수 2700에서 전산 해석을 수행하였으며, 교량 단면 주위의 압력계수와 박리현상 그리고 후류에서의 Vortex shedding 현상이 모두 적절한 분포가 나타나는 것을 확인할 수 있었다. 따라서 가상경계법을 이용하여 각진 구조물에 대한 주위 유동장해석에 대한 가능성을 확인하였으며, 풍동실험과의 결과비교를 통하여 가상경계법을 이용하여 교량 단면 주위의 유동장 해석 결과를 정량적으로 비교할 것이다.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.18 no.2
• /
• pp.55-62
• /
• 2004

In this paper, the diagnostic expert system for silicone oil-filled transformer is developed using dissolved gas analysis(DGA). There are many diagnostic methods for diagnostic oil-immersed transformer. But DGA is used to the proposed expert system since it has been verified that DGA is very efficient diagnostic method for transformer. In addition, it is resonable that fuzzy rule, degree of inclusion and fuzzy measure must be considered to handle the uncertainty nature of gas boundary and rules. The proposed expert system consists of knowledge base module, inference engine module and human-machine interface(HMI) module. The knowledge base module consists of the knowledge using the rule. The inference engine module is used to the fuzzy rule. The history of the transformer gas data is managed by the database. the effect of the proposed expert system is verified by case studies.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.577-589
• /
• 2005

The present study has numerically investigated two-dimensional laminar flow over a steadily rotating circular cylinder with a uniform planar shear, where the free-stream velocity varies linearly across the cylinder. It aims to find the combined effect of rotation and shear on the flow. Numerical simulations using the immersed boundary method are performed for the ranges of $-2.5{\le}\alpha{\le}2.5$ and $0{\le}K{\le}0.2$ at a fixed Reynolds number of Re=100, where a and K are respectively the dimensionless rotational speed and velocity gradient. Results show that the positive shear, with the upper side having the higher free-stream velocity than the lower one, favors the effect of the counter-clockwise rotation $(\alpha<0)$ but countervails that of the clockwise rotation $(\alpha>0)$. Accordingly, the absolute critical rotational speed, below which vortex shedding occurs, decreases with increasing K for $(\alpha>0)$, but increases for $\alpha>0$. The vortex shedding frequency increases with increasing \alpha (including the negative) and the variation becomes steeper with increasing K. The mean lift slightly decreases with increasing K regardless of the rotational direction. However, the mean drag and the amplitudes of the lift- and drag-fluctuations strongly depend on the direction. They all decrease with increasing K for $\alpha>0$, but increase for $\alpha<0$. Flow statistics as well as instantaneous flow folds are presented to identify the characteristics of the flow and then to understand the underlying mechanism.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.8
• /
• pp.739-747
• /
• 2010

Numerical calculations are carried out for evaluating the natural convection induced by the temperature difference between a hot inner circular cylinder and a cold outer square enclosure. A two-dimensional solution for unsteady natural convection is obtained by using the finite volume method to model an inner circular cylinder that was designed by using the immersed boundary method (IBM) for a Rayleigh number of $10^7$. In this study, we investigate the effect of the location ($\delta$) of the inner cylinder, which is located along the vertical central axis of the outer enclosure, on the heat transfer and fluid flow. The natural convection changes from unsteady to steady state depending on the $\delta$. The two critical lower bound and upper bound positions are ${\delta}_{C,L}$ = 0.05 and ${\delta}_{C,U}$ = 0.18, respectively. Within these defined bounds, the thermal and flow fields are in steady state.