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 469 records · Page 26

Light propagation in inhomogeneous universes - The ray-shooting method

The propagation of light in a clumpy universe is studied using the ray-shooting method. It is found that 'empty cones' in a clumpy universe are rare, which points out the limitation of frequently used distance measures in an inhomogeneous universe. The results show some qualitatively new features of multiple-gravitational-lens geometry, and a likely explanation for these features is presented. A statistical analysis of the simulations give amplification probability distributions for extended sources, as well as the distribution for the ratio of amplifications of compact sources to amplifications of more extended sources. Sample light curves of individual sources are obtained. Finally, the results support the validity of recent investigations on the influence of gravitational light deflection on source counts of compact extragalactic objects.

Schneider, Peter↗

Triplets of galaxies: Their dynamics, evolution, and the origin of chaos in them

Recently Karachentsev's group at The Smithsonian Astrophysical Observatory (SAO) (6-meter Telescope Observatory) published a list of 84 triple systems of galaxies with their distances, radial (line of sight) velocities, and angular sizes (Karachentseva et al., 1988). This gives a new ground for studies of the dark matter problem which fills the gap between the large cosmic scales (White, 1987; Dekel and Rees, 1987, and Einasto et al., 1977) and the scale of individual galaxies (Erickson et al., 1987). The data on the typical velocity dispersions and linear dimension of the triplets indicate that they contain considerable amounts of dark matter (see also earlier work of Karachentseva et al., 1979). Numerical simulations show that the statistical characteristics of the Karachentsev triplets can be imitated by model ensembles of triple systems with dark matter masses M sub d = (1-3 x 10(exp 12) M sub O, which is almost ten times greater than the typical mass of stellar galaxies estimated by the standard mass-to-luminosity ration (Kiseleva and Chernin, 1988). Here, the authors report that important information can be drawn from the data on the visible configurations of these systems. The statistics of configurations provide an independent evidence for dark matter in the triplets; moreover, it enables one to argue that dark matter seems to be distributed over the whole volume of the typical triplet forming its common corona rather than concentrated within individual coronae (or haloes) of the member galaxies.

Chernin, Arthur D.↗

The ideal Kolmogorov inertial range and constant

The energy transfer statistics measured in numerically simulated flows are found to be nearly self-similar for wavenumbers in the inertial range. Using the measured self-similar form, an 'ideal' energy transfer function and the corresponding energy flux rate were deduced. From this flux rate, the Kolmogorov constant was calculated to be 1.5, in excellent agreement with experiments.

Zhou, YE↗

Contact Analog/Compressed Symbology Heading Tape Assessment

A simulation assessed the performance, handling qualities and workload associated with a contact-analog, world-referenced heading tape as implemented on the Comanche Helmet Integrated Display Sight System (HIDSS) when compared with a screen-fixed, compressed heading tape. Six pilots, four active duty Army Aviators and two civilians flew three ADS-33 maneuvers and a traffic pattern in the Ames Vertical Motion Simulation facility. Small, but statistically significant advantages were found for the compressed symbology for handling qualities, workload, and some of the performance measures. It should be noted however that the level of performance and handling qualities for both symbology sets fell within the acceptable tolerance levels. Both symbology sets yield satisfactory handling qualities and performance in velocity stabilization mode and adequate handling qualities in the automatic flight control mode. Pilot comments about the contact analog symbology highlighted the lack of useful rate of change information in the heading tape and "blurring" due to the rapid movement of the heading tape. These issues warrant further study. Care must be taken in interpreting the operational significance of these results. The symbology sets yielded categorically similar data, i.e., acceptable handling qualities and adequate performance, so while the results point to the need for further study, their operational significance has yet to be determined.

Shively, R. Jay↗

Aerothermal Heating Predictions for Mars Microprobe

A combination of computational predictions and experimental measurements of the aerothermal heating expected on the two Mars Microprobes during their entry to Mars are presented. The maximum, non-ablating, heating rate at the vehicle's stagnation point (at alpha = 0 degrees) is predicted for an undershoot trajectory to be 194 Watts per square centimeters with associated stagnation point pressure of 0.064 atm. Maximum stagnation point pressure occurs later during the undershoot trajectory and is 0.094 atm. From computations at seven overshoot-trajectory points, the maximum heat load expected at the stagnation point is near 8800 Joules per square centimeter. Heat rates and heat loads on the vehicle's afterbody are much lower than the forebody. At zero degree angle-of-attack, heating over much of the hemi-spherical afterbody is predicted to be less than 2 percent of the stagnation point value. Good qualitative agreement is demonstrated for forebody and afterbody heating between CFD calculations at Mars entry conditions and experimental thermographic phosphor measurements from the Langley 20-Inch Mach 6 Air Tunnel. A novel approach which incorporates six degree-of-freedom trajectory simulations to perform a statistical estimate of the effect of angle-of-attack, and other off-nominal conditions, on heating is included.

Mitcheltree, R. A.↗

Numerical experiments with model monophyletic and paraphyletic taxa

The problem of how accurately paraphyletic taxa versus monophyletic (i.e., holophyletic) groups (clades) capture underlying species patterns of diversity and extinction is explored with Monte Carlo simulations. Phylogenies are modeled as stochastic trees. Paraphyletic taxa are defined in an arbitrary manner by randomly choosing progenitors and clustering all descendants not belonging to other taxa. These taxa are then examined to determine which are clades, and the remaining paraphyletic groups are dissected to discover monophyletic subgroups. Comparisons of diversity patterns and extinction rates between modeled taxa and lineages indicate that paraphyletic groups can adequately capture lineage information under a variety of conditions of diversification and mass extinction. This suggests that these groups constitute more than mere "taxonomic noise" in this context. But, strictly monophyletic groups perform somewhat better, especially with regard to mass extinctions. However, when low levels of paleontologic sampling are simulated, the veracity of clades deteriorates, especially with respect to diversity, and modeled paraphyletic taxa often capture more information about underlying lineages. Thus, for studies of diversity and taxic evolution in the fossil record, traditional paleontologic genera and families need not be rejected in favor of cladistically-defined taxa.

Non-NASA Center↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

Autonomous Navigation Performance During The Hartley 2 Comet Flyby

On November 4, 2010, the EPOXI spacecraft performed a 700-km flyby of the comet Hartley 2 as follow-on to the successful 2005 Deep Impact prime mission. EPOXI, an extended mission for the Deep Impact Flyby spacecraft, returned a wealth of visual and infrared data from Hartley 2, marking the fifth time that high-resolution images of a cometary nucleus have been captured by a spacecraft. The highest resolution science return, captured at closest approach to the comet nucleus, was enabled by use of an onboard autonomous navigation system called AutoNav. AutoNav estimates the comet-relative spacecraft trajectory using optical measurements from the Medium Resolution Imager (MRI) and provides this relative position information to the Attitude Determination and Control System (ADCS) for maintaining instrument pointing on the comet. For the EPOXI mission, AutoNav was tasked to enable continuous tracking of a smaller, more active Hartley 2, as compared to Tempel 1, through the full encounter while traveling at a higher velocity. To meet the mission goal of capturing the comet in all MRI science images, position knowledge accuracies of +/- 3.5 km (3-?) cross track and +/- 0.3 seconds (3-?) time of flight were required. A flight-code-in-the-loop Monte Carlo simulation assessed AutoNav's statistical performance under the Hartley 2 flyby dynamics and determined optimal configuration. The AutoNav performance at Hartley 2 was successful, capturing the comet in all of the MRI images. The maximum residual between observed and predicted comet locations was 20 MRI pixels, primarily influenced by the center of brightness offset from the center of mass in the observations and attitude knowledge errors. This paper discusses the Monte Carlo-based analysis that led to the final AutoNav configuration and a comparison of the predicted performance with the flyby performance.

flyby↗

NASA Tech Briefs, June 2014

Topics include: Real-Time Minimization of Tracking Error for Aircraft Systems; Detecting an Extreme Minority Class in Hyperspectral Data Using Machine Learning; KSC Spaceport Weather Data Archive; Visualizing Acquisition, Processing, and Network Statistics Through Database Queries; Simulating Data Flow via Multiple Secure Connections; Systems and Services for Near-Real-Time Web Access to NPP Data; CCSDS Telemetry Decoder VHDL Core; Thermal Response of a High-Power Switch to Short Pulses; Solar Panel and System Design to Reduce Heating and Optimize Corridors for Lower-Risk Planetary Aerobraking; Low-Cost, Very Large Diamond-Turned Metal Mirror; Very-High-Load-Capacity Air Bearing Spindle for Large Diamond Turning Machines; Elevated-Temperature, Highly Emissive Coating for Energy Dissipation of Large Surfaces; Catalyst for Treatment and Control of Post-Combustion Emissions; Thermally Activated Crack Healing Mechanism for Metallic Materials; Subsurface Imaging of Nanocomposites; Self-Healing Glass Sealants for Solid Oxide Fuel Cells and Electrolyzer Cells; Micromachined Thermopile Arrays with Novel Thermo - electric Materials; Low-Cost, High-Performance MMOD Shielding; Head-Mounted Display Latency Measurement Rig; Workspace-Safe Operation of a Force- or Impedance-Controlled Robot; Cryogenic Mixing Pump with No Moving Parts; Seal Design Feature for Redundancy Verification; Dexterous Humanoid Robot; Tethered Vehicle Control and Tracking System; Lunar Organic Waste Reformer; Digital Laser Frequency Stabilization via Cavity Locking Employing Low-Frequency Direct Modulation; Deep UV Discharge Lamps in Capillary Quartz Tubes with Light Output Coupled to an Optical Fiber; Speech Acquisition and Automatic Speech Recognition for Integrated Spacesuit Audio Systems, Version II; Advanced Sensor Technology for Algal Biotechnology; High-Speed Spectral Mapper; "Ascent - Commemorating Shuttle" - A NASA Film and Multimedia Project DVD; High-Pressure, Reduced-Kinetics Mechanism for N-Hexadecane Oxidation; Method of Error Floor Mitigation in Low-Density Parity-Check Codes; X-Ray Flaw Size Parameter for POD Studies; Large Eddy Simulation Composition Equations for Two-Phase Fully Multicomponent Turbulent Flows; Scheduling Targeted and Mapping Observations with State, Resource, and Timing Constraints;

Source record↗

Simulating the Diffusion of Hydrogen in Amorphous Silicates: A ‘Jumping’ Migration Process and its Implications for Solar Wind Implanted Lunar Volatiles

We use molecular dynamics (MD) simulations to better explain the movement of atomic hydrogen in amorphous silica and quantify the planetary science implications of these findings. Previous MD simulations had a large range of predicted values and did not agree well with experiment. Our simulations sample atomic motion for a longer duration and consider a wider range of temperatures than previous simulations. In contrast to constant atomic motion, the hydrogen atoms were shown to undergo random intermittent jumps from one oxygen atom to another, the number of which increase with temperature. Predicted diffusion coefficients had a better agreement to experimental values than previous MD simulations, suggesting the importance of longer simulation durations for better statistics. The low activation energy and jumps observed at lunar temperatures do not support the theory of diurnal variations in OH content for an undamaged amorphous silica surface. Instead, we conclude that energetic solar wind impacts can induce two competing atomic hydrogen motion processes in the exposed surface: A prompt effect that induces jumps in the temperature spike volume, but also a long term effect of damage in the structure that traps atomic hydrogen. We then use SDTrimSP to quantify the damage created during exposure and MD to demonstrate the H retention and trapping near these defects. Damage was shown to be dependent on impact energy, with defects easily retaining implanted hydrogen. MD results like those presented herein on unweathered surfaces are therefore most relevant to magnetic anomalies. As a result, we demonstrate the importance of lunar volatile models to account for the damage state of the substrate when modelling hydrogen diffusion, retention, and subsequent OH/water production.

Liam S. Morrissey↗

Uniform Distribution Technique for Neutrino Beam Scan Simulation

In Fermilab's neutrino facilities such as the Neutrinos at the Main Injector (NuMI) and the upcoming Long Baseline Neutrino Facility (LBNF), a proton beam strikes high-power target, producing positively and negatively charged pions and kaons. There is a need for detailed simulations in order to capture all particle interactions and beam propagation from protons on target to short-lived mesons decaying into muons and neutrinos. The generation of individual beam simulations is a resource-intensive and time-consuming process. In this paper, we describe a method through which many simulation samples with high statistics can be generated to study the effects of beam scan across a target for given beam configurations.

43 PARTICLE ACCELERATORS↗

Size Up or Size Down? National Analysis of Heat Pump Sizing and Impacts

Electrification of on-site fossil fuel combustion in buildings is recognized as a key component of achieving global greenhouse gas emissions targets. Air-source heat pumps are an efficient approach for space heating electrification. However, a major barrier to heat pump adoption is the high installation cost relative to furnaces and air conditioners. Because the installation cost of heat pumps - especially cold climate models - is dependent on their size, it is important to understand the distribution of heat pump sizes using different sizing methods and factors that impact heat pump sizing. In this study, we use sub-hourly physics simulations of 550,000 statistically representative dwelling units to analyze the distribution of heat pump size in the U.S. for three different efficiency levels of heat pumps, with and without building insulation and air sealing upgrades. We will present results from the national analysis of different factors affecting heat pump sizing along with the impact of sizing decisions on upgrade costs and operating costs.

electrification↗

Stress Rupture Life Reliability Measures for Composite Overwrapped Pressure Vessels

Composite Overwrapped Pressure Vessels (COPVs) are often used for storing pressurant gases onboard spacecraft. Kevlar (DuPont), glass, carbon and other more recent fibers have all been used as overwraps. Due to the fact that overwraps are subjected to sustained loads for an extended period during a mission, stress rupture failure is a major concern. It is therefore important to ascertain the reliability of these vessels by analysis, since the testing of each flight design cannot be completed on a practical time scale. The present paper examines specifically a Weibull statistics based stress rupture model and considers the various uncertainties associated with the model parameters. The paper also examines several reliability estimate measures that would be of use for the purpose of recertification and for qualifying flight worthiness of these vessels. Specifically, deterministic values for a point estimate, mean estimate and 90/95 percent confidence estimates of the reliability are all examined for a typical flight quality vessel under constant stress. The mean and the 90/95 percent confidence estimates are computed using Monte-Carlo simulation techniques by assuming distribution statistics of model parameters based also on simulation and on the available data, especially the sample sizes represented in the data. The data for the stress rupture model are obtained from the Lawrence Livermore National Laboratories (LLNL) stress rupture testing program, carried out for the past 35 years. Deterministic as well as probabilistic sensitivities are examined.

Murthy, Pappu L. N.↗

Fast Gaussian Process Estimation for Large-Scale In Situ Inference using Convolutional Neural Networks

Exascale computing will bring with it significant I/O limitations. One foreseeable consequence of such restrictions is that the user can save only a small fraction of complex simulation data to disk for subsequent analysis. An alternative is to fit statistical models to data in situ, that is, inside the simulation as it runs. This option requires extremely fast statistical estimation to avoid slowing down the simulation. Gaussian processes (GPs) have state-of-the-art predictive performance for modeling spatial data. However, standard estimation methods for GPs scale quite poorly to large data sets as parameter estimation requires inverting a covariance matrix to the size of the data set. In the presented work, we use a convolutional neural network (CNN) to predict the GP parameters for a spatial data set, from a simulation or otherwise, rather than optimize the parameters directly. Here, our presented case study models spatial data from E3SM, the Department of Energy’s Exascale climate model. The CNN is trained on synthetic data simulated from GP models with known parameters and then applied to data from the climate simulation. In the presented examples, the neural network scheme produces parameter estimates that compare well with standard methods such as maximum likelihood estimation in predictive performance but is obtained four orders of magnitude faster.

big data↗

Parametric description of intermittent probability distribution functions in solar wind and magnetohydrodynamic turbulence

ABSTRACT In this work, we find empirical evidence that the scale-dependent statistical properties of solar wind and magnetohydrodynamic (MHD) turbulence can be described in terms of a family of parametric probability distribution functions (PDFs) known as Normal Inverse Gaussian (NIG). Understanding these PDFs is one of the most important goals in turbulence theory, as they are inherently connected to the intermittent properties of solar wind turbulence. We investigate the properties of PDFs of Elsasser increments based on a large statistical sample from solar wind observations and high-resolution numerical simulations of MHD turbulence. In order to measure the PDFs and their corresponding properties, three experiments are presented: fast and slow solar wind for experimental data and a simulation of reduced MHD (RMHD) turbulence. Conditional statistics on a 23-yr-long sample of WIND data near 1 au and high-resolution pseudo-spectral simulation of steadily driven RMHD turbulence on a $2048^3$ mesh are used to construct scale-dependent PDFs. The empirical PDFs are fitted to NIG distributions, which depend on four free parameters. Our analysis shows that NIG distributions accurately capture the evolution of the PDFs, with scale-dependent parameters, from large scales characterized by a Gaussian distribution, turning to exponential tails within the inertial range and stretched exponentials at dissipative scales. We also show that empirically-measured NIG parameters exhibit well-defined scaling properties that are similar across the three empirical data sets, which may be indicative of universal behaviour.

Astronomy & Astrophysics↗

Automated Simulation Updates based on Flight Data

A statistically-based method for using flight data to update aerodynamic data tables used in flight simulators is explained and demonstrated. A simplified wind-tunnel aerodynamic database for the F/A-18 aircraft is used as a starting point. Flight data from the NASA F-18 High Alpha Research Vehicle (HARV) is then used to update the data tables so that the resulting aerodynamic model characterizes the aerodynamics of the F-18 HARV. Prediction cases are used to show the effectiveness of the automated method, which requires no ad hoc adjustments by the analyst.

Morelli, Eugene A.↗

The Atacama Cosmology Telescope: map-based noise simulations for DR6

Abstract The increasing statistical power of cosmic microwave background (CMB) datasets requires a commensurate effort in understanding their noise properties. The noise in maps from ground-based instruments is dominated by large-scale correlations, which poses a modeling challenge. This paper develops novel models of the complex noise covariance structure in the Atacama Cosmology Telescope Data Release 6 (ACT DR6) maps. We first enumerate the noise properties that arise from the combination of the atmosphere and the ACT scan strategy. We then prescribe a class of Gaussian, map-based noise models, including a new wavelet-based approach that uses directional wavelet kernels for modeling correlated instrumental noise. The models are empirical, whose only inputs are a small number of independent realizations of the same region of sky. We evaluate the performance of these models against the ACT DR6 data by drawing ensembles of noise realizations. Applying these simulations to the ACT DR6 power spectrum pipeline reveals a ∼ 20% excess in the covariance matrix diagonal when compared to an analytic expression that assumes noise properties are uniquely described by their power spectrum. Along with our public code,mnms, this work establishes a necessary element in the science pipelines of both ACT DR6 and future ground-based CMB experiments such as the Simons Observatory (SO).

Astronomy & Astrophysics↗