• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.389-392
• /
• 2002

In many industrial applications, multiphase flow analysis is the norm rather than an exception as compared to more-conventional single-phase investigation. This paper describes the implementation of the multiphase flow simulation capability in the general purpose CFD software AVL FIRE/SWIFT. The governing equations are discretized based on a finite volume method (FVM) suitable fur very complex geometry, The pressure field is obtained using the SIMPLE algorithm. Depending on the characteristics of the multiphase flow to be examined, the user can choose either the two-fluid model or an explicit interface-tracking model based on the Volume-of-Fluid approach. For truly 'multi'-phase flow problems, it is also possible to apply a hybrid model where certain phases are explicitly tracked while the other phases are handled by the two fluid model. In order to demonstrate the capability of the method, applications to the Taylor bubble flow simulations are presented.

• Nuclear Engineering and Technology
• /
• v.52 no.4
• /
• pp.708-720
• /
• 2020

Core-wide temperature distribution in sodium-cooled fast reactor plays a key role in its decay heat removal process, however the prediction for temperature distribution is quite complex due to the conjugate heat transfer between the assembly flow and the inter-wrapper flow. Hybrid medium model has been proposed for conjugate heat transfer modeling in the core. The core is modeled with a Realistic modeled inter-wrapper flow and hybrid medium modeled assembly flow. To validate present model, simulations for a three-assembly model were performed with Realistic modeling, traditional porous medium model and hybrid medium model, respectively. The influences of Uniform/Non-Uniform power distribution among assemblies and the Peclet number within the assembly flow have been considered. Compared to traditional porous medium model, present model shows a better agreement with in Realistic modeling prediction of the temperature distribution and the radial heat transfer between the inter-wrapper flow and the assembly flow.

• Journal of Ocean Engineering and Technology
• /
• v.35 no.1
• /
• pp.50-58
• /
• 2021

The demand for eco-friendly energy is expected to increase due to the recently strengthened environmental regulations. In particular, the flow inside the pipe used in a cargo handling system (CHS) or fuel gas supply system (FGSS) of hydrogen transport ships and hydrogen-powered ships exhibits a very complex pattern of multiphase-thermal flow, including the boiling phenomenon and high accuracy analysis is required concerning safety. In this study, a feasibility study applying the boiling model was conducted to analyze the multiphase-thermal flow in the pipe considering the phase change. Two types of boiling models were employed and compared to implement the subcooled boiling phenomenon in nucleate boiling numerically. One was the "Rohsenow boiling model", which is the most commonly used one among the VOF (Volume-of-Fluid) boiling models under the Eulerian-Eulerian framework. The other was the "wall boiling model", which is suitable for nucleate boiling among the Eulerian multiphase models. Moreover, a comparative study was conducted by combining the nucleate site density and bubble departure diameter model that could influence the accuracy of the wall boiling model. A comparison of the Rohsenow boiling and the wall boiling models showed that the wall boiling model relatively well represented the process of bubble formation and development, even though more computation time was consumed. Among the combination of models used in the wall boiling model, the simulation results were affected significantly by the bubble departure diameter model, which had a very close relationship with the grid size. The present results are expected to provide useful information for identifying the characteristics of various parameters of the boiling model used in CFD simulations of multiphase-thermalflow, including phase change and selecting the appropriate parameters.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.1
• /
• pp.27-41
• /
• 2014

In this paper, we present an efficient numerical method for multiphase image segmentation using a multiphase-field model. The method combines the vector-valued Allen-Cahn phase-field equation with initial data fitting terms containing prescribed interface width and fidelity constants. An efficient numerical solution is achieved using the recently developed hybrid operator splitting method for the vector-valued Allen-Cahn phase-field equation. We split the modified vector-valued Allen-Cahn equation into a nonlinear equation and a linear diffusion equation with a source term. The linear diffusion equation is discretized using an implicit scheme and the resulting implicit discrete system of equations is solved by a multigrid method. The nonlinear equation is solved semi-analytically using a closed-form solution. And by treating the source term of the linear diffusion equation explicitly, we solve the modified vector-valued Allen-Cahn equation in a decoupled way. By decoupling the governing equation, we can speed up the segmentation process with multiple phases. We perform some characteristic numerical experiments for multiphase image segmentation.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2728-2735
• /
• 2022

The research on the flow and heat transfer characteristics of lead bismuth(LBE) is significant for the thermal-hydraulic calculation, safety analysis and practical application of lead-based fast reactors(LFR). In this paper, a new CFD model is proposed to solve the thermal-hydraulic analysis of LBE. The model includes two parts: turbulent model and turbulent Prandtl, which are the important factors for LBE. In order to find the best model, the experiment data and design of 19-pin hexagonal rod bundle with spacer grid, undertaken at the Karlsruhe Liquid Metal Laboratory (KALLA) are used for CFD calculation. Furthermore, the turbulent model includes SST k - 𝜔 and k - 𝜀; the turbulent Prandtl includes Cheng-Tak and constant (Pr_t =1.5,2.0,2.5,3.0). Among them, the combination between SST k - 𝜔 and Cheng-Tak is more suitable for the experiment. But in the low Pe region, the deviation between the experiment data and CFD result is too much. The reason may be the inlet-effect and when Pe is in a low level, the number of molecular thermal diffusion occupies an absolute advantage, and the buoyancy will enhance. In order to test and verify versatility of the model, the NCCL performed by the Nuclear Thermal-hydraulic Laboratory (Nuthel) of Xi'an Jiao tong University is used for CFD to calculate. This paper provides two verification examples for the new universal model.

• Journal of Korea Society of Waste Management
• /
• v.34 no.7
• /
• pp.685-696
• /
• 2017

This study focuses on computational particle fluid dynamics (CPFD) modeling for the fast pyrolysis of biomass in a conical spouted bed reactor. The CPFD simulation was conducted to understand the hydrodynamics, heat transfer, and biomass fast pyrolysis reaction of the conical spouted bed reactor and the multiphase-particle in cell (MP-PIC) model was used to investigate the fast pyrolysis of biomass in a conical spouted bed reactor. A two-stage semi-global kinetics model was applied to model the fast pyrolysis reaction of biomass and the commercial code (Barracuda) was used in simulations. The temperature of solid particles in a conical spouted bed reactor showed a uniform temperature distribution along the reactor height. The yield of fast pyrolysis products from the simulation was compared with the experimental data; the yield of fast pyrolysis products was 74.1wt.% tar, 17.4wt.% gas, and 8.5wt.% char. The comparison of experimental measurements and model predictions shows the model's accuracy. The CPFD simulation results had great potential to aid the future design and optimization of the fast pyrolysis process for biomass.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.209-219
• /
• 2013

The present paper deals with the extension of AUSMPW+ scheme into two-fluid model for multiphase flow. AUSMPW+ scheme is the improvement of a single-phase AUSM+ scheme by designing pressure-based weighting functions to prevent oscillations near a wall and shock instability after a strong shock. Recently, Kitamura and Liou assessed a family of AUSM-type schemes with two-fluid model governing equations [K. Kitamura and M.-S. Liou, Comparative study of AUSM-Family schemes in compressible multi-phase flow simulations, ICCFD7-3702 (2012)]. It was observed that the direct application of the single-phase AUSMPW+ did not provide satisfactory results for most of numerical test cases, which motivates the current study. It turns out that, by designing pressure-based weighting functions, which play a key role in controlling numerical diffusion for two-fluid model, problems reported in can be overcome. Various numerical experiments validate the proposed modification of AUSMPW+ scheme is accurate and robust to solve multiphase flow within the framework of two-fluid model.

• Nuclear Engineering and Technology
• /
• v.52 no.10
• /
• pp.2196-2203
• /
• 2020

A dryout mechanism model for rectangular narrow channels at high pressure conditions is developed by assuming that the Kelvin-Helmholtz instability triggered the occurrence of dryout. This model combines the advantages of theoretical analysis and empirical correlation. The unknown coefficients in the theoretical derivation are supported by the experimental data. Meanwhile, the decisive restriction of the experimental conditions on the applicability of the empirical correlation is avoided. The expression of vapor phase velocity at the time of dryout is derived, and the empirical correlation of liquid film thickness is introduced. Since the CHF value obtained from the liquid film thickness should be the same as the value obtained from the Kelvin-Helmholtz critical stability under the same condition, the convergent CHF value is obtained by iteratively calculating. Comparing with the experimental data under the pressure of 6.89-13.79 MPa, the average error of the model is -15.4% with the 95% confidence interval [-20.5%, -10.4%]. And the pressure has a decisive influence on the prediction accuracy of this model. Compared with the existing dryout code, the calculation speed of this model is faster, and the calculation accuracy is improved. This model, with great portability, could be applied to different objects and working conditions by changing the expression of the vapor phase velocity when the dryout phenomenon is triggered and the calculation formula of the liquid film.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1995.05a
• /
• pp.238-244
• /
• 1995

• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.446-454
• /
• 2021

Critical heat flux (CHF) has traditionally been evaluated using look-up tables or empirical correlations for nuclear power plants. However, under complex moving condition, it is necessary to reconsider the CHF characteristics since the conventional CHF prediction methods would no longer be applicable. In this paper, the additional forces caused by motions have been added to the annular film dryout (AFD) mechanistic model to investigate the effect of moving condition on CHF. Moreover, a theoretical model of the natural circulation loop with additional forces is established to reflect the natural circulation characteristics of the loop system. By coupling the system loop with the AFD mechanistic model, a CHF prediction program called NACOM for natural circulation loop under moving condition is developed. The effects of three operating conditions, namely stationary, inclination and rolling, on the CHF of the loop are then analyzed. It can be clearly seen that the moving condition has an adverse effect on the CHF in the natural circulation system. For the calculation parameters in this paper, the CHF can be reduced by 25% compared with the static value, which indicates that it is important to consider the effects of moving condition to retain adequate safety margin in subsequent thermal-hydraulic designs.