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

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.5
• /
• pp.535-549
• /
• 1998

A numerical study for performance analysis of a counterflow type forced draft tower and natural draft cooling tower has been performed based on the method using the finite volume method with non-orthogonal body fitted and non-staggered grid system. For solving the coupling problem between water and air, air enthalpy balance, moisture fraction balance, water enthalpy balance, and water mass balance equations are solved with Navier-Stoke’s equations simultaneously. For the effect of turbulence, the standard k-$\varepsilon$ turbulent model is implied in this analysis. The predicted result of the present analysis is compared with the experimental data and the commercial software result to validate the present study, The predicted results show good agreement with the experimental data and the commercial software result. To investigate the influence of the cooling tower design parameters such as approach, range and wet bulb temperature, parametric studies are also peformed.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.8 no.1
• /
• pp.1-13
• /
• 2000

A numerical study for performance analysis of a crossflow-type forced draft cooling tower has been performed based on the finite volume method with non-orthogonal body fitted, and non-staggered grid system. For solving the coupling problem between water and air, air enthalpy, moisture fraction, water enthalpy, and water mass balance equations are solved with Navier-Stoke's equations simultaneously. For the effect of turbulence, the standard k-$\varepsilon$ turbulent model is implied in this analysis. The predicted result of the present analysis is compared with the experimental data and the commercial software result to validate the present study. The predicted results show good agreement with the experimental data and the commercial software result. To investigate the influence of the cooling tower design parameters such as approach, range and wet bulb temperature, parametric studies are also performed.

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

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.6
• /
• pp.766-773
• /
• 1999

A numerical analysis has been performed solidification problem using the fixed grid-enthalpy method with enthalpy-porosity relation. A modified standard $k-\varepsilon$ model was applied to describe the influence of turbulent flow. Computational procedures are based on the finite volume method and the non-staggered grid system. Comparisons with the different three experimental results show that applying a modified standard $k-\varepsilon$model in mushyzone is better than the previous computation results. This paper includes another EMBR's influences such as change of velocity field, Increasement of temperature and dispersion of flow out of nozzle into the flow field. These EMBR's influences are compared to case without EMBR.

• Journal of Ship and Ocean Technology
• /
• v.5 no.1
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.317-328
• /
• 2002

A numerical investigation on the turbulent flows over a three-dimensional steep hill is presented. The numerical model developed for the present work is based on the finite volume method and the SIMPLE algorithm with a non-staggered grid system. Standard $k-{\varepsilon}$ model and Shih's non-linear model are tested for the validation of the prediction accuracy in the 3D separated flow. Comparisons of the mean velocity and turbulence profiles between the numerical predictions and the measurements show good agreement. The Shih's non-linear model is found to predict mean flow and turbulence better than the Standard $k-{\varepsilon}$. Flow patterns have also been examined to explain the difference in the cavity zone between 2D and 3D hills.

• Nuclear Engineering and Technology
• /
• v.32 no.2
• /
• pp.169-179
• /
• 2000

A detailed numerical analysis of the evolution of thermal stratification in a curved piping system in a nuclear power plant is performed. A finite volume based thermal-hydraulic computer code has been developed employing a body-fitted, non-orthogonal curvilinear coordinate for this purpose. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of momentum interpolation method. The SIMPLE algorithm is employed for the pressure and velocity coupling, and the convection terms are approximated by a higher-order bounded scheme. The thermal-hydraulic computer code developed in the present study has been applied to the analysis of thermal stratification in a curved duct and some of the predicted results are compared with the available experimental data. It is shown that the predicted results agree fairly well with the experimental measurements and the transient formation of thermal stratification in a curved duct is also well predicted.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.5 s.149
• /
• pp.538-549
• /
• 2006

A code is developed to simulate a viscous flow field around a deformable body using the hybrid Cartesian/immersed boundary method. In this method, the immersed boundary(IB) nodes are defined near the body boundary then velocities at the IB nodes are reconstructed based on the interpolation along the normal direction to the body surface. A new method is suggested to define the IB nodes so that a closed fluid domain is guaranteed by a set of IB nodes and the method is applicable to a zero-thickness body such as a sail. To validate the developed code, the vorticity fields are compared with other recent calculations where a cylinder orbits and moves into its own wake. It is shown the code can handle a sharp trailing edge at Reynolds number of $10^5$ under moderate requirements on girds. Finally the developed code is applied to simulate the vortex shedding behind a deforming foil with flapping tail like a fish. It is shown that the acceleration of fluids near the flapping tail contributes to the generation of the thrust for propulsion.