• Proceedings of the Korea Institute of Fire Science and Engineering Conference
• /
• 2012.04a
• /
• pp.270-273
• /
• 2012

오픈 소스인 FDS(fire dynamic simulator)는 건물, 터널내의 화재나 연기, 열기류의 거동을 연구하기 위하여 국내외적으로 광범위하게 이용되고 있다. 체계적인 연구와 확장이 가능하도록 소스코드와 프로그램구조, 각종 메뉴얼을 갖추고 있으며 향후 개발 방향을 온라인을 통해서 소스코드와 함께 공개하고 있다. 비압축성 비정상해석을 근간으로 하고 있으며 난류유동을 해석할 수 있도록 DNS와 LES모델을 가지고 있다. 화재, 연소, 스프링클러, 화재 확산 등의 모델링을 제공하고 있다. 이러한 모델을 바탕으로 다양한 시나리오의 재난, 피난에 적용할 수 있다. 향후 이러한 기본 모델을 바탕으로 새로운 재난 시나리오에 따라 새로운 알고리즘의 적용하기 위해서는 FDS 기본적인 구조와 모델, 그리고 한계점을 이해할 필요가 있다. 본 연구에서는 이러한 FDS모델을 더욱 확장하기 위한 일환으로 FDS(V5.5.3)의 기본적인 구조을 파악하고 몇 가지 검증모델(verification)에 적용하였다. 또한 이를 향후 FDS의 소스코드를 확장할 수 있는 근간으로 삼고자 한다.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.20 no.2
• /
• pp.121-127
• /
• 2008

This paper investigates the effects of outdoor pressure fluctuations on natural ventilation through an opening on a building envelope. The ventilation airflow rate depends on the magnitude and the period of the pressure fluctuations, the size of the opening relative to the space volume, and the resistance characteristics of the opening. Non-dimensional parameters have been derived, which determine indoor pressure responses due to outdoor pressure fluctuations. The flow regions are categorized into (1) synchronized region, (2) opening resistance region, and (3) transition region depending on the non-dimensional parameter derived. Pressure fluctuations and flow characteristics are investigated numerically using the 4th order Runge-Kutta method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.7
• /
• pp.820-830
• /
• 1999

The coupled turbulent flow and solidification is considered in a typical slab continuous easting process using commercial program FIDAP. Standard $k-{\varepsilon}$ turbulence model is modified to decay turbulent viscosity in the mushy zone and laminar viscosity is set to a sufficiently large value at the solid region. This coupled turbulent flow and solidification model also contains thermal contact resistance due to the mold powder and air gap between the strand and mold using an effective thermal conductivity. From the computed flow pattern, the trajectory of inclusion particles was calculated. The comparison between the predicted and experimental solidified shell thickness shows a good agreement.

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

파랑의 전파와 변형에 대한 연구에는 수심방향으로 적분한 2차원방정식인 완경사방정식과 Boussinesq 방정식을 기반으로 한 수치모형을 이용한 연구가 최근까지 가장 활발하게 진행되어 오고 있다. 그러나 실제 구조물의 설계에는 2차원 수치모형에서 고려할 수 없는 수심방향 유속에 기인한 정확도의 문제로 인해 구조물의 형상과 재원을 설계하기 위한 정교한 수치모형실험이 어려워 주로 수리모형실험에 의존해 왔다. 수리모형실험은 실제 현상을 가장 잘 재현해낼 수 있어 신뢰성이 매우 높지만 다양한 실험을 수행하기가 어렵고 많은 시간과 비용이 소요되는 단점이 있다. 이에 따라 최근 수심방향으로 완전한 운동방정식인 Navier-Stokes 방정식을 푸는 3차원 수치모형에 대한 연구가 활발히 진행되고 있다. 이론적으로 매우 우수한 모형이긴 하나 정확도 높은 결과를 얻기 위해서는 매우 조밀한 격자를 필요로 하기 때문에 아직까지 막대한 계산시간이 필요하다는 단점이 있으나 컴퓨터 기술이 급격한 속도로 발전하고 있어 Navier-Stokes 방정식 모형의 적용 가능성은 계속 높아지고 있다. 파랑변형을 다루는 수치모형실험을 수행할 때 외부조파를 사용할 경우 구조물이나 지형에 의해 반사되어 나온 파랑이 조파지점에 도달할 때 실험영역으로 재 반사되는 문제가 발생한다. 이를 해결하기 위해 내부조파기법의 개발에 대한 연구가 필수적이었으며, 자유수면변위를 변수로 사용하는 모형의 경우 그 연구가 매우 활발하게 진행되어 왔다. 한편, Navier-Stokes 방정식 모형의 경우 자유수면변위를 변수로 사용하는 2차원 모형에 비해 상대적으로 연구가 미흡하였다. 본 연구에서는 기존의 연직 2차원 Navier-Stokes 방정식 모형에 사용된 연속방정식에 질량 원천항을 추가하는 내부조파기법을 도입하여 3차원 수치모형에서 고립파를 내부조파하고, 급경사에서의 고립파의 처오름 및 처내림 현상을 수리모형 실험결과와 비교 및 분석하였다. 수치모형은 Navier-Stokes 방정식을 엇갈림 격자체계에서 계산하는 동수압 모형으로서, Two-step projection 기법을 사용하는 유한차분모형을 사용하였다. 본 수치모형은 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 계산하는 LES(large-eddy-simulation) 기반의 수치모형으로, 난류 모델링을 위해 Smagorinsky LES 모형을 사용한다. 또한, 압력장의 계산을 위해 Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하였으며, 자유수면 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 수치모형실험이 전체적으로 수리모형실험에서 관측한 파랑의 처오름 및 처내림 현상을 잘 재현하고 있는 것으로 나타났으며, 정량적인 비교를 통해 수치모형의 성능을 검증하였다.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.4
• /
• pp.241-248
• /
• 2010

This paper deals with the performance analysis and design of the Darrieus-type vertical axis turbine to evaluate the effect of key design parameters such as number of blade, blade chord, pitch and camber. The commercial CFD software FLUENT was employed as an unsteady Reynolds-Averaged Navier-Stokes (RANS) solver with k-e turbulent model. Grid system was modelled by GAMBIT. Basic numerical methodology of the present study is appeared in Jung et al. (2009). Two-dimensional analysis was mostly adopted to avoid the barrier of massive calculation required for parametric study. It was found that the highly efficient turbine model could be designed through the optimization of design parametrrs.

• Journal of the Korean Institute of Gas
• /
• v.7 no.4 s.21
• /
• pp.24-29
• /
• 2003

Numerical analysis was conducted to predict the damage of indoor gas explosion for the propagation of explosion flame. Indoor gas diffusion distribution due to gas leakage was obtained by diffusion equation that adopted initial conditions from reference. Enthalpy of each gas-mixture ratio and reduced mechanism was applied to calculate flame temperature, and laminar combustion velocities with the variant of each gas concentration from reference were applied to the gas mixture. Turbulent combustion velocity was modeled by coupling of turbulent energy and laminar combustion velocity in k-$\epsilon$ model. For the analysis of flame propagation cartesian and cylindrical coordinate were used to indoor position and flame propagation respectively. The study analyzes the cause of pressure rising with the variation of flame propagation by glass damage, and the result shows that indoor pressure rising with ignition position varies window dimension.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.24 no.2
• /
• pp.187-207
• /
• 2019

In order to develop a high performance ocean model, we used Julia, a Just-In-Time compile language, and to obtain the solution of the momentum equation, we made the code to solve the Poisson equation by the Successive Over-Relaxation method. And then we made two models to test Julia calculation codes. First, a simple channel form is modeled to test constant source/sink conditions. Second, the simplified Yellow Sea was modeled to test tidal forcing, Coriolis forces, and the effect of vertical eddy diffusivity coefficients. The model has been tested with a total of eight cases in the two scenarios. As a result of the test, the depth-averaged current speed of the three cases in Scenario 1 converged perfectly to the theoretical value, and that showed well a vertical flow velocity gradient due to the bottom friction. Also, the result of Scenario 2 represented well the amphidromic points of Yellow Sea and the tidal characteristics of mid-western and southwestern coast of Korea. Therefore, it is considered that the ocean model using Julia language has developed successfully, this suggests that the ocean model has come to the stage of successful transition from a classical compile language to a Just-In-Time compile language.

• Journal of Wetlands Research
• /
• v.13 no.1
• /
• pp.79-94
• /
• 2011

As demonstrated in study for non-submerged groynes, the flow field is predominantly two-dimensional, with mainly horizontal eddies. The eddies shed form the tips of the groynes and migrate in the flow direction. These eddies have horizontal dimensions in the order of tens of meters and time-scales in the order of minutes. In the standard flow simulations, these motions are usually not resolved, due to a too coarse grid, too large time steps and, more importantly, the use of inadequate turbulence modelling. using for example a k-${\varepsilon}$ model, it is necessary to introduce substantial modifications. Therefore simulation resolved in this study, were carried out using the DELFT-3D-MOR programme, which is part of the DELFT3D software package of WL/Delft Hydraulics and In this study, apply a two-dimensional depth-averaged model, taking an horizontal large eddy simulation(HLES). The bed morphology computed when using HLES, as well as the associated time-scale, is similar to what has been obseved in a field case. When using a mean-flow model with-out HELS, the bed morphology is less realistic and the morphological time-scale is much larger. This slow development is the result of neglecting(or averaging). the strong velocity fluctuations associated with the time-varying eddy formation.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.1
• /
• pp.64-72
• /
• 1997

An analytical study has been performed to investigate the unsteady ignition characteristics of pulse combustion. In many combustion applications, strain rate of the flow can significantly affect the combustion features; ignition, extinction, and reignition. In the pulse combustion, two jets (hot combustion gases and fresh mixtures) coming from the opposite side of the combustor will collide in the combustor forming a stagnation region where the chemical reaction is suppressed by the strain rate until this becomes below the critical value. In this research, the method of large activation energy asymptotic is adopted with one step irreversible kinetics to examine the ignition response to the periodic variation of the strain rate of flow. The results show the variation of the maximum value of strain rate can determine whether the ignition or extinction occur.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.8
• /
• pp.99-106
• /
• 2003

Geometrical design variables are numerically investigated to improve aerodynamic performance of the supersonic impulse turbine of a turbopump in a liquid rocket engine. Aerodynamic redesign was performed for maximization of the blade power. Four design variables considered are blade angle, blade thickness and radii of upper and lower arc blade with appropriate constraints. A fast Navier-Stokes solver was developed and Chien's k-$\varepsilon$ turbulence modelling was used for turbulence closure. In initial shape, a flow separation was found in the middle of blade chord. However, it disappeared in final shape via its geometrical design variable change. About 3.2 percent of blade power was increased from this research.