Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model reduction”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Model-based Analysis of Mixed Uranium(VI) Reduction by Biotic and Abiotic Pathways During in Situ Bioremediation

Uranium bioremediation has emerged as a potential strategy of cleanup of radionuclear contamination worldwide. An integrated geochemical & microbial community model is a promising approach to predict and provide insights into the bioremediation of a complicated natural subsurface. In this study, an integrated column-scale model of uranium bioremediation was developed, taking into account long-term interactions between biotic and abiotic processes. It is also combined with a comprehensive thermodynamic analysis to track the fate and cycling of biogenic species. As compared with other bioremediation models, the model increases the resolution of the connection of microbial community to geochemistry and establishes direct quantitative correlation between overall community evolution and geochemical variation, thereby accurately predicting the community dynamics under different sedimentary conditions. The thermodynamic analysis examined a recently identified homogeneous reduction of U(VI) by Fe(II) under dynamic sedimentary conditions across time and space. It shows that the biogenic Fe(II) from Geobacter metabolism can be removed rapidly by the biogenic sulphide from sulfate reducer metabolism, hence constituting one of the reasons that make the abiotic U(VI) reduction thermodynamically infeasible in the subsurface. Further analysis indicates that much higher influent concentrations of both Fe(II) and U(VI) than normal are required to for abiotic U(VI) reduction to be thermodynamically feasible, suggesting that the abiotic reduction cannot be an alternative to the biotic reduction in the remediation of uranium contaminated groundwater.

bioremediation↗
DOE OSTI2021

Improved MPACT energy deposition and explicit heat generation coupling with CTF

The CASL neutronics code MPACT assumes local energy deposition with equilibrium delayed energy for steady-state and transient calculations. These approximations limit the accurate representation of the heat generation in space and its variations over time, both of which are essential for power distribution and thermal–hydraulic coupling analyses. This paper presents an explicit energy deposition model considering the neutron and gamma heating of all regions and the time-dependent delayed energy. The MPACT/CTF interfaces are updated to transfer the explicit heat sources. The new model agrees well with MCNP for problems without thermal–hydraulic feedback. The MPACT/CTF coupled calculation for a hot-full-power quarter-core case using the new model exhibits a reduction of peak pin power by 2.3% and fuel centerline temperature by 17 K. The new model also shows a keff difference up to -100 pcm in depletion and a 30 K peak fuel temperature reduction in a reactivity-initiated-accident problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

How Noninnocent Spectator Species Improve the Oxygen Reduction Activity of Single-Atom Catalysts: Microkinetic Models from First-Principles Calculations

Graphene-based single-atom catalysts are promising alternatives to platinum-based catalysts for fuel cell applications. Different transition metals have been screened using electronic structure methods by estimating onset potentials from the most endergonic elementary reaction step. Here, we calculate onset potentials for the oxygen reduction reaction on metal atoms embedded in N-substituted graphene di-vacancies by virtue of first-principles-informed microkinetic analysis. We find that for more oxophilic metals (Cr, Fe, Mn, and Ru), purely thermodynamic models systematically underestimate onset potentials. Furthermore, the oxophilic metals (Cr, Fe, Mn, and Ru) are oxidized under reaction conditions, leading to an increase in activity compared to their reduced state. Importantly, coadsorbed O m H n species actively participate in the reaction, which requires a dynamic treatment of spectator species. These findings highlight the limitations of thermodynamic analyses for electrocatalytic processes, which commonly assume the same oxidation state for each metal, and show that deviations between computational and experimental onset potentials cannot be solely attributed to the shortcomings of the electronic structure methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Assessment of Technoeconomic Opportunities in Automation for Nuclear Microreactors

Achieving full decarbonization of all economic sectors remains a challenge, especially in niche markets. For example, remote communities and industrial or mining activities detached from the main electric grid heavily rely on fossil fuels, similar to urban and industrial microgrids with combined heat and power needs. A combination of renewables and energy storage is often not suitable due to cost, reliability, intermittency, and large storage requirements. Small nuclear reactors with a flexible purpose could serve these applications. Microreactors (MR) are a class of reactors that are compact, factory manufactured, transportable, and self-regulating. Typically, they generate much less power than their large reactor counterparts. The main advantages of microreactors include the versatile nature of the energy produced, the reliability of supply, and freedom from having to transport and store large quantities of fuels on-site, coupled with the absence of dependence on an electrical grid. A strong business case is needed to move from the microreactor prototype to the commercialization phase. In fact, fossil fuels are still relatively inexpensive, and in the near term, carbon credits will be available to virtually compensate for emissions. For microreactors, one of the main costs in operation and maintenance (O&M) is their staffing levels. In this study, we investigate how to optimize the number (and thus the cost) of workers, moving from a traditional, fully manned, on-site personnel approach to an unmanned, remote personnel approach. We examine four different staffing models that can be implemented as the technology matures and evolves. We estimate the staffing needs of each model and build a business case to justify the substitution of on-site personnel with adequate technologies. To do so, we propose a cost model to quantify potential cost reductions from automating O&M activities. The model accounts for both the reduction in cost derived from the reduced number of full-time-equivalent (FTE) employees and the increase in cost derived from the need to buy new control hardware as needed. Applying the cost model that we created to different scenarios, an on-site O&M cost reduction exceeding 80% can be expected. Additionally, we found that it is more impactful to focus on automating routine O&M tasks rather than attempting to automate transient management (shutdowns, restarts, monitoring condition deviations). In fact, transients typically account for less than 1% of the total FTE time spent on the reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Elucidating Mass Transport Regimes in Gas Diffusion Electrodes for CO 2 Electroreduction

Gas diffusion electrodes (GDEs) have shown promising performance for the electrochemical reduction of CO 2 (CO 2 R). In this study, a resolved, pore scale model of electrochemical reduction of CO 2 within a liquid-filled catalyst layer is developed. Three CO 2 mass transport regimes are identified in which the CO 2 penetration depth is controlled by CO 2 consumption in the electrolyte, CO 2 conversion along the solid-electrolyte double-phase boundaries (DPBs), and CO 2 conversion concentrated around the gas–solid–electrolyte triple-phase boundaries (TPBs). While it is possible for CO 2 R to be localized around the TPBs, in systems with submicron pore radii operating at <1 A cm –2 CO 2 R will be distributed across the DPBs within the catalyst layer. This validates the assumption of pore-scale uniformity implicit in popular, volume-averaged GDE models. Overall, the CO 2 conversion efficiency depends strongly on the governing mass transport regime, and operational-phase diagrams are constructed to guide the catalyst layer design.

10 SYNTHETIC FUELS↗
DOE OSTI2023

Experimental Validation of Models of a Hull-Based Tuned Mass Damper System for a Semisubmersible Floating Offshore Wind Turbine Platform

Floating offshore wind turbine designs can be further optimized if the controller and remaining systems are designed together, known as control co-design. To effectively perform control co-design, modeling tools predicting the influence of the control systems on the response of a system must be validated. This article presents an experimental validation that utilizes a scale model of a semisubmersible platform for an offshore wind turbine that is fitted with adjustable tuned mass dampers. These dampers can be tuned to attenuate either the hull-pitch resonance or the tower-bending resonance. The data from the experiment are used to validate state-of-the-art modeling tools. It is shown that the models capture the overall effects of the tuned mass dampers; however, some models overpredict the reduction in platform pitch motion when the dampers are tuned to the pitch resonance. The relative reduction in the tower-base bending moment is more consistently captured by the models when the dampers are tuned to the tower-bending resonance. However, there are significant differences in the absolute level of tower-base bending moment between the models and the experiment. Much of the differences observed are a consequence of the challenge with accurately predicting the baseline platform resonance motion and the tower-bending moment, which should be addressed in future modeling efforts.

17 WIND ENERGY↗
NASA NTRS2016

Probabilistic Prognosis of Non-Planar Fatigue Crack Growth

Quantifying the uncertainty in model parameters for the purpose of damage prognosis can be accomplished utilizing Bayesian inference and damage diagnosis data from sources such as non-destructive evaluation or structural health monitoring. The number of samples required to solve the Bayesian inverse problem through common sampling techniques (e.g., Markov chain Monte Carlo) renders high-fidelity finite element-based damage growth models unusable due to prohibitive computation times. However, these types of models are often the only option when attempting to model complex damage growth in real-world structures. Here, a recently developed high-fidelity crack growth model is used which, when compared to finite element-based modeling, has demonstrated reductions in computation times of three orders of magnitude through the use of surrogate models and machine learning. The model is flexible in that only the expensive computation of the crack driving forces is replaced by the surrogate models, leaving the remaining parameters accessible for uncertainty quantification. A probabilistic prognosis framework incorporating this model is developed and demonstrated for non-planar crack growth in a modified, edge-notched, aluminum tensile specimen. Predictions of remaining useful life are made over time for five updates of the damage diagnosis data, and prognostic metrics are utilized to evaluate the performance of the prognostic framework. Challenges specific to the probabilistic prognosis of non-planar fatigue crack growth are highlighted and discussed in the context of the experimental results.

Leser, Patrick E.↗
DOE OSTI2026

Optimal Operation of Residential High Performance Water Heater for Reduction of Electricity Cost and Peak Demand Through Field Validation

Water heating accounts for about 18% of a typical US home’s energy use. Modern water heaters have enabled control options through APIs, offering customers the opportunity to reduce their energy cost and peak demand by dynamically adjusting settings. A water heater’s capacity to store energy using its storage tank makes it an asset for peak demand reduction and energy cost savings. For this reason, a mixed-integer linear programming model is proposed to minimize the energy cost of a high-performance water heater while also reducing the peak demand of the residential household under a time-of-use utility rate by dynamically changing the water heater’s running mode. Specifically, a multi-objective optimization model is formulated to determine the mode settings of the water heater considering hot water use, time-of-use rate, and peak demand limit of the residential household. The mode settings are associated with different dead bands of water temperature for triggering on/off action of the heat pump and heating element. A 66-gal hybrid electric high performance water heater was used for numerical simulation and practical experiments. The simulation results were well aligned with measurements of practical experiments, validating the soundness of the thermodynamic model. In addition, reductions of energy cost, enabling affordability, and reducing peak demand are demonstrated. The research team also developed a software framework with dashboards to automatically and continuously monitor and manage devices.

Liu, Guodong [ORNL] (ORCID:0000000213498608)↗
DOE OSTI2022

Reinforcement learning for online adaptation of model predictive controllers: Application to a selective catalytic reduction unit

Here we present a novel application of reinforcement learning (RL) for online dynamic tuning of model predictive controllers (MPC). Applying a state-action-reward-state-action (SARSA) algorithm for temporal difference learning with a control-specific reward function improves the error tracking performance of a standard MPC formulation. The proposed RL approach is also readily adaptable to other MPCs, or entirely different control approaches. Practical details for the implementation of the RL-MPC algorithm are also presented. The proposed algorithm is applied to a case study of controlling nitrogen oxide (NO x ) emissions in an industrial selective catalytic reduction (SCR) unit, a control problem characterized by significant nonlinearity and time delay. Along with an RL-MPC formulation for NOx control, another MPC is proposed to mitigate ammonia slip and decrease ammonia consumption in the SCR. Results showing the efficacy of the RL-MPC for NO x control through learning and implementation on the nonlinear SCR dynamic model are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2013

Effect of Two Advanced Noise Reduction Technologies on the Aerodynamic Performance of an Ultra High Bypass Ratio Fan

A wind tunnel experiment was conducted in the NASA Glenn Research Center anechoic 9- by 15-Foot Low-Speed Wind Tunnel to investigate two new advanced noise reduction technologies in support of the NASA Fundamental Aeronautics Program Subsonic Fixed Wing Project. The goal of the experiment was to demonstrate the noise reduction potential and effect on fan model performance of the two noise reduction technologies in a scale model Ultra-High Bypass turbofan at simulated takeoff and approach aircraft flight speeds. The two novel noise reduction technologies are called Over-the-Rotor acoustic treatment and Soft Vanes. Both technologies were aimed at modifying the local noise source mechanisms of the fan tip vortex/fan case interaction and the rotor wake-stator interaction. For the Over-the-Rotor acoustic treatment, two noise reduction configurations were investigated. The results showed that the two noise reduction technologies, Over-the-Rotor and Soft Vanes, were able to reduce the noise level of the fan model, but the Over-the-Rotor configurations had a significant negative impact on the fan aerodynamic performance; the loss in fan aerodynamic efficiency was between 2.75 to 8.75 percent, depending on configuration, compared to the conventional solid baseline fan case rubstrip also tested. Performance results with the Soft Vanes showed that there was no measurable change in the corrected fan thrust and a 1.8 percent loss in corrected stator vane thrust, which resulted in a total net thrust loss of approximately 0.5 percent compared with the baseline reference stator vane set.

Hughes, Christoper E.↗
DOE OSTI2019

Investigation of LSM-YSZ Composite Cathode Performance Degradation with a Multistep Charge Transfer Model

The long-term performance degradation of LSM/YSZ composite cathodes is investigated via calibrated multiphysics simulation with a multistep oxygen reduction reaction (ORR) model and structural coarsening data from a phase field study. Multiphysics simulations, including a multistep ORR mechanism, are first developed for better understanding of the limiting processes in solid oxide fuel cells. The multistep ORR mechanism includes two parallel pathways, surface pathway and bulk pathway, which consist of elementary reaction steps. The multiphysics model is simultaneously calibrated with experimental polarization curves and impedance behavior for various air/fuel supply conditions. To our knowledge, this is the first time that a multistep ORR model is calibrated with such datasets. Next, the calibrated simulations are utilized to simulate a 2D half-cell constructed with measured microstructural data and random heterogeneity. Finally, the long-term performance degradation of the half-cell is predicted by the calibrated multiphysics model coupled with structural coarsening trends simulated using a phase field-based coarsening model. Degradation of both polarization curves and impedance behavior is investigated. Thorough analyses, including changes of contributions from different pathways, the resistance components, and overall reaction order, are performed to provide more insights into cathode performance degradation due to grain coarsening phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Global and Regional Radiative Forcing from 20 Reductions in BC, OC and SO4 an HTAP2 Multi-Model Study

In the Hemispheric Transport of Air Pollution Phase 2 (HTAP2) exercise, a range of global atmospheric general circulation and chemical transport models performed coordinated perturbation experiments with 20% reductions in emissions of anthropogenic aerosols, or aerosol precursors, in a number of source regions. Here, we compare the resulting changes in the atmospheric load and vertically resolved profiles of black carbon (BC), organic aerosols (OA) and sulfate (SO4/ from 10 models that include treatment of aerosols. We use a set of temporally, horizontally and vertically resolved profiles of aerosol forcing efficiency (AFE) to estimate the impact of emission changes in six major source regions on global radiative forcing (RF) pertaining to the direct aerosol effect, finding values between. 51.9 and 210.8mW/sq m/Tg for BC, between -2.4 and -17.9mW/sq m/Tg for OA and between -3.6 and -10.3W/sq m/Tg for SO4. In most cases, the local influence dominates, but results show that mitigations in south and east Asia have substantial impacts on the radiative budget in all investigated receptor regions, especially for BC. In Russia and the Middle East, more than 80 % of the forcing for BC and OA is due to extra-regional emission reductions. Similarly, for North America, BC emissions control in east Asia is found to be more important than domestic mitigations, which is consistent with previous findings. Comparing fully resolved RF calculations to RF estimates based on vertically averaged AFE profiles allows us to quantify the importance of vertical resolution to RF estimates. We find that locally in the source regions, a 20% emission reduction strengthens the radiative forcing associated with SO4 by 25% when including the vertical dimension, as the AFE for SO4 is strongest near the surface. Conversely, the local RF from BC weakens by 37% since BC AFE is low close to the ground. The fraction of BC direct effect forcing attributable to intercontinental transport, on the other hand, is enhanced by one-third when accounting for the vertical aspect, because long-range transport primarily leads to aerosol changes at high altitudes, where the BC AFE is strong. While the surface temperature response may vary with the altitude of aerosol change, the analysis in the present study is not extended to estimates of temperature or precipitation changes.

profiles of aerosol forcing efficiency (AFE)↗
DOE OSTI2026

Multi-Scale Modeling and Prototype Development for Electrochemical CO2 Reduction (CRADA Final Report)

In this CRADA project, Lawrence Livermore National Laboratory, Stanford University, SLAC National Laboratory, and TotalEnergies collaboratively executed a multidisciplinary investigation of electrochemical reduction of CO2 to produce sustainable fuels and chemicals. Overall, the project led to an increased understanding of the fundamental processes involved in CO2 electrolysis, from the atomistic scale to the full electrolyzer device scale, ultimately leading to design guidelines for CO2 electrolyzers that will help in their future commercialization. As the model systems, Ag- and Cu-based catalysts were investigated in various forms depending on the electrochemical platform that was utilized to study the activity, selectivity, and durability towards electrochemical CO2 reduction. By employing experimental, theoretical, and computational techniques, the project team experimentally validated multi-physics models, evaluated the experimental levers that lead to increased electrolyzer reaction selectivity and energy efficiency, and used computational optimization to design higher performance electrodes. The learnings of this project were extensively documented in publicly available peer-reviewed journal publications and conference presentations, which serve as a foundation for further work to build from.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Open‐source photovoltaic model pipeline validation against well‐characterized system data

Abstract All freely available plane‐of‐array (POA) transposition models and photovoltaic (PV) temperature and performance models in pvlib‐python and pvpltools‐python were examined against multiyear field data from Albuquerque, New Mexico. The data include different PV systems composed of crystalline silicon modules that vary in cell type, module construction, and materials. These systems have been characterized via IEC 61853‐1 and 61853‐2 testing, and the input data for each model were sourced from these system‐specific test results, rather than considering any generic input data (e.g., manufacturer's specification [spec] sheets or generic Panneau Solaire [PAN] files). Six POA transposition models, 7 temperature models, and 12 performance models are included in this comparative analysis. These freely available models were proven effective across many different types of technologies. The POA transposition models exhibited average normalized mean bias errors (NMBEs) within ±3%. Most PV temperature models underestimated temperature exhibiting mean and median residuals ranging from −6.5°C to 2.7°C; all temperature models saw a reduction in root mean square error when using transient assumptions over steady state. The performance models demonstrated similar behavior with a first and third interquartile NMBEs within ±4.2% and an overall average NMBE within ±2.3%. Although differences among models were observed at different times of the day/year, this study shows that the availability of system‐specific input data is more important than model selection. For example, using spec sheet or generic PAN file data with a complex PV performance model does not guarantee a better accuracy than a simpler PV performance model that uses system‐specific data.

14 SOLAR ENERGY↗
DOE OSTI2019

Modeling anaerobic soil organic carbon decomposition in Arctic polygon tundra: insights into soil geochemical influences on carbon mineralization

Abstract. Rapid warming of Arctic ecosystems exposes soil organic matter (SOM) to accelerated microbial decomposition, potentially leading to increased emissions of carbon dioxide ( CO 2 ) and methane ( CH 4 ) that have a positive feedback on global warming. Current estimates of the magnitude and form of carbon emissions from Earth system models include significant uncertainties, partially due to the oversimplified representation of geochemical constraints on microbial decomposition. Here, we coupled modeling principles developed in different disciplines, including a thermodynamically based microbial growth model for methanogenesis and iron reduction, a pool-based model to represent upstream carbon transformations, and a humic ion-binding model for dynamic pH simulation to build a more versatile carbon decomposition model framework that can be applied to soils under varying redox conditions. This new model framework was parameterized and validated using synthesized anaerobic incubation data from permafrost-affected soils along a gradient of fine-scale thermal and hydrological variabilities across Arctic polygonal tundra. The model accurately simulated anaerobic CO 2 production and its temperature sensitivity using data on labile carbon pools and fermentation rates as model constraints. CH 4 production is strongly influenced by water content, pH, methanogen biomass, and presence of competing electron acceptors, resulting in high variability in its temperature sensitivity. This work provides new insights into the interactions of SOM pools, temperature increase, soil geochemical feedbacks, and resulting CO 2 and CH 4 production. The proposed anaerobic carbon decomposition framework presented here builds a mechanistic link between soil geochemistry and carbon mineralization, making it applicable over a wide range of soils under different environmental settings.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

An Observation‐Driven Approach to Improve Vegetation Phenology in a Global Land Surface Model

An empirical model calibration approach is presented that aims to approximate missing biosphere processes in a global land surface model without the need for substantial model structural changes. The strategy is implemented here using the NASA Catchment‐CN land surface model and Moderate Resolution Imaging Spectroradiometer (MODIS) observations of the fraction of absorbed photosynthetically active radiation (FPAR). Existing plant functional types (PFTs) of the Catchment‐CN model are divided into three subtypes, based on the bias between the model‐simulated and MODIS‐observed FPAR. Separate sets of vegetation parameters for each subtype are then calibrated at a small number of grid cells with homogeneous, single‐PFT land cover, using MODIS FPAR reference observations from 2003 to 2009. The effectiveness of the empirical approach at improving the realism of modeled vegetation dynamics is investigated with two global model simulations for the period 2010–2016, one using the newly calibrated parameter values and the other using the original values. Globally, the calibrated parameters reduce the root mean square error (RMSE) of the modeled FPAR with respect to MODIS by 0.029 (∼10%) on average. In some regions, substantially larger RMSE reductions are achieved. RMSE reductions are primarily driven by model bias reductions, with neutral effects on the temporal correlation skill. While the empirical approach is suitable for achieving consistent model improvements, it is shown to be sensitive to the characteristics of the model error, specifically a dominance of the bias component in the case of Catchment‐CN. Ultimately, more fundamental model structural changes may be required to achieve better improvements.

J. Kolassa↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗