Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation codes”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The MCNP6 Analytic Criticality Benchmark Suite

Analytical benchmarks provide an invaluable tool for verifying computer codes used to simulate neutron transport. Several collections of analytical benchmark problems [1-4] are used routinely in the verification of production Monte Carlo codes such as MCNP® [5,6]. Verification of a computer code is a necessary prerequisite to the more complex validation process. The verification process confirms that a code performs its intended functions correctly. The validation process involves determining the absolute accuracy of code results vs. nature. In typical validations, results are computed for a set of benchmark experiments using a particular methodology (code, cross-section data with uncertainties, and modeling) and compared to the measured results from the set of benchmark experiments. The validation process determines bias, bias uncertainty, and possibly additional margins. Verification is generally performed by the code developers, while validation is generally performed by code users for a particular application space. The VERIFICATION_KEFF suite of criticality problems [1,2] was originally a set of 75 criticality problems found in the literature for which exact analytical solutions are available. Even though the spatial and energy detail is necessarily limited in analytical benchmarks, typically to a few regions or energy groups, the exact solutions obtained can be used to verify that the basic algorithms, mathematics, and methods used in complex production codes perform correctly. The present work has focused on revisiting this benchmark suite. A thorough review of the problems resulted in discarding some of them as not suitable for MCNP benchmarking. For the remaining problems, many of them were reformulated to permit execution in either multigroup mode or in the normal continuous-energy mode for MCNP. Execution of the benchmarks in continuous-energy mode provides a significant advance to MCNP verification methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Arc Parameter Modeling (Arc Breakdown Experiments and Modeling: Confined Geometries and Alternate Gap Materials)

Sandia Nationa.l Laboratories (SNL) and The Boeing Company (Boeing) collaborated on this PTS to develop the capability to predict thermal and electrical characteristics of short arcs that may appear in equipment subjected to transient currents of hundreds or thousands of amperes. The final tasks of this CRADA included experiments and computer simulations used to investigate the effect of an arc breakdown between two planar electrodes and confined by various potting materials. The arc is generated by a current waveform with a rise of 8 μs and a 20 μs fall time. Voltage, current and light emissions are recorded. The developing shock wave is observed with an optical fiber embedded in the potting material and recorded using a photonic Doppler velocimetry (PDV) diagnostic system. The experiments were modeled using the ALEGRA-MHD simulation code and comparisons are made between the experiments and simulations. This product contains Protected CRADA Information which was produced December 2015 under CRADA No. 1651. l 5 and is not to be further disclosed for a period of 5 years from the date it was produced except as expressly provided for in the CRADA. Further dissemination only as authorized by the originating partners: the Boeing Company and Sandia National Laboratories Radiation Sciences Center, or higher U.S. Department of Energy programmatic authority. Please note that the project leads from both SNL and Boeing have all retired or left for other reasons at the time of the creation of this document. The final project manager collected information for this document from the final proposal and the final report.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Flexible Tool for Beam Induced Background Simulations at a Muon Collider

A Muon Collider represents a very interesting possibility for a future machine to explore the energy frontier in particle physics.However, to reach the needed luminosity, beam intensities of the order of $10^{9-11}$ muons per bunch are needed. In this context, the Beam Induced Background must be taken into account for its effects on the magnets and on the detector. Several optimizations can be conceived with the aim to mitigate them. In this view, it is of crucial importance to develop a flexible tool that allows to easily reconstruct the machine geometry in a Monte Carlo code, allowing to simulate in detail the interaction of muon decay products in the machine, while being able to change the machine optics itself to find the best configuration.In this contribution, a possible approach to such a purpose is presented, based on FLUKA for the Monte Carlo simulation and on LineBuilder for the geometry reconstruction.First results for the 1.5TeV machine optics developed by the MAP collaboration are discussed in this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1974

Geothermal reservoir simulation

The prediction of long-term geothermal reservoir performance and the environmental impact of exploiting this resource are two important problems associated with the utilization of geothermal energy for power production. Our research effort addresses these problems through numerical simulation. Computer codes based on the solution of partial-differential equations using finite-element techniques are being prepared to simulate multiphase energy transport, energy transport in fractured porous reservoirs, well bore phenomena, and subsidence.

Mercer, J. W., Jr.↗
DOE OSTI2016

Uncertainty quantification and sensitivity analysis applications to fuel performance modeling

Best-estimate fuel performance codes such as BISON currently under development at the Idaho National Laboratory, utilize empirical and mechanistic lower-length-scale informed correlations to predict fuel behavior under normal operating and accident reactor conditions. Traditionally, best-estimate results are presented using the correlations with no quantification of the uncertainty in the output metrics of interest. However, there are associated uncertainties in the input parameters and correlations used to determine the behavior of the fuel and cladding under irradiation. Therefore, it is important to perform uncertainty quantification and include confidence bounds on the output metrics that take into account the uncertainties in the inputs. In addition, sensitivity analyses can be performed to determine which input parameters have the greatest influence on the outputs. In this paper we couple the BISON fuel performance code to the DAKOTA uncertainty analysis software to analyze a representative fuel performance problem. The case studied in this paper is based upon rod 1 from the IFA-432 integral experiment performed at the Halden Reactor in Norway. The rodlet is representative of a BWR fuel rod. The input parameters uncertainties are broken into three separate categories including boundary condition uncertainties (e.g., power, coolant flow rate), manufacturing uncertainties (e.g., pellet diameter, cladding thickness), and model uncertainties (e.g., fuel thermal conductivity, fuel swelling). Utilizing DAKOTA, a variety of statistical analysis techniques are applied to quantify the uncertainty and sensitivity of the output metrics of interest. Specifically, we demonstrate the use of sampling methods, polynomial chaos expansions, surrogate models, and variance-based decomposition. The output metrics investigated in this study are the fuel centerline temperature, cladding surface temperature, fission gas released, and fuel rod diameter. The results highlighted the importance of the thermal conductivity of the fuel, the thermal conductivity of the cladding, the gap thickness, the fuel thermal expansion, and the system power factor as influential parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Validation of the Modelling of Corrosion Driven Expansion in Grouted Waste-forms - 16125

One of the current preferred solutions for the disposal of Intermediate Level radioactive waste in the UK is grout encapsulation inside steel boxes for long term storage. These boxes are then to be entombed in a Geological Disposal Facility once one is available. However, both prior to and during geological disposal, waste contents can react with the water present in the encapsulant grout and begin to corrode. This corrosion product occupies a larger volume (by several hundred percent) than the original material and the corrosion occurs over the course of decades and centuries, long after the encapsulant grout has set and hardened. The expanding waste then creates stresses in the encapsulant causing both it and the surrounding box to deform and, in extreme circumstances, and after long times, for the structural integrity of the box to fail. This behaviour needs to be quantified before the box can be approved for use as any excessive deformation could lead to a loss of containment or, since the boxes are stacked on top of each other, to a significant accident in the stores. Due to the long timescales involved, physical testing of the phenomenon in full scale is impossible. Instead, a combination of chemical modelling and finite element analysis is used to predict the extent of the corrosion in the waste and its effect of the integrity of the box. There are a number of difficulties in modelling the corrosion-driven expansion and subsequent deformation of the boxes. Although the corrosion occurs over a very long time scale, the event is highly non-linear and dynamic. All of the material models, including the steel and grout are non-linear. The deformation is in the order of 10s of centimetres, well beyond the static elastic range. Finally, the components must be allowed to move relative to each other which involve a lot of part to part contacts. This means that simple static modelling, which would solve quite quickly, cannot be used. Instead a dynamic analysis code, LS-DYNA, has been used to expand the product over time to simulate the corrosion-driven expansion. This software is generally used for modelling car crashes or scenarios where waste containers have been dropped. The analyses therefore include a significant inertial component. This paper will detail the NNL's approach to modelling this long-term, quasi-static phenomenon in the dynamic, impact oriented software code LS-DYNA. It will cover the problems inherent in this approach, the errors that it creates as well as the tests and solutions that NNL worked on in order to overcome them and be confident that the solution could be validated. This work will ultimately be used to determine the maximum inventory which can be stored in each box based on the waste content and the potential extent of expansion as well as any modification that could be made to the box design or the processing streams that may improve the storage solution. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Dynamic load balancing with enhanced shared-memory parallelism for particle-in-cell codes

Furthering our understanding of many of today’s interesting problems in plasma physics – including plasma based acceleration and magnetic reconnection with pair production due to quantum electrodynamic effects – requires large-scale kinetic simulations using particle-in-cell (PIC) codes. However, these simulations are extremely demanding, requiring that contemporary PIC codes be designed to efficiently use a new fleet of exascale computing architectures. To this end, the key issue of parallel load balance across computational nodes must be addressed. We discuss the implementation of dynamic load balancing by dividing the simulation space into many small, self-contained regions or ‘‘tiles,’’ along with shared-memory (e.g., OpenMP) parallelism both over many tiles and within single tiles. The load balancing algorithm can be used with three different topologies, including two space-filling curves. Here, we tested this implementation in the code Osiris and show low overhead and improved scalability with OpenMP thread number on simulations with both uniform load and severe load imbalance. Compared to other load-balancing techniques, our algorithm gives order-of-magnitude improvement in parallel scalability for simulations with severe load imbalance issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Simulation of Gas Production from Multilayered Hydrate-Bearing Media with Fully Coupled Flow, Thermal, Chemical and Geomechanical Processes Using TOUGH+Millstone. Part 3: Production Simulation Results

The TOUGH+Millstone simulator has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH 4 - /sub>-hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH+Hydrate simulator, v2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the evolution of hydrate-bearing systems and includes the most recent physical properties relationships, coupled seamlessly with (b) Millstone v1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of the TOUGH+Millstone code are demonstrated in the simulation and analysis of the system flow, thermal, and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the third paper of this series, we apply the simulators described in parts 1 and 2 to a problem of gas production from a complex, multilayered system of hydrate-bearing sediments in an oceanic environment. We perform flow simulations of constant-pressure production via a vertical well and compare those results to a coupled flow-geomechanical simulation of the same process. Finally, the results demonstrate the importance of fully coupled geomechanics when modeling the evolution of reservoir properties during production.

58 GEOSCIENCES↗
DOE OSTI2016

PONDEROMOTIVE ACCELERATION IN CORONAL LOOPS

Ponderomotive acceleration has been asserted to be a cause of the first ionization potential (FIP) effect, the well-known enhancement in abundance by a factor of 3–4 over photospheric values of elements in the solar corona with FIP less than about 10 eV. It is shown here by means of numerical simulations that ponderomotive acceleration occurs in solar coronal loops, with the appropriate magnitude and direction, as a “by-product” of coronal heating. The numerical simulations are performed with the HYPERION code, which solves the fully compressible three-dimensional magnetohydrodynamic equations including nonlinear thermal conduction and optically thin radiation. Numerical simulations of coronal loops with an axial magnetic field from 0.005 to 0.02 T and lengths from 25,000 to 75,000 km are presented. In the simulations the footpoints of the axial loop magnetic field are convected by random, large-scale motions. There is a continuous formation and dissipation of field-aligned current sheets, which act to heat the loop. As a consequence of coronal magnetic reconnection, small-scale, high-speed jets form. The familiar vortex quadrupoles form at reconnection sites. Between the magnetic footpoints and the corona the reconnection flow merges with the boundary flow. It is in this region that the ponderomotive acceleration occurs. Mirroring the character of the coronal reconnection, the ponderomotive acceleration is also found to be intermittent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Multigroup Scattering in Monte Carlo Radiation Transport Codes [Slides]

Legendre truncations to multigroup scattering distributions are not amenable to Monte Carlo sampling due to negative values. We have implemented two moment-preserving methods in MGMC that capture the shape of the truncation, are non-negative over [-1, 1], and can be efficiently sampled on CPUs and GPUs. MGMC can now simulate neutrons with anisotropic scattering mechanics and MGMC shows good agreement with LANL production codes PARTISN.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Estimation of Fission Product Transport Parameters for Cesium in the AGR-3/4 TRISO Fuel Experiment

A one-dimensional (1D) finite-element model of fission product transport in the AGR-3/4 experiment has been developed using the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and implemented in the fuel performance code, BISON. The model resolves capsule-specific geometries, materials, and temperature histories and simulates radial migration of fission products from the fuel compact through the inner ring, outer ring, and into the sink ring. Model parameters governing diffusion and sorption were estimated for key fission products – cesium (Cs), and europium (Eu) – by simultaneously fitting modeled isotopic concentration profiles and total ring inventories to a post-irradiation experimental measurement. These data include gamma scanning, liquid scintillation for Sr-90, radial deconsolidation leach-burn-leach analysis, tomographic reconstructions, and destructive physical sampling. A mortar-based interfacial sorption framework was implemented to enforce physically consistent mass transfer and flux conservation across gas gaps. Two classes of parameter sets were derived: a least-squares best-fit, and a safety-oriented conservative-fit, what applies strong penalties for underprediction of sink inventories. Across all twelve capsules, the model successfully reproduces the dominant radial transport trends for Cs, Sr, with decreasing concentrations from the compact outward through successive rings. Cs behavior is captured most consistently, while strontium predictions reveal systematic trade-offs between compact accuracy and conservative sink-ring bounding. The results demonstrate that sink ring weighted calibration provides conservative, safety relevant bounds on low temperature fission product transport, but at the cost of underpredicting compact inventories for Sr isotopes. These discrepancies highlight the need for additional physics, including fast-slow diffusion model, incorporating trapping mechanism in the transport behavior. Overall, this work establishes a robust, capsule-specific modeling framework for AGR-3/4 fission product transport and provides a defensible basis for parameter selection in source-term and fuel performance analyses for high temperature gas-cooled reactors.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Enhanced Performance Fast Reactors with Engineered Passive Safety Systems

This report summarizes the work and results of NEUP Project #15-8251, "Enhanced Performance Fast Reactors with Engineered Passive Safety System," performed from 2015-2018 at UC Berkeley and Argonne National Laboratory. This project concerns the incorporation of a novel engineered passive safety system, called the Autonomous Reactivity Control (ARC) system, into Sodium-cooled Fast Reactors (SFRs) to improve their transient and economic performance. The first part of this project consisted of fleshing out the basic principles and the manufacturing process for an SFR fuel assembly with the ARC system installed, with the details provided in high-quality, step-by-step CAD drawings. Once the detailed design is laid out, high-fidelity simulations were performed to optimize the dimensions of important ARC system components and to quantify the standalone time-response of ARC systems in a typical SFR environment. These high-fidelity simulations were used to create reduced-order models of ARC system response which could be incorporated into the SAS4A/SASSYS-1 whole-core transient models. In order to accurately account for the time-response of ARC systems in whole-core transient simulations, modifications were made to the SAS4A/SASSYS-1 transient code. These new capabilities were tested and incorporated into the production version of the code, and have been used extensively in this analysis. Before performing whole-core transient calculations, the two reference cores selected for this study to cover a wide design space were simulated and characterized. The two reference cores, (1) a medium-sized oxide-fueled Advanced Burner Reactor and (2) a large metal-fueled Breed-and-Burn core, challenge conventional means of achieving passive safety. Because the two cores vary substantially, simulating the ARC system in both allows for conclusions to be drawn about the applicability of ARC systems in different SFRs. Once both cores were fully characterized, transient calculations were performed to understand the complicated interplay between the ARC system and other inherent reactivity feedback mechanisms. By performing parametric studies of the transient behavior, ARC system designs that best improve transient performance in each core were identified. Once transient improvements were established, the transient calculations were extended to improve reactor economics through design modifications. Upon establishing the potential benefits of ARC system inclusion, the study examines the negative neutronic implications of incorporating ARC systems into the reference cores. The study concludes by discussing a number of knowledge gaps that exist which may have important implications for ARC system implementation. Ultimately, the results of this project indicate that ARC systems may provide an effective means to improve transient performance in certain types of systems. The transient simulations in the oxide Advanced Burner Reactor showed that significant margin to coolant boiling and fuel melting could be gained through the use of a properly designed ARC system, thus improving inherent safety. While oscillatory behavior was induced in the transient response for certain combinations of ARC system parameters, the oscillation-free design space that can be effectively used to improve performance was found to be large, and the additional margin was demonstrated to allow economical improvements to the core design. Conversely, difficulties were encountered in the Breed-and-Burn core due to fundamental aspects of the core design. The combination of a particularly strong positive coolant void feedback and the use of metallic fuel allowed for diverging oscillations to manifest in response to certain transient initiators, and in the end it was found that the current ARC system design could not simultaneously enable passive safety in all transients examined. In order for complete passive safety to be achieved, a mechanism for slowing the reverse actuation of ARC systems is necessary. Regardless, the standard ARC system was able to improve the performance of two out of the three examined transients by increasing the margin to boiling in Loss of Heat Sink transients and avoiding boiling in Transient Overpower accidents. Additionally, no significant barriers to implementation were found, although the issue of tritium production should be closely monitored in any core incorporating ARC systems. Still, a number of technology gaps exist, and these are discussed at the conclusion of the report to recommend directions for future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Warthog: Coupling Nek5000 Thermal Hydraulics to BISON Fuel Performance through the Giraffe Interface

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program under the US Department of Energy (DOE) consists of three product lines. The Fuels Product Line (FPL), hosted at Idaho National Laboratory (INL), has produced a multiphysics framework with the name Multiphysics Object Oriented Simulation Environment (MOOSE) (Gaston et al. 2015). The FPL has also produced several high-fidelity physics codes including the fuel performance code Bison (Hales et al. 2016). The Reactors Product Line (RPL), hosted at Argonne National Laboratory (Argonne), has developed the high-order spectral-element Computational Fluid Dynamics (CFD) code Nek5000 (Fischer et al. 2015). The Integration Product Line (IPL), hosted at Oak Ridge National Laboratory (ORNL), is tasked with ensuring tools from the other product lines can be coupled together so that users may perform high-fidelity multiphysics simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Simulations of particle interactions in a high-current RFQ

High-intensity beams of composite ions are susceptible to particle interactions that do not affect single ions, such as stripping of electrons and charge exchange. Beam dynamics simulation techniques used in accelerator design do not often include these particle interactions, and so the transmission of the real beam can be significantly different from the design simulations. We have modified the IMPACT-T code to include particle interactions as part of the beam dynamics simulations, including space-charge effects and beam losses. We have carried out validation simulations of interactions in a molecular hydrogen ion beam through an existing high-intensity deuteron RFQ design, as beam commissioning is often carried out using molecular hydrogen beams, to avoid deuteron-deuteron interactions in the accelerator leading to neutron activation. The new code framework successfully models the interactions within the beam, tracks the different source and product ions through the accelerator, and calculates the collective effects such as space charge from all particles together. Finally, we also discuss other simulation codes that could be modified to include similar interaction effects, and possible applications to other types of accelerators.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Data From: Influence of Hydrological Perturbations and Riverbed Sediment Characteristics on Hyporheic Zone Respiration of CO2 and N-2, Journal of Geophysical Research-Biogeosciences

This data package contains pumping data (.txt), parameter matrices, and R code (.R, .RData) to perform bootstrapping for parameter selection for the bioclogging model development. The pumping data were collected from the Russian River Riverbank Filtration site located in Sonoma County, California from 2010-2017 from three riverbank collection wells located alongside the study site. The pumping data is directly correlated with water table oscillations, so the code performs these correlations and simulates stochastic versions of water table oscillations. See Metadata Description.pdf for full details on dataset production. This dataset must be used with the R programming language. This dataset and R code is associated with the publication "Influence of Hydrological Perturbations and Riverbed Sediment Characteristics on Hyporheic Zone Respiration of CO2 and N-2"This research was supported by the Jane Lewis Fellowship from the University of California, Berkeley, the Sonoma County Water Agency (SCWA), the Roy G. Post Foundation Scholarship, the U.S. Department of Energy, Office of Science Graduate Student Research (SCGSR) Program, U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under award DE-AC02-05CH11231, and the UFZ-Helmholtz Centre for Environmental Research, Leipzig, Germany.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Assessment of buffer-IPyC thermomechanical debonding behavior using new experimental strength data in BISON

TRIstructural ISOtropic (TRISO) fuel is a nuclear fuel commonly used in High Temperature Gas-cooled Reactors (HTGRs). A single sub-millimeter-diameter TRISO fuel particle consists of a spherical fuel kernel surrounded by four coating layers: a low-density pyrocarbon buffer layer, an inner pyrolytic carbon (IPyC) layer, a silicon carbide (SiC) layer, and an outer pyrolytic carbon (OPyC) layer. The kernel is commonly made of UO2 or a mixture of uranium carbide and uranium oxide (UCO). During reactor operation, the TRISO coating layers are subjected to irradiation-induced dimensional changes and the associated thermomechanical behavior of each layer. One of the observed behaviors is gap formation between the buffer and IPyC layer due to the porous buffer’s irradiation-induced shrinkage exceeding that of the IPyC layer. Not all irradiated particles will experience buffer-IPyC gap formation. The debonding may be partial, or it may be nearly total. However, from post-irradiation examination of UCO TRISO fuels irradiated as part of the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program, it was concluded that partial buffer-IPyC debonding was the most common type of buffer-IPyC interaction. To predict TRISO thermomechanical performance, multi-physics models have been built that are being continually updated and refined. The BISON code is a finite element-based nuclear fuel performance code that may be used for 1D, 2D, and 3D TRISO particle simulations. This code is used to calculate fuel temperature, kernel swelling, buffer densification, thermal and irradiation creep, fracture, and fission gas production and release. One of the recent additions to the BISON code is the ability to model the process of layer debonding. This paper will focus on the simulation results of the improved BISON debonding model that will utilize updated strengths measured from irradiated AGR TRISO fuel particles. The new experimental strength data from micromechanical tests of irradiated TRISO fuel samples were exercised in the BISON simulations and compared to baseline strength data to assess their applicability in the models. This also includes updated buffer-IPyC bond strengths to simulate layer delamination. Based on current experimental observations it is noted that the buffer-IPyC separation occurs not exactly at the junction of these two layers, but more on the side of the buffer layer. That observation is also implemented in the TRISO interface debonding model. This improved modeling approach using experimental strength data to characterize buffer-IPyC debonding and its potential subsequent cracking will be presented in the paper along with comparisons to available experimental observations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Comparison of equilibrium ohmic and nonequilibrium swarm models for monitoring conduction electron evolution in high-altitude EMP calculations

Here, atmospheric electromagnetic pulse (EMP) events are important physical phenomena that occur through both man-made and natural processes. Radiation-induced currents and voltages in EMP can couple with electrical systems, such as those found in satellites, and cause significant damage. Due to the disruptive nature of EMP, it is important to accurately predict EMP evolution and propagation with computational models. CHAP-LA (Compton High Altitude Pulse-Los Alamos) is a state-of-the-art EMP code that solves Maxwell inline images equations for gamma source-induced electromagnetic fields in the atmosphere. In EMP, low-energy, conduction electrons constitute a conduction current that limits the EMP by opposing the Compton current. CHAP-LA calculates the conduction current using an equilibrium ohmic model. The equilibrium model works well at low altitudes, where the electron energy equilibration time is short compared to the rise time or duration of the EMP. At high altitudes, the equilibration time increases beyond the EMP rise time and the predicted equilibrium ionization rate becomes very large. The ohmic model predicts an unphysically large production of conduction electrons which prematurely and abruptly shorts the EMP in the simulation code. An electron swarm model, which implicitly accounts for the time evolution of the conduction electron energy distribution, can be used to overcome the limitations exhibited by the equilibrium ohmic model. We have developed and validated an electron swarm model previously in Pusateri et al. (2015). Here we demonstrate EMP damping behavior caused by the ohmic model at high altitudes and show improvements on high-altitude, upward EMP modeling obtained by integrating a swarm model into CHAP-LA.

54 ENVIRONMENTAL SCIENCES↗