• 에너지공학
• /
• 제11권3호
• /
• pp.237-246
• /
• 2002

한국형차세대원자로 APR-1400의 안전감압계통이 작동하면 물, 공기 및 증기가 sparger를 통해 격납건물 내 핵연료재장전 수조로 차례로 방출된다. 방출 과정 중 생기는 여러 현상 중에서 수조 내의 공기 기포군은 저주파, 고진폭의 진동 하중을 발생하며, 주파수가 침수 구조물의 고유 주파수와 거의 같은 경우에는 구조물에 심각한 영향을 줄 수 있다. 이러한 현상은 복잡하기 때문에 주파수와 하중에 대한 규명은 주로 실험에 의존해 왔으며 수치해석적 연구는 이루어지지 않았다. 본 연구에서는 sparger를 통해 수조 내로 방출되는 공기 기포군의 거동에 대한 수치해석을 상용 열수력 해석 코드인 FLUENT Version 4.5를 사용하여 수행하였다. 다상유동 해석모델중 VOF(Volume Of Fluid)모델을 사용하여 물, 공기 및 증기 등의 다상유동을 모의하였다. 해석결과를 sparger 개발을 위해 ABB-Atom이 수행하였던 실험결과와 비교하여 만족할만한 결과를 얻었다.

• 에너지공학
• /
• 제12권4호
• /
• pp.259-266
• /
• 2003

안전감압계통 작동시 수조에서 발생하는 공기 기포군의 진동 하중을 줄이기 위해 ABB-Atom sparger에 는 하중저감 링이 설치되어 있다. 하중저감 링이 압력장에 미치는 영향을 알아보기 위해, 본 연구에서는 하중저감 링이 없는 ABB-Atom sparger를 통해 수조 내로 방출되는 공기 기포군의 진동에 대한 수치해석을 상용 열수력 해석 코드인 FLUENT 4.5를 사용하여 수행하였다. 코드에 내재된 다상유동 모델 중 VOF(Volume Of Fluid)모델을 사용하여 물, 공기 및 증기 유동을 모의하였다. 해석결과를 기존의 해석결과와 비교하여 하중저감 링은 벽면 압력 하중을 줄이는 효과가 있음을 확인하였다. 아울러 배관내의 공기량과 배관 입구 조건이 벽면 압력 진동에 미치는 영향도 분석하였다.

• 대한기계학회논문집B
• /
• 제35권7호
• /
• pp.723-733
• /
• 2011

본 연구에서는 비정렬격자계와 체적포착법을 사용하여 표면장력이 지배적인 다상유동의 수치해석 방법을 제시하였다. 먼저 표면장력에 대한 CSF(Continuum Surface Force) 모델을 비정렬격자계에 적용할 수 있도록 수치해석 방법을 확립시켜 Myong(2009)이 개발한 비정렬격자계와 체적포착법을 사용한 수치 해석코드에 삽입하였다. 테스트 문제로 오직 표면장력만이 존재하는 평형상태의 정적(static) 액적 및 비평형상태의 동적(dynamic) 액적 문제에 적용하여, 이 해석방법의 유용성과 정확도를 평가하였다. 연구결과, 매끄러운 곡률 계산을 위해 필요한 필터로 본 연구에서 제안한 Laplacian 필터와 함께 CSF 모델로는 밀도보정(density-scaled)한 CSF 모델이 예측성능이 우수한 것으로 나타났다. 또한 표면장력 계산을 위한 이 모델을 채용한 본 수치해석방법은 표면장력이 지배적인 다상유동인 평형상태의 정적 액적 및 비평형상태의 동적 액적 문제 모두에 대해 정확성과 유용성이 입증되었다.

• 한국해양공학회지
• /
• 제26권6호
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

• 대한조선학회논문집
• /
• 제51권1호
• /
• pp.58-66
• /
• 2014

Numerical simulations for the flows containing free surface remain difficult problems because the drastic differences of physical properties of water and air, The difference of densities makes the solution instable in particular. For the stabilities of the solutions, the most typical methods to simulate free surface flows, such as Volume Of Fluid(VOF) and Level-Set(LS) methods, impose transient zones where the physical prosperities are continuously distributed. The thickness of the transient zone is the source of the numerical errors. The other side, marker-density method does not use such a transient zone. In the traditional marker-density method, however, the air velocities of free surface cells are extrapolated from the water velocity, and the pressures on the free surface are extrapolated from the air pressures for the stability of the solution. In this study, the marker-density method is modified for the decrease of such numerical errors. That is, the pressure on the free surface is determined to coincide with the pressure gradient terms of the governing equations, and the velocity of free surface cells are calculated with the governing equations. Two-dimensional steady spilling breakers behind of a submersed hydrofoil and three-dimensional spilling breaker near a wedge shaped ship model are simulated using INHAWAVE-II including the modified marker-density(MMD) method. The results are compared with the results of Fluent V6.3 including VOF method and several published research results.

• 대한토목학회논문집
• /
• 제28권6B호
• /
• pp.627-634
• /
• 2008

개수로의 난류흐름을 VOF(Volume of Fluid)기법을 채용한 RANS(Reynolds averaged Navier-Stokes) 방정식 모형을 사용하여 수치모의할 때 벽면함수의 거칠기를 산정하기 위해 고려해야 하는 점들을 연구하였다. 거친 벽면상의 흐름을 위한 벽면함수의 거칠기 상수(roughness constant)는 관수로 흐름의 실험을 통하여 얻어진 값을 사용한다. 그러나 개수로 흐름에서는 이 거칠기 상수가 Froude 수에 따라 변화하므로 이를 고려할 수 있어야 하며, 개수로에서 광범위하게 사용되는 Manning 조도계수에 상응하는 벽면 거칠기 높이(roughness height)를 산정하여 사용할 수 있어야 있다. 본 연구에서는 모형에 입력되는 벽면함수의 거칠기 높이와 Manning 조도계수 사이의 관계를 분석하였다. 이를 바탕으로 수치모형의 특성이 고려되고 Manning 조도계수의 함수로 표현되는 유효 거칠기 높이 산정식을 제안하였다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제8권6호
• /
• pp.589-601
• /
• 2016

The hybrid grid was adopted and numerical prediction analysis of propeller unsteady bearing force considering free surface was performed for mode and full-scale KCS hull-propeller-rudder system by employing RANS method and VOF model. In order to obtain the propeller velocity under self-propulsion point, firstly, the numerical simulation for self-propulsion test of full-scale ship is carried out. The results show that the scale effect of velocity at self-propulsion point and wake fraction is obvious. Then, the transient two-phase flow calculations are performed for model and full-scale KCS hull-propeller-rudder systems. According to the monitoring data, it is found that the propeller unsteady bearing force is fluctuating periodically over time and full-scale propeller's time-average value is smaller than model-scale's. The frequency spectrum curves are also provided after fast Fourier transform. By analyzing the frequency spectrum data, it is easy to summarize that each component of the propeller bearing force have the same fluctuation frequency and the peak in BFP is maximum. What's more, each component of full-scale bearing force's fluctuation value is bigger than model-scale's except the bending moment coefficient about the Y-axis.

• 한국전산유체공학회지
• /
• 제15권1호
• /
• pp.81-87
• /
• 2010

In this study, we conducted resistance and propulsion performance test of ship composed of the Resistance Test, Propeller Open Water Test and Self Propulsion Test using the CFD(Computational Fluid Dynamics). We used commercial RANS(Reynolds Averaged Navier Stokes equation) solver, as a calculating tool. The unstructured grids were used in a bow and stern of ship, having complex shape, for a convenience of generating grids, and the structured grids were adopted in a central hull and rest of hull having a relatively simple shape which is called hybrid grid method. In addition, The sliding mesh method was adopted to rotate a propeller directly in the Propeller Open Water and Self Propulsion Test. The Resistance Test and Self Propulsion Test were calculated using Volume of Fluid (VOF) model and considering a free surface. And all The three cases were applied realizable k-epsilon model as the turbulence model. The results of calculations were verified for the suitability of calculations by comparing MOERI's EFD results.

• 대한기계학회논문집B
• /
• 제30권8호
• /
• pp.780-787
• /
• 2006

Several types of different header designs are numerically studied to have uniform distribution of two phase flow in the evaporator header having multi-channels. The different geometries include the inlet tube position into the header and the width of header. In the numerical calculation, two types of two-phase model such as homogeneous model and VOF(Volume Of Fluid) model are employed. In this study, the mal-distribution number, $M_d$, is newly defined to evaluate the averaged level of the flow distribution in the whole passes of the evaporator. As results, two phase flow in the header can be visualized using post-processing of numerical results. Furthermore, the optimum position of the inlet tube into the header and the width of header can be proposed for the better distribution of refrigerant(R-134a) flow.

• 해양환경안전학회지
• /
• 제20권3호
• /
• pp.312-317
• /
• 2014

In this paper, the wave run-up height and depression depth around air-water interface-piercing circular cylinder have been numerically studied. The Reynolds Averaged Navier-Stokes equations (RANS) and continuity equations are solved with Reynolds Stress model (RSM) and volume of fluid (VOF) method as turbulence model and free surface modeling, respectively. A commercial Computational Fluid Dynamics (CFD) software "Star-CCM+" has been used for the current simulations. Various Froude numbers ranged from 0.2 to 1.6 are used to investigate the change of air-water interface structures around the cylinder and experimental data and theoretical values by Bernoulli are compared. The present results showed a good agreement with other studies. Kelvin waves behind the cylinder were generated and its wave lengths are longer as Froude numbers increase and they have good agreement with theoretical values. And its angles are smaller with the increase of Froude numbers.