Stabilized FE simulation of prototype thermalhydraulics problems with integrated adjointbased 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 adjointbased 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 fullycoupled stabilized FE methods. Initial results are presented that show the promise of these computational techniques in the context of nuclear reactor relevant prototype thermalhydraulics problems.
 Authors:
 Sandia National Laboratories, Computational Mathematics Department (United States)
 (United States)
 Sandia National Laboratories, Multiphysics Applications Department (United States)
 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; NAVIERSTOKES EQUATIONS; REACTORS; REYNOLDS NUMBER; SENSITIVITY; SIMULATION; THERMAL HYDRAULICS
Citation Formats
Shadid, J.N., Email: 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 thermalhydraulics problems with integrated adjointbased capabilities. United States: N. p., 2016.
Web. doi:10.1016/J.JCP.2016.04.062.
Shadid, J.N., Email: 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 thermalhydraulics problems with integrated adjointbased capabilities. United States. doi:10.1016/J.JCP.2016.04.062.
Shadid, J.N., Email: 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 thermalhydraulics problems with integrated adjointbased capabilities". United States.
doi:10.1016/J.JCP.2016.04.062.
@article{osti_22572348,
title = {Stabilized FE simulation of prototype thermalhydraulics problems with integrated adjointbased capabilities},
author = {Shadid, J.N., Email: 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 adjointbased 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 fullycoupled stabilized FE methods. Initial results are presented that show the promise of these computational techniques in the context of nuclear reactor relevant prototype thermalhydraulics problems.},
doi = {10.1016/J.JCP.2016.04.062},
journal = {Journal of Computational Physics},
number = ,
volume = 321,
place = {United States},
year = 2016,
month = 9
}

Stabilized FE simulation of prototype thermalhydraulics problems with integrated adjointbased capabilities
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 adjointbasedmore »Cited by 1 
Thermalhydraulics modeling for prototype testing of the W7X high heat flux scraper element
The longpulse operation of the Wendelstein 7X (W7X) stellarator experiment is scheduled to begin in 2020. This operational phase will be equipped with watercooled 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 steadystate 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 watercooled carbon fiber reinforced carbon composite monoblock units. Multiple fullscale prototypes have beenmore » 
Monte Carlo simulation methods in momentbased scalebridging algorithms for thermal radiativetransfer problems
We present a momentbased acceleration algorithm applied to Monte Carlo simulation of thermal radiativetransfer problems. Our acceleration algorithm employs a continuum system of moments to accelerate convergence of stiff absorption–emission physics. The combination of energyconserving 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 linearizationbased method of Fleck and Cummings [1]. A factor of 40 reduction in total computational time is achieved with the newmore »