Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “range prediction”

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

A physics-informed deep learning model of the hot tail runaway electron seed

A challenging aspect of the description of a tokamak disruption is evaluating the hot tail runaway electron seed that emerges during the thermal quench. This problem is made challenging due to the requirement of describing a strongly non-thermal electron distribution, together with the need to incorporate a diverse range of multiphysics processes, including magnetohydrodynamic instabilities, impurity transport, and radiative losses. Here this work develops a physics-informed neural network (PINN) tailored to the solution of the hot tail seed during an axisymmetric thermal quench. Here, a PINN is developed to identify solutions to the adjoint relativistic Fokker–Planck equation in the presence of a rapid quench of the plasma's thermal energy. It is shown that the PINN is able to accurately predict the hot tail seed across a range of parameters, including the thermal quench timescale, initial plasma temperature, and local current density, in the absence of experimental or simulation data. The hot tail PINN is verified by comparison with direct Monte Carlo simulations, with excellent agreement found across a broad range of thermal quench conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data analytics approach to predict high-temperature cyclic oxidation kinetics of NiCr-based Alloys

Abstract Although of practical importance, there is no established modeling framework to accurately predict high-temperature cyclic oxidation kinetics of multi-component alloys due to the inherent complexity. We present a data analytics approach to predict the oxidation rate constant of NiCr-based alloys as a function of composition and temperature with a highly consistent and well-curated experimental dataset. Two characteristic oxidation models, i.e., a simple parabolic law and a statistical cyclic oxidation model, have been chosen to numerically represent the high-temperature oxidation kinetics of commercial and model NiCr-based alloys. We have successfully trained machine learning (ML) models using highly ranked key input features identified by correlation analysis to accurately predict experimental parabolic rate constants ( k p ). This study demonstrates the potential of ML approaches to predict oxidation kinetics of alloys over wide composition and temperature ranges. This approach can also serve as a basis for introducing more physically meaningful ML input features to predict the comprehensive cyclic oxidation behavior of multi-component high-temperature alloys with proper constraints based on the known underlying mechanisms.

36 MATERIALS SCIENCE↗

Multiscale Land-Atmosphere Interaction in Tropical Ecosystems

A leading source of uncertainty in the magnitude of future climate change arises from feedback between the physical climate system and the carbon cycle. Output from 15 Earth System Models (ESMs) in the Coupled Model Intercomparison Project (CMIP5) shows that when forced by identical emissions, predicted atmospheric CO 2 concentrations in 2100 ranged from 780 ppm to 1120 ppm due to variations in carbon cycle feedback among the models. The resulting difference in radiative forcing of climate is almost 2 W m -2 , comparable to the total uncertainty in the physical climate system for the representative concentration pathway (RCP) 8.5 scenario in 2100. Scaling future carbon sources and sinks to match contemporary atmospheric CO 2 reduced this uncertainty by a factor of five, suggesting that contemporary observations might dramatically constrain future predictions, but their scaling had no mechanistic basis, so it provides little insight on carbon cycle processes.

54 ENVIRONMENTAL SCIENCES↗

Predictability of the early summer surface air temperature over Western South Asia

Variability of the Surface Air Temperature (SAT) over the Western South Asia (WSA) region leads to frequent heatwaves during the early summer (May-June) season. The present study uses the European Centre for Medium-Range Weather Forecast’s fifth-generation seasonal prediction system, SEAS5, from 1981 to 2022 based on April initial conditions (1-month lead) to assess the SAT predictability during early summer season. The goal is to evaluate the SEAS5’s ability to predict the El Niño-Southern Oscillation (ENSO) related interannual variability and predictability of the SAT over WSA, which is mediated through upper-level (200-hPa) geopotential height anomalies. This teleconnection leads to anomalously warm surface conditions over the region during the negative ENSO phase, as observed in the reanalysis and SEAS5. We evaluate SEAS5 prediction skill against two observations and three reanalyses datasets. The SEAS5 SAT prediction skill is higher with high spatial resolution observations and reanalysis datasets compared to the ones with low-resolution. Overall, SEAS5 shows reasonable skill in predicting SAT and its variability over the WSA region. Moreover, the predictability of SAT during La Niña is comparable to El Niño years over the WSA region.

54 ENVIRONMENTAL SCIENCES↗

Kibble–Zurek mechanism of Ising domains

The formation of topological defects after a symmetry-breaking phase transition is an overarching phenomenon that encodes the underlying dynamics. The Kibble-Zurek mechanism (KZM) describes these non-equilibrium dynamics of second-order phase transitions and predicts a power-law relationship between the cooling rates and the density of topological defects. It has been verified as a successful model in a wide variety of physical systems, including structure formation in the early Universe and condensed-matter materials. However, it is uncertain if the KZM mechanism is also valid for topologically trivial Ising domains, one of the most common and fundamental types of domain in condensed-matter systems. Here we show that the cooling rate dependence of Ising domain density follows the KZM power law in two different three-dimensional structural Ising domains: ferro-rotation domains in NiTiO 3 and polar domains in BiTeI. However, although the KZM slope of NiTiO 3 agrees with the prediction of the 3D Ising model, the KZM slope of BiTeI exceeds the theoretical limit, providing an example of steepening KZM slope with long-range dipolar interactions. Finally, our results demonstrate the validity of KZM for Ising domains and reveal an enhancement of the power-law exponent for transitions of non-topological quantities with long-range interactions. The Kibble-Zurek mechanism is shown to apply to structural Ising domains in three-dimensional materials. Long-range interactions modify the critical exponents away from theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cross-scale covariance for material property prediction

A simulation can stand its ground against an experiment only if its prediction uncertainty is known. The unknown accuracy of interatomic potentials (IPs) is a major source of prediction uncertainty, severely limiting the use of large-scale classical atomistic simulations in a wide range of scientific and engineering applications. Here we explore covariance between predictions of metal plasticity, from 178 large-scale (~10 8 atoms) molecular dynamics (MD) simulations, and a variety of indicator properties computed at small-scales (≤10 2 atoms). All simulations use the same 178 IPs. In a manner similar to statistical studies in public health, we analyze correlations of strength with indicators, identify the best predictor properties, and build a cross-scale “strength-on-predictors” regression model. This model is then used to estimate regression error over the statistical pool of IPs. Small-scale predictors found to be highly covariant with strength are computed using expensive quantum-accurate calculations and used to predict flow strength, within the statistical error bounds established in our study.

36 MATERIALS SCIENCE↗

Accurate and efficient predictions of keyhole dynamics in laser materials processing using machine learning-aided simulations

The keyhole phenomenon has been widely observed in laser materials processing, including laser welding, remelting, cladding, drilling, and additive manufacturing. Keyhole-induced defects, primarily pores, dramatically affect the performance of final products, impeding the broad use of these laser-based technologies. The formation of these pores is typically associated with the dynamic behavior of the keyhole. So far, the accurate characterization and prediction of keyhole features, particularly keyhole depth, as a function of time, has been a challenging task. In situ characterization of keyhole dynamic behavior using the synchrotron X-ray technique is informative but complicated and expensive. Current simulations are generally hindered by their poor accuracy and generalization abilities in predicting keyhole depths due to the lack of accurate laser absorptance data. In this study, we develop a machine learning-aided simulation method that accurately predicts keyhole dynamics, especially in keyhole depth fluctuations, over a wide range of processing parameters. In two case studies involving titanium and aluminum alloys, we achieve keyhole depth prediction with a mean absolute percentage error of 10 %, surpassing those simulated using the ray-tracing method with an error margin of 30 %, while also reducing computational time. This exceptional fidelity and efficiency empower our model to serve as a cost-effective alternative to synchrotron experiments. Our machine learning-aided simulation method is affordable and readily deployable for a large variety of materials, opening new doors to eliminate or reduce defects for a wide range of laser materials processing techniques.

Computational fluid dynamics↗

Collective dynamics and long-range order in thermal neuristor networks

Abstract In the pursuit of scalable and energy-efficient neuromorphic devices, recent research has unveiled a novel category of spiking oscillators, termed “thermal neuristors.” These devices function via thermal interactions among neighboring vanadium dioxide resistive memories, emulating biological neuronal behavior. Here, we show that the collective dynamical behavior of networks of these neurons showcases a rich phase structure, tunable by adjusting the thermal coupling and input voltage. Notably, we identify phases exhibiting long-range order that, however, does not arise from criticality, but rather from the time non-local response of the system. In addition, we show that these thermal neuristor arrays achieve high accuracy in image recognition and time series prediction through reservoir computing, without leveraging long-range order. Our findings highlight a crucial aspect of neuromorphic computing with possible implications on the functioning of the brain: criticality may not be necessary for the efficient performance of neuromorphic systems in certain computational tasks.

Science & Technology - Other Topics↗

A Nonergodic Ground-Motion Model for the San Francisco Bay Area for Small-Magnitude Earthquakes

ABSTRACT Recently, generative models have become a computationally efficient alternative to physics-based numerical simulations of ground motions. Neural networks can learn from existing ground-motion data to generate unobserved ground-motion data at new source and site locations. A key challenge with generative models is ensuring that predicted ground motions remain within a physically realistic range. For this purpose, we developed an empirical, nonergodic ground-motion model (GMM) for small-magnitude earthquakes in the San Francisco Bay area based on about 5000 recordings per component for Mw ≤ 4 earthquakes. The nonergodic GMM predicts spatially varying median source, site, and path effects for both the Fourier amplitude spectrum (FAS) and the Fourier phase derivative (a proxy for duration), as well as the corresponding epistemic uncertainty for each term. For FAS, our model shows above-average source and site effects in the western part of the region and below-average effects in the eastern part, with regional effects exhibiting larger spatial correlation lengths with increasing frequency. For duration, the source term is negligible for small-magnitude earthquakes, and the site term leads to site-specific variations up to 5 s. Path effects for FAS and duration depend on the source–site pair and are extrapolated spatially using recent methods for path-effect modeling. The aleatory variability of the within-site within-path residuals is similar to the variability found in previous studies for other regions. The nonergodic model provides two key contributions: first, median adjustment terms that are transferable to larger magnitude earthquakes, further reducing aleatory variability in probabilistic seismic hazard analysis; second, region-specific criteria for validating machine learning-based ground-motion generators to evaluate whether synthetic ground motions exhibit physically realistic source, site, and path effects.

Lacour, Maxime↗

Deep learning-based model for progress variable dissipation rate in turbulent premixed flames

A deep neural network (DNN) based large eddy simulation (LES) model for progress variable dissipation rate in turbulent premixed flames is presented. The DNN model is trained using filtered data from direct numerical simulations (DNS) of statistically planar turbulent premixed flames with n-heptane as fuel. Training data was comprised of flames with varying turbulence levels leading to a range of Karlovitz numbers. Through a-priori tests the DNN model is shown to predict the subfilter contribution to progress variable dissipation rate accurately over a range of filter widths and for all Karlovitz numbers examined in this study. Superior performance of the DNN model relative to an established physics-based model is also demonstrated. Additionally, transferability of the DNN model is highlighted by a-priori evaluation of the model using filtered DNS data from multiple cases with different Karlovitz numbers and fuel species than those that were used for training the model.

33 ADVANCED PROPULSION SYSTEMS↗

Comparative study of Minnesota functionals performance on ferroelectric BaTiO 3 and PbTiO 3

Density functional theory based simulations are the leading tool for the computational investigation of ferroelectrics and the parametrization of their classical potentials. However, the predictions often depend strongly on the exchange correlation functional. The most popular choices, LDA and GGA, tend to underestimate or overestimate some structural, electric, and energy properties. These impede development of highly accurate classical potentials that extend the reach of first-principles simulations to finite temperatures and realistic sizes. In this work, we investigate the performance of recently developed Minnesota exchange correlation functionals on prototypical ferroelectrics BaTiO 3 , and PbTiO 3 , in comparison with some popular ones. We find that there exists strong correlation between predictions for some properties (tetragonality, phase energy difference, and polarization) by different functionals. Along the correlation line, we find a range of functionals (including some from the Minnesota suite) whose predictions fall in between those of LDA and GGA and, therefore, offer an improvement. A way to relatively rank functional performance with respect to chosen benchmarks is proposed and applied to identify the top performers for BaTiO 3 and PbTiO 3 . Here, the performance is found to be material dependent. Therefore, we propose that the performance assessment carried out in this work is employed to "screen" functionals prior to their use on ferroelectrics. Quick and computationally inexpensive, it is likely to lead to improved descriptions, especially for classical potential parametrization. The Minnesota functional suite is found to be suitable for this task.

36 MATERIALS SCIENCE↗

Parameterized anomalous transport model for current-carrying collisionless plasmas in pulsed power inertial confinement fusion

Current delivery in pulsed power inertial confinement fusion is influenced by collisionless current-carrying microturbulent plasmas, which are sourced from electrode surfaces. In this setting, the lower hybrid drift instability—triggered by plasma acceleration—is a leading candidate driver of difficult-to-predict momentum and energy transport. To characterize the nonlinear state of the microturbulent plasma, a parameterized anomalous transport model is developed for the instability, with analytic formulas for anomalous collision frequency, resistivity, and species heating rates. The formulas are expressed in terms of linear-theory variables and four dimensionless parameters that characterize the macroscopic plasma state. The model is built on linear theory analysis, power law analysis, and quasilinear theory analysis, and is validated using a series of nonlinear continuum kinetic Vlasov–Poisson simulations. The theoretical and computational investigation demonstrates that the anomalous collision frequency associated with the instability can be reliably approximated, within about a factor of five or better, by the unscaled linear theory growth rate of the fastest-growing wavenumber mode. This finding enables efficient calculation of anomalous resistivity and species heating rates over a wide range of plasma conditions, resulting in improved predictive capabilities.

Complex functions↗

A Universal Augmentation Framework for Long-Range Electrostatics in Machine Learning Interatomic Potentials

Most current machine learning interatomic potentials (MLIPs) rely on short-range approximations, without explicit treatment of long-range electrostatics. To address this, we recently developed the Latent Ewald Summation (LES) method, which infers electrostatic interactions, polarization, and Born effective charges (BECs), just by learning from energy and force training data. Here, in this study, we present LES as a standalone library, compatible with any short-range MLIP, and demonstrate its integration with methods such as MACE, NequIP, Allegro, CACE, CHGNet, and UMA. We benchmark LES-enhanced models on distinct systems, including bulk water, polar dipeptides, and gold dimer adsorption on defective substrates, and show that LES not only captures correct electrostatics but also improves accuracy. Additionally, we scale LES to large and chemically diverse data by training MACELES-OFF on the SPICE set containing molecules and clusters, making a universal MLIP with electrostatics for organic systems, including biomolecules. MACELES-OFF is more accurate than its short-range counterpart (MACE-OFF) trained on the same data set, predicts dipoles and BECs reliably, and has better descriptions of bulk liquids. By enabling efficient long-range electrostatics without directly training on electrical properties, LES paves the way for electrostatic foundation MLIPs.

Kim, Dongjin [University of California, Berkeley, ↗

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of the 230 Th( p ,2n)Pa229 and 230 Th( p ,3n)Pa228 reaction cross sections from 14.1 to 16.9 MeV

Actinium-225 is of interest for medical isotope production and there is on-going research into methods of producing Ac 225 , either directly or via the decay of its parent isotopes ( Th 229 , Pa 229 , and Ra 225 ). One method that has been suggested is the Th 230 ( p , 2 n ) Pa 229 reaction. However, there is no available cross-section data for this reaction in the literature. Purpose: Measure the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reaction cross sections in the energy range where the ( p , 2 n ) reaction is predicted to peak to determine the feasibility of Ac 225 production via the Th 230 ( p , 2 n ) reaction. Methods: Targets naturally enriched in Th 230 were irradiated at the Center for Accelerator Mass Spectrometry at Lawrence Livermore National Laboratory with energies ranging from 14.1 to 16.9 MeV. Furthermore, chemical processing was used to separate the protactinium activation products, followed by γ -ray spectroscopy to measure the activities of Pa 228 , 229 , 230 , 232 produced in the irradiation. Results: We find that excitation functions are reported for the first time in the literature for the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reactions in this energy range. The peak measured value of the Th 230 ( p , 2 n ) reaction was found to be 182 ± 12 mb at 14.4 ± 0.1 MeV. The Th 232 ( p , n ) Pa 232 reaction was used to verify the experimental conditions, the measured values are reported and are comparable to the existing literature values. From the γ -ray spectrometry data, the half-life of Pa 229 was measured as 1.5 ± 0.1 days, which is within the error of the half-life reported in the evaluated nuclear data as well as in the recent measurements, and the half-life of Pa 228 was measured as 19.5 ± 0.4 hours. Conclusions: Overall, the Th 230 ( p , 2 n ) Pa 229 reaction could reasonably be used for Ac 225 isotope production, although significant amounts of relatively isotopically pure Th 230 would be needed for significant production because the low alpha-decay branching ratio of Pa 229 and long half-life of Th 229 inhibit the in-growth of significant amounts of Ac 225 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Field validation of isotropic analytical models for simulating fabric shades

Fabric roller shades are common shading materials used in commercial and residential buildings. Accurately characterizing and modeling shades helps practitioners select the appropriate product and its control strategy based on climate and occupants' priorities, such as visual comfort and view to outdoors. Previous studies established a generalized method for modeling complex fenestration systems using data-driven tabulated bidirectional scattering distribution functions. However, deploying such a method at scale to all fabric shading products on the market is too costly and time-consuming. Analytical models that are based on a limited set of measurements (e.g., normal-normal and normal-hemispherical visible transmittance and reflectance, and directional cut-off angles) can be used to model the wide variety of shading products on the market. This study evaluates the performance of two isotropic analytical models, Roos-Wienold and Modified-Kotey, for modeling fabric roller shades, with a focus on the model's ability to predict occupant visual comfort. The performance evaluation was conducted through laboratory and field measurements and simulations. The results showed that both models are sufficient for predicting vertical illuminance at seated eye-level. Roos-Wienold model was able to predict binary visual comfort classification (glare/no-glare) under a wide range of luminance conditions, while Modified-Kotey model did not perform as well under high-contrast low-adaptation conditions. Both models are insufficient in predicting visual comfort at a four-point scale (e.g., imperceptible, perceptible, disturbing, intolerable). The two isotropic models become less accurate when the fabric exhibits high anisotropy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

EVALUATION OF OPERATION TEMPERATURES UNDER NATURAL CONVECTION HELIUM FLOW IN A CONFINED CAVITY

The Spallation Neutron Source (SNS) is a high-power accelerator-based pulsed neutron source led by Oak Ridge National Laboratory (ORNL) to achieve high fluxes of neutrons for scientific experiments. Active and passive cooling of the systems and parts forming the SNS have been considered to warrant the safe operation of the facility. The diverse cooling systems make use of conjugated heat transfer mechanisms to provide a stable operation temperature for all components in the machine. Thermal power deposited into stainless-steel piping lines due to particle radiation may reach values of up to 1.2 W/cc in the regions located closer to the center of the lower IRP. These energy deposition levels, in not actively cooled components, such as the transfer-line-outer-vacuum-layer may increase the temperature of the component beyond design requirement limitations. The evaluation of the operation temperatures for the former components relies in the assumption that a low-pressure helium atmosphere provides enough heat removal capacity based on natural convection phenomena. In this work, the evaluation of steady state temperatures in components such as the CMS transfer lines has been evaluated using computational fluid dynamics (CFD), analytical correlations and experimental measurements. The companion experiments were conducted in a closed helium system at pressures varying from 1.1 to 1.5 bar. A copper rod was affixed horizontally between viewing windows and heated at constant power, and measurements were made of both the rod temperature and ambient temperature via a system fiberoptic distributed temperature sensors and RTDs. It was found that the measured heat transfer coefficients agree well with the predictions of Churchill and Chu correlations across the range of cases considered. Additionally, the ambient helium volume above the rod was imaged via background oriented schlieren (BOS), and these data was used to determine the line-averaged density gradients in this region. These gradients were compared to simulation data to validate the predictions of natural convection simulations.

Dominguez-Ontiveros, Elvis [ORNL] (ORCID:000000018↗