• Nuclear Engineering and Technology
• /
• 제42권6호
• /
• pp.636-655
• /
• 2010

A thermal-hydraulic code, named CUPID, has been being developed for the realistic analysis of transient two-phase flows in nuclear reactor components. The CUPID code development was motivated from very practical needs, including the analyses of a downcomer boiling, a two-phase flow mixing in a pool, and a two-phase flow in a direct vessel injection system. The CUPID code adopts a two-fluid, three-field model for two-phase flows, and the governing equations are solved over unstructured grids with a semi-implicit two-step method. This paper presents an overview of the CUPID code development and assessment strategy. It also presents the code couplings with a system code, MARS, and, a three-dimensional reactor kinetics code, MASTER.

• 한국전산유체공학회지
• /
• 제9권3호
• /
• pp.73-80
• /
• 2004

A LS(Level Set) formulation is developed for computing two-phase flows on non- orthogonal meshes. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. The LS formulation is implemented into a general purpose program for 3-D flows and verified through several test problems.

• 한국전산유체공학회지
• /
• 제12권2호
• /
• pp.21-25
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Nuclear Engineering and Technology
• /
• 제48권1호
• /
• pp.64-71
• /
• 2016

Void fraction is an important parameter in the oil industry. This quantity is necessary for volume rate measurement in multiphase flows. In this study, the void fraction percentage was estimated precisely, independent of the flow regime in gas-liquid two-phase flows by using ${\gamma}-ray$ attenuation and a multilayer perceptron neural network. In all previous studies that implemented a multibeam ${\gamma}-ray$ attenuation technique to determine void fraction independent of the flow regime in two-phase flows, three or more detectors were used while in this study just two NaI detectors were used. Using fewer detectors is of advantage in industrial nuclear gauges because of reduced expense and improved simplicity. In this work, an artificial neural network is also implemented to predict the void fraction percentage independent of the flow regime. To do this, a multilayer perceptron neural network is used for developing the artificial neural network model in MATLAB. The required data for training and testing the network in three different regimes (annular, stratified, and bubbly) were obtained using an experimental setup. Using the technique developed in this work, void fraction percentages were predicted with mean relative error of <1.4%.

• 한국유체기계학회 논문집
• /
• 제11권3호
• /
• pp.7-13
• /
• 2008

The cavitating flow simulation is of practical importance for many engineering systems, such as pump, turbine, nozzle, Infector, etc. In the present work, a solver for two-phase flows has been developed and applied to simulate the cavitating flows past hydrofoils. The governing equation is the two-phase Navier-Stokes equation, comprised of the continuity equation of liquid and vapor phase. The momentum and energy equation is in the mixture phase. The solver employs an implicit, dual time, preconditioned algorithm using finite difference scheme in curvilinear coordinates. An experimental data and other numerical data were compared with the present results to validate the present solver. It is concluded that the present numerical code has successfully accounted for two-phase Navier-Stokes model of cavitation flow.

• 한국항공우주학회지
• /
• 제37권1호
• /
• pp.1-16
• /
• 2009

최근 공학적으로 응용 범위가 넓은 액체-기체의 2상 유동장에 대한 전산 해석이 주목을 받고 있다. 본 연구에서는 밀도차가 큰 압축성 2상 유동장을 해석하기에 적합하도록 기존에 개발된 RoeM과 AUSMPW+ 공간 이산화 기법을 확장하였다. RoeM과 AUSMPW+에서 사용되는 충격파 포착항을 2상 혼합류의 상태방정식으로부터 새롭게 정의하여, 밀도와 음속이 상이한 두 유체에서도 정확성을 보장하면서 충격파 불안정성을 제거할 수 있도록 하였다. 개발된 2상 유동 RoeM과 AUSMPW+를 몇 가지 예제 문제에 적용하여 검증하였으며, 해석 결과는 두 수치기법이 일반적인 2상 유동에서도 충격파 안정적이며 정확한 특성을 갖는 것을 보여주었다.

• Nuclear Engineering and Technology
• /
• 제47권4호
• /
• pp.416-423
• /
• 2015

It had been disputed how to apply wall drag to the dispersed phase in the framework of the conventional two-fluid model for two-phase flows. Recently, Kim et al. [1] introduced the volume-averaged momentum equation based on the equation of a solid/fluid particle motion. They showed theoretically that for dispersed two-phase flows, the overall two-phase pressure drop by wall friction must be apportioned to each phase, in proportion to each phase fraction. In this study, the validity of the proposed wall drag model is demonstrated though one-dimensional (1D) simulations. In addition, it is shown that the existing form loss model incorrectly predicts the motion of the dispersed phase. A new form loss model is proposed to overcome that problem. The newly proposed form loss model is tested in the region covering the lower plenum and the core in a nuclear power plant. As a result, it is shown that the new models can correctly predict the relative velocity of the dispersed phase to the surrounding fluid velocity in the core with spacer grids.

• Industrial Engineering and Management Systems
• /
• 제7권2호
• /
• pp.84-92
• /
• 2008

Scheduling is one of the most important issues in the planning and operation of production systems. This paper investigates a general job shop scheduling problem with reentrant work flows and sequence dependent setup times. The disjunctive graph representation is used to capture the interactions between machines in job shop. Based on this representation, four two-phase heuristic procedures are proposed to obtain near optimal solutions for this problem. The obtained solutions in the first phase are substantially improved by reversing the direction of some critical disjunctive arcs of the graph in the second phase. A comparative study is conducted to examine the performance of these proposed algorithms.

• 한국전산유체공학회지
• /
• 제13권4호
• /
• pp.86-95
• /
• 2008

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

• 대한수학회지
• /
• 제51권1호
• /
• pp.163-187
• /
• 2014

We present a multi-stage Roe-type numerical scheme for a model of two-phase flows arisen from the modeling of deflagration-to-detonation transition in granular materials. The first stage in the construction of the scheme computes the volume fraction at every time step. The second stage deals with the nonconservative terms in the governing equations which produces states on both side of the contact wave at each node. In the third stage, a Roe matrix for the two-phase is used to apply on the states obtained from the second stage. This scheme is shown to capture stationary waves and preserves the positivity of the volume fractions. Finally, we present numerical tests which all indicate that the proposed scheme can give very good approximations to the exact solution.