Search NASA⌕ Search

SEARCH · Search NASA

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 487 records · Page 27

Smart-PGSim: Using Neural Network to Accelerate AC-OPF Power Grid Simulation

In this work we address the problem of accelerating complex power-grid simulation through machine learning ( ML). Specifically, we develop a framework, Smart-PGSim,which generates multitask-learning (MTL) neural network (NN)models to predict the initial values of variables critical to the problem convergence. MTL models allow information sharing when predicting multiple dependent variables while including customized layers to predict individual variables. We show that,to achieve the required accuracy, it is paramount to embed domain-specific constraints derived from the specific power-grid components in the MTL model. Smart-PGSim then employs the predicted initial values as a high-quality initial condition for the power-grid numerical solver (warm start), resulting in both higher performance compared to state-of-the-art solutions while maintaining the required accuracy. Smart-PGSim brings 2.60×speedup on average (up to 3.28×) computed over 10,000 problems, without losing solution optimality.

machine learning, neural networks↗

A new soil mixing layer model for simulating conservative solute loss from initially saturated soil to surface runoff

Surface water pollution due to solute loss from soil to surface runoff is a serious threat to the environment, and it is of great significance to predict the solute loss. While soil mixing layer theory has been widely used to predict the solute loss, it has been observed that the conventional soil mixing layer models with constant mixing layer depth underestimate solute concentration in surface runoff, and thus underestimate the solute loss after surface runoff starts. This study introduced a new soil mixing layer model with a new concept that the mixing layer depth increases downward with time after surface runoff occurs. This new concept was implemented by adding a mixing front between the soil mixing layer and its underlying soil layer, which is conceptually similar to the wetting front of soil infiltration given by the Green-Ampt equation. This is equivalent to treating the soil mixing layer as a deforming control volume. A mechanistic explanation for the new model was presented. This study also developed a closed form solution of the new model. The new model was evaluated by simulating five laboratory experiments conducted by two groups of researchers. The first two experiments have ponding water, and consider different drainage and runoff conditions at both early and late runoff times. The other three experiments do not have ponding water, and consider different infiltration conditions at the bottom of soil boxes. Simulation results of the five experiments show that the new model resolves the underestimation problem of conventional soil mixing layer models. The new model is subject to a total of six limits such as that the new model only considers conservative solute. Three limits are related to the assumptions used to derive the closed form solution, and they may be alleviated in a future research as briefly discussed at the end of the paper.

54 ENVIRONMENTAL SCIENCES↗

A Note on the Probability of Initiation Problem

The purpose of this note is to draw attention to Pál for his derivation of the Pál-Bell equation which is used by LLNL and LANL to model Probability of Initiation (POI) problems. Both laboratories credit the equation to Bell and Lee. Although the equation formulated by Bell and Lee and the equation derived by Pál appear to be the same, the equation by Bell and Lee is flawed. Scattering, which was absent in Bell’s original formulation of a POI problem, was later introduced by Bell and Lee by replacing a fission term with a scattering term. Such a replacement, however, leads to an incorrect average number $\overline{v}$ of neutrons that is produced in a fission event. By approaching a POI problem from probabilistic point of view, Pál formulates the problem as a branching process that i) describes the growth of a population of neutrons in a fissile medium, and ii) predicts $\overline{v}$ correctly. There are five reasons for writing this report. The first is to clarify the failing in Bell and Lee’s derivation. The second is to establish the origin of the equation used by the laboratories. The third is to elucidate the fixed point technique formulated by Mokhtar-Kharroubi to solve the Pál-Bell equation. The fourth is to derive the α eigenvalue that is associated with a super-critical situation. The fifth is to describe in the Appendix a second order positivity preserving technique that accelerates the convergence of the first order fixed point method developed by Mokhtar-Kharroubi.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Complementing the CCS Class VI Well Permit Process with DOE-NETL's SMART Initiative Tools and Workflows

This is a presentation on model explorer developed under SMART initiative Task 2. Our team will present the current progress of the model explorer in using machine learning models to accelerate CCS project at GWPC meeting. Model explorer bring new capabilities, (fast, Realtime, and accurate) that can help CCS stakeholders including regulatory agencies, public and site operators make faster decisions and process information and data.

Hosseini, Seyyed↗

Effects of finite element discretization and model simplification on calculations of ductile failure initiation

The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.

97 MATHEMATICS AND COMPUTING↗

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗

Status Update: Deposition Modeling For SNF Canister CISCC

This report fulfills the M3 milestones M3SF-21PN010207025 & M3SF-20PN0102070412. During fiscal year (FY) 2020, Pacific Northwest National Laboratory (PNNL) worked to further develop the FY 2019 deposition and particle tracking models. This status report outlines these efforts and presents the progress made so far. Model development work is ongoing and is planned to continue in FY 2021. The FY 2020 model development included work on: Wind Effects. Model development and sensitivity studies, investigated how wind direction and speed affect deposition; Brownian Motion. Implementing Brownian Motion into existing models; Particle Size Variability. Depending on the particle composition, the diameter of the particle may vary with changes in relative humidity. Models were developed to analyze this; Multiphase and Fluid Film Modeling. Investigating canister surface wetting and drying, and how this effects overall deposition. Models were developed to analyze this; Difusophoresis. Performing initial work to implement diffusiophoresis into the existing models; Turbophoresis. Performing initial work to implement turbophoresis into the existing models. Much of this work will continue into FY21. The authors present initial results and discus current and future work.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimating nucleon substructure properties in a unified model of p-Pb and Pb-Pb collisions

We apply a well tested hybrid transport model, which couples viscous hydrodynamics to a hadronic afterburner, to describe bulk observables in proton-lead and lead-lead collisions at $\sqrt{s}$$_ {NN}$ = 5.02 TeV. The quark-gluon plasma (QGP) initial conditions are modeled using the parametric TRENTo model with additional nucleon substructure parameters to vary the number and size of hot spots inside each nucleon, followed by a pre-equilibrium free streaming stage to match the full energy-momentum tensor of the initial state onto viscous hydrodynamics. Initial condition and QGP medium parameters, such as the temperature dependence of the QGP shear and bulk viscosities, are then calibrated using Bayesian parameter estimation to describe charged particle yields, mean p T and anisotropic flow harmonics of both collision systems in a single self-consistent framework. We find that the hybrid model provides a compelling, simultaneous description of both collision systems using appropriately chosen model parameters, and present new posterior estimates for the size and shape of the nucleon and temperature dependence of QGP shear and bulk viscosities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Inclusive jet and hadron suppression in a multistage approach

We present a new study of jet interactions in the quark-gluon plasma created in high-energy heavy-ion collisions, using a multistage event generator within the jetscape framework. We focus on medium-induced modifications in the rate of inclusive jets and high transverse momentum (high-p T ) hadrons. Scattering-induced jet energy loss is calculated in two stages: a high virtuality stage based on the matter model, in which scattering of highly virtual partons modifies the vacuum radiation pattern, and a second stage at lower jet virtuality based on the LBT model, in which leading partons gain and lose virtuality by scattering and radiation. Coherence effects that reduce the medium-induced emission rate in the matter phase are also included. The TRENTO model is used for initial conditions, and the (2 + 1)dimensional VISHNU model is used for viscous hydrodynamic evolution. Jet interactions with the medium are modeled via 2-to-2 scattering with Debye screened potentials, in which the recoiling partons are tracked, hadronized, and included in the jet clustering. Holes left in the medium are also tracked and subtracted to conserve transverse momentum. Calculations of the nuclear modification factor (R AA ) for inclusive jets and high-p T hadrons are compared to experimental measurements at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC). Within this framework, we find that with one extra parameter which codifies the transition between stages of jet modification—along with the typical parameters such as the coupling in the medium, the start and stop criteria, etc.—we can describe these data at all energies for central and semicentral collisions without a rescaling of the jet transport coefficient q.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Home charging for all: Techno-economic and life cycle assessment of multi-unit dwelling electric vehicle charging hubs

Ubiquitous electric vehicle adoption can drastically reduce greenhouse gas emissions (GHG) but will require equitable home charging infrastructure for all residences. Unlike single-family homes, Multi-Unit Dwellings (MUD) currently lack market share and access to nearby charging infrastructure partly due to expensive capital costs for reluctant residential property-owners. Further, this work evaluates the levelized cost of charging (LCOC) for Battery Electric Vehicles (BEVs) at MUD community charging hubs through a techno-economic analysis (TEA) that leverages real-world charging data and costs. Three different MUD charging types are investigated for a baseline and optimistic case: 1.9-kW Level 1 (L1), 6.6-kW Level 2 (L2), and 50-kW Direct Current Fast Charging (DCFC) stations under three different ownership models: resident, utility, and private company. Results demonstrate L1 and L2 chargers to be less expensive than the gasoline equivalent across the United States. Further, utility and private company ownership models, which avoid initial costs for residential property owners, result in a large LCOC premium for L1 and moderate LCOC premium for L2 relative to the resident ownership model. In contrast, DCFC is shown to be expensive for baseline scenarios but economical for optimistic scenarios especially under private company ownership. This work also performs a cradle to grave (C2G) life cycle assessment (LCA) of an average passenger BEV and gasoline conventional vehicle (CV) using yearly (grid mix & vehicle parameters), hourly (grid mix), and state-level (grid mix) resolution. Results show BEVs to have lower GHG emissions (-86% to-10%) than gasoline CVs in the contiguous U.S. Next, the system boundary of the TEA is extended to the total cost of ownership (TCO) of BEVs. The TCO is then coupled with the C2G GHG emissions to calculate the cost of GHG emissions reduction. Ultimately, the cost of GHG emissions reduction from MUD BEVs relative to gasoline CVs is shown to be negative for every scenario except baseline DCFC, meaning MUD BEV charging infrastructure can be a cost-effective endeavor to reduce GHG emissions.

33 ADVANCED PROPULSION SYSTEMS↗

DEPLOYING FAST CHARGING INFRASTRUCTURE FOR ELECTRIC VEHICLES IN URBAN NETWORKS: AN ACTIVITY-BASED APPROACH

This paper explores an important problem under the domain of network modeling, the optimal configuration of charging infrastructure for electric vehicles (EVs) in urban networks considering EV users' daily activities and charging behavior. This study proposes a charging behavior simulation model considering different initial state of charge (SOC), travel distance, availability of home chargers, and the daily schedule of trips for each traveler. The proposed charging behavior simulation model examines the complete chain of trips for EV users as well as the interdependency of trips traveled by each driver. The problem of finding the optimum charging configuration is then formulated as a mixed-integer nonlinear programming problem that considers the dynamics of travel time and travel distance, the interdependency of trips made by each driver, limited range of EVs, remaining battery capacity for recharging, waiting time in queue, and detour to access a charging station. This problem is solved using a metaheuristic approach for a large-scale case network. A series of examples are presented to demonstrate the model efficacy and explore the impact of energy consumption on the final SOC and the optimum charging infrastructure.

Chain of Trips↗

Machine Learning for LBNF Beam Diagnostics

This paper focuses on developing a machine learning model for predicting initial beam parameters for the Long Baseline Neutrino Facility (LBNF) beamline using downstream muon monitor data. Parameters such as proton beam position on target, sigma on target, focusing horn current, and focusing horn tilt are parameters we anticipate to be predictable based on the muon monitors. Uncertainty in initial beam condition measurements are a major contributor to uncertainty in downstream flux, and over operation time beam misalignment can occur [1]. A machine learning model has promise to detect anomalies along the beamline based on discrepancies between predicted configurations and measured configurations, and thus can expedite error detection and handling. A PyTorch neural network is defined, trained, and tested. The developed model currently does not provide reliable predictions, with the lowest loss being 0.09.. Further steps to improve the model’s accuracy are discussed, as well as future plans to detect anomalous beam conditions using a digital twin.

O'Brien, Bridget [Fermilab]↗

An experimental, theoretical, and kinetic modeling study of post-flame oxidation of ammonia

The post -flame oxidation rate of ammonia was investigated in a novel atmospheric pressure flow reactor at temperatures of 1280 ± 16 K and as a function of residence time and mixture composition (1-10% O 2 , dry and moist). The experimental results, as well as selected data from literature, were analyzed using an updated detailed chemical kinetic model. The medium temperature, very lean conditions enhance the importance of reactions of the nitroxyl (HNO) intermediate. High-level theory was used to calculate the rate constant for HNO + NH 2 , indicating that this step is significantly faster than values used in literature. Furthermore, a trajectory based approach was used to determine collision efficiencies for selected bath gases for HNO + M. The experimental results show that the NH 3 oxidation rate increases with temperature and O 2 concentration, while the presence of water vapor slightly inhibits reaction. Formation of NO and N 2 O was strongly promoted at higher levels of O 2 . Modeling results agreed well with the measurements, except at the lowest level of O 2 . The predicted oxidation rate of NH 3 was shown to result from a delicate balance between chain branching and terminating steps involving NH 2 , H 2 NO, and HNO. Recent theoretical work on reactions of these species by Klippenstein and coworkers and Stagni et al. was instrumental in improving modeling predictions. After initiation, NO reached a pseudo -steady-state level, where the pathways to NO were largely balanced by the NH 2 + NO reaction. Nitric oxide was partly oxidized to NO 2 , with the NH 2 + NO 2 reaction responsible for most of the N 2 O formation. Novelty and significance statement: This study provides the first detailed kinetic analysis of the lean postflame oxidation of ammonia, based on time -resolved flow reactor data in a novel reactor. In addition to the post -flame oxidation rate of ammonia, data for formation of NO and N 2 O were compared with modeling predictions. The medium temperature, very lean conditions enhance the importance of reactions of the HNO and H 2 NO intermediates. Inclusion in the model of results from recent high-level theoretical work, including present calculations for HNO + NH 2 and HNO + M, was crucial for capturing the observed behavior. It is argued that the post -initiation steady-state NH 3 oxidation rates constitute important data for model validation, along with ignition delays and laminar flame speeds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Climate Response to Emissions Reductions Due to COVID-19: Initial Results From CovidMIP

Many nations responded to the COVID-19 pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO 2 , other greenhouse gases and ozone and aerosol precursors. We perform a coordinated Intercomparison, CovidMIP, of Earth System model simulations to assess the impact on climate of these emissions reductions. Eleven models performed multiple initial-condition ensembles to produce over 280 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over East Asia) and associated increases in surface shortwave radiation levels. However, any impact on near-surface temperature or rainfall during 2020-2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID-19-related emission reductions on near-term climate.

54 ENVIRONMENTAL SCIENCES↗

Assignment and calibration of relative permeability by hydrostratigraphic units for multiphase flow analysis, case study: CO 2 -EOR operations at the Farnsworth Unit, Texas

Included in the most critical factors for geological CO 2 storage site screening, selection, and operation is effective simulations of multiphase flow and transport. Relative permeability is probably the greatest source of potential uncertainty in multiphase flow simulation, second only to intrinsic permeability heterogeneity. The specific relative permeability relationship assigned greatly impacts forecasts of CO 2 trapping mechanisms, phase behavior, and long-term plume movement. A primary goal of this study is to evaluate the impacts and implications of different methods of assigning relative permeability relationships for CO 2 -EOR model forecasts. Most simulation studies published in the literature base selection of relative permeability functions on the geologic formation or rock type alone. In this study, we initially implemented reservoir model grids with previously-identified hydrostratigraphic units based on porosity and permeability relationship of the Morrow ‘B’ Sandstone, then assigned relative permeability functions for those hydrostratigraphic units. Specific, constrained relative permeability relationships were created and assigned to each hydrostratigraphic unit using petrophysical data and Mercury Intrusion Capillary Pressure (MICP) measurements, from core samples of each hydrostratigraphic unit. Results of forward simulations with the newly-calibrated models will be compared to those of previous models as well as to simulation results for a range of different relative permeability relationships. The study site is the Farnsworth Unit (FWU) in the northeast Texas Panhandle, an active CO 2 -EOR operation. The target formation is the Morrow ‘B’ Sandstone, a clastic formation composed of medium to course sands.

58 GEOSCIENCES↗

Engineering infinite-range SU($\textit{n}$) interactions with spin-orbit-coupled fermions in an optical lattice

Here, we study multilevel fermions in an optical lattice described by the Hubbard model with on-site SU($\textit{n}$) -symmetric interactions. We show that in an appropriate parameter regime this system can be mapped onto a spin model with all-to-all SU($\textit{n}$)-symmetric couplings. Raman pulses that address internal spin states modify the atomic dispersion relation and induce spin-orbit coupling, which can act as a synthetic inhomogeneous magnetic field that competes with the SU($\textit{n}$) exchange interactions. We investigate the mean-field dynamical phase diagram of the resulting model as a function of $\textit{n}$ and different initial configurations that are accessible with Raman pulses. Consistent with previous studies for $\textit{n}$ = 2 , we find that for some initial states the spin model exhibits two distinct dynamical phases that obey simple scaling relations with $\textit{n}$. Moreover, for $\textit{n}$ > 2 we find that dynamical behavior can be highly sensitive to initial intraspin coherences. Our predictions are readily testable in current experiments with ultracold alkaline-earth-metal(-like) atoms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High speed PLIF study of the Richtmyer–Meshkov instability upon re-shock

The Richtmyer–Meshkov instability (RMI) of a twice-shocked gas interface is researched using high-speed planar-laser induced fluorescence (PLIF) in the Wisconsin Shock Tube Laboratory’s vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium-acetone mixture and argon. This IC is accelerated by a shock of nominal strength M = 1.8, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An estimate of the light gas mole fraction is extracted from high-speed imaging using an iterative process that accounts for the nonlinear temperature dependence of the acetone’s fluorescence quantum yield (FQY) and absorption cross-section. A vorticity deposition model for the initial growth rate after reshock is compared with the Mikaelian model for re-shock. Previously used in literature, the number of generations is shown to naturally arise from a normalisation of the scalar transport equation. A self-similar analysis is then further performed using the mole fraction data to explore the evolution of the RMI after reshock and the higher order moments of the light gas mole fraction are compared with a proposed model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-driven causal model discovery and personalized prediction in Alzheimer's disease

Abstract With the explosive growth of biomarker data in Alzheimer’s disease (AD) clinical trials, numerous mathematical models have been developed to characterize disease-relevant biomarker trajectories over time. While some of these models are purely empiric, others are causal, built upon various hypotheses of AD pathophysiology, a complex and incompletely understood area of research. One of the most challenging problems in computational causal modeling is using a purely data-driven approach to derive the model’s parameters and the mathematical model itself, without any prior hypothesis bias. In this paper, we develop an innovative data-driven modeling approach to build and parameterize a causal model to characterize the trajectories of AD biomarkers. This approach integrates causal model learning, population parameterization, parameter sensitivity analysis, and personalized prediction. By applying this integrated approach to a large multicenter database of AD biomarkers, the Alzheimer’s Disease Neuroimaging Initiative, several causal models for different AD stages are revealed. In addition, personalized models for each subject are calibrated and provide accurate predictions of future cognitive status.

Zheng, Haoyang (ORCID:0000000168358242)↗