• 에너지공학
• /
• 제17권4호
• /
• pp.218-226
• /
• 2008

가압 경수로의 주요 기기에서 발생할 수 있는 과도 2상 유동(Two-phase flow) 현상에 대한 해석을 수행하기 위해 원자로 기기 열수력 해석 코드를 개발 중에 있다. 개발중인 기기 열수력 해석 코드는 지배 방정식으로 Two-phase, three-field model을 사용하고 있으며, 복잡한 기하학적 형상의 원자로 기기를 모사하기 위해 비정렬 격자계(Unstructured grid)를 활용하고 있다. 수치해석 기법으로는, 원자로 계통 해석코드 REIAP5가 사용 중이며 대부분의 원자로 내 2상 유동 조건에서 안정적이며 정확하다고 알려진 Semi-implicit 방법을 적용하였다. 그러나 기존의 Semi-implicit 방법은 1차원, 엇갈림격자(Staggered grid)에 대해 개발되었기 때문에 이를 다차원, 비정렬, 비엇갈림 격자(Non-staggered grid)에 적용하기 위해 기존의 Semi-implicit 방법을 수정하였다. 본 논문에서는 수정된 Semi-implicit 방법을 소개하고 이를 이용해 수행한 예비 계산결과를 수록하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• 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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• 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.

• 대한기계학회논문집
• /
• 제19권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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 춘계학술대회논문집
• /
• pp.290-297
• /
• 2009

가압 경수로의 주요 기기에서 발생할 수 있는 과도 2상 유동(Two-phase flow) 현상에 대한 해석을 수행하기 위해 원자로 기기 열수력 해석 코드를 개발 중에 있다. 개발 중인 기기 열수력 해석 코드는 지배 방정식으로 Two-phase, three-field model을 사용하고 있으며, 복잡한 기하학적 형상의 원자로 기기를 모사하기 위해 비정렬 격자계(Unstructured grid)를 활용하고 있다. 수치해석 기법으로는, 원자로 계통 해석코드 RELAP5가 사용 중이며 대부분의 원자로 내 2상 유동 조건에서 안정적이며 정확하다고 알려진 Semi-implicit 방법을 적용하였다. 그러나 기존의 Semi-implicit 방법은 1차원, 엇갈림격자(Staggered grid)에 대해 개발되었기 때문에, 이를 다차원, 비정렬, 비엇갈림 격자(Non-staggered grid)에 적용하기 위해 기존의 Semi-implicit 방법을 수정하였다. 본 논문에서는 Semi-implicit 방법의 대류항을 이차정확도를 갖도록 확장하였으며, 이차정확도에 의한 수치확산의 감소를 평가하기 위해 수행된 수치시험의 결과를 기술하였다. 이차정확도 및 일차정확도로 계산된 값을 해석해 또는 격자 수렴성 시험을 통해 평가해 본 결과, 이차정확도 계산시 수치 확산의 감소 확인하였다.

• 한국전산유체공학회지
• /
• 제15권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.

• 한국항공우주학회지
• /
• 제30권5호
• /
• pp.101-108
• /
• 2002

본 연구는 비압축성 열린 공동 유동에 대한 수치적 모사이다. 2차원 Navier-Stokes 방적식을 제어체적에 대해 엇갈림 격자계를 이용하여 공간에 대해서는 C-QUICK을 시간에 대해서는 내재적 기법을 이용하여 적분하였다. 압력장은 SIMPLE-C 알고리즘에 의하여 계산 되었다. 정상 모드에서는 나타나지 않지만 전단층 모드의 경우에 나타나는 세 번째 소용돌이가 공동 유동의 안정성에 중요한 역할을 하는 것을 알 수 있었다. 이를 바탕으로 공동 앞전 아래 벽면에 Jet blowing을 정상 상태와 비정상 상태로 가하여 그 영향을 알아보았다. 주기적인 blowing 인 경우 가진 주기와 위상차 그리고 속도 크기가 중요한 변수이며 이 변수들에 의한 영향을 연구하였다.

• 설비공학논문집
• /
• 제8권3호
• /
• pp.375-385
• /
• 1996

Turbulent flow and heat transfer characteristics around staggered tube banks were studied using the 3-D Navier-Stokes equations and energy equation governing a steady incompressible flow, which were reformulated in a non-orthogonal coordinate system with cartesian velocity components and discretized by the finite volume method with a non-staggered variable arrangement. The predicted turbulent kinetic energy using RNG $k-{\varepsilon}$ model was lower than that of standard $k-{\varepsilon}$ model but showed same result for mean flow field quantities. The prediction of the skin friction coefficient using RNG $k-{\varepsilon}$ model showed better trend with experimental data than standard $k-{\varepsilon}$ model result. The inclined flow showed higher velocity and skin friction coefficient than transverse flow because of extra strain rate ($\frac{{\partial}w}{{\partial}y}$). Also, this was why the inclined flow showed higher local heat transfer coefficient than the transverse flow.