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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A Generic Modeling and Development Approach for WECC Composite Load Model

The composite load model (CLM) proposed by Western Electricity Coordinating Council (WECC) gains increasing traction in the industry, particularly in North America. At the same time, it is recognized that further improvements in the structure, initialization and aggregation methods are needed. However, lack of an open implementation of WECC CLM becomes a roadblock for many researchers to work on improving it. This paper presents the first open reference implementation of WECC CLM. Individual load components and the CLM are first developed and tested in Matlab, then translated to the high performance computing (HPC), parallel simulation platform GridPACK. Main contributions of the paper include: 1) presenting some important yet undocumented details of modeling and initializing the CLM, particularly for a parallel simulation platform like GridPACK; 2) implementation details of the load components such as the single-phase air conditioner motor; 3) implementing the CLM in a modular and extensible manner. The implementation has been tested at both the component as well as system levels and bench-marked against commercial simulation programs.

Huang, Qiuhua↗
DOE OSTI2022

Bimolecular Peroxy Radical (RO 2 ) Reactions and Their Relevance in Radical Initiated Oxidation of Hydrocarbons

The kinetics of peroxy radical (RO 2 ) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH 3 O 2 + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO 2 + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO 2 + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO 2 + OH reaction is a significant sink of RO 2 radicals at the early stage of autoignition, affecting fuel oxidation through RO 2 ↔ QOOH, RO 2 ↔ alkene + HO 2 , or RO 2 + HO 2 ↔ ROOH + O 2 . Additionally, the kinetic stability of the trioxide formed from RO 2 + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). In conclusion, the present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

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

Short-time polymerization and convolution kinetics in the analysis of the initiation step of isoprene polymerization over neodymium catalyst

The paper considers a convolution model for polymerization of 1,3-dienes in the presence of organometallic catalysts. It is shown that initial steps of polymerization can be represented as a set of shifted parallel reactions of the chain propagation, in which shift value are governed by the kinetics of the initiation step. In this case, polymer accumulation curve is described by an equation containing convolution of functions which allow one to solve an inverse problem of restoring the function of active site accumulation. Based on the example of isoprene polymerization over neodymium catalyst, it was demonstrated that the results of the inverse problem solution are in good agreement with the experimental data on concentration of the active sites. These data were obtained by a method of controlled introduction of inhibitor in short-time polymerization.

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

Higgs-mass predictions in the MSSM and beyond

Abstract Predictions for the Higgs masses are a distinctive feature of supersymmetric extensions of the Standard Model, where they play a crucial role in constraining the parameter space. The discovery of a Higgs boson and the remarkably precise measurement of its mass at the LHC have spurred new efforts aimed at improving the accuracy of the theoretical predictions for the Higgs masses in supersymmetric models. The “ Precision SUSY Higgs Mass Calculation Initiative ” (KUTS) was launched in 2014 to provide a forum for discussions between the different groups involved in these efforts. This report aims to present a comprehensive overview of the current status of Higgs-mass calculations in supersymmetric models, to document the many advances that were achieved in recent years and were discussed during the KUTS meetings, and to outline the prospects for future improvements in these calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

A Recipe for implementing the Arrhenius-Shock-Temperature State Sensitive WSD (AWSD) model, with parameters for PBX 9502

A reactive ow model for the tri-amino-tri-nitro-benzene (TATB) based plastic bonded explosive PBX 9502 is presented. This newly devised model is based primarily on the shock temperature of the material, along with local pressure, and accurately models a broader range of detonation and initiation scenarios. The equation of state for the reactants and products, as well as the thermodynamic closure of pressure and temperature equilibration are carried over from the Wescott-Stewart-Davis (WSD) model7,8. Thus, modifying an existing WSD model in a hydrocode should be rather straightforward.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2018

Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution

Hydrodynamics, oxygen transfer, and oxygen uptake in bubble-column and airlift bioreactors were studied using multiphase Euler-Euler computational fluid dynamics (CFD) simulations. Interphase mass transfer of oxygen was modeled with a liquid mass-transfer coefficient and sustained driving force provided by a phenomenological model for microbial oxygen uptake rate (OUR). Laboratory-scale reactor simulations of a bubble-column (0.15 m diameter and 0.75 m initial liquid height) showed reasonable agreement with gas holdup and mass transfer experiments reported in the literature. Commercial-scale bubble-column and airlift reactor simulations (5 m diameter and 25 m initial liquid height) were simulated using this validated model. Similar trends for the variation of overall gas holdup and oxygen concentrations with superficial velocity were observed for both of the commercial-scale reactors. However, the simulation results indicate differences in hydrodynamics, such as better recirculation in the airlift reactors that enables favorable distribution of oxygen in the reactor.

09 BIOMASS FUELS↗
DOE OSTI2023

QuDPy: A Python-based tool for computing ultrafast non-linear optical responses

Nonlinear Optical Spectroscopy is a well-developed field with theoretical and experimental advances that have benefited multiple disciplines, including chemistry, biology, and physics. However, for the accurate interpretation of the corresponding multi-dimensional spectra, there is a need for precise quantum dynamical simulations based on model Hamiltonians. In this article, we present the initial release of our code, QuDPy (Quantum Dynamics in Python), which provides a robust numerical platform for performing quantum dynamics simulations based on model systems, including open quantum systems. Furthermore, a distinguishing feature of our approach is the ability to specify various high-order optical response pathways in the form of double-sided Feynman diagrams through a straightforward input syntax. This syntax outlines the time-ordering of ket-sided or bra-sided optical interactions acting on the time-evolving density matrix of the system. We utilize the quantum dynamics capabilities of QuTip to simulate the spectral response of complex systems, allowing us to compute virtually any $n$-th order optical response of the model system. To illustrate the utility of our approach, we provide a series of example calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Beam-energy dependence of correlations between mean transverse momentum and anisotropic flow of charged particles in Au+Au collisions at RHIC

The correlation between the mean transverse momentum, [ p T ], and the squared anisotropic flow, $v^2_n$, on an event-by-event basis has been suggested to be influenced by the initial conditions in heavy-ion collisions. We present measurements of the variances and covariance of [ p T ] and $v^2_n$, along with their dimensionless ratio, for Au+Au collisions at various beam energies: $\sqrt{s_{NN}}$ = 14.6, 19.6, 27, 54.4, and 200 GeV. Our measurements reveal a distinct energy-dependent behavior in the variances and covariances. In addition, the dimensionless ratio displays a similar behavior across different beam energies. We compare our measurements with hydrodynamic models and similar measurements from Pb+Pb collisions at the Large Hadron Collider (LHC). These findings provide valuable insights into the beam energy dependence of the specific shear viscosity ( η/s ) and initial-state effects, allowing for differentiating between different initial-state models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Characterizing the initial conditions of heavy-ion collisions at the LHC with mean transverse momentum and anisotropic flow correlations

Correlations between mean transverse momentum [pT] and anisotropic flow coefficients v 2 or v 3 are measured as a function of centrality in Pb–Pb and Xe–Xe collisions at $\sqrt{s{NN}}$ = 5.02 TeV and 5.44 TeV, respectively, with ALICE. In addition, the recently proposed higher-order correlation between [p T ], v 2 , and v 3 is measured for the first time, which shows an anticorrelation for the presented centrality ranges. These measurements are compared with hydrodynamic calculations using IP-Glasma and T R ENTo initialstate shapes, the former based on the Color Glass Condensate effective theory with gluon saturation, and the latter a parameterized model with nucleons as the relevant degrees of freedom. The data are better described by the IP-Glasma rather than the T R ENTo based calculations. In particular, Trajectum and JETSCAPE predictions, both based on the T R ENTo initial state model but with different parameter settings, fail to describe the measurements. As the correlations between [p T ] and v n are mainly driven by the correlations of the size and the shape of the system in the initial state, these new studies pave a novel way to characterize the initial state and help pin down the uncertainty of the extracted properties of the quark–gluon plasma recreated in relativistic heavy-ion collisions

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

R-parity violating decays of Wino chargino and wino neutralino LSPs and NLSPs at the LHC

The R-parity violating decays of both Wino chargino and Wino neutralino LSPs are analyzed within the context of the B − L MSSM “heterotic standard model”. These LSPs correspond to statistically determined initial soft supersymmetry breaking parameters which, when evolved using the renormalization group equations, lead to an effective theory satisfying all phenomenological requirements; including the observed electroweak vector boson and Higgs masses. The explicit decay channels of these LSPs into standard model particles, the analytic and numerical decay rates and the associated branching ratios are presented. The decay lengths of these RPV interactions are discussed. It is shown that the vast majority of these decays are “prompt”, although a small, but calculable, number correspond to “displaced vertices” of various lengths. It is demonstrated that for a Wino chargino LSP, the NLSP is the Wino neutralino with a mass only slightly higher than the LSP — and vice-versa. As a consequence, we show that both the Wino chargino and Wino neutralino LSP/NLSP R-parity violating decays should be simultaneously observable at the CERN LHC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

GrainGNN: A dynamic graph neural network for predicting 3D grain microstructure

We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Nuclear theory and event generators for charge-changing neutrino reactions

Semi-inclusive $CC\nu$ cross sections based on factorized cross sections are studied for a selection of spectral function models with the objective of facilitating the choice of models for use as input into event generators. The basic formalism for such cross sections is presented along with an introduction to constructing spectral functions for simple models based on the independent particle shell model (IPSM), the relativistic Fermi gas model (RFG) and a local density approximation (LDA) based on the RFG. Spectral functions for these models are shown for $^{16}$O along with a more sophisticated model which includes nucleon-nucleon interactions.. Inclusive and semi-inclusive cross sections are calculated for these models. Although the inclusive cross sections are all of similar size and shape, the semi-inclusive cross sections are subtantially different depending upon whether the spectral functions contain some features associated with the nuclear shell model or are based on the RFG and LDA models. Here, calculations of average values and standard deviations of the initial neutrino energy using the semi-inclusive cross section for the various models are presented and indicate that there may be simple kinematical descriptions of the average neutrino energy which is common to all of these models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Capturing Infrastructure Interdependencies for Power Outages Prediction During Extreme Events

As extreme weather events such as hurricanes, severe thunderstorms, and floods grow in frequency and intensity, the disruption of power grid systems poses significant challenges, including widespread electrical outages, economic losses, and threats to public safety. This paper presents a forward-looking approach that leverages geographical graph-based machine learning models to predict county-level maximum power outages during such events. By capturing the intricate interdependencies within power system networks, our approach aims to provide precise and actionable predictions that can optimize emergency response efforts and enhance grid resilience. Through the integration of real-world data, including hurricane advisories and power outage records, we have trained and benchmarked multiple machine learning models, demonstrating the feasibility and potential of this method. While our initial results are promising, this paper also charts a course for advancing these models, addressing the remaining challenges, and ultimately transforming how we anticipate and respond to the impacts of extreme weather on power systems.

Lee, Sangkeun (Matt) [ORNL] (ORCID:000000021317511↗
DOE OSTI2017

GeneratorSE: A Sizing Tool for Variable-Speed Wind Turbine Generators

This report documents a set of analytical models employed by the optimization algorithms within the GeneratorSE framework. The initial values and boundary conditions employed for the generation of the various designs and initial estimates for basic design dimensions, masses, and efficiency for the four different models of generators are presented and compared with empirical data collected from previous studies and some existing commercial turbines. These models include designs applicable for variable-speed, high-torque application featuring direct-drive synchronous generators and low-torque application featuring induction generators. In all of the four models presented, the main focus of optimization is electromagnetic design with the exception of permanent-magnet and wire-wound synchronous generators, wherein the structural design is also optimized. Thermal design is accommodated in GeneratorSE as a secondary attribute by limiting the winding current densities to acceptable limits. A preliminary validation of electromagnetic design was carried out by comparing the optimized magnetic loading against those predicted by numerical simulation in FEMM4.2, a finite-element software for analyzing electromagnetic and thermal physics problems for electrical machines. For direct-drive synchronous generators, the analytical models for the structural design are validated by static structural analysis in ANSYS.

17 WIND ENERGY↗
DOE OSTI2018

Bison Fuel Performance Modeling on RIA in a PWR

This report is prepared to fulfill a Level 2 milestone for CASL on “deliver and document RIA fuel performance modeling using Bison fuel performance code.” The work performed is focused on modeling the full length commercial fuel rod response under a design basis Reactivity Initiated Accident (RIA), in particular, a hypothetical Control Rod Ejection Accident (CREA) in a Pressurized Water Reactor (PWR). Fuel rods with different burnups, 0 (fresh fuel), 20, 40, and 60 GWd/tU have been modeled with a reactivity worth insertion of $1.24 at Hot Zero Power (HZP) condition. Bison failure modeling method, and a brief summary of validations on selected RIA cases from CABRI and NSRR testing facilities are also described. The work performed has demonstrated the successful application of Bison code in the modeling of thermal and mechanical responses in RIAs for PWRs. In particular, Bison code has a more mature capability in the modeling of the frictional contact between pellet and cladding of the high burnup fuels, which is critical for modeling the Pellet Clad Mechanical Interaction (PCMI). The modeling of RIA using Bison code shows fuels in PWR environment could have different responses than in the simulated RIA tests, and results indicate that a transition from PCMI to Departure from Nucleate Boiling (DNB) could happen at relatively lower enthalpies: a better heat transfer from the fuel to coolant could increase cladding temperature and ductility, but at the same time also reduces the margin to DNB, which could cause high temperature oxidation and ballooning type of failures. The precise determination of the transition into DNB may rely on a more accurate input of the power histories in the simulated RIAs from VERA, the core simulator of CASL. The limiting conditions in PWRs could be quite different than those tests with non-typical testing conditions, and the improvement of the modeling of thermal responses, such as, clad coolant heat transfer, gap heat transfer, clad thermal inertia, and radial power profiles could be important considerations of the RIAs. In modeling the fuel responses in a fast transient, a number of variables could affect the initial burnup states, and some material models, e.g., fuel mechanical models, fission gas release, and corrosion models, and parameters can affect modeling results as well. Those uncertainties and advanced modeling capabilities need to be investigated in the future study of RIAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Cloud, Aerosol, and Complex Terrain Interactions (CACTI) Field Campaign Report

General circulation models and downscaled regional models exhibit persistent biases in deep convective initiation location and timing, cloud top height, stratiform area and precipitation fraction, and anvil coverage (e.g., Del Genio 2012, Del Genio et al. 2012, Hohenegger and Stevens 2013, Song et al. 2013). Despite important impacts on the distribution of atmospheric heating, moistening, momentum, and precipitation (e.g., Hartmann et al. 1984, Fritsch et al. 1986, Houze 1989, 2004, Donner et al. 2001, Del Genio and Kovari 2002, Schumacher et al. 2004, Nesbitt et al. 2006, Storelvmo 2012), nearly all climate models fail to represent mesoscale convective organization (Del Genio 2012), while system evolution is not represented at all (Ovchinnikov et al. 2006). Recent advances in cumulus parameterization coupled with increasing model resolution have improved predictions, but even relatively higher-resolution models without parameterized deep convection have some persistent kinematic and microphysical biases (e.g., Blossey et al. 2007, Matsui et al. 2009, Luo et al. 2010, Lang et al. 2011, Varble et al. 2011, Fridlind et al. 2012, Hagos et al. 2014, Varble et al. 2014a-b, Fan et al. 2017, Stanford et al. 2017, Han et al. 2019). To improve representation of convective systems in models requires adequate characterization of their predictability as a function of environmental conditions. Because of the significant sensitivities of deep convective initiation, intensity, lifetime, propagation, and mesoscale convective organization to many factors including multi-scale atmospheric circulations, ambient environmental stability, humidity, wind distributions, and cloud microphysical processes, this characterization relies on comprehensively observing many cases of convective initiation, non-initiation, organization, and non-organization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Droplet and Particle Methods to Investigate Turbulent Particle Laden Jets

The SARS-CoV-2 pandemic has heightened the interest in particle-laden turbulent jets generated by breathing, talking, coughing and sneezing, and how these can contribute to disease transmission. We present quantitative measurement methods for such flows, while exploring and offering improvements for common shortcomings. We generate jets consisting of either liquid droplets or solid particles in an isothermal, quiescent and electrically isopotential experimental chamber that was constructed to control the effects of ambient forcing on jet behavior. For liquid droplets, we find promise in surface deposition analysis based on fluorescent tracer use. For particles, we explore the performance of commercially available adhesive sampling strips and develop conductive grounded carbon tape based sampling strips. We explore ways in which the smallest of thermal gradients or electrostatic charge issues can affect particle dispersion, and suggest practical methods to address these issues. The developed methods are applied to study the simultaneous deposition of (Formula presented.) 25, 50 and 200 μm solid particles from a particle laden turbulent jet with a mean velocity of 33.2 m/s. The deposition location as a function of particle size was compared to results from a simple numerical RANS model, and illustrates ways in which imprecise initial or boundary conditions can lead to a notable deviation from experimental results. The differences in deposition pattern seen in experimental and numerical results despite a carefully controlled environment and characterized particle ejection indicate the need for a more stringent numerical model validation, especially when studying fate and transport of mid-range (neither purely aerosol or ballistic) sized particles.

42 ENGINEERING↗
DOE OSTI2021

Synthetic observables for electron-capture supernovae and low-mass core collapse supernovae

Stars in the mass range from 8 M ⊙ to 10 M⊙ are expected to produce one of two types of supernovae (SNe), either electron-capture supernovae (ECSNe) or core-collapse supernovae (CCSNe), depending on their previous evolution. Either of the associated progenitors retain extended and massive hydrogen-rich envelopes and the observables of these SNe are, therefore, expected to be similar. In this study, we explore the differences in these two types of SNe. Specifically, we investigate three different progenitor models: a solar-metallicity ECSN progenitor with an initial mass of 8.8 M ⊙ , a zero-metallicity progenitor with 9.6 M ⊙ , and a solar-metallicity progenitor with 9 M ⊙ , carrying out radiative transfer simulations for these progenitors. Here, we present the resulting light curves for these models. The models exhibit very low photospheric velocity variations of about 2000 km s -1 ; therefore, this may serve as a convenient indicator of low-mass SNe. The ECSN has very unique light curves in broad-bands, especially the U band, and does not resemble any currently observed SN. This ECSN progenitor being part of a binary will lose its envelope for which reason the light curve becomes short and undetectable. The SN from the 9.6 M ⊙ progenitor exhibits also quite an unusual light curve, explained by the absence of metals in the initial composition. The artificially iron-polluted 9.6 M ⊙ model demonstrates light curves closer to normal SNe IIP. The SN from the 9 M ⊙ progenitor remains the best candidate for so-called low-luminosity SNe IIP like SN 1999br and SN 2005cs.

79 ASTRONOMY AND ASTROPHYSICS↗