DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Interface capturing simulations of bubble population effects in PWR subchannels

Abstract

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fullymore » resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.« less

Authors:
 [1]; ORCiD logo [2]; ORCiD logo [1]
  1. North Carolina State Univ., Raleigh, NC (United States)
  2. Argonne National Lab. (ANL), Argonne, IL (United States). Nuclear Engineering Division
Publication Date:
Research Org.:
Argonne National Laboratory (ANL), Argonne, IL (United States)
Sponsoring Org.:
USDOE Office of Nuclear Energy (NE); USDOE Office of Science (SC)
OSTI Identifier:
1770592
Alternate Identifier(s):
OSTI ID: 1691969
Grant/Contract Number:  
AC02-06CH11357; AC05-00OR22725
Resource Type:
Accepted Manuscript
Journal Name:
Nuclear Engineering and Design
Additional Journal Information:
Journal Volume: 365; Journal ID: ISSN 0029-5493
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
22 GENERAL STUDIES OF NUCLEAR REACTORS; DNS; interface capturing; PWR subchannel

Citation Formats

Cambareri, Joseph J., Fang, Jun, and Bolotnov, Igor A. Interface capturing simulations of bubble population effects in PWR subchannels. United States: N. p., 2020. Web. doi:10.1016/j.nucengdes.2020.110709.
Cambareri, Joseph J., Fang, Jun, & Bolotnov, Igor A. Interface capturing simulations of bubble population effects in PWR subchannels. United States. https://doi.org/10.1016/j.nucengdes.2020.110709
Cambareri, Joseph J., Fang, Jun, and Bolotnov, Igor A. Tue . "Interface capturing simulations of bubble population effects in PWR subchannels". United States. https://doi.org/10.1016/j.nucengdes.2020.110709. https://www.osti.gov/servlets/purl/1770592.
@article{osti_1770592,
title = {Interface capturing simulations of bubble population effects in PWR subchannels},
author = {Cambareri, Joseph J. and Fang, Jun and Bolotnov, Igor A.},
abstractNote = {As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.},
doi = {10.1016/j.nucengdes.2020.110709},
journal = {Nuclear Engineering and Design},
number = ,
volume = 365,
place = {United States},
year = {Tue Jun 09 00:00:00 EDT 2020},
month = {Tue Jun 09 00:00:00 EDT 2020}
}

Works referenced in this record:

Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations
journal, November 1988


Scalable Implicit Flow Solver for Realistic Wing Simulations with Flow Control
journal, November 2014

  • Rasquin, Michel; Smith, Cameron; Chitale, Kedar
  • Computing in Science & Engineering, Vol. 16, Issue 6
  • DOI: 10.1109/MCSE.2014.75

An Adaptive Level Set Approach for Incompressible Two-Phase Flows
journal, January 1999

  • Sussman, Mark; Almgren, Ann S.; Bell, John B.
  • Journal of Computational Physics, Vol. 148, Issue 1
  • DOI: 10.1006/jcph.1998.6106

Prediction of turbulence suppression and turbulence modeling in bubbly two-phase flow
journal, June 1993


Development of Swirling Flow in a Rod Bundle Subchannel
journal, August 2002

  • McClusky, Heather L.; Holloway, Mary V.; Beasley, Donald E.
  • Journal of Fluids Engineering, Vol. 124, Issue 3
  • DOI: 10.1115/1.1478066

A stabilized finite element method for computing turbulence
journal, May 1999


Direct numerical simulation of turbulent channel flow up to
journal, June 2015


Mapping of the Lateral Flow Field in Typical Subchannels of a Support Grid With Vanes
journal, November 2003

  • McClusky, Heather L.; Holloway, Mary V.; Conover, Timothy A.
  • Journal of Fluids Engineering, Vol. 125, Issue 6
  • DOI: 10.1115/1.1625688

DNS of turbulent flow with hemispherical wall roughness
journal, January 2015


Turbulence suppression in bubbly two-phase flow
journal, September 1990


A continuum method for modeling surface tension
journal, June 1992


Interface tracking simulations of bubbly flows in PWR relevant geometries
journal, February 2017


Bubble tracking analysis of PWR two-phase flow simulations based on the level set method
journal, November 2017


Direct numerical simulation of reactor two-phase flows enabled by high-performance computing
journal, April 2018


A stabilized finite element method for the incompressible Navier-Stokes equations using a hierarchical basis
journal, January 2001


Interface Tracking Investigation of Geometric Effects on the Bubbly Flow in PWR Subchannels
journal, August 2018


Influence of Bubbles on the Turbulence Anisotropy
journal, April 2013

  • Bolotnov, Igor A.
  • Journal of Fluids Engineering, Vol. 135, Issue 5
  • DOI: 10.1115/1.4023651

A parallel adaptive mesh method for the numerical simulation of multiphase flows
journal, October 2013


Validation of the sub-channel code F-COBRA-TF
journal, September 2008


Computation of incompressible bubble dynamics with a stabilized finite element level set method
journal, October 2005

  • Nagrath, Sunitha; Jansen, Kenneth E.; Lahey, Richard T.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 194, Issue 42-44
  • DOI: 10.1016/j.cma.2004.11.012