Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical simulation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2025

Direct Numerical Simulation of Flow Past a Gaussian Bump at a High Reynolds Number

We present the results from a direct numerical simulation of a spanwise-periodic turbulent flow past a Gaussian bump. The problem setup is designed to investigate the interaction of an incoming turbulent boundary layer with the strong favorable and adverse pressure gradients generated by the Gaussian bump as the flow passes over it at a Reynolds number of 340000 based on the bump height, or 4 million based on the bump length. The statistical results from the present simulation are compared against our earlier results at a Reynolds number of 2 million. An internal layer, which forms beneath the strongly accelerated boundary layer over the windward side of the bump, is found to generate its near-wall turbulence stress peaks in closer proximity of the wall in the higher Reynolds-number case. Furthermore, the logarithmic layer of the higher Reynolds-number boundary layer appears more resistant to changes induced by strong acceleration and surface curvature effects over the same region. Despite a nearly identical flow separation point in the two flows, the detached shear layer grows at a faster rate and subsequently reattaches at an earlier point in the higher Reynolds-number flow. The surface pressure and skin-friction distributions over the attached flow region compare well against the corresponding experimental data for both flows. However, some differences appear in the separated flow region, which are attributed to the three-dimensionality of the experimental model setup that is not included in the simulation owing to the spanwise periodic assumption. Comparisons with the stereoscopic particle image velocimetry measurements on the central plane of the experimental model over the windward side of the bump show reasonable overall agreement in the mean velocity components, but the turbulence stress components do not agree well at some streamwise locations. Comparisons over the leeward side of the bump show that the mean separated shear layer in the simulation is tilted significantly more toward the wall than the experimental shear layer on the central plane. This mismatch in the mean shear layer orientation is due to the experimental model three-dimensionality and tunnel end-wall effects, which are not modeled in the present spanwise-periodic simulation.

Favorable and Adverse Pressure Gradients↗
DOE OSTI2025

MPACT Safeguards Modeling: FY25 Update

Sandia National Laboratories develops and maintains several open-source software packages to support material accountancy analyses. This includes the Material Accountancy Performance Indicator Toolkit (MAPIT), the Fissile Facility Flow Modeler (F3M) and the Separation and Safeguards Performance Model Library (SSPM-L). MAPIT is responsible for performing statistical safeguards analyses on bulk and itemized data from nuclear fuel cycle facilities and can operate on real or synthetic data. MAPIT is the only open-source software for such analyses. F3M is a library of modules, built in MATLAB Simulink, that contain pre made blocks to represent different generic fuel cycle processes. These blocks can be used together in a modular fashion to represent and simulate nuclear fuel cycle processes with the goal of improving facility-level accountancy during the design phase. F3M is also an open-source library. Finally, the SSPM-L library is a series of completed models built from F3M. The library includes facility models such as a generic PUREX facility and a fuel fabrication facility. The SSPM-L library is not open source, but is available to collaborators with a relevant use case. These tools include modeling and simulation pipelines to simulate nuclear fuel cycle facilities and the underlying software needed to simulate measurement uncertainty and perform statistical analyses. Together, these tools can perform end-to-end nuclear material accountancy analyses. This report documents the various improvements made to these tools in FY25. Specifically, we added new statistical test, new statistical modeling capabilities, new fuel cycle facility models, and launched a new open-source model component library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimizing LSST observing strategy for weak lensing systematics

ABSTRACT The Legacy Survey of Space and Time (LSST) survey will provide unprecedented statistical power for measurements of dark energy. Consequently, controlling systematic uncertainties is becoming more important than ever. The LSST observing strategy will affect the statistical uncertainty and systematics control for many science cases; here, we focus on weak lensing (WL) systematics. The fact that the LSST observing strategy involves hundreds of visits to the same sky area provides new opportunities for systematics mitigation. We explore these opportunities by testing how different dithering strategies (pointing offsets and rotational angle of the camera in different exposures) affect additive WL shear systematics on a baseline operational simulation, using the ρ-statistics formalism. Some dithering strategies improve systematics control at the end of the survey by a factor of up to ∼3–4 better than others. We find that a random translational dithering strategy, applied with random rotational dithering at every filter change, is the most effective of those strategies tested in this work at averaging down systematics. Adopting this dithering algorithm, we explore the effect of varying the area of the survey footprint, exposure time, number of exposures in a visit, and exposure to the Galactic plane. We find that any change that increases the average number of exposures (in filters relevant to WL) reduces the additive shear systematics. Some ways to achieve this increase may not be favorable for the WL statistical constraining power or for other probes, and we explore the relative trade-offs between these options given constraints on the overall survey parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

A Fast-Time Study of Aircraft Reordering in Arrival Sequencing and Scheduling

In order to ensure that the safe capacity of the terminal area is not exceeded, Air Traffic Management ATM often places restrictions on arriving flights transitioning from en route airspace to terminal airspace. This restriction of arrival traffic is commonly referred to as arrival flow management, and includes techniques such as metering, vectoring, fix-load balancing, and the imposition of miles-in-trail separations. These restrictions are enacted without regard for the relative priority which airlines may be placing on individual flights based on factors such as crew criticality, passenger connectivity, critical turn times, gate availability, on-time performance, fuel status, or runway preference. The development of new arrival flow management techniques which take into consideration priorities expressed by air carriers will likely reduce the economic impact of ATM restrictions on the airlines and lead to increased airline economic efficiency by allowing airlines to have greater control over their individual arrival banks of aircraft. NASA and the Federal Aviation Administration (FAA) have designed and developed a suite of software decision support tools (DSTs) collectively known as the Center TRACON Automation System (CTAS). One of these tools, the Traffic Management Advisor (TMA) is currently being used at the Fort Worth Air Route Traffic Control Center to perform arrival flow management of traffic into the Dallas/Fort Worth airport (DFW). The TMA is a time-based strategic planning tool that assists Traffic Management Coordinators (TMCs) and En Route Air Traffic Controllers in efficiently balancing arrival demand with airport capacity. The primary algorithm in the TMA is a real-time scheduler which generates efficient landing sequences and landing times for arrivals within about 200 no a. from touchdown. This scheduler will sequence aircraft so that they arrive in a first- come - first-served (FCFS) order. While FCFS sequencing establishes a fair order based on estimated times of arrival, it does not take into account individual airline priorities among incoming flights. NASA is exploring the possibility of allowing airlines to express relative arrival priorities to air traffic management through the development of new CTAS scheduling algorithms which take into consideration airline arrival preferences. The accommodation of airline priorities in arrival sequencing and scheduling would under most circumstances result in a deviation from a "natural" or FCFS arrival order. As a First step toward developing airline influenced sequencing algorithms, an investigation was conducted to determine the feasibility of reordering arrival traffic from a strict FCFS sequence. A fast-time simulation has been developed which allows statistical evaluation of sequencing and scheduling algorithms for arrival traffic at the Dallas/Fort Worth Airport. In contrast to real-time simulation or field tests, which would require on the order of ninety minutes to examine a single traffic rush period, the fast-time simulation allows examination of multiple rush periods in a matter of seconds.

Carr, Greg↗
NASA NTRSDate not supplied

Objective Measurement Assessment of Departure Advisories for Ramp Controllers from a Human-In-The-Loop Simulation

NASA has developed and demonstrated an integrated arrival, departure, and surface concept and technology for efficient air traffic management in busy terminal environments under the Airspace Technology Demonstration 2 (ATD-2) subproject. In April/May 2019, NASA conducted a human-in-the-loop (HITL) simulation with ramp and tower controllers having experience at Dallas/Fort Worth International Airport (DFW). The purpose of this HITL was to evaluate the impacts of various surface metering goals in ramp operations at DFW and test new features of the ATD-2 ramp controller decision support tools. This paper evaluates the quantitative metrics from the simulation results related to airport performance and surface metering given various departure scheduling advisories for ramp controllers. The objective measurements compared include the controller’s compliance with the target times for pushback and spot arrival, the number of metered flights, aircraft gate hold and taxi time, runway throughput, and the number of aircraft on the surface. The simulation results show that there were no statistically significant benefits from surface metering at DFW for airport performance due to several simulation artifacts notably different from real operations. Considering controller workload and situation awareness, following either gate pushback or spot arrival advisory could be a better metering option than utilizing both advisories when surface metering is applied at DFW.

airport ramp traffic control decision support↗
NASA NTRS2014

Reply to Comment by Laprise on 'the Added Value to Global Model Projections of Climate Change by Dynamical Downscaling: a Case Study over the Continental U.S. Using the GISS-ModelE2 and WRF Models'

In his comment, Laprise raises several points that we agree merit consideration. His primary critique is that our study [Racherla et al., 2012] tested the ability of the WRF regional climate model to reproduce historical temperature and precipitation change relative to the driving global climate model (GCM) using only a single simulation rather than an ensemble. He asserts that the observed changes are smaller than the internal variability in the climate system (i.e., not statistically significant) and that thus a single simulation should not necessarily be able to capture the observations. Laprise points out that the statistical signal is reduced for a multi-decadal trend such as the one we analyzed in comparison with mean climatology and cites two studies showing that for particular climate parameters it can take any years for a signal to be discerned over internal variability. He states that The results of theexperiment as designed were strongly influenced by the presence of internal variability and sampling errors,which masked the rather small climate changes that may have occurred as a consequence of changes inforcing during the period considered. While Laprise discusses statistics in general terms at some length, for the actual climate trends examined in our study, he offers no evidence that the forced signal was smallcompared with internal variability. The two studies he cites [de Ela et al., 2013; Maraun, 2013] do not provide convincing evidence as they concern climate variables averaged over different times and areas. One in fact examines extreme precipitation events, which by definition are rare and thus have a lower significance level. We accept the general point that it is important to consider internal variability, and as noted in our paper we agree that an ensemble of simulations is in principle an optimal, though computationally expensive, approach. While we did not present the statistical significance of the observations in our original paper, we have now evaluated those for the regional temperature trends used in our study to evaluate the added value of WRF and thus can analyze data as to the magnitude of the trends with respect to internal variability.

climate↗
NASA NTRS1994

A statistical inference approach for the retrieval of the atmospheric ozone profile from simulated satellite measurements of solar backscattered ultraviolet radiation

NASA's Mission to Planet Earth (MTPE) will address important interdisciplinary and environmental issues such as global warming, ozone depletion, deforestation, acid rain, and the like with its long term satellite observations of the Earth and with its comprehensive Data and Information System. Extensive sets of satellite observations supporting MTPE will be provided by the Earth Observing System (EOS), while more specific process related observations will be provided by smaller Earth Probes. MTPE will use data from ground and airborne scientific investigations to supplement and validate the global observations obtained from satellite imagery, while the EOS satellites will support interdisciplinary research and model development. This is important for understanding the processes that control the global environment and for improving the prediction of events. In this paper we illustrate the potential for powerful artificial intelligence (AI) techniques when used in the analysis of the formidable problems that exist in the NASA Earth Science programs and of those to be encountered in the future MTPE and EOS programs. These techniques, based on the logical and probabilistic reasoning aspects of plausible inference, strongly emphasize the synergetic relation between data and information. As such, they are ideally suited for the analysis of the massive data streams to be provided by both MTPE and EOS. To demonstrate this, we address both the satellite imagery and model enhancement issues for the problem of ozone profile retrieval through a method based on plausible scientific inferencing. Since in the retrieval problem, the atmospheric ozone profile that is consistent with a given set of measured radiances may not be unique, an optimum statistical method is used to estimate a 'best' profile solution from the radiances and from additional a priori information.

Bonavito, N. L.↗
DOE OSTI2025

Sensitivity of an integrated experiment to uncertainty in the high explosive equations of state

Traditionally, hydrodynamics simulations are performed with a single equation-of-state (EOS) to describe each material. These EOSs typically have a physics-informed functional form with adjustable parameters that are calibrated in order to replicate small-scale data. However, because the calibration data have uncertainty and there are typically inherent degeneracies in fitting the EOS, there are actually multiple EOSs that might be consistent with calibration data. In this work, we perform uncertainty quantification (UQ) for the reactant and product equations of state for the high explosive PBX 9501 to yield an ensemble of EOSs that match the uncertain small-scale calibration data. We then simulate an experiment of an explosively formed penetrator repeatedly with different EOSs to both validate the UQ analysis and determine the effects of EOS uncertainty on the prediction of quantities of interest in the experiment. In general, we find good agreement between the simulation predictions and the experimental measurements, and we identify an EOS variable that contributes most directly to the spread in the predictions as the EOSs are varied.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Geometric Sensitivity of Residual Fields in Metal Additive Manufacturing

The sensitivity of geometrical parameters on the residual stress and distortion as computed by a metal additive manufacturing simulation is presented. A nozzle-like geometry was chosen with its length, wall thickness, and outer radius varied. A Latin hypercube sampling of 300 specimens produced the combinations of the geometrical parameters. Full thermomechanical simulations predicted the residual fields. Statistical analysis including regression curves shows that there is a strong correlation between these residual fields and the geometry, hinting that it may be possible to predictably influence the residual displacements and residual stresses by varying the overall geometry.

Multiscale Modeling↗
NASA NTRS1995

Towards Principled Experimental Study of Autonomous Mobile Robots

We review the current state of research in autonomous mobile robots and conclude that there is an inadequate basis for predicting the reliability and behavior of robots operating in unengineered environments. We present a new approach to the study of autonomous mobile robot performance based on formal statistical analysis of independently reproducible experiments conducted on real robots. Simulators serve as models rather than experimental surrogates. We demonstrate three new results: 1) Two commonly used performance metrics (time and distance) are not as well correlated as is often tacitly assumed. 2) The probability distributions of these performance metrics are exponential rather than normal, and 3) a modular, object-oriented simulation accurately predicts the behavior of the real robot in a statistically significant manner.

Object-Oriented Simulation↗
NASA NTRS2003

Simulation of Wind Profile Perturbations for Launch Vehicle Ascent Flight Systems Design Assessments

Ideally, a statistically representative sample of measured high-resolution wind profiles with wavelengths as small as tens of meters is required for assessment of launch vehicle ascent flight systems component capability and vehicle operability for a selected launch site. At most potential launch sites a sample of high-resolution measured wind profiles may not exist. Representative samples of Rawinsonde wind profiles are more likely to be available because of the extensive network of measurement sites established for routine measurements at 12-hr intervals in support of national and international weather observing and forecasting activity. Such a sample, although large enough to statistically represent relatively large wavelength perturbations, would be inadequate for launch system design assessment applications because the Rawinsonde system can accurately measure wind perturbations with wavelengths no smaller than 2000m (1000m altitude increment). Wavelengths less than 2000m in the raw Rawinsonde data, which tend to be dominated by un-damped spurious balloon motion and radar tracking system noise, are filtered within the data processing scheme. The Kennedy Space Center (KSC) Jimsphere wind profiles (150/month and seasonal pairs) are the only adequate high resolution (approximately 150 to 300m effective resolution, but over-sampled at 25m intervals) data that have been used extensively in launch vehicle design, operability assessments and operational protection of vehicle systems for wind perturbation uncertainty. Jimsphere wind profiles have been measured at a few other potential launch sites but the number of profiles is relatively small and the samples are not statistically representative of the site dependent wind profile variability. A simulation process has been developed for enhancement of measured low-resolution Rawinsonde profiles that are more likely to be available in the vicinity of potential launch sites and are a statistically representative sample of wind profile perturbation wavelengths greater than 2000m. The enhancement produces perturbed wind profiles with wavelengths as small as desired for application in launch vehicle ascent flight simulations and design assessments.

Adelfang, S. I.↗
DOE OSTI2022

FlameNEST: explicit profile likelihoods with the Noble Element Simulation Technique

We present FlameNEST, a framework providing explicit likelihood evaluations in noble element particle detectors using data-driven models from the Noble Element Simulation Technique. FlameNEST provides a way to perform statistical analyses on real data with no dependence on large, computationally expensive Monte Carlo simulations by evaluating the likelihood on an event-by-event basis using analytic probability elements convolved together in a single TensorFlow multiplication. Furthermore, this robust framework creates opportunities for simple inter-collaboration analyses which will be fundamental for the future of experimental dark matter physics.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Extensions to models of the galaxy–halo connection

ABSTRACT We explore two widely used empirical models for the galaxy–halo connection, subhalo abundance matching (SHAM) and the halo occupation distribution (HOD), and compare them with the hydrodynamical simulation IllustrisTNG (TNG) for multiple statistics quantifying the galaxy distribution at $n_{\rm gal}\approx 1.3\times 10^{-3}\, ({\rm Mpc}\,h^{-1})^{-3}$. We observe that in their most straightforward implementations, both models fail to reproduce the two-point clustering measured in TNG. We find that SHAM models that use the relaxation velocity, Vrelax, and the peak velocity, Vpeak, perform best, and match the clustering reasonably well, although neither captures adequately the one-halo clustering. Splitting the total sample into sub-populations, we discover that SHAM overpredicts the clustering of high-mass, blue, star-forming, and late-forming galaxies and underpredicts that of low-mass, red, quiescent, and early-forming galaxies. We also study various baryonic effects, finding that subhaloes in the dark-matter-only simulation have consistently higher values of their SHAM-proxy properties than their full-physics counterparts. We then consider a 2D implementation of the HOD model augmented with a secondary parameter (environment, velocity anisotropy, σ2Rhalf-mass, and total potential) tuned so as to match the two-point clustering of the IllustrisTNG galaxies on large scales. We analyse these galaxy populations adopting alternative statistical tools such as galaxy–galaxy lensing, void–galaxy cross-correlations, and cumulants of the density field, finding that the hydrodynamical galaxy distribution disfavours σ2Rhalf-mass and the total potential as secondary parameters, while the environment and velocity anisotropy samples are consistent with full physics across all statistical probes examined. Our results demonstrate the power of examining multiple statistics for determining the secondary parameters that are vital for understanding the galaxy–halo connection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Aemulus $\nu$: precise predictions for matter and biased tracer power spectra in the presence of neutrinos

Here we present the Aemulusν simulations: a suite of 150 (1.05 h -1 Gpc) 3 N-body simulations with a mass resolution of 3.51 × 10 10 Ω cb /0.3 h -1 M ⊙ in a wνCDM cosmological parameter space. The simulations have been explicitly designed to span a broad range in σ8 to facilitate investigations of tension between large scale structure and cosmic microwave background cosmological probes. Neutrinos are treated as a second particle species to ensure accuracy to 0.5 eV, the maximum neutrino mass that we have simulated. By employing Zel'dovich control variates, we increase the effective volume of our simulations by factors of 10-10 5 depending on the statistic in question. As a first application of these simulations, we build new hybrid effective field theory and matter power spectrum surrogate models, demonstrating that they achieve ≤ 1% accuracy for k ≤ 1 hMpc -1 and 0 ≤ z ≤ 3, and ≤ 2% accuracy for k ≤ 4 hMpc -1 for the matter power spectrum. We publicly release the trained surrogate models, and estimates of the surrogate model errors in the hope that they will be broadly applicable to a range of cosmological analyses for many years to come.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Interacting scales and energy transfer in isotropic turbulence

The dependence of the energy transfer process on the disparity of the interacting scales is investigated in the inertial and far-dissipation ranges of isotropic turbulence. The strategy for generating the simulated flow fields and the choice of a disparity parameter to characterize the scaling of the interactions is discussed. The inertial range is found to be dominated by relatively local interactions, in agreement with the Kolmogorov assumption. The far-dissipation is found to be dominated by relatively non-local interactions, supporting the classical notion that the far-dissipation range is slaved to the Kolmogorov scales. The measured energy transfer is compared with the classical models of Heisenberg, Obukhov, and the more detailed analysis of Tennekes and Lumley. The energy transfer statistics measured in the numerically simulated flows are found to be nearly self-similar for wave numbers in the inertial range. Using the self-similar form measured within the limited scale range of the simulation, an 'ideal' energy transfer function and the corresponding energy flux rate for an inertial range of infinite extent are constructed. From this flux rate, the Kolmogorov constant is calculated to be 1.5, in excellent agreement with experiments.

Zhou, YE↗
DOE OSTI2025

The Impacts of Rotational Mixing on the Precipitation Simulated by a Convection Permitting Model

With increased availability of computational resources, regional and global scale convection-permitting model (CPM, Δx ~ 1–10 km) simulations are becoming more common. CPMs have improved accuracy in their representation of deep convection and mesoscale convective systems (MCSs) compared to coarser resolution models. However, CPMs still exhibit convective cloud and precipitation biases relative to observations, notably a lesser frequency of light precipitation rates and greater frequency of heavy precipitation rates. In this work we hypothesize that these CPM biases are related to under-resolved mixing between convective updrafts and their surrounding environment. To test this hypothesis, we introduce a parameterization to the Weather Research and Forecasting model (WRF) that adds a small angular rotation of the grid-scale flow about the axis perpendicular to the plane of convective drafts. This rotated flow is then allowed to alter advection of moisture and hydrometeors. The effects of such mixing on precipitation characteristics are evaluated in month-long 4-km grid spacing simulations over the Amazon. The enhanced mixing transports moisture and condensate from convective cores to other areas including downdrafts. This increases the frequency of low-precipitable water and light precipitation. It also decreases the frequency of intense precipitation from isolated deep convection and MCSs, increases cloud top temperatures, reduces radar echo-top heights, and increases overall precipitation by altering the relationship of precipitation with precipitable water, in better agreement with observations. The results suggest when optimized using multiple observations, such an approach may provide a path toward more accurate representation of convection and precipitation statistics in convection-permitting simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulating Z 2 lattice gauge theory on a quantum computer

The utility of quantum computers for simulating lattice gauge theories is currently limited by the noisiness of the physical hardware. Various quantum error mitigation strategies exist to reduce the statistical and systematic uncertainties in quantum simulations via improved algorithms and analysis strategies. We perform quantum simulations of Z 2 gauge theory with matter to study the efficacy and interplay of different error mitigation methods: readout error mitigation, randomized compiling, rescaling, and dynamical decoupling. We compute Minkowski correlation functions in this confining gauge theory and extract the mass of the lightest spin-1 state from fits to their time dependence. Quantum error mitigation extends the range of times over which our correlation function calculations are accurate by a factor of 6 and is therefore essential for obtaining reliable masses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2021

Mars 2020 Navigation Performance

On February 18, 2021, the Perseverance rover and Ingenuity helicopter demonstration landed at Jezero Crater. The Entry, Descent, and Landing (EDL) architecture, largely the same used to land Curiosity at Gale Crater on August 6, 2012, required high-fidelity flight dynamics simulation with two independent tools to verify performance. The process for creating the EDL simulation using the Dynamics Simulator for Entry, Descent and Surface landing (DSENDS) tool will be discussed, along with its use for for independent verification of the EDL statistical analysis results and reference trajectory simulation. Analysis and usage details both in development and cruise, along with post-landing assessment of the prediction performance of the simulation, will also be discussed.

Wagner, Sean↗