Numerical Simulation of Boiling Two-Phase Flow in Tight-Lattice Rod Bundle by 3-Dimensional Two-Fluid Model Code ACE-3D

Author(s):  
Hiroyuki Yoshida ◽  
Takeharu Misawa ◽  
Kazuyuki Takase

Two-fluid model can simulate two phase flow less computational cost than inter-face tracking method and particle interaction method. Therefore, two-fluid model is useful for thermal hydraulic analysis in large-scale domain such as a rod bundle. Japan Atomic Energy Agency (JAEA) develops three dimensional two-fluid model analysis code ACE-3D, which adopts boundary fitted coordinate system in order to simulate complex shape channel flow. In this paper, boiling two-phase flow analysis in a tight lattice rod bundle is performed by ACE-3D code. The parallel computation using 126CPUs is applied to this analysis. In the results, the void fraction, which distributes in outermost region of rod bundle, is lower than that in center region of rod bundle. At height z = 0.5 m, void fraction in the gap region is higher in comparison with that in center region of the subchannel. However, at height of z = 1.1m, higher void fraction distribution exists in center region of the subchannel in comparison with the gap region. The tendency of void fraction to concentrate in the gap region at vicinity of boiling starting point, and to move into subchannel as water goes through rod bundle, is qualitatively agreement with the measurement results by neutron radiography. To evaluate effects of two-phase flow model used in ACE-3D code, numerical simulation of boiling two-phase in tight lattice rod bundle with no lift force model (neglecting lift force acting on bubbles) is also performed. From the comparison of numerical results, it is concluded that the effects of lift force model are not so large on overall void fraction distribution in tight lattice rod bundle. However, higher void fraction distribution in center region of the subchannel was not observed in this simulation. It is concluded that the lift force model is important for local void fraction distribution in rod bundles.

Author(s):  
Shimo Yu ◽  
Xiao Yan ◽  
Junyi Zhang

Two-phase flow is an important and common phenomenon in nuclear reactor systems, and the characteristics of two-phase flow such as heat transfer and pressure drop strongly depend on the radial distribution of void fraction. This paper is presenting the CFD simulation for void fraction radial distribution of mono- and poly-disperse air-water two phase flow using Euler-Euler two-fluid model. Interfacial forces including transverse forces such as lift, wall and turbulent dispersion forces are taken into account, Furthermore the bubble size distribution and bubble break-up and coalescence processes are taken into account in case of a poly-disperse flow by using the S-Gamma model. The sauter mean diameter and interfacial area concentration (IAC) distribution can also be obtained. The simulation results are compared to an experimental database of MT-LOOP test facility (FZD, Germany)[1].


2016 ◽  
Vol 139 (1) ◽  
Author(s):  
Avinash Vaidheeswaran ◽  
Deoras Prabhudharwadkar ◽  
Paul Guilbert ◽  
John R. Buchanan ◽  
Martin Lopez de Bertodano

A new two-fluid model averaging in the near-wall region is proposed to ensure consistent matching of the two-phase k–ε turbulence model with the two-phase logarithmic law of the wall (Marie J. L., Moursali, E., and Tran-Cong, S., 1997, “Similarity Law and Turbulence Intensity Profiles in a Bubbly Boundary Layer,” Int. J. Multiphase Flow, 23(2), pp. 227–247). The void fraction distribution obtained with the averaging procedure is seen to conform to the two-phase wall function approach which is based on a double step function void fraction distribution. In particular, the proposed averaging technique is shown to achieve grid convergence in the near-wall region, which could not be obtained otherwise. Computational fluid dynamics (CFD) results with the proposed technique are in good agreement with experiments on upward bubbly flows over a flat plate, and upward and downward flows in pipes. An additional advantage of the proposed technique is that it replaces the wall force model, which has a significant degree of uncertainty in turbulent flow modeling, with a simpler geometric constraint.


2003 ◽  
Vol 2003 (0) ◽  
pp. 189
Author(s):  
Masaki Misawa ◽  
Akihiko Minato ◽  
Akio Suzuki ◽  
Naoki Ichikawa ◽  
Masaharu Kuroda

2021 ◽  
Vol 33 (3) ◽  
pp. 033324
Author(s):  
Alejandro Clausse ◽  
Martín López de Bertodano

Author(s):  
David Heinze ◽  
Thomas Schulenberg ◽  
Lars Behnke

A simulation model for the direct contact condensation of steam in subcooled water is presented that allows determination of major parameters of the process, such as the jet penetration length. Entrainment of water by the steam jet is modeled based on the Kelvin–Helmholtz and Rayleigh–Taylor instability theories. Primary atomization due to acceleration of interfacial waves and secondary atomization due to aerodynamic forces account for the initial size of entrained droplets. The resulting steam-water two-phase flow is simulated based on a one-dimensional two-fluid model. An interfacial area transport equation is used to track changes of the interfacial area density due to droplet entrainment and steam condensation. Interfacial heat and mass transfer rates during condensation are calculated using the two-resistance model. The resulting two-phase flow equations constitute a system of ordinary differential equations, which is solved by means of the explicit Runge–Kutta–Fehlberg algorithm. The simulation results are in good qualitative agreement with published experimental data over a wide range of pool temperatures and mass flow rates.


Author(s):  
Wei Yao ◽  
Christophe Morel

In this paper, a multidimensional two-fluid model with additional turbulence k–ε equations is used to predict the two-phase parameters distribution in freon R12 boiling flow. The 3D module of the CATHARE code is used for numerical calculation. The DEBORA experiment has been chosen to evaluate our models. The radial profiles of the outlet parameters were measured by means of an optical probe. The comparison of the radial profiles of void fraction, liquid temperature, gas velocity and volumetric interfacial area at the end of the heated section shows that the multidimensional two-fluid model with proper constitutive relations can yield reasonably predicted results in boiling conditions. Sensitivity tests show that the turbulent dispersion force, which involves the void fraction gradient, plays an important role in determining the void fraction distribution; and the turbulence eddy viscosity is a significant factor to influence the liquid temperature distribution.


Sign in / Sign up

Export Citation Format

Share Document