• 한국안전학회지
• /
• 제24권1호
• /
• pp.14-20
• /
• 2009

Direct numerical simulations(DNS) were performed for the prediction of transient buoyant jet diffusion flames where the Froude numbers(Fr) are 5 and 160, respectively. The thermodynamic and transport properties were evaluated using CHEMKIN package to enhance the prediction performance of the DNS code. The simulated buoyant jet diffusion flame of Fr=5 and 160 showed the transient, dynamic motion well. It was identified that the buoyant jet flames were flickered periodically, and the simulated flickering frequency of the jet diffusion flame of Fr=5 was 12.5Hz, which was in good agreement with the experimental results. The flamelet structures of the buoyant jet diffusion flames could be well understood by comparing the scalar dissipation rates(SDR) and the heat release rates(HRR) of the flames. It was found that the SDR was strongly coupled with the HRR in the buoyant jet diffusion flames.

• Ocean Systems Engineering
• /
• 제3권4호
• /
• pp.257-273
• /
• 2013

Top tensioned riser (TTR) is often used in a floating oil/gas production system deployed in deep water for oil/gas transport. This study focuses on the extension of the existing numerical code, known as CABLE3D, to allow for static and dynamic simulation of a TTR connected to a floating structure through a tensioner system or buoyancy can, and restrained by riser guides at different elevations. A tensioner system usually consists of three to six cylindrical tensioners. Although the stiffness of individual tensioner is assumed to be linear, the resultant stiffness of a tensioner system may be nonlinear. The vertical friction between a TTR and the hull at its riser guide is neglected assuming rollers are installed there. Near the water surface, a TTR is forced to move horizontally due to the motion of the upper deck of a floating structure as well as related riser guides. The extended CABLE3D is then integrated into a numerical code, known as COUPLE, for the simulation of the dynamic interaction among the hull of a floating structure, such as spar or TLP, its mooring system and riser system under the impact of wind, current and waves. To demonstrate the application of the extended CABLE3D and its integration with COUPLE, the numerical simulation is made for a truss spar under the impact of Hurricane "Ike". The mooring system of the spar consists of nine mooring lines and the riser system consists of six TTRs and two steel catenary risers (SCRs).

• 대한조선학회논문집
• /
• 제52권5호
• /
• pp.387-394
• /
• 2015

Before salvaging a wrecked ship, the physics-based simulation is needed to predict lifting force before real operation by floating crane or barge. Procedures affecting lifting force for the salvage can be divided into three stages. At the first stage, the bottom breakout force for the wrecked ship to escape from seabed sediment should be calculated. At the second step, the current force acting on the wrecked ship while lifting from the seabed to near sea surface should be considered. Finally, buoyancy change near at the sea surface when the wrecked ship start to escape from the water should be considered. In the previous studies, only the breakout force at the first stage was calculated based on simple assumption of embedment depth and contact area of the wrecked ship. Therefore, we develop a program for salvage simulation including whole stages. It is composed of four modules such as the equations of motion, time integration, force calculation, and visualization. As a result, it is applied to simulate lifting the wrecked ship according to various environmental loads including seabed sediments.

• 환경영향평가
• /
• 제32권6호
• /
• pp.473-487
• /
• 2023

1970~1980년대 준공된 우리나라의 다목적댐 저수지에는 발전방류를 위한 고정식 취수구가 심층에 설치되어 있어 일부 댐 하류에서는 냉수 방류에 의한 농작물 냉해, 안개 일수 증가 등의 문제점이 제기된 바 있다. 본 연구에서는 고정식 취수구를 통해 심층 취수가 이루어지고 있는 합천호를 대상으로 취수 수심이 저수지의 수온 성층 구조와 방류 수온에 미치는 영향을 분석하였다. 3차원 수리수질모형인 AEM3D를 이용하여 합천호의 연직 수온 분포를 재현하고 계절별 수온성층 구조를 분석하였으며, 수문 조건에 따른 수온성층 변화를 비교 분석하기 위하여 풍수해와 갈수해를 대상으로 모델링 하였다. 또한 취수심 변경 시나리오를 적용함으로써 취수 수심이 수온성층 구조에 미치는 영향을 분석하였다. 모의 결과 심층 취수를 표층 취수로 변경할 경우 수온약층의 형성 위치가 풍수해 6.5 m, 갈수해 6.8 m 감소하여 더 얕은 수심에 형성될 것으로 분석되었다. 또한 수체 안정도 지수인 Schmidt Stability Index (SSI)와 Buoyancy frequency (N²)가 증가하여 수온성층 강도가 증가하는 것으로 나타났다. 표층 취수시 연평균 방류수온이 풍수해 3.5℃, 갈수해 5.0℃ 증가하여 하류하천의 영향은 감소하나, 호내의 저수온층 수체적과 수온성층 강도가 증가하므로 추후 수질관리를 위해 취수심을 댐 운영의 주요인자로 고려해야 할 것으로 판단된다.

• International Journal of Air-Conditioning and Refrigeration
• /
• 제10권3호
• /
• pp.117-128
• /
• 2002

This paper describes the smoke filling process of a fire field model based on a self-deve-loped SMEP (Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy production term. Also it solves the radiation equation using the discrete ordinates method. Compressibility is assumed and the perfect gas law is used. Comparison of the calculated upper-layer average tempera-ture and smoke layer clear height with the zone models has shown reasonable agreement. The zone models used are the CFAST and the NBTC one-room. For atrium fires with ceiling glass the ceiling heat flux by solar heat causes a high smoke temperature near the ceiling. However, it has no effect on the smoke movement such as the smoke layer clear heights that are important in fire safety. In conclusion, the smoke layer clear heights that are important in evacuation activity except the early of a fire were not as sensitive as the smoke layer tem-perature to the nature of ceiling heat flux condition. Thus, a fire sensor in atrium with ceiling glass has to consider these phenomena.

• 한국연소학회지
• /
• 제13권4호
• /
• pp.26-36
• /
• 2008

Experimental study is conducted to identify the existence of a shrinking flame disk and to clarify its flame characteristics through the inspection of critical mole fraction at flame extinction and edge flame oscillation at low strain rate flames. Experiments are made as varying global strain rate, velocity ratio, and burner distance. The transition from a shrinking flame disk to a flame hole is verified through gradient measurements of maximum flame temperature. The evidence of edge flame oscillation in flame disk is also provided through numerical simulation in microgravity. It is found at low strain rate flame disks in normal gravity that buoyancy effects are importantly contributing to lateral heat loss to burner rim, and is proven through critical mole fraction at flame extinction, edge flame oscillation, and measurements of flame temperature gradient along flame disk surface.

• 한국환경과학회지
• /
• 제23권7호
• /
• pp.1223-1232
• /
• 2014

Density currents have been easily observed in environmental flows, for instance turbidity currents and pollutant plumes in the oceans and rivers. In this study, we explored the propagation dynamics of density currents using the FLOW-3D computational fluid dynamics code. The renormalization group (RNG) $k-{\varepsilon}$ scheme, a turbulence numerical technique, is employed in a Reynold-averaged Navier-Stokes framework (RANS). The numerical simulations focused on two different types of intrusive density flows: (1) propagating into a two-layer ambient fluid; (2) propagating into a linearly stratified fluid. In the study of intrusive density flows into a two-layer ambient fluid, intrusive speeds were compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting propagation speed of the density currents. We also numerically reproduced the effect of the ratio of current depth to the overall depth of fluid. The numerical model provided excellent agreement with the analytical values. It was also clearly demonstrated that RNG $k-{\varepsilon}$ scheme within RANS framework is able to accurately simulate the dynamics of density currents. Simulations intruding into a continuously stratified fluid with the various buoyancy frequencies are carried out. These simulations demonstrate that three different propagation patterns can be developed according to the value of $h_n/H$ : (1) underflows developed with $h_n/H=0$ ; (2) overflows developed when $h_n/H=1$ ; (3) intrusive interflow occurred with the condition of 0 < $h_n/H$ < 1.

• 대한기계학회논문집
• /
• 제8권6호
• /
• pp.536-543
• /
• 1984

본 연구에서는 유동 조건을 표시하는 무차원 수로 널리 이용되고 있는 Richa- rdson수(R$_{i}$)를 변수로 하여 R$_{i}$의 변화에 따른 온도 분포, 표면 유동 속도와 확산율 등을 연구하였다. 아울러 이러한 열 구배에 의한 성층화 유동(stratifield flow)에서 R$_{i}$가 유동 특성과 어떤 상관 관계를 가지는 가를 알아보기 위하여 같 은 R$_{i}$값을 갖는 다른 조건들의 실험을 수행하였다. 비슷한 R$_{i}$값에서의 온 도 분포를 비교 분석함으로써 아직까지 정립되고 있지 않는 열오염 현상 해석에서의 지배 변수로서의 R$_{i}$의 위치를 확인하고, 부력 효과에 의한 흐름의 난류 구조가 어떻게 영향 받는가 하는 문제에도 관심을 두었다.제에도 관심을 두었다.

• 한국태양에너지학회 논문집
• /
• 제29권5호
• /
• pp.28-34
• /
• 2009

The atrium in the building has many advantages including its aesthetic and functional effect. But the upper part of the atrium has the thermal problem of overheating due to insolation through the window. But natural ventilation aided buoyancy effect can be a solution to make comfortable indoor environment. Proper design of openings is very important to improve thermal environment in the atrium. In this study, thermal evaluations were performed to improve thermal environment in the atrium. Indoor thermal environment of an atrium at Seoul was measured in the field and simulated with Computational Fluid Dynamics( CFD) code. The turbulent flow model adopted is $K-{\varepsilon}$ model. The results of computer simulations are compared with the measurements at the point in the atrium. In order to improve the indoor ventilation environment of the atrium, thermal environment evaluations of six alternatives were conducted. After evaluations of the results, the design guidelines of an atrium are suggested.

• Ocean Systems Engineering
• /
• 제6권4호
• /
• pp.345-362
• /
• 2016

This study aims to develop the numerical method to estimate level ice impact load and investigate the dynamic interaction between an arctic Spar with sloped surface and drifting level ice. When the level ice approaches the downward sloped structure, the interaction can be decomposed into three sequential phases: the breaking phase, when ice contacts the structure and is bent by bending moment; the rotating phase, when the broken ice is submerged and rotated underneath the structure; and the sliding phase, when the submerged broken ice becomes parallel to the sloping surface causing buoyancy-induced fictional forces. In each phase, the analytical formulas are constructed to account for the relevant physics and the results are compared to other existing methods or standards. The time-dependent ice load is coupled with hull-riser-mooring coupled dynamic analysis program. Then, the fully coupled program is applied to a moored arctic Spar with sloped surface with drifting level ice. The occurrence of dynamic resonance between ice load and spar motion causing large mooring tension is demonstrated.