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

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Prediction of The TREAT Facility Power and Temperature of Sirius-2c Nuclear Propulsion Fuel Experiment

This paper presents simulation results of the Sirius-2c full power experiment that was performed at the TREAT facility to support experiment modeling and transient analysis of the NASA-sponsored Sirius series of experiments for nuclear thermal propulsion applications. We presents a novel approach to predict the control rod motion and power of the reactor and the temperature of the specimen during transient experiments using the Griffin and Bison packages of the Multiphysics Object Oriented Simulation Environment (MOOSE). First, the model predicts the transient rod motion that results in the desired reactor power shape. Then, the temperature of the specimen is predicted based on the power deposited in the specimen. The prediction results are in good agreement with the measured reactor power and control rod position and over-predicted the specimen temperature which is mainly related to specimen surface emissivities. The reactor reaches a peak power of 80 MW and a specimen temperature of 2500 K and considering the right specimen emissivities the predicted temperature is in an excellent agreement with the measured values. This model can be utilized in helping and supporting the design and optimization a new experiments that are going to be irradiated in the TREAT facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advanced Reactors-Intermediate Heat Exchanger (IHX) Coupling: Theoretical Modeling and Experimental Validation

The overall goal of the research project was to model the behavior of the advanced reactorintermediate heat exchange system and to develop advanced control techniques for off-normal conditions. The specific objectives defined for the project were: 1. To develop the steady-state thermal hydraulic design of the intermediate heat exchanger (IHX); 2. To develop mathematical models to describe the advanced nuclear reactor-IHX-chemical process/power generation coupling during normal and off-normal operations, and to simulate models using multiphysics software; 3. To develop control strategies using genetic algorithm or neural network techniques and couple these techniques with the multiphysics software; 4. To validate the models experimentally The project objectives were accomplished by defining and executing four different tasks corresponding to these specific objectives. The first task involved selection of IHX candidates and developing steady state designs for those. The second task involved modeling of the transient and offnormal operation of the reactor-IHX system. The subsequent task dealt with the development of control strategies and involved algorithm development and simulation. The last task involved experimental validation of the thermal hydraulic performances of the two prototype heat exchangers designed and fabricated for the project at steady state and transient conditions to simulate the coupling of the reactor- IHX-process plant system. The experimental work utilized the two test facilities at The Ohio State University (OSU) including one existing High-Temperature Helium Test Facility (HTHF) and the newly developed high-temperature molten salt facility.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Conservative tightly-coupled simulations of stochastic multiscale systems

Multiphysics problems often involve components whose macroscopic dynamics is driven by microscopic random fluctuations. The fidelity of simulations of such systems depends on their ability to propagate these random fluctuations throughout a computational domain, including subdomains represented by deterministic solvers. When the constituent processes take place in nonoverlapping subdomains, system behavior can be modeled via a domain-decomposition approach that couples separate components at the interfaces between these subdomains. Its coupling algorithm has to maintain a stable and efficient numerical time integration even at high noise strength. We propose a conservative domain-decomposition algorithm in which tight coupling is achieved by employing either Picard's or Newton's iterative method. Coupled diffusion equations, one of which has a Gaussian white-noise source term, provide a computational testbed for analysis of these two coupling strategies. Fully-converged (“implicit”) coupling with Newton's method typically outperforms its Picard counterpart, especially at high noise levels. This is because the number of Newton iterations scales linearly with the amplitude of the Gaussian noise, while the number of Picard iterations can scale superlinearly. At large time intervals between two subsequent inter-solver communications, the solution error for single-iteration (“explicit”) Picard's coupling can be several orders of magnitude higher than that for implicit coupling. Increasing the explicit coupling's communication frequency reduces this difference, but the resulting increase in computational cost can make it less efficient than implicit coupling at similar levels of solution error, depending on the communication frequency of the latter and the noise strength. This trend carries over into higher dimensions, although at high noise strength explicit coupling may be the only computationally viable option.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Polarity effects and apparent ion recombination in microionization chambers

Purpose: Microchambers demonstrate anomalous voltage-dependent polarity effects. Existing polarity and ion recombination correction factors do not account for these effects. As a result, many commercial microchamber models do not meet the specification of a reference-class ionization chamber as defined by the American Association of Physicists in Medicine. The purpose of this investigation is to determine the cause of these voltage-dependent polarity effects. Methods: A series of microchamber prototypes were produced to isolate the source of the voltage-dependent polarity effects. Parameters including ionization-chamber collecting-volume size, stem and cable irradiation, chamber assembly, contaminants, high-Z materials, and individual chamber components were investigated. Measurements were performed with electrodes coated with graphite to isolate electrode conductivity. Chamber response was measured as the potential bias of the guard electrode was altered with respect to the collecting electrode, through the integration of additional power supplies. Ionization chamber models were also simulated using COMSOL Multiphysics software to investigate the effect of a potential difference between electrodes on electric field lines and collecting volume definition. Results: Investigations with microchamber prototypes demonstrated that the significant source of the voltage-dependent polarity effects was a potential difference between the guard and collecting electrodes of the chambers. The voltage-dependent polarity effects for each prototype were primarily isolated to either the guard or collecting electrode. Polarity effects were reduced by coating the isolated electrode with a conductive layer of graphite. Polarity effects were increased by introducing a potential difference between the electrodes. COMSOL simulations further demonstrated that for a given potential difference between electrodes, the collecting volume of the chamber changed as the applied voltage was altered, producing voltage-dependent polarity effects in the chamber response. Ionization chamber measurements and COMSOL simulations demonstrated an inverse relationship between the chamber collecting volume size and the severity of voltage-dependent polarity effects on chamber response. The effect of a given potential difference on chamber polarity effects was roughly ten times greater for microchambers as compared to Farmer-type chambers. Stem and cable irradiations, chamber assembly, contaminants, and high-Z materials were not found to be a significant source of the voltage-dependent polarity effects. Conclusions: A potential difference between the guard and collecting electrodes was found to be the primary source of the voltage-dependent polarity effects demonstrated by microchambers. For a given potential difference between electrodes, the relative change in the collecting volume is smaller for larger-volume chambers, illustrating why these polarity effects are not seen in larger-volume chambers with similar guard and collecting electrode designs. Thus, for small-volume chambers, it is necessary to reduce the potential difference between the guard and collecting electrodes in order to reduce polarity effects for reference dosimetry measurements.

60 APPLIED LIFE SCIENCES↗

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Designing, Modeling, Constructing, and Testing a Flat Panel Speaker and Sound Diffuser for a Simulator

The goal of this project was to design, model, build, and test a flat panel speaker and frame for a spherical dome structure being made into a simulator. The simulator will be a test bed for evaluating an immersive environment for human interfaces. This project focused on the loud speakers and a sound diffuser for the dome. The rest of the team worked on an Ambisonics 3D sound system, video projection system, and multi-direction treadmill to create the most realistic scene possible. The main programs utilized in this project, were Pro-E and COMSOL. Pro-E was used for creating detailed figures for the fabrication of a frame that held a flat panel loud speaker. The loud speaker was made from a thin sheet of Plexiglas and 4 acoustic exciters. COMSOL, a multiphysics finite analysis simulator, was used to model and evaluate all stages of the loud speaker, frame, and sound diffuser. Acoustical testing measurements were utilized to create polar plots from the working prototype which were then compared to the COMSOL simulations to select the optimal design for the dome. The final goal of the project was to install the flat panel loud speaker design in addition to a sound diffuser on to the wall of the dome. After running tests in COMSOL on various speaker configurations, including a warped Plexiglas version, the optimal speaker design included a flat piece of Plexiglas with a rounded frame to match the curvature of the dome. Eight of these loud speakers will be mounted into an inch and a half of high performance acoustic insulation, or Thinsulate, that will cover the inside of the dome. The following technical paper discusses these projects and explains the engineering processes used, knowledge gained, and the projected future goals of this project

Dillon, Christina↗

Initial Neutronics and Thermal-Hydraulic Coupling for Spent Nuclear Fuel Canister

This report documents work performed supporting the US Department of Energy (DOE) Nuclear Energy Spent Fuel and Waste Disposition, Spent Fuel and Waste Science and Technology, under work breakdown structure element 1.08.01.03.05, “Direct Disposal of Dual Purpose Canisters.” In particular, this appendix fulfills the M3 milestone, M3SF-19OR010305015, “Multiphysics criticality consequence analysis capability development status report,” as Revision 1 to M3SF-19OR010305016, “Initial neutronic and thermal hydraulic coupling for waste package” within work package SF-19OR01030501, “Direct Disposal of Dual Purpose Canisters–ORNL.” This report presents the initial development status of a multiphysics criticality consequence simulation framework, Terrenus. Terrenus currently couples two physics codes: the Monte Carlo radiation transport code Shift and the thermal-hydraulics code COBRA-SFS. The coupling capability has been demonstrated using a simplified 3 x 3 fuel pin lattice. Terrenus will be developed further to include (1) general geometry package support, allowing modeling of a spent nuclear fuel (SNF) cask with all structural details, (2) a depletion solver to determine the change in nuclide composition at the end of each time step, and (3) a mechanics code to determine any structural impact due to a criticality event. Future research and development works include identifying/developing a modern two-phase thermal-hydraulics code and developing an approach to solve fast neutronic transients. The goal of this multiphysics framework is to determine the feasibility of direct disposal of currently loaded SNF canisters by including or excluding a criticality event from a repository performance analysis framework in terms of consequences. The revision 1 to this report includes extending Terrenus to a 17 x 17 pressurized water reactor (PWR) fuel assembly. We also investigated the phenomena necessary to quench a critical reaction, and investigated the water level necessary to achieve criticality in a typical canister. Finally, we describe the initial implementation of a critical power search capability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Lessons Learned from Open-Source Software Training Toward Expanding the Fusion Workforce

Within the United States (U.S.), the Fusion Innovation Research Engine (FIRE) Collaboratives seek to accelerate fusion technology development through wide-ranging community-driven research activities that bring together industry, laboratories, research institutes, and academia. Increasing the maturity of key components and systems for future fusion power plants (FPPs) toward commercialization, will, by necessity, increase the need for knowledgeable engineers, designers, and researchers in industry capable of integrating and further improving these technologies within industry FPP concepts. Further, various modeling and simulation packages support these programs and are under-development within them to drive design iteration and the development of FPP digital twins. To derive the greatest benefit from model output and capabilities, training these same specialists will be vital. Thus, the development of effective and accessible software onboarding and professional development opportunities focused on fusion will be key to keeping the pace of growth high in the coming decade and beyond. A similar need exists within the advanced fission reactor community, driven by an ever-growing need for carbon-free baseload power for industrial and data center applications. Within the U.S. and around the world, several open-source packages and frameworks exist to support this endeavor, and two we will highlight here are the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and OpenMC. These packages, notably, are also being utilized in the FIRE Collaboratives program. In this presentation, we will highlight the lessons learned from over 30 years of combined software development and training experience, focused on developing strong technical software foundations within the nuclear workforce, from students to professional engineers & scientists. We will connect this experience to present and emerging needs in the fusion energy community, and, finally, will outline possible paths forward to develop a large, robust, global fusion workforce.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

High-to-Low Activities across VVI, AMA and PHI

This milestone report documents work performed in CASL that was broadly focused on the use of validated high-fidelity simulation codes to inform, improve and calibrate low-fidelity codes to facilitate multiphysics coupling, design, performance analysis, and risk assessment. This high-to-low (Hi2Lo) initiative has the following objectives: employ physics encapsulated in high-fidelity codes to improve the functionality of low-fidelity codes, quantify biases in low-fidelity codes, develop and calibrate physical models and closure relations so that low-fidelity codes more accurately incorporate requisite physics, and develop low-fidelity codes that can be effectively coupled to simulate multiphysics behavior. In combination, these efforts are designed to impact both CIPS and DNB. For example, the development and experimental design-based Bayesian calibration of dynamic gas conductance models in CTFFuel will provide the capability to better model crud-induced power shifts. Whereas the present emphasis has focused on the high-fidelity codes STAR-CCM+ and BISON and low-fidelity codes CTF and CTFFuel, future emphasis will be expanded to include MAMBA. During the reporting period, Hi2Lo efforts focused on two broad components: demonstration of the experimental design-based Bayesian calibration framework to calibrate CTF using both STAR-generated synthetic data and experimental data, and the use of STAR to inform and improve the turbulence modeling and grid heat transfer modeling capabilities in CTF. A complementary milestone focused on the implementation and verification of fuel models in CTFFuel, for which we demonstrated initial calibration results. Future Hi2Lo work will focus on the use of the experimental design-based Bayesian calibration framework, demonstrated this year for STAR and CTF, to calibrate parameters in the fuel swelling, fuel irradiation-induced densification, and fuel relocation models in CTFFuel using both experimental and BISON-generated data. A second component will focus on the continued use of STAR to improve grid modeling capabilities in CTF.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessing the Effect of Parametric Uncertainty on the FRAPCON-3.4 Thermal Conductivity Model

In the absence of a well-accepted, mathematical description of how units of material respond to one or more external actions, scientists must often rely upon empirically-derived models to make predictions. These models are often employed to model very complex phenomena, with an unfortunate side-effect that intuition, as to the values of particular model parameters, may not serve to guide their use. Given our increased reliance on modeling and simulation to make predictions in the absence of experimental data, it befits the scientific and engineering communities to explore and report upon uncertainty quantification techniques applied to previously adopted (or at least well-accepted) empirical models that derive from or pertain to substantial experimental data sets. This report represents a collection of three methodologies aimed at assessing the predictive capability of the empirical thermal conductivity model adopted by the nuclear fuel performance code, FRAPCON-3.4. Each of these methodologies considers the effect of uncertain parameters – a plausible reality in the context of empirically-derived models – on the ability of the model to predict uranium dioxide conductivity data from open literature sources. The results lead the authors to question the predictive capability of the FRAPCON model for predicting the thermal conductivities associated with irradiated fuel samples. The report concludes with a preliminary examination of Idaho National Laboratory’s (INL) nuclear fuel performance code, BISON, developed under INL’s Multiphysics Object Oriented Simulation Environment (MOOSE).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parallel Algorithms and Software for Nuclear, Energy, and Environmental Applications Part I: Multiphysics Algorithms

There is a growing trend within energy and environmental simulation to consider tightly coupled solutions to multiphysics problems. This can be seen in nuclear reactor analysis where analysts are interested in coupled flow, heat transfer and neutronics, and in nuclear fuel performance simulation where analysts are interested in thermomechanics with contact coupled to species transport and chemistry. In energy and environmental applications, energy extraction involves geomechanics, flow through porous media and fractured formations, adding heat transport for enhanced oil recovery and geothermal applications, and adding reactive transport in the case of applications modeling the underground flow of contaminants. These more ambitious simulations usually motivate some level of parallel computing. Many of the physics coupling efforts to date utilize simple code coupling or first-order operator splitting, often referred to as loose coupling. While these approaches can produce answers, they usually leave questions of accuracy and stability unanswered. Additionally, the different physics often reside on distinct meshes and data are coupled via simple interpolation, again leaving open questions of stability and accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High fidelity multiphysics tightly coupled model for a lead cooled fast reactor concept and application to statistical calculation of hot channel factors

A tightly coupled multiphysics code system is established using the MOOSE framework for hot channel factor (HCF) evaluation on a Lead Fast Reactor (LFR) concept. The coupled system is driven by the Griffin multiphysics coupling capability under which the MOOSE Heat Transfer module and NekRS computational fluid dynamics solver are coupled for conjugate heat transfer using the Cardinal application. The coupled capability is demonstrated on an LFR assembly model based on materials and geometry of a prototypical lead-cooled fast reactor design by Westinghouse Electric Company, LLC. Moreover, the work integrates the Multiphysics Object Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM) to perform calculations for statistical analysis of HCF. Furthermore, the coupling strategy and workflow demonstrated in this paper is not only useful for predicting accurate hot channel factors for different kinds of advanced reactors but also for other engineering applications such as control rod worth assessment, generation of high-fidelity database for Artificial intelligence (AI)/machine learning (ML) training, design optimization and multi-resolution modeling.

Cardinal↗

Modeling quantum cascade lasers: Coupled electron and phonon transport far from equilibrium and across disparate spatial scales

Quantum cascade lasers (QCLs) are high-power coherent light sources in the midinfrared and terahertz parts of the electromagnetic spectrum. They are devices in which the electronic and lattice systems are far from equilibrium, strongly coupled to one another, and the problem bridges disparate spatial scales. Here, we present our ongoing work on the multiphysics and multiscale simulation of far-from-equilibrium transport of charge and heat in midinfrared QCLs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗

Multiphysics modeling of two-phase film boiling within porous corrosion deposits

Porous corrosion deposits on nuclear fuel cladding, known as CRUD, can cause multiple operational problems in light water reactors (LWRs). CRUD can cause accelerated corrosion of the fuel cladding, increase radiation fields and hence greater exposure risk to plant workers once activated, and induce a downward axial power shift causing an imbalance in core power distribution. In order to facilitate a better understanding of CRUD's effects, such as localized high cladding surface temperatures related to accelerated corrosion rates, we describe an improved, fully-coupled, multiphysics model to simulate heat transfer, chemical reactions and transport, and two-phase fluid flow within these deposits. Our new model features a reformed assumption of 2D, two-phase film boiling within the CRUD, correcting earlier models' assumptions of single-phase coolant flow with wick boiling under high heat fluxes. This model helps to better explain observed experimental values of the effective CRUD thermal conductivity. Finally, we propose a more complete set of boiling regimes, or a more detailed mechanism, to explain recent CRUD deposition experiments by suggesting the new concept of double dryout specifically in thick porous media with boiling chimneys. - Highlights: • A two-phase model of CRUD's effects on fuel cladding is developed and improved. • This model eliminates the formerly erroneous assumption of wick boiling. • Higher fuel cladding temperatures are predicted when accounting for two-phase flow. • Double-peaks in thermal conductivity vs. heat flux in experiments are explained. • A “double dryout” mechanism in CRUD is proposed based on the model and experiments.

42 ENGINEERING↗

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↗