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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Data assimilation for combustion ignition delay time simulation using schlieren image velocimetry

This study sought to improve the accuracy of simulating spray penetration and combustion ignition delay by means of data assimilation (DA). The simulations were conducted using the Reynolds-averaged Navier–Stokes (RANS) equations and assimilating the schlieren image data. In DA, an ensemble square root filter (EnSRF) was used to build the statistical model, making the simulation results more accurate without any change in the governing equations. Recognizing that the spray-cone injection angle has a large effect on penetration, we created ensemble members with different injection angles. And we applied the two-component velocity distribution calculated via SIV and updated both velocity and temperature by using a DA statistical model derived from RANS ensemble simulations. The ignition delay time is generally known to vary even under the same experimental conditions because it is influenced by many factors. In this study, we attempted the transient DA-assisted RANS simulation to predict the ignition delay time even when the temporal resolution and accuracy of the observation data ware insufficient. Our trials offer an example of how a combination of techniques can be effectively used to assimilate experimental data obtained under restricted conditions.

Combustion simulation↗
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 OSTI2021

Angular Momentum Removal by Neutron and γ-Ray Emissions during Fission Fragment Decays

In this work, we investigate the angular momentum removal from fission fragments (FFs) through neutron and γ-ray emission, finding that about half the neutrons are emitted with angular momenta ≥ 1.5 ℏ and that the change in angular momentum after the emission of neutrons and statistical γ rays is significant, contradicting usual assumptions. Per fission event, in our simulations, the neutron and statistical γ-ray emissions change the spin of the fragment by 3.5 – 5 ℏ , with a large standard deviation comparable to the average value. Such wide angular momentum removal distributions can hide any underlying correlations in the fission fragment initial spin values. Within our model, we reproduce data on spin measurements from discrete transitions after neutron emissions, especially in the case of light FFs. The agreement further improves for the heavy fragments if one removes from the analysis the events that would produce isomeric states. Finally, we show that while in our model the initial FF spins do not follow a sawtoothlike behavior observed in recent measurements, the average FF spin computed after neutron and statistical γ emissions exhibits a shape that resembles a sawtooth. This suggests that the average FF spin measured after statistical emissions is not necessarily connected with the scission mechanism as previously implied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Prospects for silvicultural enhancement of fire resistance in mesic westside forests of the Pacific Northwest

Increasing wildfire activity in mesic, temperate Pacific Northwest forests west of the Cascade Range crest has stimulated interest in understanding whether alternative forest management practices could reduce risk of stand-replacing fire. To explore how management can enhance fire resistance in these forests and assess tradeoffs among resistance enhancement, carbon sequestration and storage, and economic returns, we conducted 40-year simulations of stand development with BioSum, a framework for conducting landscape analysis with the Forest Vegetation Simulator (FVS), utilizing a statistically representative and spatially balanced sample of Forest Inventory and Analysis (FIA) plots. Simulation outcomes under business-as-usual silviculture were contrasted with fire-aware silviculture, and treatment optimization logic was developed and applied to represent landscape-scale outcomes under business-as-usual and fire-focused management scenarios. Simulation results indicate that fire-aware prescriptions and fire-focused management can meaningfully enhance stand- and landscape-scale fire resistance of westside forests under less than extreme fire weather, but at the cost of lower economic returns and reduced net carbon storage and sequestration over the 40-year analysis window. Shifting from business-as-usual regeneration harvests with short rotations to fire-aware, episodic selection harvest improved fire resistance the most, especially in young privately-owned forests, and with only modest tradeoffs in carbon and economic outcomes. While fire-aware treatments generally reduced net present value from forest operations over business-as-usual, most treatments still generated positive net present value and could be implemented without subsidy. Fire-aware prescriptions that removed and utilized non-merchantable harvest residues instead of burning them, via either pile or broadcast burning, partially mitigated carbon emissions associated with fire-aware treatments, with about the same improvement in fire resistance. Given the currently limited institutional and financial capacity to implement fire resistance enhancing treatments at scale, the insights from this analysis may aid managers seeking to elevate fire resistance to prioritize where and how to manage.

Science & Technology - Other Topics↗
NASA NTRS1996

Guidance and Control Design for High-Speed Rollout and Turnoff (ROTO)

A ROTO architecture, braking and steering control law and display designs for a research high speed Rollout and Turnoff (ROTO) system applicable to transport class aircraft are described herein. Minimum surface friction and FMS database requirements are also documented. The control law designs were developed with the aid of a non-real time simulation program incorporating airframe and gear dynamics as well as steering and braking guidance algorithms. An attainable objective of this ROTO system, as seen from the results of this study, is to assure that the studied aircraft can land with runway occupancy times less then 53 seconds. Runway occupancy time is measured from the time the aircraft crosses the runway threshold until its wing tip clears the near side of the runway. Turnoff ground speeds of 70 knots onto 30 degree exits are allowed with dry and wet surface conditions. Simulation time history and statistical data are documented herein. Parameters which were treated as variables in the simulation study include aircraft touchdown weight/speed/location, aircraft CG, runway friction, sensor noise and winds. After further design and development of the ROTO control system beyond the system developed earlier, aft CG MD-11 aircraft no longer require auto-asymmetric braking (steering) and fly-by-wire nose gear steering. However, the auto ROTO nose gear hysteresis must be less than 2 degrees. The 2 sigma dispersion certified for MD-11 CATIIIB is acceptable. Using this longitudinal dispersion, three ROTO exits are recommended at 3300, 4950 and 6750 feet past the runway threshold. The 3300 foot exit is required for MD-81 class aircraft. Designs documented in this report are valid for the assumptions/models used in this simulation. It is believed that the results will apply to the general class of transport aircraft; however further effort is required to validate this assumption for the general case.

Goldthorpe, S. H.↗
NASA NTRS2012

Comparison of Statistical Estimation Techniques for Mars Entry, Descent, and Landing Reconstruction from MEDLI-like Data Sources

Flight data from an entry, descent, and landing (EDL) sequence can be used to reconstruct the vehicle's trajectory, aerodynamic coefficients and the atmospheric profile experienced by the vehicle. Past Mars missions have contained instruments that do not provide direct measurement of the freestream atmospheric conditions. Thus, the uncertainties in the atmospheric reconstruction and the aerodynamic database knowledge could not be separated. The upcoming Mars Science Laboratory (MSL) will take measurements of the pressure distribution on the aeroshell forebody during entry and will allow freestream atmospheric conditions to be partially observable. This data provides a mean to separate atmospheric and aerodynamic uncertainties and is part of the MSL EDL Instrumentation (MEDLI) project. Methods to estimate the flight performance statistically using on-board measurements are demonstrated here through the use of simulated Mars data. Different statistical estimators are used to demonstrate which estimator best quantifies the uncertainties in the flight parameters. The techniques demonstrated herein are planned for application to the MSL flight dataset after the spacecraft lands on Mars in August 2012.

Dutta, Soumyo↗
NASA NTRS2012

Broadband Noise Prediction When Turbulence Simulation Is Available - Derivation of Formulation 2B and Its Statistical Analysis

We show that a simple modification of Formulation 1 of Farassat results in a new analytic expression that is highly suitable for broadband noise prediction when extensive turbulence simulation is available. This result satisfies all the stringent requirements, such as permitting the use of the exact geometry and kinematics of the moving body, that we have set as our goal in the derivation of useful acoustic formulas for the prediction of rotating blade and airframe noise. We also derive a simple analytic expression for the autocorrelation of the acoustic pressure that is valid in the near and far fields. Our analysis is based on the time integral of the acoustic pressure that can easily be obtained at any resolution for any observer time interval and digitally analyzed for broadband noise prediction. We have named this result as Formulation 2B of Farassat. One significant consequence of Formulation 2B is the derivation of the acoustic velocity potential for the thickness and loading terms of the Ffowcs Williams-Hawkings (FW-H) equation. This will greatly enhance the usefulness of the Fast Scattering Code (FSC) by providing a high fidelity boundary condition input for scattering predictions.

Farassat, Fereidoun↗
NASA NTRS1972

Simulation of superposed coherent and chaotic radiation of arbitrary spectral shape.

A simulation technique is described and utilized to generate superposed coherent and chaotic (thermal) radiation of arbitrary spectral shape. The statistical properties of a simulated radiation field with a Lorentzian spectral density are investigated with a photoelectron counting experiment. The experimental photocount distribution and normalized mth-order factorial moments are compared to theory and verify that the simulated radiation field, an appropriately modulated laser beam, has the expected statistical properties. The concept described indicates that one may, in principle, generate a thermal source of comparable intensity to that of a laser, arbitrary spectral shape, and bandwidth.

Ruggieri, N. F.↗
NASA NTRS1992

Statistics of the cosmic Mach number from numerical simulations of a cold dark matter universe

Results are presented of an analysis of the cosmic Mach number, M, the ratio of the streaming velocity, v, to the random velocity dispersion, sigma, of galaxies in a given patch of the universe, which was performed on the basis of hydrodynamical simulations of the cold dark matter scenario. Galaxy formation is modeled by application of detailed physical processes rather than by the ad hoc assumption of 'bias' between dark matter and galaxy fluctuations. The correlation between M and sigma is found to be very weak for both components. No evidence is found for a physical 'velocity bias' in the quantities which appear in the definition of M. Standard cold-dark-matter-dominated universes are in conflict, at a statistically significant level, with the available observation, in that they predict a Mach number considerably lower than is observed.

Suto, Yasushi↗
DOE OSTI2020

Perturbation theory for modeling galaxy bias: validation with simulations of the Dark Energy Survey

We describe perturbation theory (PT) models of galaxy bias for applications to photometric galaxy surveys. We model the galaxy-galaxy and galaxy-matter correlation functions in configuration space and validate against measurements from mock catalogs designed for the Dark Energy Survey (DES). We find that an effective PT model with five galaxy bias parameters provides a good description of the 3D correlation functions above scales of 4 Mpc/$h$ and $z < 1$. Our tests show that at the projected precision of the DES-Year 3 analysis, two of the non-linear bias parameters can be fixed to their co-evolution values, and a third (the $k^2$ term for higher derivative bias) set to zero. The agreement is typically at the 2 percent level over scales of interest, which is the statistical uncertainty of our simulation measurements. To achieve this level of agreement, our {\it fiducial} model requires using the full non-linear matter power spectrum (rather than the 1-loop PT one). We also measure the relationship between the non-linear and linear bias parameters and compare them to their expected co-evolution values. We use these tests to motivate the galaxy bias model and scale cuts for the cosmological analysis of the Dark Energy Survey; our conclusions are generally applicable to all photometric surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Solving inverse problems in stochastic models using deep neural networks and adversarial training

Inverse problems associated with stochastic models constitute a significant portion of scientific and engineering applications. In such cases the unknown quantities are distributions. The applicability of traditional methods is limited because of their demanding assumptions or prohibitive computational consumption; for example, maximum likelihood methods require closed-form density functions, and Markov Chain Monte Carlo needs a large number of simulations. We propose a new method that estimates the unknown distribution by matching the statistical properties between observed and simulated random processes. We leverage the expressive power of neural networks to approximate the unknown distribution and use a discriminative neural network for computing the statistical discrepancies between the observed and simulated random processes. Here we demonstrated numerically that the proposed methods can estimate both the model parameters and learn complicated unknown distributions.

42 ENGINEERING↗
DOE OSTI2021

Quasi-brittle porous material: Simulated effect of stochastic air void structure on compressive strength

Highlights: • A simulation procedure to study the effect of porosity on the compressive strength of a quasi-brittle material is introduced • Based on the simulated results, the relationship between compressive strength and air void fraction (porosity) is presented. • A linear relationship was achieved between the cubic root of porosity and the simulated compressive strength. The effect of porosity comprised of spherical air voids on the compressive strength of quasi-brittle material was studied via simulations. The simulated porous structures were based on pore size distributions of two mortar samples measured by X-ray microtomography. While the simulation method set practical limits on the size of sample, the base of the statistics was established by simulating 128 small structures generated by sampling from pore structures of two mortars. By studying the application of the classical strength-porosity formulas to the simulated data, a new simple model was formed. A linear relationship was achieved between the cubic root of air void fraction (porosity) and the simulated compressive strength. The reasons for scattering of simulated strength around fitted trend remained unresolved in this study; no clear dependence on pore number or other distribution properties was observed. With the presented simulation approach, the dependence of compressive strength on porosity is achieved independently of disturbances that occur in experimental studies creating understanding of compressional behavior of low porosity materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Prototype Courthouse Building Energy Model: Building Characteristics and Energy Simulation Results

To coordinate commercial building energy research activities, the US Department of Energy uses a suite of commercial prototype building models, which includes 16 building types and covers 80% of US commercial floorspace. Efforts are underway to expand this suite by developing prototype models for additional building types, including a courthouse. Informed by building design guides, databases, documented projects, and personal communication with courthouse design experts, a 3-story, 4-courtroom, 69,324 ft 2 building was considered as the prototype courthouse to represent an average-sized courthouse in the United States. This report documents the building and system characteristics of the prototype courthouse model combined with climate-specific and construction vintage–specific requirements for the building envelope and systems from ASHRAE Standard 90.1. The energy simulation results are presented for 4 construction vintages (2004, 2007, 2010, and 2013) and 15 ASHRAE climate zones. To verify the energy simulation results, energy use statistics of existing courthouses from different building database are also presented. Based on the simulated energy use, the site energy use intensity of the prototype courthouse ranged from 45 to 159 kBtu/ft 2 compared with the average energy use intensity of 94.7 kBtu/ft 2 for the courthouse building type in the 2012 Commercial Building Energy Consumption Survey (CBECS) Data.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1989

Research in LMSS propagation

The Virginia Tech Satellite Communications Group has participated in the Land Mobile Satellite System (LMSS) program through JPL sponsorship since 1985. Involvement has mainly been in modeling and simulation of propagation characteristics and effects. Models developed to predict cummulative fade distributions for fading LMSS signals include LMSSMOD and the Simple Models which approximate LMSSMOD. Models to predict the mean and standard deviation of signal attenuation through roadside vegetation, namely the Average Path Model, were developed. In the area of simulation, efforts have centered around the development of a software simulator that uses data bases derived from experimental data to generate simulated data with arbitrary statistical behavior. This work has progressed to the development of an integrated analysis and simulation package, LIPS. The basic theory and results for the models and simulator have been previously documented in reports and papers. All LMSS activities are summarized and details of this year's efforts are given.

Barts, R. M.↗
DOE OSTI2022

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗
NASA NTRS1997

Factor Analysis Applied the VFY-218 RCS Data

Present statistical factor analysis of computer simulations and measurement data for the VFY-218 configuration. Factor analysis try to quantify the statistical grouping of measurements and simulations.

Woo, Alex↗
NASA NTRS1991

Monte Carlo simulation of launchsite winds at Kennedy Space Center

This paper develops and validates an easily implemented model for simulating random horizontal wind profiles over the Kennedy Space Center (KSC) at Cape Canaveral, Florida. The model is intended for use in Monte Carlo launch vehicle simulations of the type employed in mission planning, where the large number of profiles needed for statistical fidelity of such simulation experiments makes the use of actual wind measurements impractical. The model is based on measurements made at KSC and represents vertical correlations by a decaying exponential model which is parameterized via least-squares parameter optimization against the sample data. The validity of the model is evaluated by comparing two Monte Carlo simulations of an asymmetric, heavy-lift launch vehicle. In the first simulation, the measured wind profiles are used, while in the second, the wind profiles are generated using the stochastic model. The simulations indicate that the use of either the measured or simulated wind field results in similar launch vehicle performance.

Queen, Eric M.↗
DOE OSTI2025

Experimental and Computational Characterization of a Modified Sioutas Cascade Impactor for Respirable Radioactive Aerosols

Oak Ridge National Laboratory is collecting and characterizing aerosols released when spent nuclear fuel (SNF) rods are fractured in bending. An aerosol collection system was designed and tested to collect respirable sized (<10 μm aerodynamic diameter [AED]) particulates inside a hot cell facility. The setup is a modified version of the commercially available Sioutas cascade impactor, to which additional stages were added to expand the aerosol collection range from 2.5 to ~15 μm AED. To accommodate the additional stages and specific test conditions, the operating flow rate for aerosol collection was reduced, and testing was conducted by using pressure drop measurements, surrogate dust collection, and particle size characterization. The fluid flow distribution within the cascade and its stages was simulated in STAR-CCM+, and the stage-wise pressure drops obtained using the computational fluid dynamics model were then compared to experimental data. Lagrangian particle simulations were also performed, and stage-wise collection statistics were obtained from the simulation for comparison with the experimental data obtained using SNF-surrogate dust particles. The results provide valuable insights into the stage-wise particle collection characteristics of the modified cascade impactor and can also be used to improve the prediction accuracy of the manufacturer-determined analytical correlations.

aerosol modeling↗