• Journal of Advanced Marine Engineering and Technology
• /
• v.18 no.3
• /
• pp.94-101
• /
• 1994

To develope the new simulator for the analysis of fluid flow information, the influence of various convective difference schemes were evaluated. General curvilinear coordinate system with nonorthogonal grids was adopted for the successful analysis of various complex geometries. Computation results show that if we can not obtain full grid numbers within available computational environment, we need to use higher order finite difference schemes to keep the prediction accuracy.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.64-64
• /
• 2016

하천은 인간에게 용수의 이용 및 하천호안의 휴식처로써의 이용을 통해 직접적인 영향을 주고, 하천구조물의 심미적 영향, 랜드마크로써의 역할을 통해 간접적인 영향을 준다. 또한, 하천은 하천생태계에 서식하는 동 식물에게 영향을 준다. 그러나 하천유사로 인해 통수능이 감소하고, 하천구조물 주변에 침식을 야기할 뿐만 아니라, 댐과 저수지에 유사의 퇴적으로 저수용량의 감소시킨다. 그러므로 이를 예측하는 것은 경제적, 환경적으로 중요하다. 하상변동의 모의를 위해 기존의 2차원 모형은 만곡흐름에서 유동의 helical flow를 고려하지 않아 예측이 부정확하였다. 본 연구에서는 천수방정식을 이용한 하상변동 수치모의에 helical flow의 영향을 고려하였다. 하천과 같은 천수영역에서의 흐름 및 하상변동을 해석하기 위하여 수심평균 된 Navier-Stokes equations인 천수방정식을 이용하였다. 지배방정식은 곡선 좌표계에서 유한체적법으로 차분하였고, 비엇갈림격자를 사용하였다. 지배방정식의 닫힘 문제를 해결하기위해 0-방정식 난류모형을 사용하였고, "time marching" 기법의 적용을 위해 계산단계분할 방법을 이용하였다. 비엇갈림격자의 사용으로 인해 검사체적의 면에서의 유속이 필요하여 pressure-velocity coupling을 사용하여 유속의 진동을 줄였다. 또한, 만곡부의 helical flow를 모의하기위해 helical flow intensity model을 도입하였다. 앞에서 계산한 흐름을 바탕으로 유사량 산정공식과 Exner 방정식을 이용하여 하상변동을 모의하였다. 흐름의 검증, helical flow의 영향에 대한 확인, 하상변동의 적용을 위해 선행연구의 실험이 사용되었다.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.2
• /
• pp.22-31
• /
• 1997

We developed a general purpose program for the analysis of flows in the combustor with recirculating flow regime and simulated the flows. The program uses non-staggered grids based on finite volume method and the primitive variables are cartesian velocities. The combustion model is irreversible one step reaction with infinite chemistry The Favre averaged governing equations are considered and the clipped gaussian distribution is considered as a probability density function of the conserved scalar. We calculated turbulent diffusion flame with recirculating flow regime. Simulation shows two recirculating regions like experimental results. Velocity, turbulent kinetic energy, temperature and concentration distribution in simulation agree well with experimental data.

• Proceedings of the Korean Institute of Industrial Safety Conference
• /
• 1997.11a
• /
• pp.43-48
• /
• 1997

The smoke filling process for the atrium space containing a fire source is simulated using two types of deterministic fire modus: Zone model and Field model. The zone mode used is the CFAST(version 1.6) mode developed at the Building and Fire Research laboratories, NIST in the USA. The lied model is a self-developed fire field model based on Computational Fluid Dynamics(CFD) theories. This article is focused on finding out the smoke movement and temperature distribution in atrium space which is cubic in shape. A computational procedure for predicting velocity and temperature distribution in fro-induced flow is based on the solution, in finite volume method and non-staggered grid system, of 3-dimensional equations for the conservation of mass, momentum, energy, species and so forth. The fire model i. e. Zone model and Field model predicted similar results for Ire clear height and the smoke layer temperature.

• Fire Science and Engineering
• /
• v.10 no.3
• /
• pp.10-18
• /
• 1996

In foreign country such as USA and Japan, considerable research has been done regarding the spread of smoke in room of fire involvement by using computer. but, in our country it has not been, so, this paper presents a detailed qualitative description of phenomena which occure during typical fire scenarios through numberical analysis. The governing equations are solved by using FVM method with non-staggered grid. The SIMPLE method for pressure-velocity couple and power-law scheme for convention terms are used. It shows that a plume is formed, hot plume gases impinge on the ceiling and they spread across it. then, it eventually reaches the bounding walls of the enclosure. It takes 20s for smoke to fill the enclosure.

• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.218-226
• /
• 2008

A component-scale two-phase analysis code has been developed for a realistic simulation of two-phase flow transients in a light water nuclear reactor component. In the code, a two-fluid three-field model is adopted and the governing equations are solved on an unstructured mesh. For the numerical solution scheme, the semi-implicit method used in the RELAP5 code was selected, which has been proved to be very stable and accurate for most of practical applications. However, some modifications were needed for its application to an unstructured non-staggered grid. This paper presents the modified semi-implicit numerical method for unstructured grid and the preliminary results of the calculations.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.109-114
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1821-1826
• /
• 2003

In this study, we discuss CIP method, which can treat compressible/incompressible problem and multi-phase problem. We can apply this method to the general equations by using CCUP. In this paper, non-staggered grid arrangement and predictor-corrector method are used to enhance later development and the solution accuracy and convergence performance. To validate the numerical algorithm proposed in this paper, the two-dimensional unsteady flow in the driven cavity is simulated. The numerical results of this subject using the present algorithm are compared with other numerical results. As a result, it is prived that the present scheme gives stable and improved solutions, and the results even coarse grid are in good agreement with other result using a fine grid arrangement.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.1
• /
• pp.212-218
• /
• 1995

Turbulent flow and heat transfer characteristics around staggered tube bundles were studied using a non-orthogonal boundary fitted coordinate system and the low Reynolds .kappa. - .epsilon. turbulence model suggested by Lam and Bremhorst. The predicted flow characteristics for two tube pitches and tube arrangement showed good agreement with the experimental data except the strongly curved region. The predicted Nusselt number was compared with measurements obtained in the staggered rough bundles and it revealed the similar trend to measurements, but the location of the maximum and minimum heat transfer differed somewhat from the measurements.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.