Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiphysics simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Continuous Integration, In-Code Documentation, and Automation for Nuclear Quality Assurance Conformance

The Multiphysics Object Oriented Simulation Environment (MOOSE) is an open-source, finite element framework for solving highly coupled sets of nonlinear equations. The development of the framework and applications occurs concurrently using an agile, continuous-integration software package. Included in the framework is an in-code, extensible documentation system. Using these two tools in union with the repository management tools GitHub and GitLab, a software quality plan was created and followed such that MOOSE and a MOOSE-based application (BISON) have been shown to meet the American Society of Mechanical Engineers’ Nuclear Quality Assurance-1 standard. The approach relies heavily on automation for both testing and documentation. The resulting effort demonstrates that a rigorous software quality plan may be implemented that incurs a minimal impact on day-to-day development of the software, satisfying the stringent guidelines necessary to operate the software in a safety function within a nuclear facility.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A MOOSE-Based Model for Fission Product Transport and Source Term Estimation for High-Temperature Gas-Cooled Reactors

Thanks to fuel elements containing tristructural isotropic (TRISO) particles combined with a low core power density and passive feedback mechanisms leading to modest temperature rises in the event of accidental events, high-temperature gas-cooled reactors (HTGRs) offer a high degree of reliability in terms of fission product retention. While the anticipated source term for HTGRs is expected to be very low, it is important to provide a quantitative estimate of radiological releases during nominal and accidental conditions. Here, we propose a computationally efficient mechanistic source term methodology relying on the Multiphysics Object Oriented Simulation Environment (MOOSE) for tracking fission product transport from TRISO particles up to the coolant pressure boundary, as well as modeling the transport and potential deposition of these nuclides inside the reactor coolant loop. The proposed computational scheme is applied to estimate source term inventories for a representative 10-MW(thermal) prismatic high-temperature microreactor and is qualitatively compared against known release fractions. In addition to providing an alternate analysis tool, this MOOSE model can help reactor designers quantify the influence of key design parameters relevant for studies of radiological dose consequences.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Cryogenic System Preliminary Design for a 0.5m-Long, Conduction-Cooled Nb3Sn Undulator Magnet Prototype

There are several NbTi superconducting undulator (SCU) magnets currently in operation at the Advanced Photon Source (APS) at Argonne National Laboratory (ANL). The development on Nb3Sn-based superconducting undulator magnets at APS is underway due to the potential to further enhance the performance of the SCUs. Superconducting undulator magnets need to keep temperature gradients minimized in order to retain thermal and operating current margin. We have designed the thermal links for efficient heat conduction using 3D finite element analysis (FEA) simulation in COMSOL Multiphysics software, which was later used for the calculation of the temperature distribution across a 0.5 m long, conduction-cooled, Nb3Sn undulator magnet prototype that includes both conductive heat transfer and radiative heating components. We have modelled the evolution of the thermal properties of the magnet winding as well as other cold parts during cool-down from ambient temperature with the operation of a SHI RDK-415D cryocooler to examine the estimated time that is needed for the cooldown as well as the baseline temperature we could achieve. A key result was that a maximum coil delta T = 0.11 K along the designed thermal links and a temperature range of 3.42 K-3.53 K for the winding were predicted at steady state. A pair of G10 support rods for carrying the magnet as well as a pair of current leads made of metals were also designed and optimized for the testing system, targeting for an operation temperature of 4 K and a coil current of 950 A DC.

Zhang, Danlu↗
DOE OSTI2020

Tomographic Reconstruction of the Neutron Time-Energy Spectrum from a Dense Plasma Focus

A dense plasma focus (DPF) nuclear fusion device is an attractive pulsed neutron source in many applications, due to its relatively large neutron yield, produced in a short time duration. In order to design a DPF that generates neutrons within a specified time profile, or generates a neutron energy spectrum with specific properties, it is necessary to be able to characterize and model the results of the fusion process releasing the neutrons. The time-energy spectrum of fusion neutrons is an ideal quantity to use to validate multiphysics codes that simulate the pinch and fusion processes, because it is a quantity that requires the physics of each of the stages leading up to and ending in the fusion reaction to be simulated correctly. In particular, since DPF fusion neutrons are not monoenergetic -- and there can often be several fusion pinches creating neutrons -- a computer simulation matching high quality neutron spectrum measurements provides great confidence in the fidelity of the simulation. In order to make such a comparison, it is first necessary to have quality measurements from which to infer the spectrum. In this work we pose neutron spectroscopy as the classical tomographic inverse problem from neutron time of flight data at multiple distances, but enhanced by using an additional measurement of the time profile of the fusion pinch near its source and by using detector pairs set up in a geometry that allows for scatter background subtraction. The detector pairs enhance the quality of the time of flight measurements, and the additional constraint posed by the measured time profile allows for reconstructions discretized as finely as the time measurements and in energy as finely as 100 keV, without the problem being underdetermined. We present results from a Deuterium-fueled DPF at the U.S. Department of Energy's Nevada National Security Site and show that we can infer the time-energy spectrum from our measurements for both single and multi-pinch fusion reactions with equal fidelity.

Catenacci, Jared↗
DOE OSTI2023

Numerical Study of Solar Receiver Tube With Modified Surface Roughness for Enhanced and Selective Absorptivity in Concentrated Solar Power Tower

Abstract Concentrated solar power (CSP) is a reliable renewable energy source that is progressively lowering its cost of energy. However, the heat loss due to reflected and emitted radiation hinders the maximum achievable thermal efficiency for solar receiver tubes on the solar tower. Current solar selective coatings cannot withstand the high temperatures that come with state-of-the-art CSP towers often needing to be recoated soon after initial operation. We intend to use Inconel 718 with different additive manufacturing (AM) practices to construct surfaces that allow for more light-trapping to occur. By adjusting printing parameters, we can tailor a surface to allow for more absorption while diminishing emitted radiation heat loss. By using COMSOL Multiphysics, we can generate these theoretical surfaces to emulate a printed surface, and using the coupled Multiphysics we can simulate how the surface dictates radiation properties. Our results show that by having a rougher surface we can enhance the absorptivity of Inconel 718 (IN718) by 38.8%. We expect this work to transform how solar absorber tubes are manufactured without using selective coatings and supplement the US Department of Energy (DOE) 2030 SunShot Initiative.

14 SOLAR ENERGY↗
DOE OSTI2016

Configuration analysis and optimization on multipolar Galatea trap

Multipolar Galatea magnetic trap simulation model was established with the finite element simulation software COMSOL Multiphysics. Analyses about the magnetic section configuration show that better magnetic configuration should make more plasma stay in the weak magnetic field rather than the annular magnetic shell field. Then an optimization model was established with axial electromagnetic force, weak magnetic field area and average magnetic mirror ratio as the optimization goals and with the currents of myxines as design variables. Select appropriate weight coefficients and get optimization results by applying genetic algorithm. Results show that the superiority of the target value of typical application parameters, including the average magnetic mirror can reduce more than 5%, the weak magnetic field area can increase at least 65%, at the same time, axial electromagnetic force acting on the outer myxines can be reduced to less than 50 N. Finally, the results were proved by COMSOL Multiphysics and the results proved the optimized magnetic trap configuration with more plasma in the weak magnetic field can reduce the plasma diffusion velocity and is more conducive for the constraint of plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

CARDINAL

Cardinal is a code that automates the workflow for iteratively solving neutral particle transport via the OpenMC code and computational fluid dynamics via the NekRS code within the Multiphysics Object Oriented Simulation Environment (MOOSE), a finite element framework for solving general partial differential equations. Cardinal establishes a mapping between the Monte Carlo geometry representation used in OpenMC and the spectral element mesh used in NekRS to other physics applications built upon the MOOSE framework, allowing high-resolution particle transport and fluid dynamics physics feedback to the MOOSE application ecosystem. Heat generation rates as solved by OpenMC are used as heat sources in MOOSE, while the temperature and density fields from NekRS are used as conjugate heat transfer boundary conditions and/or source fields in MOOSE. OpenMC, Nek5000, and MOOSE are open source, community-developed simulation codes.

NOVAK, APRILJEAN RODGERS↗
DOE OSTI2017

Clustering of Observed Diurnal Cycles of Precipitation over the United States for Evaluation of a WRF Multiphysics Regional Climate Ensemble

The diurnal cycle of precipitation during the summer season over the contiguous United States is examined in eight distinct regions. These were identified using cluster analysis applied to the diurnal cycle characteristics at 2141 rainfall gauges over the 10-yr period 1991–2000. Application of the clustering technique provides a physically meaningful way of identifying regions for comparison of model results with observations. The diurnal cycle for each region is specified in terms of 1) total precipitation, 2) frequency of precipitation occurrence, and 3) intensity of precipitation per occurrence on an hourly basis averaged over the 10-yr period. The amplitude and phase of each element of the diurnal cycle was obtained from harmonic analysis and has been compared with the results of a 24-member multiphysics ensemble of simulations produced by the Weather Research and Forecast (WRF) Model on a region-by-region basis. Three cumulus schemes, two radiation schemes, two microphysics schemes, and two planetary boundary layer schemes were included in the ensemble. Simulations of total precipitation showed reasonable agreement with observations in regions where the diurnal cycle is directly influenced by solar radiation, (e.g., the U.S. Southeast), but they were less successful in regions where other factors influence the diurnal cycle (e.g., the central United States). The diurnal cycle of precipitation frequency and intensity showed substantial biases in the simulations of all eight regions, namely, overestimation of occurrences and underestimation of intensities. Simulations were sensitive to the cumulus and radiation schemes but were largely insensitive to either microphysics or planetary boundary layer schemes.

Mooney, P. A.↗
DOE OSTI2015

Off-diagonal Jacobian support for Nodal BCs

In this brief note, we describe the implementation of o-diagonal Jacobian computations for nodal boundary conditions in the Multiphysics Object Oriented Simulation Environment (MOOSE) framework. There are presently a number of applications based on the MOOSE framework that solve complicated physical systems of partial dierential equations whose boundary conditions are often highly nonlinear. Accurately computing the on- and o-diagonal Jacobian and preconditioner entries associated to these constraints is crucial for enabling ecient numerical solvers in these applications. Two key ingredients are required for properly specifying the Jacobian contributions of nonlinear nodal boundary conditions in MOOSE and nite element codes in general: 1. The ability to zero out entire Jacobian matrix rows after \normal" assembly has taken place, and eciently replace the contents of those rows with the contributions associ- ated to the nodal boundary conditions. 2. An interface for allowing users to specify both the residual and Jacobian contributions associated to their nonlinear boundary condition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Monitoring, Modeling, and Diagnosis of Alkali-Silica Reaction in Small Concrete Samples

Assessment and management of aging concrete structures in nuclear power plants require a more systematic approach than simple reliance on existing code margins of safety. Structural health monitoring of concrete structures aims to understand the current health condition of a structure based on heterogeneous measurements to produce high-confidence actionable information regarding structural integrity that supports operational and maintenance decisions. This report describes alkali-silica reaction (ASR) degradation mechanisms and factors influencing the ASR. A fully coupled thermo-hydro-mechanical-chemical model developed by Saouma and Perotti by taking into consideration the effects of stress on the reaction kinetics and anisotropic volumetric expansion is presented in this report. This model is implemented in the GRIZZLY code based on the Multiphysics Object Oriented Simulation Environment. The implemented model in the GRIZZLY code is randomly used to initiate ASR in a 2D and 3D lattice to study the percolation aspects of concrete. The percolation aspects help determine the transport properties of the material and therefore the durability and service life of concrete. This report summarizes the effort to develop small-size concrete samples with embedded glass to mimic ASR. The concrete samples were treated in water and sodium hydroxide solution at elevated temperature to study how ingress of sodium ions and hydroxide ions at elevated temperature impacts concrete samples embedded with glass. Thermal camera was used to monitor the changes in the concrete sample and results are summarized.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

BIGHORN Computational Fluid Dynamics Theory, Methodology, and Code Verification & Validation Benchmark Problems

This document presents the theoretical background for a hybrid finite-element / finite-volume fluid flow solver, namely BIGHORN, based on the Multiphysics Object Oriented Simulation Environment (MOOSE) computational framework developed at the Idaho National Laboratory (INL). An overview of the numerical methods used in BIGHORN are discussed and followed by a presentation of the formulation details. The document begins with the governing equations for the compressible fluid flow, with an outline of the requisite constitutive relations. A second-order finite volume method used for solving the compressible fluid flow problems is presented next. A Pressure-Corrected Implicit Continuous-fluid Eulerian (PCICE) formulation for time integration is also presented. The multi-fluid formulation is being developed. Although multi-fluid is not fully-developed, BIGHORN has been designed to handle multi-fluid problems. Due to the flexibility in the underlying MOOSE framework, BIGHORN is quite extensible, and can accommodate both multi-species and multi-phase formulations. This document also presents a suite of verification & validation benchmark test problems for BIGHORN. The intent for this suite of problems is to provide baseline comparison data that demonstrates the performance of the BIGHORN solution methods on problems that vary in complexity from laminar to turbulent flows. Wherever possible, some form of solution verification has been attempted to identify sensitivities in the solution methods, and suggest best practices when using BIGHORN.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

SIERRA Code Coupling Module: Arpeggio User Manual Version 4.44

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2017

MOOSE Framework Speed and Memory Optimizations

The Multiphysics Object Oriented Simulation Environment (MOOSE) framework has been developed to solve a variety of complex simulations ranging from groundwater flows to full reactor simulations. MOOSE excels at solving fully-implicit, fully-coupled, highly nonlinear models using the Jacobian Free Newton-Krylov (JFNK) method relying on years of expertise in several dependent packages including PETSc, libMesh, Hypre and many others. The modular design of MOOSE enables the rapid development of new physics models that can be coupled with relative ease within the framework for building up simulations of different levels of complexity. With the ongoing adoption of MOOSE to new application areas along with streamlining of existing simulations, the priorities must be adjusted to include code optimizations. This includes both execution speed and the memory footprint of MOOSE-based simulations. These two categories can be further broken down into ''micro-optimizations'' or source code optimizations and algorithmic optimizations where source code is restructured to accomplish the same tasks through a more efficient means. However, before any code optimizations can occur, one must have robust way of measuring these critical simulation properties to identify the best candidates for optimizations along with having a way to identify the effectiveness of proposed optimizations to the framework. This report discusses the preliminary work towards the implementation of a robust execution baseline of the framework and derivative applications under a variety of different simulations. In addition, several potential optimization target areas in both the MOOSE framework and dependent libraries have been identified. An implementation strategy and potential benefits are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

FY17 Status Report on NEAMS Neutronics Activities

Under the U.S. DOE NEAMS program, the high-fidelity neutronics code system has been developed to support the multiphysics modeling and simulation capability named SHARP. The neutronics code system includes the high-fidelity neutronics code PROTEUS, the cross section library and preprocessing tools, the multigroup cross section generation code MC2-3, the in-house meshing generation tool, the perturbation and sensitivity analysis code PERSENT, and post-processing tools. The main objectives of the NEAMS neutronics activities in FY17 are to continue development of an advanced nodal solver in PROTEUS for use in nuclear reactor design and analysis projects, implement a simplified sub-channel based thermal-hydraulic (T/H) capability into PROTEUS to efficiently compute the thermal feedback, improve the performance of PROTEUS-MOCEX using numerical acceleration and code optimization, improve the cross section generation tools including MC2-3, and continue to perform verification and validation tests for PROTEUS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

NEAMS-IPL MOOSE Midyear Framework Activities

Multiphysics Object Oriented Simulation Environement (MOOSE) is a modular, pluggable framework for building complex simulations. The ability to add new objects with custom syntax is a core capability that makes MOOSE a powerful platform for coupling multiple applications together within a single environment. Documentation is not only a required component for a complete software quality assurance plan, but is critical to the overall success and wide-spread adoption of the platform by other researchers and institutions. Therefore, the notion of a modular, pluggable framework has been extended to the creation of documentation content. This report details the methodology for defining the NQA-1 compliant requirements in a manner that allows for the automatic generation of the Software Requirement Specification (SRS), Software Design Description (SDD), and Requirement Traceabiltity Matrix (RTM) from any repository revision.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Evaluation of statistical variation of microstructural properties and temperature effects on creep fracture of Grade 91

This report describes three improvements made to a physically-based model for creep and creep rupture in Grade 91 steel developed as part of the Advanced Reactor Technologies program: 1. Progress on transitioning the model framework to the MOOSE finite element package to improve parallel scalability, increase the physical size of the simulations, and incorporate multiphysics effects into the model. 2. The development of a physically-based method of scaling the model parameters to accurately capture creep in Grade 91 for temperatures in the range of 450° to 500° C. 3. Extending the model to account for microstructural statistical variations, specifically capturing the effect of grain boundary energy on the physical parameters underling the model for grain boundary void nucleation and cavitation. Put together, these improvements result in an accurate, predictive model for the physical response of Grade 91 steel over the expected use temperature range for the material in future liquid metal cooled fast reactors. The model can be used to more accurately predict engineering properties of the material for very long service lives, which could lead to safer, more economical future advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗