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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Flight investigation of a four-dimensional terminal area guidance system for STOL aircraft

A series of flight tests and fast-time simulations were conducted, using the augmentor wing jet STOL research aircraft and the STOLAND 4D-RNAV system to add to the growing data base of 4D-RNAV system performance capabilities. To obtain statistically meaningful data a limited amount of flight data were supplemented by a statistically significant amount of data obtained from fast-time simulation. The results of these tests are reported. Included are comparisons of the 4D-RNAV estimated winds with actual winds encountered in flight, as well as data on along-track navigation and guidance errors, and time-of-arrival errors at the final approach waypoint. In addition, a slight improvement of the STOLAND 4D-RNAV system is proposed and demonstrated, using the fast-time simulation.

Neuman, F.↗
DOE OSTI2021

Quantum simulation of molecules without fermionic encoding of the wave function

Abstract Molecular simulations generally require fermionic encoding in which fermion statistics are encoded into the qubit representation of the wave function. Recent calculations suggest that fermionic encoding of the wave function can be bypassed, leading to more efficient quantum computations. Here we show that the two-electron reduced density matrix (2-RDM) can be expressed as a unique functional of the unencoded N -qubit-particle wave function without approximation, and hence, the energy can be expressed as a functional of the 2-RDM without fermionic encoding of the wave function. In contrast to current hardware-efficient methods, the derived functional has a unique, one-to-one (and onto) mapping between the qubit-particle wave functions and 2-RDMs, which avoids the over-parametrization that can lead to optimization difficulties such as barren plateaus. An application to computing the ground-state energy and 2-RDM of H 4 is presented.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

GEANT4 parameter tuning using Professor

The GEANT4 toolkit is used extensively in high energy physics to simulate the passage of particles through matter and to predict effects such as detector efficiencies and smearing. GEANT4 uses many underlying models to predict particle interaction kinematics, and uncertainty in these models leads to uncertainty in high energy physics measurements. The GEANT4 collaboration recently made free parameters in some models accessible through partnership with GEANT4 developers. We present a study of the impact of varying parameters in three GEANT4 hadronic physics models on agreement with thin target datasets and describe fits to these datasets using the Professor model tuning framework [1]. We find that varying parameters produces substantially better agreement with some datasets, but that more degrees of freedom are required for full agreement. This work is a first step towards a common framework for propagating uncertainties in GEANT4 models to high energy physics measurements, and we outline future work required to complete that goal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗
NASA NTRS1999

Mars Characterization for Future Missions

The purpose of this research project was to create statistical products which could of use to the engineering and scientific communities planning future missions to Mars. This has been accomplished. Using simulations of Mars' atmosphere under a variety of conditions, we have created statistical databases quantifying the behavior of Mars' atmosphere under a variety of conditions expected to be encountered on forthcoming missions. These data are now being incorporated into a new version of Mars-GRAM.

Bridger, Alison F. C.↗
DOE OSTI2024

Simulating cosmological evolution by quantum quench of an atomic Bose-Einstein condensate

In cosmological evolution, it is the homogeneous scalar field (inflaton) that drives the universe to expand isotropically and to generate standard model particles. However, to simulate cosmology, atomic gas research has focused on the dynamics of Bose-Einstein condensates (BEC) with continuously applied forces. In this paper we argue that a complementary approach needs also to be pursued; we thus consider the analog BEC experiments in a nondriven and naturally closed atomic system. We implement this using a BEC in an optical lattice which, after a quench, freely transitions from an unstable to a stable state. Furthermore, this dynamical evolution displays the counterpart “preheating,” “reheating,” and “thermalization” phases of cosmology. Importantly, our studies of these analog processes yield tractable analytic models. Additionally, of great utility to the cold atom community, such understanding elucidates the dynamics of nonadiabatic condensate preparation. Indeed, the dynamical processes discussed here are generic and in future cold atom reequilibration experiments it will be important to observe both the preheating stage, corresponding to a fragmented condensate, and the reheating stage, corresponding to a particle cloud.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2010

Failure Bounding And Sensitivity Analysis Applied To Monte Carlo Entry, Descent, And Landing Simulations

In the study of entry, descent, and landing, Monte Carlo sampling methods are often employed to study the uncertainty in the designed trajectory. The large number of uncertain inputs and outputs, coupled with complicated non-linear models, can make interpretation of the results difficult. Three methods that provide statistical insights are applied to an entry, descent, and landing simulation. The advantages and disadvantages of each method are discussed in terms of the insights gained versus the computational cost. The first method investigated was failure domain bounding which aims to reduce the computational cost of assessing the failure probability. Next a variance-based sensitivity analysis was studied for the ability to identify which input variable uncertainty has the greatest impact on the uncertainty of an output. Finally, probabilistic sensitivity analysis is used to calculate certain sensitivities at a reduced computational cost. These methods produce valuable information that identifies critical mission parameters and needs for new technology, but generally at a significant computational cost.

Gaebler, John A.↗
NASA NTRS2001

Simulated Performance of the Orbiting Wide-angle Light Collectors (OWL) Experiment

The Orbiting Wide-angle Light collectors (OWL) experiment is in NASA's mid-term strategic plan and will stereoscopically image, from equatorial orbit, the air fluorescence signal generated by airshowers induced by the ultrahigh energy (E greater than few x 10(exp 19) eV) component of the cosmic radiation. The use of a space-based platform enables an extremely large event acceptance aperture and thus will allow a high statistics measurement of these rare events. Detailed Monte Carlo simulations are required to quantify the physics potential of the mission as well as optimize the instrumental parameters. This paper reports on the results of the GSFC Monte Carlo simulation for two different, OWL instrument baseline designs. These results indicate that, assuming a continuation of the cosmic ray spectrum (theta approximately E(exp -2.75), OWL could have an event rate of 4000 events/year with E greater than or equal to 10(exp 20) eV. Preliminary results, based upon these Monte Carlo simulations, indicate that events can be accurately reconstructed in the detector focal plane arrays for the OWL instrument baseline designs under consideration.

Krizmanic, J. F.↗
NASA NTRS1994

TDRS orbit determination by radio interferometry

In support of a NASA study on the application of radio interferometry to satellite orbit determination, MITRE developed a simulation tool for assessing interferometry tracking accuracy. The Orbit Determination Accuracy Estimator (ODAE) models the general batch maximum likelihood orbit determination algorithms of the Goddard Trajectory Determination System (GTDS) with the group and phase delay measurements from radio interferometry. ODAE models the statistical properties of tracking error sources, including inherent observable imprecision, atmospheric delays, clock offsets, station location uncertainty, and measurement biases, and through Monte Carlo simulation, ODAE calculates the statistical properties of errors in the predicted satellites state vector. This paper presents results from ODAE application to orbit determination of the Tracking and Data Relay Satellite (TDRS) by radio interferometry. Conclusions about optimal ground station locations for interferometric tracking of TDRS are presented, along with a discussion of operational advantages of radio interferometry.

Pavloff, Michael S.↗
NASA NTRS2000

A statistical model of the human core-temperature circadian rhythm

We formulate a statistical model of the human core-temperature circadian rhythm in which the circadian signal is modeled as a van der Pol oscillator, the thermoregulatory response is represented as a first-order autoregressive process, and the evoked effect of activity is modeled with a function specific for each circadian protocol. The new model directly links differential equation-based simulation models and harmonic regression analysis methods and permits statistical analysis of both static and dynamical properties of the circadian pacemaker from experimental data. We estimate the model parameters by using numerically efficient maximum likelihood algorithms and analyze human core-temperature data from forced desynchrony, free-run, and constant-routine protocols. By representing explicitly the dynamical effects of ambient light input to the human circadian pacemaker, the new model can estimate with high precision the correct intrinsic period of this oscillator ( approximately 24 h) from both free-run and forced desynchrony studies. Although the van der Pol model approximates well the dynamical features of the circadian pacemaker, the optimal dynamical model of the human biological clock may have a harmonic structure different from that of the van der Pol oscillator.

NASA Discipline Regulatory Physiology↗
NASA NTRS2012

Dipole Alignment in Rotating MHD Turbulence

We present numerical results from long-term CPU and GPU simulations of rotating, homogeneous, magnetohydrodynamic (MHD) turbulence, and discuss their connection to the spherically bounded case. We compare our numerical results with a statistical theory of geodynamo action that has evolved from the absolute equilibrium ensemble theory of ideal MHD turbulence, which is based on the ideal MHD invariants are energy, cross helicity and magnetic helicity. However, for rotating MHD turbulence, the cross helicity is no longer an exact invariant, although rms cross helicity becomes quasistationary during an ideal MHD simulation. This and the anisotropy imposed by rotation suggests an ansatz in which an effective, nonzero value of cross helicity is assigned to axisymmetric modes and zero cross helicity to non-axisymmetric modes. This hybrid statistics predicts a large-scale quasistationary magnetic field due to broken ergodicity , as well as dipole vector alignment with the rotation axis, both of which are observed numerically. We find that only a relatively small value of effective cross helicity leads to the prediction of a dipole moment vector that is closely aligned (less than 10 degrees) with the rotation axis. We also discuss the effect of initial conditions, dissipation and grid size on the numerical simulations and statistical theory.

Shebalin, John V.↗
DOE OSTI2022

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Turbulence statistics in fully developed channel flow at low Reynolds number

A direct numerical simulation of a turbulent channel flow is performed. The unsteady Navier-Stokes equations are solved numerically at a Reynolds number of 3300, based on the mean centerline velocity and channel half-width, with about 4 million grid points. All essential turbulence scales are resolved on the computational grid and no subgrid model is used. A large number of turbulence statistics are computed and compared with the existing experimental data at comparable Reynolds numbers. Agreements as well as discrepancies are discussed in detail. Particular attention is given to the behavior of turbulence correlations near the wall. A number of statistical correlations which are complementary to the existing experimental data are reported for the first time.

Kim, John↗
NASA NTRS1990

The wet tropospheric range correction - Product intercomparisons and the simulated effect for tropical Pacific altimeter retrievals

Four data sets used to correct satellite altimetry for the effects of tropospheric water vapor over the tropical Pacific Ocean are evaluated statistically, with a focus on the period 1979-1983. Satellite altimeter measurements are simulated by degrading model sea-level solutions with selected water-vapor data, and seasonal and interannual variations are investigated. The results are presented in extensive contour maps and characterized in detail. The standard deviation between two satellite-based data sets (Nimbus-7 SMMR and TOVS) was found to be 2-4 cm, while both the rms difference between seasonal cycles of the two sets and the difference between the mean fields were 1-2 cm. Larger differences (4-8 cm) were found between these data sets and those based on operational analyses (ECMWF and Fleet Numerical Oceanography Center). The maximum interannual anomaly in the water-vapor correction was 12 cm, for the eastern Pacific at the peak of the 1982-1983 ENSO event.

Zimbelman, Darrell F.↗
NASA NTRS2020

Local and Remote Mean and Extreme Temperature Response to Regional Aerosol Emissions Reductions

The climatic implications of regional aerosol and precursor emissions reductions implemented to protect human health are poorly understood. We investigate the mean and extreme temperature response to regional changes in aerosol emissions using three coupled chemistry–climate models: NOAA GFDL CM3, NCAR CESM1, and NASA GISS-E2. Our approach contrasts a long present-day control simulation from each model (up to 400 years with perpetual year 2000 or 2005 emissions) with 14 individual aerosol emissions perturbation simulations (160–240 years each). We perturb emissions of sulfur dioxide (SO2) and/or carbonaceous aerosol within six world regions and assess the statistical significance of mean and extreme temperature responses relative to internal variability determined by the control simulation and across the models. In all models, the global mean surface temperature response (perturbation minus control) to SO2 and/or carbonaceous aerosol is mostly positive (warming) and statistically significant and ranges from +0.17 K (Europe SO2) to -0.06 K (US BC). The warming response to SO2 reductions is strongest in the US and Europe perturbation simulations, both globally and regionally, with Arctic warming up to 1 K due to a removal of European anthropogenic SO2 emissions alone; however, even emissions from regions remote to the Arctic, such as SO2 from India, significantly warm the Arctic by up to 0.5 K. Arctic warming is the most robust response across each model and several aerosol emissions perturbations. The temperature response in the Northern Hemisphere midlatitudes is most sensitive to emissions perturbations within that region. In the tropics, however, the temperature response to emissions perturbations is roughly the same in magnitude as emissions perturbations either within or outside of the tropics. We find that climate sensitivity to regional aerosol perturbations ranges from 0.5 to 1.0 K (W m(exp -2))(exp -1) depending on the region and aerosol composition and is larger than the climate sensitivity to a doubling of CO2 in two of three models. We update previous estimates of regional temperature potential (RTP), a metric for estimating the regional temperature responses to a regional emissions perturbation that can facilitate assessment of climate impacts with integrated assessment models without requiring computationally demanding coupled climate model simulations. These calculations indicate a robust regional response to aerosol forcing within the Northern Hemisphere midlatitudes, regardless of where the aerosol forcing is located longitudinally. We show that regional aerosol perturbations can significantly increase extreme temperatures on the regional scale. Except in the Arctic in the summer, extreme temperature responses largely mirror mean temperature responses to regional aerosol perturbations through a shift of the temperature distributions and are mostly dominated by local rather than remote aerosol forcing.

Westervelt, Daniel M.↗
NASA NTRS1978

Experiments in monthly mean simulation of the atmosphere with a coarse-mesh general circulation model

The Hansen atmospheric model was used to compute five monthly forecasts (October 1976 through February 1977). The comparison is based on an energetics analysis, meridional and vertical profiles, error statistics, and prognostic and observed mean maps. The monthly mean model simulations suffer from several defects. There is, in general, no skill in the simulation of the monthly mean sea-level pressure field, and only marginal skill is indicated for the 850 mb temperatures and 500 mb heights. The coarse-mesh model appears to generate a less satisfactory monthly mean simulation than the finer mesh GISS model.

Lutz, R. J.↗
NASA NTRS1990

Patterns in simulated turbulent channel flow

The surface stress is a key diagnostic in wall bounded turbulent flows. Large fluctuations in the stress are believed to be associated with intermittent 'bursting' events during which a large proportion of the turbulence production takes place. If this is so, then a detailed investigation of the structure of the surface stress and its spatial relationship to events within the flow could have wide application in drag reduction and other aspects of flow control. The initial phase of this project, therefore, concentrated on the surface stress field. The first objective is to carry out a statistical analysis of the instantaneous surface stress in a simulated turbulent channel flow, including comparison with multipoint experimental data from a zero pressure gradient turbulent boundary layer. The second objective is to apply a simple pattern recognition procedure to educe the characteristic spatial structure of various flow diagnostics. The final objective is to extend the pattern recognition analysis to examine the whole three dimensional structure of the flow.

Stretch, D.↗
NASA NTRS1992

Modeling and simulation of mobile satellite propagation

Mobile satellite systems are subject to severe fading due to blockage of the line-of-sight (LOS) path by roadside vegetation. A thorough understanding of the fading effects is necessary for the design of a reliable land mobile satellite system (LMSS). Analytical and empirical models are presented for predicting fade statistics for vegetative shadowing of mobile satellite terminals. A software simulator for generating simulated fade data is also presented. A physical model relating physical path parameters to propagation model parameters is presented, and results using the model are shown.

Barts, R. M.↗