Techno-economic uncertainty analysis of wet waste-to-biocrude via hydrothermal liquefaction
Not Available
Li, Shuyun↗
SEARCH · Search NASA
Results for “Coupled process modeling”
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.
Exploring the energy landscape of RBMs: reciprocal space insights into bosons, hierarchical learning and symmetry breaking
Deep generative models have become ubiquitous due to their ability to learn and sample from complex distributions. Despite the proliferation of various frameworks, the relationships among these models remain largely unexplored, a gap that hinders the development of a unified theory of AI learning. In this work, we address two central challenges: clarifying the connections between different deep generative models and deepening our understanding of their learning mechanisms. We focus on Restricted Boltzmann Machines (RBMs), a class of generative models known for their universal approximation capabilities for discrete distributions. By introducing a reciprocal space formulation for RBMs, we reveal a connection between these models, diffusion processes, and systems of coupled bosons. Our analysis shows that at initialization, the RBM operates at a saddle point, where the local curvature is determined by the singular values of the weight matrix, whose distribution follows the Marc̆enko-Pastur law and exhibits rotational symmetry. During training, this rotational symmetry is broken due to hierarchical learning, where different degrees of freedom progressively capture features at multiple levels of abstraction. This leads to a symmetry breaking in the energy landscape, reminiscent of Landau’s theory. This symmetry breaking in the energy landscape is characterized by the singular values and the weight matrix eigenvector matrix. We derive the corresponding free energy in a mean-field approximation. We show that in the limit of infinite size RBM, the reciprocal variables are Gaussian distributed. Our findings indicate that in this regime, there will be some modes for which the diffusion process will not converge to the Boltzmann distribution. To illustrate our results, we trained replicas of RBMs with different hidden layer sizes using the MNIST dataset. Our findings not only bridge the gap between disparate generative frameworks but also shed light on the fundamental processes underpinning learning in deep generative models.
97 MATHEMATICS AND COMPUTING↗
PFLOTRAN-SIP: A PFLOTRAN Module for Simulating Spectral-Induced Polarization of Electrical Impedance Data
Spectral induced polarization (SIP) is a non-intrusive geophysical method that collects chargeability information (the ability of a material to retain charge) in the time domain or its phase shift in the frequency domain. Although SIP is a temporal method, it cannot measure the dynamics of flow and solute/species transport in the subsurface over long times (i.e., 10–100 s of years). Data collected with the SIP technique need to be coupled with fluid flow and reactive-transport models in order to capture long-term dynamics. To address this challenge, PFLOTRAN-SIP was built to couple SIP data to fluid flow and solute transport processes. Specifically, this framework couples the subsurface flow and transport simulator PFLOTRAN and geoelectrical simulator E4D without sacrificing computational performance. PFLOTRAN solves the coupled flow and solute-transport process models in order to estimate solute concentrations, which were used in Archie’s model to compute bulk electrical conductivities at near-zero frequency. These bulk electrical conductivities were modified while using the Cole–Cole model to account for frequency dependence. Using the estimated frequency-dependent bulk conductivities, E4D simulated the real and complex electrical potential signals for selected frequencies for SIP. These frequency-dependent bulk conductivities contain information that is relevant to geochemical changes in the system. This study demonstrated that the PFLOTRAN-SIP framework is able to detect the presence of a tracer in the subsurface. SIP offers a significant benefit over ERT in the form of greater information content. It provided multiple datasets at different frequencies that better constrained the tracer distribution in the subsurface. Consequently, this framework allows for practitioners of environmental hydrogeophysics and biogeophysics to monitor the subsurface with improved resolution.
54 ENVIRONMENTAL SCIENCES↗
Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System
The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.
54 ENVIRONMENTAL SCIENCES↗
Sensitivity evaluation of the Kudryavtsev permafrost model
Modeling is an important way to assess current and future permafrost spatial distribution and dynamics, especially in data poor areas like the Arctic region. Here, we evaluate a physics-based analytical model, Kudryavtsev's active layer model, which is widely used because it has relatively few data requirements. This model was recently incorporated into a component modeling toolbox, allowing for coupled modeling of permafrost and geomorphic processes over geological timescales. However, systematic quantitative assessment of the influence of its controlling parameters on permafrost temperature and active layer thickness predictions has not been undertaken before. In this work, we investigate the sensitivity of the Kudryavtsev's active layer model by Monte Carlo simulations to generate probability distributions for input parameters and compare predictions with a comprehensive benchmark dataset of in-situ permafrost observations over entire Alaska. Predicted permafrost surface temperature is highly dependent on mean annual air temperature (r = 0.78 on average), annual temperature amplitude (–0.41), and winter-averaged snow thickness (0.30). Uncertainty of predicted permafrost temperature is relatively small (RMSE = 1 °C), when air temperature and snow depth are well constrained. Similarly, RMSE between simulated and observed ALT at stations is ~0.08 m. However, under given air temperature and snow conditions, soil water content bias can significantly affect modeled active layer thickness (RMSE = 0.1 m or 40% of the observed active layer thickness). If soil water content has a large bias, improvements in other parameters may not significantly improve the active layer predictions of the Kudryavtsev's model.
54 ENVIRONMENTAL SCIENCES↗
Multiphysics analysis system for heat pipe cooled micro-reactors employing PRAGMA-OpenFOAM-ANLHTP
A multiphysics analysis system for neutronics/thermo-mechanical/heat-pipe analysis of heat pipe cooled micro-reactors was developed using the PRAGMA code as the neutronics engine. PRAGMA, which used to be a GPU-based continuous-energy MC code for power reactor applications, now has an extended geometry package to handle geometries with unstructured meshes generated by the ANSYS Design-Modeler and Meshing. The NVIDIA ray tracing engine OptiX was exploited for efficient neutron transport on unstructured geometry. On the multiphysics side, the open-source CFD tool OpenFOAM and one-dimensional heat pipe analysis code ANLHTP were adopted. The manager-worker system based on the MPI dynamic process management (DPM) model enables efficient coupling of codes employing different parallelization schemes. With all the features, the multiphysics analysis of the one-sixth symmetrical MegaPower 2D core was performed and it demonstrated the benefit of the tight integration of the three-way coupling system and one-to-one geometry coupling strategy. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Understanding Microstructure Variability in Vapor-Deposited Energetic Materials by Using Phase-Field Methods
Critical components, such as detonators, in Sandia's stockpile contain heterogeneous materials whose performance and reliability depend on accurate, predictive models of coupled, complex phenomena to predict their synthesis, processing, and operation. Ongoing research in energetic materials has shown that microstructural properties, such as density, pore-size, morphology, and specific surface area are correlated to their initiation threshold and detonation behavior. However, experiments to study these specific characteristics of energetic materials are challenging and time consuming. Therefore, in this work, we turn to mesoscale modeling methods that may be capable of reproducing some observed phenomena to refine and predict outcomes beforehand. Even so, we have no physics-based modeling capability to predict how the microstructure of an energetic material will evolve over near- and long-term time scales. Thus, the goal of this work is to (i) identify any knowledge gaps in how the underlying microstructure forms and evolves during the synthesis process, and (ii) develop and test a mesoscale phase-field model for vapor deposition to capture critical mechanisms of microstructure formation, evolution, and variability in vapor-deposited energetic materials, such as processing conditions, material properties, and substrate interactions. Based on this work, the phase-field method is shown to be a valuable tool for developing the necessary models containing coupled, complex phenomena to investigate and understand the synthesis and processing of energetic materials.
36 MATERIALS SCIENCE↗
Natural Fracture Pattern Development Using 3D Coupled Numerical Model: Mechanics and Diagenesis
In the earth's crust, while tectonic processes can mechanically create natural fractures, cementation can simultaneously grow on the surfaces of the fractures. This crystallization can affect the growth of the fractures. To more accurately address the modeling of natural fracture development, both mechanical deformation and chemical diagenesis should be carried out in an integrated manner. Here, we achieve this with a three-dimensional (3D) simulator that couples mechanical and chemical processes for natural fracture pattern development in layer confined situations. Tectonic loading over geologic time is mimicked by displacement driven boundary conditions in a 3D displacement discontinuity method (DDM) scheme. Quartz diagenesis competes with mechanical fracture opening, resulting in three different cementation conditions - open, bridging or fully infilled – which define fracture element stiffness that is used as feedback within the 3D-DDM code during propagation. An important milestone of this work is that this coupled structural diagenetic model can reproduce power-law aperture population statistics which cannot be generated using mechanics alone. More accurately estimated mechanical fracture apertures and their cementation condition is a critical step toward better permeability estimation and fluid path analysis for subsurface reservoirs.
58 GEOSCIENCES↗
Computational Tools for Additive Manufacture of Tailored Microstructure and Properties
Additive manufacturing has the potential to revolutionize industrial hardware and unlock efficiency gains through the fabrication of geometries and architectures not possible by conventional processing. Currently most additive builds use a single set of process parameters (e.g. laser power and scan speed) which results in a part with a homogenous microstructure that provides a singular performance level. To move beyond this state, Raytheon Technologies Research Center worked to create a set of computational tools to track material evolution through each step of the additive process. Computational fluid dynamics and phase field models for microstructure evolution as a function of processing parameters, and crystal plasticity models fully coupling microstructure and mechanical properties for performance predictions were leveraged to establish a connection between additive parameters and the final microstructure. This framework was utilized to tailor spatially-varying mechanical properties in a part by appropriately controlling the microstructure evolution during the additive process. Specifically, a turbine blade was 3D printed from nickel superalloy IN718 using laser powder bed fusion with coarse grains in the airfoil section which experiences the highest temperatures and is creep limited while finer grains were printed in the root of the blade which experience higher stresses but at lower temperatures and is therefore fatigue limited. The benefit of being able to intentionally insert coarse grains in the high temperature region of the blade was showcased with a microstructure sensitive creep model that indicates longer creep life for coarser grains.
20 FOSSIL-FUELED POWER PLANTS↗
Omics-to-Reactive-Transport (ORT): A workflow linking genome-scale metabolic models with reactive transport codes
Motivation: Nutrient and contaminant behavior in the subsurface are governed by multiple coupled hydrobiogeochemical processes which occur across different temporal and spatial scales. Accurate description of macroscopic system behavior requires accounting for the effects of microscopic and especially microbial processes. Microbial processes mediate precipitation and dissolution and change aqueous geochemistry, all of which impacts macroscopic system behavior. As `omics data describing microbial processes is increasingly affordable and available, novel methods for using this data quickly and effectively for improved ecosystem models are needed. Results: We propose a workflow (`Omics to Reactive Transport – ORT) for utilizing metagenomic and environmental data to describe the effect of microbiological processes in macroscopic reactive transport models. This workflow utilizes and couples two open-source software packages: KBase (a software platform for systems biology) and PFLOTRAN (a reactive transport modeling code). We describe the architecture of ORT and demonstrate an implementation using metagenomic and geochemical data from a river system. Our demonstration uses microbiological drivers of nitrification and denitrification to predict nitrogen cycling patterns which agree with those provided with generalized stoichiometries. While our example uses data from a single measurement, our workflow can be applied to spatiotemporal metagenomic datasets to allow for iterative coupling between KBASE and PFLOTRAN. Live, interactive models, which incorporate the results from this narrative into a PFLOTRAN simulation, are available (without login) at https://pflotranmodeling.paf.subsurfaceinsights.com/pflotran-simple-model/.
Rubinstein, Rebecca L↗
Multilevel Analysis, Design, and Modeling of Coupling Advanced Nuclear Reactors and Thermal Energy Storage in an Integrated Energy System
This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (A LWR), a high-temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state and transient-state operation with minimum disruption to the nuclear power cycle. This first part of the study presents steady-state models developed using Aspen HYSYS®, with the thermal energy bypass for NPP-TES coupling being varied at up to 50%. The various components were sized using the Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The TES-coupled A-NPP steady-state models and cost functions then provided a baseline for dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools. A stochastic optimization of the various energy storage systems coupled to the A-NPPs was then performed using the Risk Analysis Virtual Environment (RAVEN) and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON). The signal processing and synthetic history capabilities of RAVEN were used to account for the unpredictable behavior of electricity markets. An autoregressive moving average (ARMA) model was used to analyze price signals from the Pennsylvania-New Jersey-Maryland (PJM) market and were applied to the HERON analysis in order to optimize a system with the best economics. Transient modeling evaluation was then performed using Modelica models within the HYBRID repository, which was developed at INL for the Department of Energy Integrated Energy Systems program for the characterization of dynamic integrated system behavior and feedback. This includes evaluation of the TES-coupled A-LWR systems’ impact on physical and thermal system response during imposed system demands. Additional TES-coupled reactor types, coupling approaches, markets, and TES technologies will be evaluated in future work.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DECOVALEX-2023: Task F2-Salt Final Report
The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2-(salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
ORT: a workflow linking genome-scale metabolic models with reactive transport codes
Abstract Motivation Nutrient and contaminant behavior in the subsurface are governed by multiple coupled hydrobiogeochemical processes which occur across different temporal and spatial scales. Accurate description of macroscopic system behavior requires accounting for the effects of microscopic and especially microbial processes. Microbial processes mediate precipitation and dissolution and change aqueous geochemistry, all of which impacts macroscopic system behavior. As ‘omics data describing microbial processes is increasingly affordable and available, novel methods for using this data quickly and effectively for improved ecosystem models are needed. Results We propose a workflow (‘Omics to Reactive Transport—ORT) for utilizing metagenomic and environmental data to describe the effect of microbiological processes in macroscopic reactive transport models. This workflow utilizes and couples two open-source software packages: KBase (a software platform for systems biology) and PFLOTRAN (a reactive transport modeling code). We describe the architecture of ORT and demonstrate an implementation using metagenomic and geochemical data from a river system. Our demonstration uses microbiological drivers of nitrification and denitrification to predict nitrogen cycling patterns which agree with those provided with generalized stoichiometries. While our example uses data from a single measurement, our workflow can be applied to spatiotemporal metagenomic datasets to allow for iterative coupling between KBase and PFLOTRAN. Availability and implementation Interactive models available at https://pflotranmodeling.paf.subsurfaceinsights.com/pflotran-simple-model/. Microbiological data available at NCBI via BioProject ID PRJNA576070. ORT Python code available at https://github.com/subsurfaceinsights/ort-kbase-to-pflotran. KBase narrative available at https://narrative.kbase.us/narrative/71260 or static narrative (no login required) at https://kbase.us/n/71260/258. Supplementary information Supplementary data are available at Bioinformatics online.
54 ENVIRONMENTAL SCIENCES↗
Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H 2 M model
Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories. This is particularly the case for bentonite, an often-preferred barrier material. Gas flow in preferential pathways is characterized by localization and spontaneous behavior, which is challenging to simulate in numerical models due to strong hydro-mechanical coupling. To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project, this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations. The proposed model utilized hydro-mechanical spatial distributions, namely Young's modulus and gas entry pressure, and elastoplasticity combined with a linear swelling model. A conceptual strain-dependent permeability approach simulated dilatancy-controlled gas flow based on hydro-mechanical coupling. To test the effectiveness of the presented approach, a gas injection test in a compacted, saturated bentonite sample was simulated using the open-source code OpenGeoSys 5.8 and compared with experimental observations. The presented methodology is capable of simulating localized gas flow in preferential pathways. The spatial distributions of Young's modulus and gas entry pressure affect the swelling pressure, relative permeability and, in combination with the strain-dependent permeability model, also the intrinsic permeability.
42 ENGINEERING↗
Multiple Degradation Mechanisms in Reinforced Concrete Structures, Modeling and Risk Analysis
The overarching goal of this project is to complement ongoing Department of Energy (DOE) Light Water Reactor Sustainability (LWRS)-funded Grizzly concrete modeling development efforts by providing improved multi-physics models for incorporation into the Grizzly and BlackBear codes. Objectives of this work identified at the outset of this project include: 1. Coupling between mechanical damage and transport processes. The mechanical damage will be characterized by nonlinear mechanical constitutive models appropriate for concrete used in nuclear power plants (NPPs). 2. Improved representation of the coupling among the transport process, building on the current moisture and thermal transport models in Grizzly. The parameters for coupling among heat and mass transport in the multi-physics model will be experimentally determined. 3. Coupling among different scale levels of concrete constituents, fine and coarse aggregates, cement paste, hydration products, and pore structure. 4. Probabilistic analysis of the random nature of the heterogeneous concrete structure and the environmental factors (ambient temperature and humidity) and their multiple effects on concrete deterioration. 5. Benchmark verification, validation and uncertainty quantification of the thermo-hygro- chemo-mechanical (THCM) concrete formulation. The end goal of this work was to provide a comprehensive and robust simulation capability that can be applied at the engineering scale for analysis of realistic deterioration scenarios in NPP concrete structures.
36 MATERIALS SCIENCE↗
Two Mechanisms Limiting the Emitted Electron Current from a Cathode to an Anode
It is known that the current of emitted electrons flowing through a plasma can saturate upon formation of a potential well adjacent to the cathode (the “space charge effect”). Here, we demonstrate another saturation mechanism that will often set a more restrictive limit on the global current. When “backflow saturation” occurs, the cathode sheath weakens to allow emitted electrons that already entered the plasma to backflow to the cathode. This effect could not be captured by studies modelling the cathode sheath by itself because its origin is coupled to processes in the interior plasma and anode sheath. By modeling a full plasma diode, we show that depending on conditions the global current can be limited in four ways; by backflow alone, by space charge alone, by both mechanisms in a stable cooperative form, or by both in a competing oscillatory form. Published by the American Physical Society 2025
Laboratory plasma↗
An Analytical Model for the Tension-Shear Coupling of Woven Fabrics with Different Weave Patterns under Large Shear Deformation
It is essential to accurately describe the large shear behavior of woven fabrics in the composite preforming process. An analytical model is proposed to describe the shear behavior of fabrics with different weave patterns, in which tension-shear coupling is considered. The coupling is involved in two parts, the friction between overlapped yarns and the in-plane transverse compression between two parallel yarns. By introducing the concept of inflection points of a yarn, the model is applicable for fabrics with different weave patterns. The analytical model is validated by biaxial tension-shear experiments. A parametric study is conducted to investigate the effects of external load, yarn geometry, and weave structure on the large shear behavior of fabrics. The developed model can reveal the physical mechanism of tension-shear coupling of woven fabrics. Moreover, the model has a high computational efficiency due to its explicit expressions, thus benefiting the material design process.
36 MATERIALS SCIENCE↗
The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles
Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.
36 MATERIALS SCIENCE↗