• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.243-248
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

• Nuclear Engineering and Technology
• /
• v.37 no.6
• /
• pp.511-524
• /
• 2005

Thermalhydraulic reactor simulation of tomorrow will require a new generation of codes combining at least three scales, the CFD scale in open medium, the component scale and the system scale. DNS will be used as a support for modelling more macroscopic models. NEPTUNE is such a new generation multi-scale platform developed jointly by CEA-DEN and EDF-R&D and also supported by IRSN and FRAMATOME-ANP. The major steps towards the next generation lie in new physical models and improved numerical methods. This paper presents the advances obtained so far in physical modelling for each scale. Macroscopic models of system and component scales include multi-field modelling, transport of interfacial area, and turbulence modelling. Two-phase CFD or CMFD was first applied to boiling bubbly flow for departure from nucleate boiling investigations and to stratified flow for pressurised thermal shock investigations. The main challenges of the project are presented, some selected results are shown for each scale, and the perspectives for future are also drawn. Direct Numerical Simulation tools with Interface Tracking Techniques are also developed for even smaller scale investigations leading to a better understanding of basic physical processes and allowing the development of closure relations for macroscopic and CFD models.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.44-53
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed for a 3D component of a nuclear system code and a component-scale analysis tool. A two-fluid three-field model is used for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. An unstructured grid is adopted for realistic simulations of the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been applied to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the modified numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing very well.

• Journal of the Korean Society of Combustion
• /
• v.8 no.4
• /
• pp.15-23
• /
• 2003

The zone conditional two-fluid equations are derived and validated against DNS database of a premixed turbulent flame. The conditional statistics of major flow variables are investigated to understand the mechanism of flame generated turbulence. The flow field in burned zone shows substantially increased turbulent kinetic energy, which is highly anisotropic due to reaction kinematics across thin f1amelets. The transverse component may be larger than the axial component for a distributed pdf of the flamelet orientation angle, while the opposite occurs due to redistribution of turbulent kinetic energy and flamelet orientation normal to the flow at the end of a flame brush. The major source or sink terms of turbulent kinetic energy are the interfacial transfer by the mean reaction rate and the work terms by fluctuating pressure and velocity on a flame surface. Ad hoc modeling of some interfacial terms may be required for further application of the two-fluid model in turbulent combustion simulations.

• Nuclear Engineering and Technology
• /
• v.56 no.5
• /
• pp.1712-1724
• /
• 2024

Spiral tubes are used in a wide range of applications and it is significant to understand the vibration introduced by two-phase flow in spiral tubes. In this paper, the numerical method is used to study the vibration induced by the gas-liquid two-phase flow in a spiral tube with different flow regimes. The pressure fluctuation characteristics at the pipe wall and the solid vibration response characteristics are obtained. The results show that the motion of small bubbles in bubbly flow leads to small pressure fluctuations with low-frequency broadband (0-50 Hz). The motion of the gas plug in the plug flow causes small amplitude periodic pressure fluctuation with a shortened low-frequency broadband (0-15 Hz) compared to the bubbly flow. The motion of the gas slug in the slug flow causes large periodic fluctuations in pressure with a significant dominant frequency (6-7 Hz). The wavy flow is very stable and has a distinct main frequency (1-2 Hz). The vibration regime in the bubbly flow and wave flow are close to the first-order mode, and the vertical vibrating component is dominant. The plug flow and slug flow excite higher-order vibration modes, and the lateral vibration component plays more important part in the vibration response.

• Journal of the Korean Society for Precision Engineering
• /
• v.30 no.3
• /
• pp.317-323
• /
• 2013

In-mold coating is a reactive fluid designed to improve the surface quality of injection molded thermoplastic substrate in functional and cosmetic properties. In this study, a mixing head for in-mold coating was designed, and mixing characteristics of two-component polyurethane flowing through runner were investigated based on flow simulations. In order to achieve uniform mixing of two components injected through straight mixing head, an impingement aftermixer was used in runner design. Semi-circular cross-section was better than circular one for runners for uniform mixing. With increasing runner length and flow rate, mixing became more uniform. In addition, the degree of mixing was more improved with decreasing viscosity of isocyanate.

• Nuclear Engineering and Technology
• /
• v.41 no.10
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.

• Nuclear Engineering and Technology
• /
• v.42 no.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.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.4
• /
• pp.369-378
• /
• 2010

The creation of the hydraulic turbine flow factor map will undoubtedly benefit its design by decreasing both the design cycle time and product cost. In this paper, the geometry and flow variables, which effectively affect the flow factor, are proposed, analyzed and determined. These flow variables are further used to create the operating condition maps by using different model approaches categorized into Response Surface Method (RSM) and Artificial Neural Network (ANN). The accuracies of models created by different approaches are compared and the performances of model approaches are analyzed. The influences of chosen variables and the combination of Principle Component Analysis (PCA) and model approaches are also studied. The comparison results between predicted and actual flow factors suggest that two-hidden-layer Feed-forward Neural Network (FFNN), and one.hidden-layer FFNN with PCA has the best performance on forming this mapping, and are accurate sufficiently for hydraulic turbine design.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.315-323
• /
• 2009

Pressure oscillations triggered by the unstable interaction of dynamic flow features of the hydraulic turbine with the hydraulic plant system - including the electrical design - can at times reach significant levels and could lead to damage of plant components or could reduce component lifetime significantly. Such a problem can arise for overload as well as for part load operation of the turbine. This paper discusses an approach to analyze the overload high pressure oscillation problem using computational fluid dynamic (CFD) modeling of the hydraulic machine combined with a network modeling technique of the hydraulic system. The key factor in this analysis is the determination of the overload vortex rope volume occurring within the turbine under the runner which is acting as an active element in the system. Two different modeling techniques to compute the flow field downstream of the runner will be presented in this paper. As a first approach, single phase flow simulations are used to evaluate the vortex rope volume before moving to more sophisticated modeling which incorporates two phase flow calculations employing cavitation modeling. The influence of these different modeling strategies on the simulated plant behavior will be discussed.