• Korean Journal of Agricultural and Forest Meteorology
• /
• v.11 no.4
• /
• pp.192-205
• /
• 2009

The unexpected wind over the Mt. Hwawang on 9 February 2009 was deadly when many spectators were watching a traditional event to burn dried grasses and the fire went out of control due to the wind. We analyzed the fatal wind based on wind flow simulations over a digitized complex terrain of the mountain with a localized heating area using a three dimensional computational fluid dynamics model, CFD_NIMR_SNU (Computational Fluid Dynamics_National Institute of Meteorological Research_Seoul National University). Three levels of fire intensity were simulated: no fire, $300^{\circ}C$ and $600^{\circ}C$ of surface temperature at the site on fire. The surface heat accelerated vertical wind speed by as much as $0.7\;m\;s^{-1}$ (for $300^{\circ}C$) and $1.1\;m\;s^{-1}$ (for $600^{\circ}C$) at the center of the fire. Turbulent kinetic energy was increased by the heat itself and by the increased mechanical force, which in turn was generated by the thermal convection. The heating together with the complex terrain and strong boundary wind induced the unexpected high wind conditions with turbulence at the mountain. The CFD_NIMR_SNU model provided valuable analysis data to understand the consequences of the fatal mountain fire. It is suggested that the place of fire was calm at the time of the fire setting due to the elevated terrain of the windward side. The suppression of wind was easily reversed when there was fire, which caused updraft of hot air by the fire and the strong boundary wind. The strong boundary wind in conjunction with the fire event caused the strong turbulence, resulting in many fire casualties. The model can be utilized in turbulence forecasting over a small area due to surface fire in conjunction with a mesoscale weather model to help fire prevention at the field.

• Korean Journal of Radiology
• /
• v.2 no.1
• /
• pp.28-36
• /
• 2001

Objective: To provide a systematic overview of the effects of various parameters on contrast enhancement within the same population, an animal experiment as well as a computer-aided simulation study was performed. Materials and Methods: In an animal experiment, single-level dynamic CT through the liver was performed at 5-second intervals just after the injection of contrast medium for 3 minutes. Combinations of three different amounts (1, 2, 3 mL/kg), concentrations (150, 200, 300 mgI/mL), and injection rates (0.5, 1, 2 mL/sec) were used. The CT number of the aorta (A), portal vein (P) and liver (L) was measured in each image, and time-attenuation curves for A, P and L were thus obtained. The degree of maximum enhancement (Imax) and time to reach peak enhancement (Tmax) of A, P and L were determined, and times to equilibrium (Teq) were analyzed. In the computed-aided simulation model, a program based on the amount, flow, and diffusion coefficient of body fluid in various compartments of the human body was designed. The input variables were the concentrations, volumes and injection rates of the contrast media used. The program generated the time-attenuation curves of A, P and L, as well as liver-to-hepatocellular carcinoma (HCC) contrast curves. On each curve, we calculated and plotted the optimal temporal window (time period above the lower threshold, which in this experiment was 10 Hounsfield units), the total area under the curve above the lower threshold, and the area within the optimal range. Results: A. Animal Experiment: At a given concentration and injection rate, an increased volume of contrast medium led to increases in Imax A, P and L. In addition, Tmax A, P, L and Teq were prolonged in parallel with increases in injection time The time-attenuation curve shifted upward and to the right. For a given volume and injection rate, an increased concentration of contrast medium increased the degree of aortic, portal and hepatic enhancement, though Tmax A, P and L remained the same. The time-attenuation curve shifted upward. For a given volume and concentration of contrast medium, changes in the injection rate had a prominent effect on aortic enhancement, and that of the portal vein and hepatic parenchyma also showed some increase, though the effect was less prominent. A increased in the rate of contrast injection led to shifting of the time enhancement curve to the left and upward. B. Computer Simulation: At a faster injection rate, there was minimal change in the degree of hepatic attenuation, though the duration of the optimal temporal window decreased. The area between 10 and 30 HU was greatest when contrast media was delivered at a rate of 2 3 mL/sec. Although the total area under the curve increased in proportion to the injection rate, most of this increase was above the upper threshould and thus the temporal window was narrow and the optimal area decreased. Conclusion: Increases in volume, concentration and injection rate all resulted in improved arterial enhancement. If cost was disregarded, increasing the injection volume was the most reliable way of obtaining good quality enhancement. The optimal way of delivering a given amount of contrast medium can be calculated using a computer-based mathematical model.

• Journal of Bio-Environment Control
• /
• v.15 no.4
• /
• pp.296-305
• /
• 2006

The heterogeneity of crop transpiration is important to clearly understand the microclimate mechanisms and to efficiently handle the water resource in greenhouses. A computational fluid dynamic program (Fluent CFD version 6.2) was developed to study the internal climate and crop transpiration distributions of greenhouses. Additionally, the global solar radiation model and a crop heat exchange model were programmed together. Those models programmed using $C^{++}$ software were connected to the CFD main module using the user define function (UDF) technology. For the developed CFD validity, a field experiment was conducted at a $17{\times}6 m^2$ plastic-covered mechanically ventilated single-span greenhouse located at Pusan in Korea. The CFD internal distributions of air temperature, relative humidity, and air velocity at 1m height were validated against the experimental results. The CFD computed results were in close agreement with the measured distributions of the air temperature, relative humidity, and air velocity along the greenhouse. The averaged errors of their CFD computed results were 2.2%,2.1%, and 7.7%, respectively.

• Korean Journal of Ecology and Environment
• /
• v.42 no.1
• /
• pp.85-94
• /
• 2009

This study examines the operational effectiveness of a Convection Current Aeration System (CCAS) in reservoir. CCAS was run from June, 2008 when the thermocline begun forming in the reservoir. This paper reviews the influence of stratification, dissolved oxygen dynamics and temperature in the lake's natural state from June to October 2008. The survey was done on a week basis. Upwelling flow effects a radius of $7{\sim}10m$ at a surface directly and was irrelevant to the strength of thermocline. On the other hand, it was affected the number of working days, and strength of thermocline at vertical profiles of the reservoir. Longer CCAS run, the deeper was the vertical direct flow area. However it didn't break the thermocline during summer season of 2008. The operating efficiency of the CCAS in the reservoir depends on hydraulics and meteological conditions. Computational Fluid Dynamics (CFD) is a very useful tool for evaluating the operating efficiency of fluid dynamics. The geometry for CFD simulation consists of a cylindrical vessel 25 m radius and 40 m height. The CCAS is located in center of domain. The non-uniform tetrahedral meshes had a bulk of the geometry. The meshes ranged from the coarse to the very fine. This is attributed to the cold water flowing into the downcomer and rising, creating a horizontal flow to the top of the CCAS. The result of CFD demonstrate a closer agreement with surveyed data for temperature and flow velocity. Theoretical dispersion volume were calculated at 8m depth, 120 m diameter working for 30 days and 10 m depth, 130 m diameter working for 50 days.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.3
• /
• pp.128-137
• /
• 2013

Due to the global warming and air pollution, interest in renewable energies has increased in recent years. In particular, the crisis of the depletion of fossil energy resources in the near future has accelerated the renewable energy technologies. Among the renewable energy resources, oceans covering almost three-fourths of earth's surface have an enormous amount of energy. For this reason, various approaches have been made to harness the tremendous energy potential. In order to achieve two purposes: to improve harbor water quality and to use wave energy, this study proposed a sea water exchange structure applying an Oscillating Water Column (OWC) wave generation system that utilizes the air flow velocity induced by the vertical motion of water column in the air chamber as a driving force of turbine. In particular, the airflow velocity in the air chamber was estimated from the time variations of water surface profile computed by using 3D-NIT model based on the 3-dimensional irregular numerical wave tank. The relationship of the frequency spectrums between the computed airflow velocities and the incident waves was analyzed. This study also discussed the characteristics of frequency spectrums in the air chamber according to the presence of the structure, wave deformations by the structure, and the power of the water and air flows were also investigated. It is found that the phase difference exists in the time series data of water level fluctuations and air flow in the air chamber and the air flow power is superior to the fluid flow power.

• Ocean and Polar Research
• /
• v.35 no.3
• /
• pp.249-258
• /
• 2013

Vertical and horizontal mixing processes in the ocean mixed layer determine sea surface temperature and temperature variability. Accordingly, simulating these processes properly is crucial in order to obtain more accurate climate simulations and more reliable future projections using an ocean general circulation model (OGCM). In this study, by using Modular Ocean Model version 4 (MOM4) developed by Geophysical Fluid Dynamics Laboratory, the upper ocean temperature and mixed layer depth were simulated with two different vertical mixing schemes that are most widely used and then compared. The resultant differences were analyzed to understand the underlying mechanism, especially in the Tropical Pacific Ocean where the differences appeared to be the greatest. One of the schemes was the so-called KPP scheme that uses K-Profile parameterization with nonlocal vertical mixing and the other was the N scheme that was rather recently developed based on a second-order turbulence closure. In the equatorial Pacific, the N scheme simulates the mixed layer at a deeper level than the KPP scheme. One of the reasons is that the total vertical diffusivity coefficient simulated with the N scheme is ten times larger, at maximum, in the surface layer compared to the KPP scheme. Another reason is that the zonal current simulated with the N scheme peaks at a deeper ocean level than the KPP scheme, which indicates that the vertical shear was simulated on a larger scale by the N scheme and it enhanced the mixed layer depth. It is notable that while the N scheme simulates a deeper mixed layer in the equatorial Pacific compared to the KPP scheme, the sea surface temperature (SST) simulated with the N scheme was cooler in the central Pacific and warmer in the eastern Pacific. We postulated that the reason for this is that in the central Pacific atmospheric forcing plays an important role in determining SST and so does a strong upwelling in the eastern Pacific. In conclusion, what determines SST is crucial in interpreting the relationship between SST and mixed layer depth.

• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.05a
• /
• pp.249-254
• /
• 2003

분석에 이용된 격자망은 동곡, 고로, 미성, 병천, 효령 및 무성지점에 각각 12개, 30개, 45개, 76개, 46개 및 1265개의 조격자를 구성하였으며 하천의 수로방향 및 경사형태를 세분화하기 위하여 각 지점에 대하여 8개, 24개, 24개, 44개, 12개 및 64개의 세격자로 분할하였다. AGNPS를 이용한 첨두유량의 모의발생 결과치가 동곡, 고로, 미성, 병천, 효령 및 무성지점에서 측정값과 비교하여 각 호우사상별로 상대오차가 1.0~25.0%, 4.0~27.0%, 7.0~29.2%, 2.0~23.9%, 3.0~25.0% 및 3.6~21.0%의 차이를 나타내었다. 분석결과에서 AMCII조건에서는 관측치와 분석결과치가 유사하게 나타났으나 AMCI조건에 대해서는 상대적으로 작은 값을 보였으며 AMCIII조건에서는 다소 큰 값으로 분석되었다. SCS방법에서 제안하는 AMC조건별 CN값을 우리 실정에 적합하도록 수정 보완하기 위한 수정 유출곡선지수 $CN_{m}$ /I과 $CN_{m}$/III을 재구성하였으며, 여기에 적용되는 수정 유출 곡선지수식의 계수 a를 추정한 결과, 기왕에 발표된 연구결과와 거의 일치된 경향을 나타내었다. 제안된 수정 CN식을 이용하여 산정한 결과치와 관측치는 거의 유사하게 나타났다. AGNPS모델에 의한 유출량 산정에 있어 수문학적 토양피복형수(CN)의 결정을 위하여 선행강우량과 토양의 공극율 및 지형인자인 각 셀마다의 유역경사를 이용하여 관계식(CN =f($X_1$, $X_2$, $X_3$))을 유도하였으며, 분석 결과에서 CN이 선행강우량과 가장 밀접한 관계가 있음을 알 수 있었으며 유역경사, 토양의 공극율 순으로 나타났다..88mg/$\ell$~의 범위로 나타났다. 무태교 지점에서의 총인의 농도는 0.52mg/$\ell$~0.99mg/$\ell$~의 범위이었다. 신천에 금호강물을 혼합한 이후에도 부유물질, 생화학적산소요구량, 암모니아태 질소, 총인 등의 농도가 개선되지 않았다. 즉 금호강물의 혼합은 신천수질환경사업소에서 배출되는 방류수에 함유되어 있을 2차 오염물질의 희석이라는 이점외의 수질개선효과는 확인되지 않았다.l years and a new type of transfer crane has been developed. Design concepts and control methods of a new crane will be introduced in this paper.and momentum balance was applied to the fluid field of bundle. while the movement of′ individual material was taken into account. The constitutive model relating the surface force and the deformation of bundle was introduced by considering a representative prodedure that stands for the bundle movement. Then a fundamental equations system could be simplified considering a steady state of the process. On the basis of the simplified model, the simulation was performed and the results co

• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.05a
• /
• pp.292-296
• /
• 2003

Autoclaving처리를 하지 않은 키토산과 autoclaving처리한 키토산과의 중금속 흡착실험을 해 본 결과 다음과 같은 결과를 도출 할 수 있었다. 1) 먼저 중금속간의 흡착능을 고찰하기 위하여 Langmuir와 Freundlich 흉착등온식에 적용시켜 본 결과 Freundlich 흡착등온식 보다는 Langmuir 흡착등온식이 보다 적합한 것으로 나타났다. 2) 시간에 따라서 autoclaving 처리한 키토산의 중금속 제거 가능성과 그 효율을 검토하기 위하여 각 중금속간의 Langmuir 흡착등온식을 이용하여 $q_{max}$를 나타냈다. 15 min > 60 min > 0 min 순으로 모든 중금속 제거 실험에서 15 min 동안 autoclaving 처리한 키토산의 중금속 흡착량이 가장 높은 것으로 나왔다. 그러므로 기존의 키토산 보다는 15 min동안 autoclaving 처리한 키토산이 중금속 제거에 더 좋은 흡착제로써 역할을 할 수 있을 것으로 판단된다. 3) 키토산을 이용한 중금속 제거에서는 $Pb^{2+}$ > $Cd^{2+}$ > $Cu^{2+}$ > $Cr^{3+}$순서로 제거가 되었다. 여러 연구자들의 실험 결과를 종합해 볼 때 Pb$^{2+}$가 중금속 중에 제거가 잘 된다는 연구 결과가 많이 발표 되었으며, $Cd^{2+}$, $Cu^{2+}$, $Cr^{3+}$의 경우에는 흡착제의 종류에 따라서 제거되는 순서가 다르다는 연구 결과들이 보고 되고 있다. 그러나 어떠한 이유로 중금속의 제거에 차이가 있는지에 대해서는 명확한 결론이 내려져 있지 않는 실정이다. 이러한 중금속간의 경쟁적인 관계에 대해 더 많은 세밀한 연구가 이루어져야 할 것 같다.는 0.52mg/$\ell$~0.99mg/$\ell$~의 범위이었다. 신천에 금호강물을 혼합한 이후에도 부유물질, 생화학적산소요구량, 암모니아태 질소, 총인 등의 농도가 개선되지 않았다. 즉 금호강물의 혼합은 신천수질환경사업소에서 배출되는 방류수에 함유되어 있을 2차 오염물질의 희석이라는 이점외의 수질개선효과는 확인되지 않았다.l years and a new type of transfer crane has been developed. Design concepts and control methods of a new crane will be introduced in this paper.and momentum balance was applied to the fluid field of bundle. while the movement of′ individual material was taken into account. The constitutive model relating the surface force and the deformation of bundle was introduced by considering a representative prodedure that stands for the bundle movement. Then a fundamental equations system could be simplified considering a steady state of the process. On the basis of the simplified model, the simulation was performed and the results could be confirmed by the experiments under various conditions.뢰, 결속 등 다차원의 개념에 대한 심도 깊은 연구와 최근 제기되고 있는 이론

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.3
• /
• pp.279-285
• /
• 2011

We perform a numerical study to determine the time of onset of natural convection in a transient hot wire (THW) device for measuring the thermal conductivity of nanofluids. The samples used in this simulation are water-based $Al_2O_3$ nanofluids with volume fractions of 1%, 4%, and 10%, and the properties are calculated by theoretical models and experimental correlations. The THW apparatus using coated wire is modeled by the control-volume-based finite difference method, and the start of natural convection is determined by observing the temperature rise of the wire under a gravity field. The onset time is 11.5 s for water and 41.6 s for water-based $Al_2O_3$ nanofluids predicted by Maxwell thermal conductivity model with a 10% volume fraction. We confirm that the onset time of natural convection of nanofluids in the cylinder increases with the nanoparticle volume fraction. We suggest a correlation for predicting the onset time on the basis of the numerical results. Finally, it is shown that the measurement error due to natural convection is negligible if the measurement using the transient hot wire method is completed before the onset of natural convection in the base fluid.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.1
• /
• pp.77-82
• /
• 2010

Our study aims to understand the flow of liquid in an open-top rectangular microchannel that can be used in micro total analysis systems ($\mu$-TAS) because it has advantages in terms of light transmission and energy efficiency. We measured the liquid velocity using particle tracking technique and conducted a simulation with computational fluid dynamics by altering the area of channel cross section and channel length for the capillary-driven flow in the open-top rectangular microchannel. When liquid water drops to an entrance of the fabricated microchannel with a height of 20 μm and a width of 20 ${\mu}m$, it flows along the microchannel by only capillary force. In the wetting behavior of the liquid, important parameters of this flow are channel size, contact angle and liquid properties such as surface tension and viscosity, which are used to control the flow of liquid in the microchannel.