• Energy Research

Abstract A novel multi-scale domain overlapping coupling methodology designed to couple a computational fluid dynamics (CFD) code with a system thermal hydraulic (STH) code was developed and its performance was investigated. The methodology has been implemented in the coupling infrastructure code Janus, developed at the University of Michigan, providing methods for the on-the-fly data transfer through memory between the commercial CFD code STAR-CCM+ and the US NRC best-estimate thermal hydraulic system code TRACE. Coupling between these two software packages is motivated by the desire to extend the range of applicability of TRACE to scenarios in which local momentum and energy transfer are important, such as three-dimensional mixing of localized slugs of deborated or cold water in the downcomer and lower plenum of a reactor pressure vessel. The intra-fluid shear forces necessary to correctly capture these effects are neglected in the TRACE equations of motion, but are readily calculated from CFD solutions. CFD/STH coupling implementations therefore have applications in reactor transients such as boron dilution scenarios, Anticipated Transient Without SCRAM (ATWS) and Main Steam Line Break (MSLB). The proposed method is based on aliasing all spatial sources and sinks of momentum in the CFD domain as frictional losses in the system code domain. The internal velocity fields and, consequently, the inertial component of the pressure field are maintained consistent between the CFD and STH domains through a complementary velocity-matching interface. In this paper, coupled simulations are performed on Cartesian and cylindrical geometry with emphasis on consistency, convergence, and stability during transient scenarios. Results show that the presented domain overlapping coupling method is capable of adjusting pressure and velocity profiles of multi-dimensional system code solutions to match CFD solutions accurately. Important characteristics of transient simulations were found to include the background flow rate, specifically the stabilizing effect of viscous forces, as well as the time derivative of the flow rate. Under certain adverse conditions, the basic coupling method is found to produce unstable behavior. A stabilization method for adjusting CFD data is laid out and found to significantly improve the method’s performance under the most challenging conditions. Recommendations are laid out for further improving the coupling via advanced time stepping methods.

A novel multiscale coupling methodology based on a domain overlapping approach has been developed to couple a computational fluid dynamics code with a best-estimate thermal hydraulic code. The methodology has been implemented in the coupling infrastructure code Janus, developed at the University of Michigan, providing methods for the online data transfer between the commercial computational fluid dynamics code STAR-CCM+ and the US NRC best-estimate thermal hydraulic system code TRACE. Coupling between these two software packages is motivated by the desire to extend the range of applicability of TRACE to scenarios in which local momentum and energy transfer are important, such as three-dimensional mixing. These types of flows are relevant, for example, in the simulation of passive safety systems including large containment pools, or for flow mixing in the reactor pressure vessel downcomer of current light water reactors and integral small modular reactors. The intrafluid shear forces neglected by TRACE equations of motion are readily calculated from computational fluid dynamics solutions. Consequently, the coupling methods used in this study are built around correcting TRACE solutions with data from a corresponding STAR-CCM+ solution. Two coupling strategies are discussed in the paper: one based on a novel domain overlapping approach specifically designed for transient operation, and a second based on the well-known domain decomposition approach. In the present paper, we discuss the application of the two coupling methods to the simulation of open and closed loops in both steady state and transient operation. The objective of this study is to examine the performance of each coupling method in terms of convergence, consistency, and numerical stability. As expected, the results produced by the two methods were found to be identical once numerical convergence is achieved and consistent with the standalone STAR-CCM+ solution in both steady state and transient cases. However, the domain overlapping method was found to achieve convergence at larger integration time steps than the domain decomposition approach and exhibited superior convergence and numerical stability characteristics in both steady state and transient scenarios.

Abstract Hydrocarbon and nuclear fuels will contribute substantially to the growing global energy portfolio for some time to come. At the same time atmospheric carbon levels and economic concerns are increasingly shaping the power generation environment. Supercritical fluid technologies are being investigated to address both concerns. For example, direct fired supercritical CO 2 power cycles offer key advantages for carbon capture and thermal efficiency. Additionally, the supercritical water reactor Generation IV nuclear reactor allows higher power densities and more efficient operation than current pressurized or boiling water reactors. As these technologies are developed, a need exists for accurate, yet cost-effective modeling and simulation solutions to assist with design and safety calculations. Supercritical fluid properties are strongly dependent on temperature which results in a number of unique flow features. Among these features are the complex buoyancy and buoyancy-turbulence interaction behaviors exhibited by supercritical flows that are not found in liquid or gas flows. Therefore, standard turbulence models with current model coefficient values are not likely to be applicable to supercritical fluid flows. The subject of this paper is efficient modeling and simulation of supercritical flows to support nascent technologies growing to maturation and operational deployment. We present a framework for developing Reynolds-Averaged Navier-Stokes turbulence models specifically equipped for the challenges of supercritical flows. A novel formulation of the algebraic heat flux model of the buoyancy production of turbulence term is used with a traditional shear stress transport model. To produce a new turbulence model for supercritical channel flows, the empirical coefficients of the resulting four equation model were calculated from data previously published by other authors regarding upward supercritical flow through heated pipes. The presented SST k t - ω t approach was validated against heated tube experiments to show predictive capabilities in moderate flow conditions, where buoyancy effects are important but not dominating of inertial effects.