skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Stabilized FE simulation of prototype thermal-hydraulics problems with integrated adjoint-based capabilities

Abstract

A critical aspect of applying modern computational solution methods to complex multiphysics systems of relevance to nuclear reactor modeling, is the assessment of the predictive capability of specific proposed mathematical models. In this respect the understanding of numerical error, the sensitivity of the solution to parameters associated with input data, boundary condition uncertainty, and mathematical models is critical. Additionally, the ability to evaluate and or approximate the model efficiently, to allow development of a reasonable level of statistical diagnostics of the mathematical model and the physical system, is of central importance. In this study we report on initial efforts to apply integrated adjoint-based computational analysis and automatic differentiation tools to begin to address these issues. The study is carried out in the context of a Reynolds averaged Navier–Stokes approximation to turbulent fluid flow and heat transfer using a particular spatial discretization based on implicit fully-coupled stabilized FE methods. Initial results are presented that show the promise of these computational techniques in the context of nuclear reactor relevant prototype thermal-hydraulics problems.

Authors:
 [1];  [2];  [3];  [1];  [4];  [3]
  1. Sandia National Laboratories, Computational Mathematics Department (United States)
  2. (United States)
  3. Sandia National Laboratories, Multiphysics Applications Department (United States)
  4. Sandia National Laboratories, Optimization and UQ Department (United States)
Publication Date:
OSTI Identifier:
22572348
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 321; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ERRORS; FLUID FLOW; HEAT TRANSFER; MATHEMATICAL MODELS; MATHEMATICAL SOLUTIONS; NAVIER-STOKES EQUATIONS; REACTORS; REYNOLDS NUMBER; SENSITIVITY; SIMULATION; THERMAL HYDRAULICS

Citation Formats

Shadid, J.N., E-mail: jnshadi@sandia.gov, Department of Mathematics and Statistics, University of New Mexico, Smith, T.M., Cyr, E.C., Wildey, T.M., and Pawlowski, R.P. Stabilized FE simulation of prototype thermal-hydraulics problems with integrated adjoint-based capabilities. United States: N. p., 2016. Web. doi:10.1016/J.JCP.2016.04.062.
Shadid, J.N., E-mail: jnshadi@sandia.gov, Department of Mathematics and Statistics, University of New Mexico, Smith, T.M., Cyr, E.C., Wildey, T.M., & Pawlowski, R.P. Stabilized FE simulation of prototype thermal-hydraulics problems with integrated adjoint-based capabilities. United States. doi:10.1016/J.JCP.2016.04.062.
Shadid, J.N., E-mail: jnshadi@sandia.gov, Department of Mathematics and Statistics, University of New Mexico, Smith, T.M., Cyr, E.C., Wildey, T.M., and Pawlowski, R.P. 2016. "Stabilized FE simulation of prototype thermal-hydraulics problems with integrated adjoint-based capabilities". United States. doi:10.1016/J.JCP.2016.04.062.
@article{osti_22572348,
title = {Stabilized FE simulation of prototype thermal-hydraulics problems with integrated adjoint-based capabilities},
author = {Shadid, J.N., E-mail: jnshadi@sandia.gov and Department of Mathematics and Statistics, University of New Mexico and Smith, T.M. and Cyr, E.C. and Wildey, T.M. and Pawlowski, R.P.},
abstractNote = {A critical aspect of applying modern computational solution methods to complex multiphysics systems of relevance to nuclear reactor modeling, is the assessment of the predictive capability of specific proposed mathematical models. In this respect the understanding of numerical error, the sensitivity of the solution to parameters associated with input data, boundary condition uncertainty, and mathematical models is critical. Additionally, the ability to evaluate and or approximate the model efficiently, to allow development of a reasonable level of statistical diagnostics of the mathematical model and the physical system, is of central importance. In this study we report on initial efforts to apply integrated adjoint-based computational analysis and automatic differentiation tools to begin to address these issues. The study is carried out in the context of a Reynolds averaged Navier–Stokes approximation to turbulent fluid flow and heat transfer using a particular spatial discretization based on implicit fully-coupled stabilized FE methods. Initial results are presented that show the promise of these computational techniques in the context of nuclear reactor relevant prototype thermal-hydraulics problems.},
doi = {10.1016/J.JCP.2016.04.062},
journal = {Journal of Computational Physics},
number = ,
volume = 321,
place = {United States},
year = 2016,
month = 9
}
  • Cited by 1
  • A critical aspect of applying modern computational solution methods to complex multiphysics systems of relevance to nuclear reactor modeling, is the assessment of the predictive capability of specific proposed mathematical models. The understanding of numerical error, the sensitivity of the solution to parameters associated with input data, boundary condition uncertainty, and mathematical models is critical. Additionally, the ability to evaluate and or approximate the model efficiently, to allow development of a reasonable level of statistical diagnostics of the mathematical model and the physical system, is of central importance. In our study we report on initial efforts to apply integrated adjoint-basedmore » computational analysis and automatic differentiation tools to begin to address these issues. The study is carried out in the context of a Reynolds averaged Navier–Stokes approximation to turbulent fluid flow and heat transfer using a particular spatial discretization based on implicit fully-coupled stabilized FE methods. We present the initial results that show the promise of these computational techniques in the context of nuclear reactor relevant prototype thermal-hydraulics problems.« less
  • The long-pulse operation of the Wendelstein 7-X (W7-X) stellarator experiment is scheduled to begin in 2020. This operational phase will be equipped with water-cooled plasma facing components to allow for longer pulse durations. Certain simulated plasma scenarios have been shown to produce heat fluxes that surpass the technological limits on the edges of the divertor target elements during steady-state operation. In order to reduce the heat load on the target elements, the addition of a “scraper element” (SE) is under investigation. The SE is composed of 24 water-cooled carbon fiber reinforced carbon composite monoblock units. Multiple full-scale prototypes have beenmore » tested in the GLADIS high heat flux test facility. Previous computational studies revealed discrepancies between the simulations and experimental measurements. In this work, single-phase thermal-hydraulics modeling was performed in ANSYS CFX to identify potential causes for such discrepancies. Possible explanations investigated were the effects of a non-uniform thermal contact resistance and a potential misalignment of the monoblock fibers. And while the difference between the experimental and computational results was not resolved by a non-uniform thermal contact resistance, the computational results provided insight into the potential performance of a W7-X monoblock unit. Circumferential temperature distributions highlighted the expected boiling regions of such a unit. Finally, simulations revealed that modest angles of fiber misalignment in the monoblocks result in asymmetries at the unit edges and provide temperature differences similar to the experimental results.« less
  • Cited by 6
  • We present a moment-based acceleration algorithm applied to Monte Carlo simulation of thermal radiative-transfer problems. Our acceleration algorithm employs a continuum system of moments to accelerate convergence of stiff absorption–emission physics. The combination of energy-conserving tallies and the use of an asymptotic approximation in optically thick regions remedy the difficulties of local energy conservation and mitigation of statistical noise in such regions. We demonstrate the efficiency and accuracy of the developed method. We also compare directly to the standard linearization-based method of Fleck and Cummings [1]. A factor of 40 reduction in total computational time is achieved with the newmore » algorithm for an equivalent (or more accurate) solution as compared with the Fleck–Cummings algorithm.« less