Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Numerical study of multiscale compaction-initiated detonation

A multiscale model of heterogeneous condensed-phase explosives is examined computationally to determine the course of transient events following the application of a piston-driven stimulus. The model is a modified version of that introduced by Gonthier (Combust Sci Technol 175(9):1679–1709, 2003. https://doi.org/10.1080/00102200302373) in which the explosive is treated as a porous, compacting medium at the macro-scale and a collection of closely packed spherical grains capable of undergoing reaction and diffusive heat transfer at the meso-scale. A separate continuum description is ascribed to each scale, and the two scales are coupled together in an energetically consistent manner. Following piston-induced compaction, localized energy deposition at the sites of intergranular contact creates hot spots where reaction begins preferentially. Reaction progress at the macro-scale is determined by the spatial average of that at the grain scale. A parametric study shows that combustion at the macro-scale produces an unsteady detonation with a cyclical character, in which the lead shock loses strength and is overtaken by a stronger secondary shock generated in the partially reacted material behind it. Finally, the secondary shock in turn becomes the new lead shock and the process repeats itself.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Initial use of Nek5000/Cardinal to improve closure models in Pronghorn

Heat transfer coefficient closure models for pebble bed reactors are built using a data-driven approach by leveraging the capabilities of an Evolutionary Algorithm entitled Particle Swarm Optimization (PSO). In the present work, the Computational Fluid Dynamics code nekRS was used in order to collect the high-fidelity flow data for a core with 1,568 pebbles. To characterize the heat transfer, multiple concentric regions were considered to extract the physical quantities of interest, e.g./ the Reynolds number. The PSO algorithm is employed as part of an inverse problem targeting determine what are the coefficients for a Nusselt number correlation to match the collected data. Such correlation should follow any given format that is defined a priori. Finally, two correlations are proposed, one with an implicit dependence on the pebbles’ wall temperatures and another expressed as a fully explicit correlation depending on the flow conditions and the position within the core. Anyway, given the generic nature of the proposed approach, correlations following different formats could be tested. Preliminary results for the high-fidelity simulation of a fast MSR core are presented. The target Reynolds number is currently 20K, with the expectation that this will increase, pending the availability of further computational resources. These simulations will be used to inform lower fidelity models, including a coarse CFD turbulence model in Pronghorn. Additionally, they will serve as a reference for the RANS models in Nek5000/NekRS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Identification of mechanisms driving heterogeneous void growth in ductile aluminum

Void growth plays a central role in ductile fracture, yet the specific mechanisms that control this remain obscure. Classical models, such as those proposed by Rice and Tracey in 1969, are able to capture average rates of void growth, but cannot capture the heterogeneity of individual void growth. Building on recent work, the present study employs laboratory-based diffraction contrast tomography and in-situ x-ray computed tomography to investigate the effect of grain structure and other microstructural factors on void growth in an Al-2219 alloy. Crystal plasticity finite element (CP-FE) modeling is used alongside experimental data to evaluate the contributions of local mechanical states, grain orientation, grain size, and neighboring microstructural features. No strong linear relationships are found with any of the considered descriptors and void growth rate. Potential complex nonlinear relationships are explored with the use of a random forest regression model, which identifies initial void volume, void aspect ratio, local normal stress state, local shear stress state, and local equivalent plastic strain (EQPS) as features that most improve void growth rate predictions. The combination of these analyses suggests that these features should be prioritized to improve models of void growth.

Diffraction contrast tomography (DCT)↗
DOE OSTI2025

Standard Library Grid-Forming Hybrid Control Inverter-Based Resource Model Specification (REGFM_C1)

This document describes a standard library grid-forming (GFM) hybrid control inverter-based resource (IBR) model. The GFM hybrid control approach implements both a typical GFM control and a typical grid-following (GFL) control inside one single inverter simultaneously, so that it can take advantage of both methods without comprising the benefits of a typical GFM. The initial version of model specification was jointly developed by Pacific Northwest National Laboratory (PNNL), Tesla Energy, and EPRI, and it was revised multiple times later to incorporate suggestions from WECC MVS members. Tesla Energy provided main control blocks to support the development of this model specification. This standard library model is developed to help the utility industry better understand the GFM technology. The model could be used to represent equipment for long-term planning studies where vendor-specific models are not available. As equipment matures and improves, standard library models will be updated to capture the new functionalities of GFMs. It is not intended that these models will always remain representative of all future GFM technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Evaluation and Validation Approaches for Simulation of Social Behavior: Challenges and Opportunities

While models of social systems have grown in scope to represent increasingly complex phenomena, efforts to validate these models against observations suffer from a number of challenges that include the lack of shared approaches, the use of single-dimensional evaluation criteria, and the difficulty in evaluating causal mechanisms. Moreover, current evaluation approaches suffer from sensitivity of the existing measurements to initial conditions and model assumptions, and failure to account for uncertainty. These effects lead to the lack of generalizability to new contexts, inability to replicate the validation approach for new models, and difficulties in interpretation of validation results. These challenges can be addressed by building a set of diverse and robust evaluation metrics and by developing best practices for comparative, reproducible validation of varied social simulation approaches against real-world data. Herein, we outline several best practices for data-drive model validation, such as defining evaluation strategies consistent with the intended use case and defining broad comparison criteria to highlight both model strengths and weaknesses. Leveraging common validation practices will enable robust analysis of simulation performance across diverse social scenarios and modeling approaches.

Saldanha, Emily G.↗
DOE OSTI2024

Machine Learning Accelerated First-Principles Study of the Hydrodeoxygenation of Propanoic Acid

The complex reaction network of catalytic biomass conversions often involves hundreds of surface intermediates and thousands of reaction steps, greatly hindering the rational design of metal catalysts for these conversions. Here, we present a framework of machine learning (ML)-accelerated first-principles studies for the hydrodeoxygenation (HDO) of propanoic acid over transition metal surfaces. The microkinetic model (MKM) is initially parametrized by ML-predicted energies and iteratively improved by identifying the rate-determining species and steps (RDS), computing their energies by density functional theory (DFT), and reparameterizing the MKM until all the RDS are computed by DFT. The Gaussian process (GP) model performs significantly better than the linear ridge regression model for predicting both the adsorption free energies and transition state free energies. Parameterized with energies from the GP model, only 5–20% of the full reaction network has to be computed by DFT for the MKM to possess DFT-level accuracy for the TOF and dominant reaction pathway. While the linear ridge regression model performs worse than the GP model, its performance is greatly improved when only transition states are predicted by the regression model and adsorption energies are computed by DFT. Overall, we find that a high accuracy in adsorption free energies is more important for a reliable MKM than a high accuracy in TS free energies. Lastly, based on the GP model with GOH and GCHCHCO as catalyst descriptors, we build two-dimensional volcano plots in activity and selectivity that can help design promising alloy catalysts for HDO reactions of organic acids.

adsorption↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2024

A Pseudoreversible Normalizing Flow for Stochastic Dynamical Systems with Various Initial Distributions

Here, we present a pseudoreversible normalizing flow method for efficiently generating samples of the state of a stochastic differential equation (SDE) with various initial distributions. The primary objective is to construct an accurate and efficient sampler that can be used as a surrogate model for computationally expensive numerical integration of SDEs, such as those employed in particle simulation. After training, the normalizing flow model can directly generate samples of the SDE’s final state without simulating trajectories. The existing normalizing flow model for SDEs depends on the initial distribution, meaning the model needs to be retrained when the initial distribution changes. The main novelty of our normalizing flow model is that it can learn the conditional distribution of the state, i.e., the distribution of the final state conditional on any initial state, such that the model only needs to be trained once and the trained model can be used to handle various initial distributions. This feature can provide a significant computational saving in studies of how the final state varies with the initial distribution. Additionally, we propose to use a pseudoreversible network architecture to define the normalizing flow model, which has sufficient expressive power and training efficiency for a variety of SDEs in science and engineering, e.g., in particle physics. We provide a rigorous convergence analysis of the pseudoreversible normalizing flow model to the target probability density function in the Kullback–Leibler divergence metric. Numerical experiments are provided to demonstrate the effectiveness of the proposed normalizing flow model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Giant Cloud Condensation Nuclei Facilitate Drizzle Formation in Stratocumulus—Insights From a Combined Observation‐Modeling Framework

The mechanism for initiating drizzle drop remains a gap in the current understanding of warm rain formation. One prevalent hypothesis suggests that the presence of Giant Cloud Condensation Nuclei (GCCN) generates drizzle‐sized drops necessary to trigger the Collision‐Coalescence (C‐C) process. Here, in this study, this hypothesis is investigated using a novel framework that integrates in situ observations, remote sensing measurements, and idealized models. Results show that GCCN can efficiently generate drizzle drops through condensation, producing a broad Droplet Size Distribution (DSD) comparable to in situ observations. The large drizzle drop and broad DSD strongly facilitate C‐C, further accelerating drizzle initiation. To compare with observation, the model‐generated DSDs are used to generate radar Doppler spectra where radar reflectivity and Doppler skewness is estimated. The simulated radar quantities correspond well with radar observations, providing critical evidence for the GCCN‐induced drizzle initiation mechanism.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering

We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.

Engineering↗
DOE OSTI2021

Cavitation model of the inflationary stage of Big Bang

In this study, we propose a model for the initial stage of the development of the universe analogous to cavitation in a liquid in a negative pressure field. It is assumed that at the stage of inflation, multiple breaks of the metric occur with the formation of areas of physical vacuum in which the generation of matter occurs. The proposed model explains the large-scale isotropy of the universe without ultrafast inflationary expansion and the emergence of a large-scale cellular (cluster) structure, as a result of the development of cavitation ruptures of a false vacuum. It is shown that the cavitation model can be considered on par with (or as an alternative to) the generally accepted inflationary multiverse model of the Big Bang.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Potential Signature of Population III Pair-instability Supernova Ejecta in the BLR Gas of the Most Distant Quasar at z = 7.54*

The search for Population III stars has fascinated and eluded astrophysicists for decades. One promising place for capturing evidence of their presence must be high-redshift objects; signatures should be recorded in their characteristic chemical abundances. We deduce the Fe and Mg abundances of the broadline region (BLR) from the intensities of ultraviolet Mg ii and Fe ii emission lines in the near-infrared spectrum of UKIDSS Large Area Survey (ULAS) J1342+0928 at z = 7.54, by advancing our novel flux-to-abundance conversion method developed for quasars up to z ~ 3. We find that the BLR of this quasar is extremely enriched, by a factor of 20 relative to the solar Fe abundance, together with a very low Mg/Fe abundance ratio: [Fe/H] = +1.36 ± 0.19 and [Mg/Fe] =-1.11 ± 0.12, only 700 million years after the Big Bang. We conclude that such an unusual abundance feature cannot be explained by the standard view of chemical evolution that considers only the contributions from canonical supernovae. While there remains uncertainty in the high-mass end of the Population III initial mass function, here we propose that the larger amount of iron in ULAS J1342+0928 was supplied by a pair-instability supernova (PISN) caused by the explosion of a massive Population III star in the high-mass end of the possible range of 150–300 M ⊙ . Chemical evolution models based on initial PISN enrichment well explain the trend in [Mg/Fe]-z all the way from z < 3 to z = 7.54. We predict that stars with very low [Mg/Fe] at all metallicities are hidden in the galaxy, and they will be efficiently discovered by ongoing new-generation photometric surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Actinide Production in the Neutron-rich Ejecta of a Neutron Star Merger

The rapid neutron-capture ("r-") process is responsible for synthesizing many of the heavy elements observed in both the solar system and Galactic metal-poor halo stars. Simulations of r-process nucleosynthesis can reproduce abundances derived from observations with varying success, but so far they fail to account for the observed overenhancement of actinides, present in about 30% of r-process-enhanced stars. In this work, we investigate actinide production in the dynamical ejecta of a neutron star merger (NSM) and explore whether varying levels of neutron-richness can reproduce the actinide boost. We also investigate the sensitivity of actinide production on nuclear physics properties: fission distribution, β-decay, and mass model. For most cases, the actinides are overproduced in our models if the initial conditions are sufficiently neutron-rich for fission cycling. Here, we find that actinide production can be so robust in the dynamical ejecta that an additional lanthanide-rich, actinide-poor component is necessary in order to match observations of actinide-boost stars. We present a simple actinide-dilution model that folds in estimated contributions from two nucleosynthetic sites within a merger event. Our study suggests that while the dynamical ejecta of an NSM are likely production sites for the formation of actinides, a significant contribution from another site or sites (e.g., the NSM accretion disk wind) is required to explain abundances of r-process-enhanced, metal-poor stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Literature Review of Documented Persistent Secondary Uranium Sources at DOE LM Sites

The U.S. Department of Energy Office of Legacy Management (LM) has the responsibility to protect human health and the environment at 101 legacy uranium sites across the nation by implementing long-term surveillance and maintenance requirements. Included in these responsibilities is to provide for the safe and environmentally sound disposal, long-term stabilization, and control of uranium mill tailings in a manner that minimizes or eliminates radiation health hazards to the public at 21 Uranium Mill Tailings Radiation Control Act (UMTRCA) Title I sites; 6 UMTRCA Title II sites; and 9 Comprehensive Environmental Response, Compensation, and Liability Act/Resource Conservation and Recovery Act sites. These sites are former uranium mills (processing sites), or sites with wastes from processing sites encapsulated in engineered disposal cells (disposal sites), or sites where components of atomic weapons were manufactured. During site characterization, groundwater contamination was often discovered and evaluated or modeled for future movement with the inclusion of uranium sorption in the groundwater plume. Much of this original work did identify uranium processing wastes as the original contaminant source. However, more recent work has identified the existence and influence of persistent secondary uranium sources that may prevent uranium in groundwater from attenuating as quickly as initial evaluations or modeling indicated. By understanding these processes, better decisions can be made to choose the appropriate remediation and closure strategies, identify long-term surveillance and maintenance activities, and more accurately estimate life-cycle costs. The purpose of this report is to summarize the literature available on documented persistence secondary uranium sources at LM sites as part of the Applied Studies and Technology Persistent Secondary Contaminant Sources Project.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗