SEARCH · Engineering Papers

Results for “simulation analytics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗

DOE OSTI2025

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Precision measurements of EFT parameters and BAO peak shifts for the Lyman- α forest

We present precision measurements of the bias parameters of the one-loop power spectrum model of the Lyman- α (Ly- α ) forest, derived within the effective field theory (EFT) of large-scale structure. We fit our model to the three-dimensional flux power spectrum measured from the ACCEL 2 hydrodynamic simulations. The EFT model fits the data with an accuracy of below 2% up to k = 2 h Mpc − 1 . Further, we analytically derive how nonlinearities in the three-dimensional clustering of the Ly- α forest introduce biases in measurements of the baryon acoustic oscillations (BAOs) scaling parameters in radial and transverse directions. From our EFT parameter measurements, we obtain a theoretical error budget of Δ α ∥ = − 0.2 % ( Δ α ⊥ = − 0.3 % ) for the radial (transverse) parameters at redshift z = 2.0 . This corresponds to a shift of − 0.3 % (0.1%) for the isotropic (anisotropic) distance measurements. We provide an estimate for the shift of the BAO peak for Ly- α -quasar cross-correlation measurements assuming analytical and simulation-based scaling relations for the nonlinear quasar bias parameters resulting in a shift of − 0.2 % ( − 0.1 % ) for the radial (transverse) dilation parameters, respectively. This analysis emphasizes the robustness of Ly- α forest BAO measurements to the theory modeling. We provide informative priors and an error budget for measuring the BAO feature—a key science driver of the currently observing Dark Energy Spectroscopic Instrument (DESI). Our work paves the way for full-shape cosmological analyses of Ly- α forest data from DESI and upcoming surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST. Published by the American Physical Society 2025

de Belsunce, Roger (ORCID:0000000336604028)↗

DOE OSTI2024

A Visual Analytic Platform for Interactive Validation of Human Mobility Simulations

Human mobility insights guide domain experts in an array of decisions, including critical infrastructure design, disaster response, epidemic modeling, national security, and policy making. Due to the inherent noise and privacy concerns in real-world individual-level mobility data, it is often preferred to leverage simulators that generate synthetic mobility data instead. However, it is critical to inspect and validate the output of such simulators to ensure the synthetic data is aligned with the characteristics of the population and the area of interest known to domain experts. While there exist many quantitative approaches for validating synthetic data, we argue it is also important to also validate such data qualitatively to capture aspects that are known to domain experts but difficult to quantify. In this work, we demonstrate a visual analytic platform that empowers domain experts to interact with their simulation outputs along spatial and temporal dimensions. By augmenting automated techniques and human skills, our visual analytic platform is a step towards interactive capabilities for model steering and quality control of mobility simulators.

Monadjemi, Shayan↗

DOE OSTI2022

Estimation of the time for steam generator trip due to cyber intrusions

The time required to trip a pressurized water reactor (PWR) by inserting malicious signals into its steam generator (SG) control system has been studied using the Generic PWR (GPWR) Simulator. A semi-analytical model is developed to approximately reproduce the simulator response and understand the dynamics of the control unit. A series of two proportional-integral controllers determines control action according to preset constants, the readings from the feedwater level sensor, and those from feedwater and steam flowrate transmitters. It is observed that the most important factor that determines whether a trip will occur is how much additional water is added to or withheld from the SG over time compared to normal operating conditions. In order to determine the effects of control action on the SG, changes in mass inventory are considered. This approach models the SG water level as a function of mass inventory and has a backward temporal memory. A Python interface is developed for the GPWR framework to automatically simulate different spoofing scenarios and post-process the related data. We observe that the trip times predominantly depend on flow mismatch and/or level errors. Controller parameters, including the integral time and gain constants, either speed up or slow down the rate of progression to a trip setpoint but do not cause a trip by themselves. The reactor can trip on a high-level signal when the reading crosses above 78%, increased from its reference level of 57%, or a low-level reading when it is below 25%. The present results show roughly how long the operators would have to respond to an attack, given a specific set of spoofing signals within the issue space analyzed. Furthermore, we have generated a simple surface by fitting a combination of exponential functions to the data obtained from the GPWR Simulator. In general, trips on a low level have been observed to occur faster than those on a high level.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1990

Computer simulation of gear tooth manufacturing processes

The use of computer graphics to simulate gear tooth manufacturing procedures is discussed. An analytical basis for the simulation is established for spur gears. The simulation itself, however, is developed not only for spur gears, but for straight bevel gears as well. The applications of the developed procedure extend from the development of finite element models of heretofore intractable geometrical forms, to exploring the fabrication of nonstandard tooth forms.

Mavriplis, Dimitri↗

NASA NTRS2012

Stagnation Point Radiative Heating Relations for Venus Entry

Improved analytic expressions for calculating the stagnation point radiative heating during entry into the atmosphere of Venus have been developed. These analytic expressions can be incorporated into entry trajectory simulation codes. Together with analytical expressions for convective heating at the stagnation point, the time-integrated total heat load at the stagnation point is used in determining the thickness of protective material required, and hence the mass of the fore body heatshield of uniform thickness.

Tauber, Michael E.↗

NASA NTRS1987

Integrated control and display research for transition and vertical flight on the NASA V/STOL Research Aircraft (VSRA)

Results of a substantial body of ground-based simulation experiments indicate that a high degree of precision of operation for recovery aboard small ships in heavy seas and low visibility with acceptable levels of effort by the pilot can be achieved by integrating the aircraft flight and propulsion controls. The availability of digital fly-by-wire controls makes it feasible to implement an integrated control design to achieve and demonstrate in flight the operational benefits promised by the simulation experience. It remains to validate these systems concepts in flight to establish their value for advanced short takeoff vertical landing (STOVL) aircraft designs. This paper summarizes analytical studies and simulation experiments which provide a basis for the flight research program that will develop and validate critical technologies for advanced STOVL aircraft through the development and evaluation of advanced, integrated control and display concepts, and lays out the plan for the flight program that will be conducted on NASA's V/STOL Research Aircraft (VSRA).

Foster, John D.↗

NASA NTRS1988

Integrated control and display research for transition and vertical flight on the NASA V/STOL Research Aircraft (VSRA)

Foster, John D.↗

NASA NTRS1974

Flight investigation of approach and flare from simulated breakout altitude of a subsonic jet transport and comparison with analytical models

Satisfactory and optimum flare windows are defined from pilot ratings and comments. Maximum flare normal accelerations, touchdown rates of sink, and total landing maneuver time increments are summarized as a function of approach airspeed margin (with respect to reference airspeed) and flare initiation altitude. The effects of two thrust management techniques are investigated. Comparisons are made with predictions from three analytical models and the results of a simulator study. The approach speed margin was found to have a greater influence on the flare initiation altitude than the absolute airspeed. The optimum airspeed was between the reference airspeed and the reference airspeed plus 10 knots. The optimum flare initiation altitude range for unrestricted landings was from 11 meters to 20 meters (36 feet to 66 feet), and the landing time in the optimum window was 8 seconds. The duration of the landing maneuver increased with increasing flare initiation altitude and with increasing speed margins on the approach.

Matheny, N. W.↗

DOE OSTI2022

Analytical and Numerical Analysis of Lithium Plating Onset in Single and Bilayer Graphite Electrodes during Fast Charging

Lithium plating is commonly observed in anodes charged at fast rates, and can lead to capacity loss and battery safety issues. The increased risk of plating has been attributed to transport limitations, and architectured electrodes may reduce plating risk. However, while theoretical studies have shown that reaction non-uniformity arises due to interplay of transport limitations, anode open circuit voltage behavior and reaction kinetics, its effect on lithium plating has not been studied. We use analytic and numerical simulations to predict onset of plating in graphite anode half-cells at high C-rates and demonstrate how anodes with layered porosities can delay plating. Simplified analytical models identify trends for plating onset and predictions are calibrated against numerical models. A calibrated numerical model of graphite demonstrates qualitative agreement with analytical model predictions. This reaction inhomogeneity mechanism occurs in the absence of lithium ion depletion, indicating that these mechanisms may contribute to capacity loss independently or simultaneously. A bilayer model of graphite exhibits delayed plating onset, and an optimization procedure is presented. This theoretical work presents quantitative and mechanistic insight on how reaction inhomogenity affects lithium metal plating onset and can be used as a guide to engineer anodes resistant to lithium plating.

Mijailovic, Aleksandar S.↗

DOE OSTI2021

A Data-Driven Global Sensitivity Analysis Framework for Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.

14 SOLAR ENERGY↗

DOE OSTI2020

Back-to-back relative-excess observable to identify the chiral magnetic effect

Background: The chiral magnetic effect (CME) is extensively studied in heavy-ion collisions at RHIC and LHC. In the commonly used reaction-plane (RP)-dependent, charge-dependent azimuthal correlator ($Δ_γ$), both the close and back-to-back pairs are included. Many backgrounds contribute to the close pairs (e.g., resonance decays, jet correlations), whereas the back-to-back pairs are relatively free of those backgrounds. Purpose: In order to reduce those backgrounds, we propose a new observable, which only focuses on the back-to-back pairs, namely, the relative back-to-back opposite-sign over same-sign pair excess ($r_{BB}$) as a function of the pair azimuthal orientation with respect to the RP ($φ_{BB}$). Methods: We use analytical calculations and toy model simulations to demonstrate the sensitivity of $r_{BB}(φ_{BB})$ to the CME and its insensitivity to backgrounds. Results: With finite CME, the $φ_{BB}$ distribution of $r_{BB}$ shows a clear characteristic modulation. Its sensitivity to background is significantly reduced compared to the previous $Δ_γ$ observable. The simulation results are consistent with our analytical calculations. Conclusions: In this work, our studies demonstrate that the $r_{BB}(φ_{BB})$ observable is sensitive to the CME signal and rather insensitive to the resonance backgrounds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1991

Technical evaluation report, AGARD Fluid Dynamics Panel Symposium on Effects of Adverse Weather on Aerodynamics

The purpose of the meeting on Effects of Adverse Weather on Aerodynamics was to provide an update of the stae-of-the-art with respect to the prediction, simulation, and measurement of the effects of icing, anti-icing fluids, and various precipitation on the aerodynamic characteristics of flight vehicles. Sessions were devoted to introductory and survey papers and icing certification issues, to analytical and experimental simulation of ice frost contamination and its effects of aerodynamics, and to the effects of heavy rain and deicing/anti-icing fluids.

Reinmann, J. J.↗

DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Simulation modeling for long duration spacecraft control systems

The use of simulation is described and it is contrasted to analytical solution techniques for evaluation of analytical reliability models. The role importance sampling plays in simulation of models of this type was also discussed. The simulator tool used for our analysis is described. Finally, the use of the simulator tool was demonstrated by applying it to evaluate the reliability of a fault tolerant hypercube multiprocessor intended for spacecraft designed for long duration missions. The reliability analysis was used to highlight the advantages and disadvantages offered by simulation over analytical solution of Markovian and non-Markovian reliability models.

Boyd, Mark A.↗

NASA NTRSDate not supplied

Early Progress Toward the Feasibility of the Centrifugal Nuclear Thermal Rocket

The Centrifugal Nuclear Thermal Rocket (CNTR) is a Nuclear Thermal Propulsion (NTP) concept designed to heat propellant directly by the reactor fuel. The primary difference between the CNTR concept and traditional NTP systems is that rather than using traditional solid fuel elements, the CNTR uses liquid fuel with the liquid contained in rotating cylinders by centrifugal force. If the concept can be successfully realized, the CNTR would have a high specific impulse (~1800 seconds) at high thrust, which may enable (i) viable near-term human Mars exploration by reducing round-trip times to 420 days and (ii) direct injection orbits for scientific missions to the Solar System outer planets and potentially Kuiper Belt objects. The CNTR could also use storable propellants such as ammonia, methane, or hydrazine at an Isp of ~900 seconds, enabling long-term in-space storage of a dormant system. Significant engineering challenges must be addressed to establish the technical viability of the CNTR. Research is presently underway to determine resolutions for these engineering challenges. In particular, research has begun on the analytical modeling and simulation of the two-phase heat transfer between the liquid metallic uranium fuel and the gaseous propellant. A paper was presented at the 2021 IAC which described these challenges and the study plan to address them. This paper describes the analytical and experimental progress to date toward resolving these challenges and establishing the engineering feasibility of the CNTR technology

High Performance Nuclear Thermal Propulsion NTP SN↗

DOE OSTI2024

Including Neutrino-driven Convection in the Force Explosion Condition to Predict Explodability of Multidimensional Core-collapse Supernovae (FEC+)

Abstract Most massive stars end their lives with core collapse. However, it is not clear which explode as a core-collapse supernova (CCSN), leaving behind a neutron star, and which collapse to a black hole, aborting the explosion. One path to predict explodability without expensive multidimensional simulations is to develop analytic explosion conditions. These analytic conditions also provide a deeper understanding of the explosion mechanism and they provide some insight into why some simulations explode and some do not. The analytic force explosion condition (FEC) reproduces the explosion conditions of spherically symmetric CCSN simulations. In this follow-up manuscript, we include the dominant multidimensional effect that aids explosion—neutrino-driven convection—in the FEC. This generalized critical condition (FEC+) is suitable for multidimensional simulations and has potential to accurately predict explosion conditions of two- and three-dimensional CCSN simulations. We show that adding neutrino-driven convection reduces the critical condition by ∼30%, which is consistent with previous multidimensional simulations.

Gogilashvili, Mariam (ORCID:0000000269448052)↗