Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conceptual models”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Techno-economic assessment for the production of algal fuels and value-added products: opportunities for high-protein microalgae conversion

Abstract Background Microalgae possess numerous advantages for use as a feedstock in producing renewable fuels and products, with techno-economic analysis (TEA) frequently used to highlight the economic potential and technical challenges of utilizing this biomass in a biorefinery context. However, many historical TEA studies have focused on the conversion of biomass with elevated levels of carbohydrates and lipids and lower levels of protein, incurring substantial burdens on the ability to achieve high cultivation productivity rates relative to nutrient-replete, high-protein biomass. Given a strong dependence of algal biomass production costs on cultivation productivity, further TEA assessment is needed to understand the economic potential for utilizing potentially lower-cost but lower-quality, high-protein microalgae for biorefinery conversion. Results In this work, we conduct rigorous TEA modeling to assess the economic viability of two conceptual technology pathways for processing proteinaceous algae into a suite of fuels and products. One approach, termed mild oxidative treatment and upgrading (MOTU), makes use of a series of thermo-catalytic operations to upgrade solubilized proteins and carbohydrates to hydrocarbon fuels, while another alternative focuses on the biological conversion of those substrates to oxygenated fuels in the form of mixed alcohols (MA). Both pathways rely on the production of polyurethanes from unsaturated fatty acids and valorization of unconverted solids for use as a material for synthesizing bioplastics. The assessment found similar, albeit slightly higher fuel yields and lower costs for the MA pathway, translating to a residual solids selling price of $899/ton for MA versus $1033/ton for MOTU as would be required to support a $2.50/gallon gasoline equivalent (GGE) fuel selling price. A variation of the MA pathway including subsequent upgrading of the mixed alcohols to hydrocarbon fuels (MAU) reflected a required solids selling price of $975/ton. Conclusion The slight advantages observed for the MA pathway are partially attributed to a boundary that stops at oxygenated fuels versus fungible drop-in hydrocarbon fuels through a more complex MOTU configuration, with more comparable results obtained for the MAU scenario. In either case, it was shown that an integrated algal biorefinery can be economical through optimal strategies to utilize and valorize all fractions of the biomass.

09 BIOMASS FUELS↗
DOE OSTI2024

Bench-Scale Development of a Novel Direct Air Capture Technology Using High-Capacity Structured Sorbents

The work performed under this project has resulted into development of a DAC technology utilizing a structured sorbent to capture CO 2 from ambient air with a key innovation of direct Joule heating of the sorbent for CO 2 desorption. A working SMA, fully integrated with an electrically resistive heating layer, high surface area support, and high CO 2 capacity sorbent coated onto a commercial ceramic monolithic substrate, was successfully developed and demonstrated over >200 adsorption-desorption cycles in a high-fidelity bench test unit directly using ambient air. A cordierite-based monolith was selected as a substrate owing to its high surface area, low bulk density, low heat capacity, and commercial availability. Reaction kinetics study conducted during this project led to development of a promoter for the base Na 2 CO 3 sorbent that could be incorporated into the sorbent to enhance to achieve higher CO 2 adsorption/desorption rates, greater working capacity, and reduced regeneration temperature. An accelerated aging study was conducted in a TGA to determine sorbent stability and no degradation in the sorbent performance was observed even after 250 adsorption-desorption cycles. An electrically resistive heating layer was developed with tunable electrical properties. The heating layer was coated onto the selected cordierite substrate. Aging studies performed showed the electrical properties and heating performance was stable after 500 heating and cooling cycles. The collective findings on the selected cordierite substrate, robust heating layer, promoter and sorbent selection were used to synthesize a full, 6”x6” SMA for bench-scale testing. The bench-scale DAC system was constructed to test full size SMAs using real ambient air for adsorption and joule heating for regeneration. After completing shakedown and commissioning of the 1 kg/day of CO 2 capacity DAC bench unit, an extended operation was performed to complete over 230 cycles with the full size SMA. This testing showed no observable degradation in sorbent performance. A detailed process model, TEA and LCA were developed for a conceptual 100,000 TPY CO 2 removal DAC facility. The LCA results show the net CO 2 e emissions from the DAC system are highly dependent on the electricity source. All other factors, including SMA manufacturing, materials for facility enclosure, etc., are minor cost contributors compared to the energy consumption required for CO 2 removal. With the successful development and validation of the SMA for the sustained performance for CO 2 removal from ambient air with joule heated regeneration in this project, a fully integrated 1 TPY bench-scale DAC system is currently in development with the support of DOE/FECM (DE-FE0032243). The project objective is to demonstrate the engineering design of the DAC system to produce a continuous, high purity CO 2 stream from ambient air. This project will address and validate key engineering features of the DAC system including the gas sealing mechanism and panels, enclosure and air contactor design, and automation sequence to achieve continuous CO 2 production.

42 ENGINEERING↗
DOE OSTI2025

Advances in Modeling Capabilities for Critical Mineral Separation Technologies: A PrOMMiS Overview

This is an oral presentation at the TechConnect conference on the work developed by PrOMMiS. PrOMMiS builds on and extends capabilities developed within the Department of Energy’s (DOE) Institute for the Design of Advanced Energy Systems (IDAES), Integrated Platform, and Water Treatment Technoeconomic Assessment Platform (WaterTAP), which have been successfully leveraged by other Department of Energy research areas. The open-source toolkit facilitates validation, reproducibility, and accountability, allowing for easy extension of the framework to other systems. This talk presents an overview of the PrOMMiS capabilities, including unit model library, advances in thermophysical properties models, and capital cost libraries for simulation and optimization of mineral processing technologies. The PrOMMiS applications include (1) conceptual design and superstructure optimization for screening different process configurations and identifying promising technologies; (2) dynamic modeling and optimization to enable the creation of digital twins; (3) surrogate modeling tools to leverage data when predictive thermodynamic models are not currently available; (4) technical risk reduction via uncertainty quantification and robust optimization to identify process designs that are robust to process variability and uncertainties; and (5) deployment of uncertainty quantification tools to maximize knowledge gained from experimental campaigns, while reducing the number of experiments required

critical minerals and materials↗
DOE OSTI2025

SAM Two-Phase Flow Model Development and Applications for Operational Transients in Advanced Reactors

As advanced nuclear technologies continue to develop, the need for the flexible operation and generation of these advanced reactors becomes necessary to maximize economic potential. As large-scale experiments are not always feasible, modeling and simulations of advanced reactors play a crucial role in design optimization and analysis. The SAM (System Analysis Module) code developed at Argonne National Laboratory is a state-of-the-art system-level thermal-hydraulic code aimed at simulating advanced reactor systems. Recent code developments have implemented two-phase flow modeling using the homogeneous equilibrium model, and a new steam generator component has been developed to utilize the two-phase flow implementation. In addition to verification tests, a load-following simulation was performed to model a realistic load-following transient in a proposed integrated system consisting of a conceptual advanced reactor known as the Advanced Burner Test Reactor (ABTR) and thermal energy storage (TES) tanks. The integrated system model uses two large TES tanks designed for sodium and a model helical coil steam generator to simulate the operational load-following transient. The flow rates of the feedwater and secondary loops are regulated to meet a prescribed steam generator load consistent with the electricity demand over a 24-h period. In conclusion, the results found the ABTR system was able to maintain stable reactor conditions and primary- and secondary-side characteristics over the course of the load-following transient.

Advanced Burner Test Reactor (ABTR)↗
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Typology for Characterizing Human Action in MultiSector Dynamics Models

The role of individual and collective human agency is increasingly recognized as a prominent and arguably paramount determinant in shaping the behavior, trajectory, and vulnerability of multisector systems. This human influence operates at multiple scales: from short-term (hourly to daily) to long-term (annually to centennial) timescales, and from the local to the global, pushing systems towards either desirable or undesirable outcomes. However, the effort to represent human systems in multisector models has been fragmented across philosophical, methodological, and disciplinary lines. To cohere insights across diverse modeling approaches, we present a new typology for classifying how human actors are represented in the broad suite of coupled human-natural system models that are applied in MultiSector Dynamics (MSD) research. The typology conceptualizes a “sector” as a system-of-systems that includes a diverse group of human actors, defined across individual to collective social levels, involved in governing, provisioning, and utilizing products, goods, or services towards some human end. We trace the salient features of modeled representations of human systems by organizing the typology around three key questions: 1) Who are the modeled actors in MSD systems? 2) What are their modeled actions? 3) How and for what purpose are these actors and actions operationalized in a computational model? The typology is used to critically examine existing models and chart the frontier of future human systems modeling for MSD research.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

A framework to implement human reliability analysis during early design stages of advanced reactors

Nuclear power plants require human actions throughout their lifecycle from design, construction, operation, and decommissioning. However, for advanced reactors (e.g., Generation IV), the reliance on human intervention in safety-related actions is expected to be reduced or completely replaced by automated actions. The Probabilistic Risk Assessment (PRA) Standard for Advanced Non-LWR Nuclear Power Plants requires that the impacts of all operator actions are captured and incorporated in the risk of the modeled plant. Moreover, the Modernization of Technical Requirements for Licensing Advanced Reactors requires human reliability analysis (HRA) to be included throughout all design and PRA development stages. However, due to the lack of details during the early design stages, HRA is often postponed until the design is mature enough. Conducting HRA in later design stages, though it may be adequate in capturing pre-, at-, and post-initiators comes short of informing the design itself in the iterative design lifecycle. Hence, this paper presents a framework to include HRA during the design's early stages, pre-conceptual or conceptual. The proposed framework provides a process for the removal of operator actions that do not contribute to the risk and the identification of all key operator actions that are critical to the safety of the design. The results of this framework are then used to inform the design of those safety-related operator actions to update the design further. Then, using information from the updated design, this framework can be reapplied to investigate the impact of the design update on human reliability. The PRA model of the X-energy's pre-conceptual Xe-100 high-temperature gas-cooled pebble-bed reactor (HTGR-PB) design is used to demonstrate the approach. In the pre-conceptual Xe-100 PRA model, also called Phase 0 PRA model, human actions were considered an integral part of analyzing the plant response to different initiating events. Hence, in this paper, all possible human actions in the Xe-100 PRA model are identified, analyzed, and removed to emulate a design relying only on the available automated control systems. The preliminary results of this assessment show how safe the Xe-100 design is even without crediting any human actions. The results also list necessary sequences in which operator actions are critical to the risk profile of the design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Mascon distribution techniques for asteroids and comets

The mass-concentration model is an approach that has been used to model the gravitational fields of irregularly shaped bodies such as asteroids and comets. By this approach, the body is treated as a collection of point masses. The method is conceptually simple, easy to program, valid down to the surface, and capable of modeling arbitrary density heterogeneities. How the mass concentrations are distributed as well as how mass is assigned to these concentrations is, however, nontrivial. These aspects significantly affect the accuracy and efficiency of the gravitational model. In this paper, we frame the distribution process in terms of numerical integration applied to finite volume meshes. We describe a new method using unstructured, curvilinear, finite volume meshes to significantly improve the accuracy of the mass-concentration model. We then compare the accuracy and efficiency of several variations of our distribution technique to those from literature using Asteroid Eros and Bennu as example bodies. Our results show that the mascon model can be as accurate as the analytic polyhedral model at the surface using an equivalent number of computational elements—i.e., mascon to surface facets. We report the improvement in the model’s performance can be mainly attributed to the volume mesh topology while mesh curving can provide modest case-dependent improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impact of Geothermal District Heating System on Flexibility of Microgrid in Tuttle, Oklahoma: Preprint

Flexibility is the capability of the power grid to maintain a balance between electricity generation and variable demand. This study presents preliminary results evaluating the impact of geothermal district heating systems on the flexibility of a conceptual microgrid in Tuttle, Oklahoma. Heating demand profiles were modeled using EnergyPlus for the district that includes two schools and 250 single-family houses. Then, geothermal energy production was modeled using GEOPHIRES to estimate how much heating demand in the district can be supplied by five different geothermal system scenarios. The results indicated that geothermal energy production varied depending on the resource temperature at different depths, system configurations, and flow rates. For the grid flexibility analysis, electricity consumptions in the five geothermal systems were estimated for pump operations to circulate water from the wells to radiators, while electricity consumption by air-source heat pump in the base case was estimated to supply the same heating load. Electricity consumption in the geothermal systems was significantly lower than those in base cases. The electricity saved by the geothermal system was then incorporated into the microgrid electrical load profiles where variable renewable electricity generation is significantly high. The results visually showed that geothermal district heating system can improve grid flexibility as a baseload during the winter season. The results also highlighted potential opportunities to save energy costs that will be further analyzed in future study.

Cambium↗
DOE OSTI2026

Model-Agnostic Signal Discovery with Machine Learning: Bridging the Gap Between Theory and Practice

Searches for new phenomena in complex scientific data are predominantly model-dependent, optimized for specific hypotheses, and therefore limited in their coverage of the space of possible signals. Recently, new AI-based model-agnostic search strategies, many of which have been pioneered in high-energy physics, have been proposed which provide a complementary paradigm, prioritizing broad exploration over tailored analyses. These techniques offer an opportunity to enhance the overall discovery potential of modern experiments, especially in regimes where theoretical guidance is scarce. In this document, we review the conceptual framework behind the main classes of AI-based model-agnostic strategies. We discuss the potential pitfalls of these methods, and strategies for their validation and interpretation. We aim for this document to serve as a useful reference both for practitioners and for researchers interested in learning more about these model-agnostic search strategies.

Amram, Oz [Fermilab] (ORCID:0000000237653123)↗
DOE OSTI2022

Single-Phase Dielectric Fluid Thermal Management for Power-Dense Automotive Power Electronics

This paper describes the design and performance of a dielectric fluid cooling concept for automotive power electronics. The concept combines a low-thermal-resistance package (which eliminates metalized ceramic substrates) with a high-performance convective cooling strategy (slot jets impinging on finned surfaces). Modeling was first used to design the cooling system to maximize thermal performance and minimize pumping power. Additionally, a prototype was then fabricated, and experiments were conducted to validate the model predictions using three fluids at various fluid flow rates (16.7 cm 3 /s [1 L/min] to 68.3 cm 3 /s [4.1 L/min]) and inlet temperatures (30C and 70C). The final design was compact (120-cm 3 total volume, including heat exchanger and conceptual power modules) and cooled 12 devices (e.g., silicon carbide). The validated model was then used to predict the junction-to-fluid thermal resistance and pumping power for various conditions, including 40C fluid temperature. The results predict thermal resistance values as low as 19 mm 2 K/W are possible using the dielectric fluid cooling approach. The dielectric fluid cooling system is predicted to provide thermal resistance and pumping power values that are approximately 56% and 90% lower, respectively, compared to an automotive power electronics cooling system.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Harmonising the land-use flux estimates of global models and national inventories for 2000–2020

As the focus of climate policy shifts from pledges to implementation, there is a growing need to track progress on climate change mitigation at the country level, particularly for the land-use sector. Despite new tools and models providing unprecedented monitoring opportunities, striking differences remain in estimations of anthropogenic land-use CO 2 fluxes between, on the one hand, the national greenhouse gas inventories (NGHGIs) used to assess compliance with national climate targets under the Paris Agreement and, on the other hand, the Global Carbon Budget and Intergovernmental Panel on Climate Change (IPCC) assessment reports, both based on global bookkeeping models (BMs). Recent studies have shown that these differences are mainly due to inconsistent definitions of anthropogenic CO 2 fluxes in managed forests. Countries assume larger areas of forest to be managed than BMs do, due to a broader definition of managed land in NGHGIs. Additionally, the fraction of the land sink caused by indirect effects of human-induced environmental change (e.g. fertilisation effect on vegetation growth due to increased atmospheric CO 2 concentration) on managed lands is treated as non-anthropogenic by BMs but as anthropogenic in most NGHGIs. We implement an approach that adds the CO 2 sink caused by environmental change in countries' managed forests (estimated by 16 dynamic global vegetation models, DGVMs) to the land-use fluxes from three BMs. This sum is conceptually more comparable to NGHGIs and is thus expected to be quantitatively more similar. Our analysis uses updated and more comprehensive data from NGHGIs than previous studies and provides model results at a greater level of disaggregation in terms of regions, countries and land categories (i.e. forest land, deforestation, organic soils, other land uses). Our results confirm a large difference (6.7 GtCO 2 yr —1 ) in global land-use CO 2 fluxes between the ensemble mean of the BMs, which estimate a source of 4.8 GtCO 2 yr —1 for the period 2000–2020, and NGHGIs, which estimate a sink of —1.9 GtCO 2 yr —1 in the same period. Most of the gap is found on forest land (3.5 GtCO 2 yr —1 ), with differences also for deforestation (2.4 GtCO 2 yr —1 ), for fluxes from other land uses (1.0 GtCO 2 yr —1 ) and to a lesser extent for fluxes from organic soils (0.2 GtCO 2 yr —1 ). By adding the DGVM ensemble mean sink arising from environmental change in managed forests (—6.4 GtCO 2 yr —1 ) to BM estimates, the gap between BMs and NGHGIs becomes substantially smaller both globally (residual gap: 0.3 GtCO 2 yr —1 ) and in most regions and countries. However, some discrepancies remain and deserve further investigation. For example, the BMs generally provide higher emissions from deforestation than NGHGIs and, when adjusted with the sink in managed forests estimated by DGVMs, yield a sink that is often greater than NGHGIs. In summary, this study provides a blueprint for harmonising the estimations of anthropogenic land-use fluxes, allowing for detailed comparisons between global models and national inventories at global, regional and country levels. This is crucial to increase confidence in land-use emissions estimates, support investments in land-based mitigation strategies and assess the countries' collective progress under the Global Stocktake of the Paris Agreement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Relationships between fault characteristics and seismic responses in a large lab-scale tri-axial injection test conducted on a faulted Castlegate Sandstone

Here to investigate mechanisms causing microseismicity (M w < 2) at a CO 2 injection site, a large-scale triaxial block experiment was carried out on a faulted (saw-cut) cubic-meter of Castlegate Sandstone. The experiment consisted of injection tests at varying differential stresses, while monitoring and recording pore pressure and acoustic emissions (AEs). During the experiment, ~33,000 AEs and ~14 mm of horizontal displacement/slip, like a strike-slip fault movement, occurred. To understand the AE responses and ascertain fault characteristics near the located AEs, we modeled the topography of the fault surface, fault aperture, and fault-gouge thickness using pre- and post-experiment laser scans of the fault surface on each half of the block. Additionally, we characterized surface roughness parallel and perpendicular to slip. Models show crushing and flattening of the fault surfaces can be linked to the spatiotemporal distribution of AEs within 50 mm of the fault surface. Approximately 65% of AEs were in areas with small aperture (≤300 μm); thicker fault gouge was observed in adjacent areas with wider aperture and shows a two-fold reduction in grain size relative to unaltered Castlegate Sandstone. This work provides a conceptual understanding on fault surface evolution, which can be applied towards modeling of seismic slip.

58 GEOSCIENCES↗
DOE OSTI2020

Non-Hermitian approach for quantum plasmonics

Herein, we examine the limits of applicability of a simple non-Hermitian model for exciton/plasmon interactions in the presence of dissipation and dephasing. The model can be used as an alternative to the more complete Lindblad density matrix approach and is computationally and conceptually simpler. We find that optical spectra in the linear regime can be adequately described by this approach. The model can fail, however, under continuous optical driving in some circumstances. In the case of two quantum dots or excitons interacting with a plasmon, the model can also describe coherences and entanglement qualitatively when both dissipation and dephasing are present and quantitatively in the limit with no dephasing. The approach, within a single excitation manifold, is also applied to assess the role of disorder for 50 quantum dots interacting with a plasmon, where we find that, on average, large enough disorder can help stabilize the ensemble average of the open quantum system toward a dark quasi-steady-state much faster than without disorder. Furthermore, while such single excitation manifold calculations in this size limit can readily be done with either the non-Hermitian or Lindblad forms, as one goes to larger Hilbert space sizes, the computational and storage advantages of the non-Hermitian approach can become more useful. Published under license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Teaching a neural network to attach and detach electrons from molecules

Abstract Interatomic potentials derived with Machine Learning algorithms such as Deep-Neural Networks (DNNs), achieve the accuracy of high-fidelity quantum mechanical (QM) methods in areas traditionally dominated by empirical force fields and allow performing massive simulations. Most DNN potentials were parametrized for neutral molecules or closed-shell ions due to architectural limitations. In this work, we propose an improved machine learning framework for simulating open-shell anions and cations. We introduce the AIMNet-NSE (Neural Spin Equilibration) architecture, which can predict molecular energies for an arbitrary combination of molecular charge and spin multiplicity with errors of about 2–3 kcal/mol and spin-charges with error errors ~0.01e for small and medium-sized organic molecules, compared to the reference QM simulations. The AIMNet-NSE model allows to fully bypass QM calculations and derive the ionization potential, electron affinity, and conceptual Density Functional Theory quantities like electronegativity, hardness, and condensed Fukui functions. We show that these descriptors, along with learned atomic representations, could be used to model chemical reactivity through an example of regioselectivity in electrophilic aromatic substitution reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗
DOE OSTI2022

Toward reducing uncertainty quantification costs in DEM models of particulate flow: Testing simple, sensitivity-based, forward uncertainty propagation techniques

The performance of two conceptually-simple uncertainty quantification techniques are tested against the rigorous nested-loop sampling technique of Roy and Oberkampf (Comput Methods Appl Mech Eng, 200: 2131–2144, 2011) (herein called full-sampling) using two very small-scale DEM-based models of particulate flow (one gas-solid flow and one granular flow). The first simplified forward uncertainty propagation technique, reduced-sampling, uses a sensitivity analysis to eliminate uncertain inputs that have little impact on the model output prior to nested-loop sampling. The second technique, boundary-sampling, uses a sensitivity analysis to inform the selection of two bounding cases for each key model output. In conclusion, the uncertainties in the model outputs obtained via the reduced- and boundary-sampling methods agree well with those from full-sampling for both the gas-solid and granular flow models while yielding computational savings of 65–75% (reduced sampling) and 94–97% (boundary sampling).

42 ENGINEERING↗
DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗