Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphysics”

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 739 records · Page 41

Modeling and experimental validation of a W f /W-fabrication by chemical vapor deposition and infiltration

Tungsten (W) has a unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall material in future fusion devices. However, its intrinsic brittleness and its further embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) with extrinsic toughening is being developed. A possible synthesis route is chemical vapor deposition (CVD) using heated W fabrics as substrate. The challenge is that the growing CVD-W can isolate domains from precursor access leading to strength-reducing pores. To deepen the process understanding and to optimize the CVD parameters, models were developed with COMSOL Multiphysics and validated experimentally. W deposition rate equations as function of the temperature and the partial pressures of the precursors H 2 and WF 6 were experimentally validated in previous work. In the present article, these equations are applied to obtain partial pressures within the CVD reactor. The results are taken as input for transient simulations in the microscale, in which W coatings, growing onto multiple adjacent W fibers, were simulated via mesh deformation and remeshing. The surface-to-surface contact of the W coatings and the corresponding potential pore formation were simulated by implementing sophisticated deposition rate stop conditions. Within the measuring uncertainties of ≃ ±1%, the models are validated successfully by experimental comparison regarding the deposition rate, pore structure, and relative densities ranging from 0.6 to 0.9.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of the oxygen reduction reaction for dense LSM thin films

In this study, the oxygen reduction reaction mechanism is investigated using numerical methods on a dense thin (La 1-x Sr x ) y MnO 3±δ film deposited on a YSZ substrate. This 1-D continuum model consists of defect chemistry and elementary oxygen reduction reaction steps coupled via reaction rates. The defect chemistry model contains eight species including cation vacancies on the A- and B-sites. The oxygen vacancy is calculated by solving species transportation equations in multiphysics simulations. Due to the simple geometry of a dense thin film, the oxygen reduction reaction was reduced to three elementary steps: surface adsorption and dissociation, incorporation on the surface, and charge transfer across the LSM/YSZ interface. The numerical simulations allow for calculation of the temperature- and oxygen partial pressure-dependent properties of LSM. The parameters of the model are calibrated with experimental impedance data for various oxygen partial pressures at different temperatures. The results indicate that surface adsorption and dissociation is the rate-determining step in the ORR of LSM thin films. With the fine-tuned parameters, further quantitative analysis is performed. The activation energy of the oxygen exchange reaction and the dependence of oxygen non-stoichiometry on oxygen partial pressure are also calculated and verified using the literature results.

36 MATERIALS SCIENCE↗

Critical Assemblies: Dragon Burst Assembly and Solution Assemblies

This work reviews the historical literature associated with the Dragon experiment and water boiler reactors operated at Los Alamos during the Manhattan Project. Frisch’s invited talk given at the American Nuclear Society’s Fast Burst Reactor Conference held at the University of New Mexico in Albuquerque, New Mexico, in 1969 is quoted. From the literature review, basic models for the Dragon experiment and for a water boiler–type assembly (aqueous homogeneous reactor) were created that can be used for conducting multiphysics simulations for criticality excursion studies. This methodology utilizes the coupled neutronic-hydrodynamic method to perform a time-dependent dynamic simulation of a criticality excursion. MCNP® was utilized to calculate important nuclear kinetic parameters that were incorporated into the models. Simulation results compare reasonably well with historic data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Structured large-pore foams improve thermal performance of LiMIT-style liquid lithium PFC

As magnetically confined fusion devices improve, the conditions at the walls become increasingly intense. Plasma facing components (PFCs) must withstand these extreme heat and particle loads without damage or degradation. Liquid lithium PFCs are known to be quite resilient, and the presence of lithium also serves to improve plasma properties. The liquid metal infused trench (LiMIT) concept is an open surface liquid lithium PFC design that has been tested extensively at the University of Illinois and in fusion devices around the world. LiMIT utilizes thermoelectric magnetohydrodynamics (TEMHD) to passively drive liquid lithium flow. This work demonstrates an extension of the LiMIT trench geometry to three dimensions. Additively manufactured large pore metallic foams maintain TEMHD drive while drastically improving heat flux handling and resistance to lithium dryout, a phenomenon where locally high TEMHD forces depresses the lithium level and exposes underlying solid structure. COMSOL multiphysics modeling of the system yields insight into the forces at play in dryout development, and shows the 3D structures can eliminate dryout. Here, low heat proof-of-concept experimental testing of the system matches computational results, and high heat flux electron beam tests more than double the proven operational range of a LiMIT-style PFC, to 6.8 MW m –2 , with no indications of dryout or impending damage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetostrictive Ultrasonic Waveguide Transducer for In-Pile Thermometry

Real-time and reliable temperature measurement on nuclear fuels is crucial for the safe operation of existing pressurized-water reactors and future advanced nuclear reactors. Magnetostrictive materials deform when subjected to a magnetic field or exhibit magnetization variation when stressed. Here, based on these properties, this study prototyped an ultrasonic thermometer (UT) consisting of a magnetostrictive waveguide, a dc coil providing appropriate magnetic biasing, and an ac coil generating an acoustic impulse and detecting the resulting acoustic echoes. By tracking the time of flight between the excitation pulse and the echoes, the UT can potentially detect nuclear fuel cladding temperature from a long distance. In this study, magnetostrictive iron–gallium alloys, or Galfenol, were selected as the waveguide due to their large magnetostriction, superior temperature survivability, excellent radiation resilience, and high mechanical robustness. A multiphysics finite-element model considering electrical, magnetic, and mechanical dynamics in the magnetostrictive UT was then developed. The model exhibited an error of 0.24% in time-of-flight simulation and, therefore, enabled computer-aided design and guided signal processing. Between room temperature and 120 °C, the new Galfenol-based UT exhibits a linear sensitivity of 162.8×10 ₋6 °C ₋1 , which is 51.7% higher than a previous magnetostrictive UT based on iron–cobalt–vanadium alloys.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Design and Analysis Progress of ITER Diagnostic Equatorial Port #09

ITER is the world's largest fusion device currently under construction in the South of France with >50 diagnostic systems to be installed inside the port plugs (PPs), the interspace (IS), or the port cell region of various diagnostic ports. The plasma facing diagnostic first wall (DFW) and its supporting diagnostic shield modules (DSM) are designed to protect front-end diagnostics from plasma neutron and plasma radiation, while providing apertures for diagnostic access to the plasma. The design of ITER port plug structures including the DFW and the DSM is largely driven by the electromagnetic loads included on the PP structural components during plasma major disruptions and the vertical displacement events (VDEs). Unlike DFW and DSM, the design of diagnostic system, however, is likely driven by the steady-state thermal loads from plasma volumetric and surface heating and the dynamic response of the in-port components attached to the port-specific DSM or closure plate under transient loads induced on the vacuum vessel and the port extension during asymmetric VDEs. Three tenant diagnostic systems are integrated into the equatorial port 09. The toroidal interferometer/polarimeter, or TIP system, is installed in the left drawer (DSM3, left looking from plasma) for measuring the plasma density so to control the fuel input. The electron cyclotron emission (ECE) system is installed in the middle drawer (DSM2) to provide the high spatial and temporal resolution measurements of electron temperature evolution and the electron thermal transport inferences. The visible/infrared wide angle viewing system is installed in the right drawer (DSM1) to provide visible and infrared viewing and temperature data of the first wall for its protection in support of machine operation. The port plug integration design and multiphysics analyses are performed following port integration requirements including the weight limit (45 tones total), shut down dose rate limits, the cooling/heating and structural integrity validation. Mass distribution for TIP and ECE DSMs has been optimized to minimize the total weight by a new design of the boron carbide shielding pocket. Furthermore, the lightened DSM maintains its front-end load distribution and the structural stiffness with minimum impact to the DFW so to better protect on-board diagnostics; while still provides sufficient front-end stiffness for its structural integrity as well as the diagnostics function requirements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Three-Dimensional, Thermodynamically and Variationally Consistent, Fully Coupled, Electro-Chemo-Thermo-Mechanical Model of Solid-State Batteries

Here we present a theoretical framework for high fidelity modeling of three-dimensional electro-chemo-thermo-mechanical coupled systems such as solid-state batteries. The framework provides a useful means of utilizing molecular level information to predict the overall performance of the system at the continuum level. The central statement of the theoretical framework is a thermodynamically consistent inf-sup problem. The total potential of the coupled system consists of the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. All governing equations derive from variational principles. We also show that the resulting initial-boundary-value problem (IBVP) is consistent with existing, widely adopted equations for modeling batteries. Moreover, the IBVP can be naturally discretized using the (variational) finite element method to simulate practical solid-state batteries with arbitrary geometries. Several examples illustrate multiphysics coupling between pairs of the physical field equations, culminating in a comprehensive three-dimensional model of charge-discharge in a solid state battery.

25 ENERGY STORAGE↗

Implementation of Creep-Damage Model for Concrete Fracture in MOOSE

Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.

36 MATERIALS SCIENCE↗

The MOOSE Thermal Hydraulics Module

The Multiphysics Object-Oriented Simulation Environment (MOOSE) is an open-source object-oriented finite element framework written in C++ (Lindsay et al., 2022). The Thermal Hydraulics Module (THM) is an optional MOOSE physics module that provides capabilities for studying thermal hydraulic systems. Its core capability lies in assembling a network of coupled components, for instance, pipes, junctions, and valves. THM provides several new systems to MOOSE to enable and facilitate thermal hydraulic simulations, most notably the Components system, which provides a higher-level syntax to MOOSE’s lower-level objects. This system is extensible by the user, but the current library primarily includes components based on a one-dimensional, single-phase, variable-area, compressible flow model, as well as heat conduction.

42 ENGINEERING↗

Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuels High Impact Problem: Coordinate Multiscale U 3 Si 2 Modeling

Since the events at the Fukushima-Daiichi nuclear power plant in March 2011 significant research has unfolded at national laboratories, universities and other institutions into alternative materials that have potential enhanced accident tolerance when compared to traditional \uo~fuel zircaloy clad fuel rods. One of the potential replacement fuels is uranium silicide (\usi) for its higher thermal conductivity and uranium density. The lower melting temperature is of potential concern during postulated accident conditions. Another disadvantage for \usi~ is the lack of experimental data under power reactor conditions. Due to the aggressive development schedule for inserting some of the potential materials into lead test assemblies or rods by 2022~\cite{bragg-sitton_2014} multiscale multiphysics modeling approaches have been used to provide insight into these materials. \\ \noindent The purpose of this letter report is to highlight the multiscale modeling effort for \usi~fuel as part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuel (ATF) High Impact Problem (HIP). This fiscal year two new models for \usi~fuel have been incorporated into the BISON fuel performance code~\cite{Williamson2012} based upon lower length scale simulations using the GRASS-SST rate theory code. The models are briefly described in the following sections.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuels High Impact Problem: Coordinate Multiscale FeCrAl Modeling

Since the events at the Fukushima-Daiichi nuclear power plant in March 2011 significant research has unfolded at national laboratories, universities and other institutions into alternative materials that have potential enhanced ac- cident tolerance when compared to traditional UO2 fuel zircaloy clad fuel rods. One of the potential replacement claddings are iron-chromium-alunimum (FeCrAl) alloys due to their increased oxidation resistance [1–4] and higher strength [1, 2]. While the oxidation characteristics of FeCrAl are a benefit for accident tolerance, the thermal neu- tron absorption cross section of FeCrAl is about ten times that of Zircaloy. This neutronic penalty necessitates thinner cladding. This allows for slightly larger pellets to give the same cold gap width in the rod. However, the slight increase in pellet diameter is not sufficient to compensate for the neutronic penalty and enriching the fuel beyond the current 5% limit appears to be necessary [5]. Current estimates indicate that this neutronic penalty will impose an increase in fuel cost of 15-35% [1, 2]. In addition to the neutronic disadvantage, it is anticipated that tritium release to the coolant will be larger because the permeability of hydrogen in FeCrAl is about 100 times higher than in Zircaloy [6]. Also, radiation-induced hardening and embrittlement of FeCrAl need to be fully characterized experimentally [7]. Due to the aggressive development schedule for inserting some of the potential materials into lead test assemblies or rods by 2022 [8] multiscale multiphysics modeling approaches have been used to provide insight into these the use of FeCrAl as a cladding material. The purpose of this letter report is to highlight the multiscale modeling effort for iron-chromium-alunimum (FeCrAl) cladding alloys as part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program through its Accident Tolerant Fuel (ATF) High Impact Problem (HIP). The approach taken throughout the HIP is to utilize lower length scale approaches (e.g., density functional theory, cluster dynamics, rate theory, phase field, and Visco-Plastic- Self-Consistent (VPSC)) to develop more physically informed models at the engineering scale for use in the BISON [9] fuel performance code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

This report details a fully coupled neutronics thermal hydraulics reference plant model for a gas-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission's Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory high performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON MOOSE-based applications. Griffin provides the reactor physics capabilities including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient solutions. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions. BISON solves thermal conduction problems for the average pebble and TRISO in the pebble-bed core thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel and moderator temperatures. We present results for the coupled steady state equilibrium core and a protected loss of flow event. Although this model is to be considered prototypical with regard to the capabilities in BlueCRAB, its results are consistent with published work by the Institute of Nuclear and New Energy Technology (INET) in China and other research entities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Useability and Optimization Improvements in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this broad use, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed or will be deployed in the coming months. The capabilities were implemented in the same order as they are covered in this report: increased support of face variables, arbitrary spatial and temporal evaluation of material properties, and the addition of a triangular meshing library in libMesh.

97 MATHEMATICS AND COMPUTING↗

NEAMS Technical Area Support in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include a factor of 10^4 improvement in dependency resolution speed, sorting of user objects, ability to compute residuals and Jacobians together, transfer fixes, support for the mortar method in finite volume discretizations, addition of generalized advection schemes for fluid simulations, 10^2 speedup in some Griffin simulations due to a new matrix-only solve type, and much more.

97 MATHEMATICS AND COMPUTING↗

Improvements in High Temperature Gas Cooled Reactor Modeling Capabilities in the Pronghorn Code

This report details the improvement of pebble bed reactor modeling capabilities in the Pronghorn code in fiscal year 2022. The following accomplishments are reported: Deployment of weakly compressible finite volume formulation to the HTR- PM reference plan model; Enable modeling of stagnant gas gaps in the finite volume formulation; Enable using all Pronghorn correlations available in the finite element version in the finite volume version; Modeling of decay heat in pebble bed reactors; Simplifying the input for multiphysics equilibrium core calculations and significant reduction of execution time; Implementation of advanced correlations developed by the Center of Excellence for Thermal-Fluids Applications in Nuclear Energy . In addition, this report includes a development plan for Pronghorn and associated NEAMS tools for prismatic gas-cooled reactors.

97 MATHEMATICS AND COMPUTING↗

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗

Griffin/Pronghorn PBMR-400 Benchmark Results [Slides]

Content Overview: 1) Pronghorn & Griffin, MOOSE based applications overview: Codes capabilities description & MOOSE ecosystem overiew; Current research activities and synergies with other projects; Different multiphysics coupling scheme description. 2) PBMR-400 Benchmark: OECD NEA Benchmark and Pronghorn/Griffin models description; Standalone NK Ex1 and TH Ex2 Results description; Coupled Steady state Ex3; Results analysis for Transiet Ex4a (100-40-100 Load Following) & Ex3 (Protected PLOFC). 3) Conclusions $ Future steps.

97 MATHEMATICS AND COMPUTING↗