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 253 records · Page 14

Uniform Beam Simulation Technique for Beam Scans and Machine Learning Studies at Fermilab

Fermilab's neutrino facilities, including NuMI and the upcoming LBNF, use proton beams to produce positively and negatively charged pions and kaons. Detailed simulations are necessary to study particle interactions and beam propagation. To efficiently analyze beam scan effects, we propose a technique to generate multiple simulation samples with high statistics. These samples can be used to develop beamline simulation based machine learning applications. In this technique, we generate a uniformly distributed single simulation data sample. We calculate Gaussian weights for each beam configurations and apply them to post-processing measurements. In this poster, we demonstrate the proposed simulation technique. This technique reduces simulation time and computing resources significantly.

Wickremasinghe, Athula↗

Orbital flight simulation utility software unit specifications, revision 1

The HP PASCAL source code defines the specifications for a Utility Software Unit (USU) designed to support orbital flight simulators such as MANHANDLE and GREAS (General Research and Engineering Analysis Simulator). Besides providing basic input/output, mathematical, matrix, quaternion, and statistical routines for such simulators, one of the primary functions of the USU is to isolate all system-dependent codes in one well-defined compartment, thereby facilitating transportation of the simulations from one computer to another. Directives are given for the PASCAL compilers of the HP-9000 Series 200 Pascal 3.0 and the HP-9000 Series 500 HP-UX 5.0 operating systems that produce a single file of relocatable code from four separate files of source code. Three of the source code files are common to both operating systems. The fourth source code file (utilspif.I) contains all of the system-dependent PASCAL code for the USU. A fifth file of source code written in C is required to interface utilspif.I with the HP-UX I/O package. The Pascal 3.0 compiler directives and the driver source code for a unit rest program and counterparts for the HP-UX 5.0 operating system are given. The major portion of the unit test program source code is common to both operating systems. Unit test results from the Pascal 3.0 operating system and results from the HP-UX operating system are given.

Wilson, S. W.↗

Effects of Motion Fidelity on Pilot Performance Using Mixed-Reality Visuals

This paper investigates the effects of simulator motion fidelity on pilot performance using mixed-reality visuals. In NASA’s Vertical Motion Simulator, eleven pilots flew a lift-plus-cruise electric vertical takeoff and landing vehicle in a vertiport approach and landing task with four different motion conditions: no motion, small hexapod motion, large hexapod motion, and full motion. The simulator motion condition statistically significantly affected: 1) the severity of cybersickness, 2) motion ratings, 3) final approach and take-off boundary crossing groundspeed, 4) longitudinal touchdown point, and 5) battery time remaining after landing. The findings suggest that motion systems are essential in simulators with virtual- or mixed-reality visuals when it is required to have pilot flying performance similar to that in real flight.

mixed reality↗

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗

Statistical analysis of atmospheric turbulence about a simulated block building

An array of towers instrumented to measure the three components of wind speed was used to study atmospheric flow about a simulated block building. Two-point spacetime correlations of the longitudinal velocity component were computed along with two-point spatial correlations. These correlations are in good agreement with fundamental concepts of fluid mechanics. The two-point spatial correlations computed directly were compared with correlations predicted by Taylor's hypothesis and excellent agreement was obtained at the higher levels which were out of the building influence. The correlations fall off significantly in the building wake but recover beyond the wake to essentially the same values in the undisturbed, higher regions.

Steely, S. L., Jr.↗

Effects of Motion Fidelity on Pilot Performance Using Mixed-Reality Visuals

This paper investigates the effects of simulator motion fidelity on pilot performance using mixed-reality visuals. In NASA’s Vertical Motion Simulator, eleven pilots flew a lift-plus-cruise electric vertical takeoff and landing vehicle in a vertiport approach and landing task under four different motion conditions: no motion, small hexapod motion, large hexapod motion, and full motion. The simulator motion condition statistically significantly affected: 1) the severity of cybersickness, 2) motion ratings, 3) final approach and take-off boundary crossing groundspeeds, 4) longitudinal touchdown point, and 5) battery time remaining after landing. These results suggest that the use of motion systems in simulators with virtual- or mixed-reality visuals is warranted.

simulation↗

Accurate low-cost methods for performance evaluation of cache memory systems

Methods of simulation based on statistical techniques are proposed to decrease the need for large trace measurements and for predicting true program behavior. Sampling techniques are applied while the address trace is collected from a workload. This drastically reduces the space and time needed to collect the trace. Simulation techniques are developed to use the sampled data not only to predict the mean miss rate of the cache, but also to provide an empirical estimate of its actual distribution. Finally, a concept of primed cache is introduced to simulate large caches by the sampling-based method.

Laha, Subhasis↗

Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. Validation on simulations

Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end simulationbased analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogs of the Dark Energy Survey (DES). Here, we implement: 1) second and third moments; 2) wavelet phase harmonics (WPH); 3) the scattering transform (ST). Our analysis is fully based on simulations, it spans a space of seven $νw$CDM cosmological parameters, and it forward models the most relevant sources of systematics of the data (masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration). We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a likelihood-free-inference approach. We validate the pipeline extensively, and we find that WPH exhibits the strongest performance when combined with second moments, followed by ST, and then by third moments. The combination of all the different statistics further enhances constraints with respect to second moments, up to 25 percent, 15 percent, and 90 percent for S 8 , Ω m , and the figure-of-merit FoM S8;Ωm , respectively. We further find that non-Gaussian statistics improve constraints on w and on the amplitude of intrinsic alignment with respect to second moments constraints. The methodological advances presented here are suitable for application to Stage IV surveys from Euclid, Rubin-LSST, and Roman with additional validation on mock catalogs for each survey. In a companion paper we present an application to DES Year 3 data.

79 ASTRONOMY AND ASTROPHYSICS↗

A Study of the Response of Deep Tropical Clouds to Mesoscale Processes: Modeling Strategies and Simulations of TOGA-COARE Convective Systems - Part 1

Interactions between deep tropical clouds over the western Pacific warm pool and the larger-scale environment are key to understanding climate change. Cloud models are an extremely useful tool in simulating and providing statistical information on heat and moisture transfer processes between cloud systems and the environment, and can therefore be utilized to substantially improve cloud parameterizations in climate models. In this paper, the Goddard Cumulus Ensemble (GCE) cloud-resolving model is used in multi-day simulations of deep tropical convective activity over the Tropical Ocean-Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). Large-scale temperature and moisture advective tendencies, and horizontal momentum from the TOGA-COARE Intensive Flux Array (IFA) region, are applied to the GCE version which incorporates cyclical boundary conditions. Sensitivity experiments show that grid domain size produces the largest response to domain-mean temperature and moisture deviations, as well as cloudiness, when compared to grid horizontal or vertical resolution, and advection scheme. It is found that a minimum grid-domain size of 500 km is needed to adequately resolve the convective cloud features. The control experiment shows that the atmospheric heating and moistening is primarily a response to cloud latent processes of condensation/evaporation, and deposition/sublimation, and to a lesser extent, melting of ice particles. Air-sea exchange of heat and moisture is found to be significant, but of secondary importance, while the radiational response is small. The simulated rainfall and atmospheric heating and moistening, agrees well with observations, and performs favorably to other models simulating this case.

Johnson, Daniel E.↗

Astrometric Detection of Extrasolar Planets: Results of a Feasibility Study with the Palomar 5 Meter Telescope

The detection of extrasolar planets around stars like the Sun remains an important goal of astronomy. We present results from Palomar 5 m observations of the open cluster NGC 2420 in which we measure some of the sources of noise that will be present in an astrometric search for extrasolar planets. This is the first time that such a large aperture has been used for high-precision astrometry. We find that the atmospheric noise is 150 micro-arcsec hr(exp 1/2) across a 90 sec field of view and that differential chromatic refraction (DCR) can be calibrated to 128 micro-arcsec for observations within 1 hr of the meridian and 45 deg of zenith. These results confirm that a model for astrometric measurements can be extrapolated to large apertures. We demonstrate, based upon these results, that a large telescope achieves the sensitivity required to perform a statistically significant search for extra solar planets. We describe an astrometric technique to detect planets, the astrometric signals expected, the role of reference stars, and the sources of measurement noise: photometric noise, atmospheric motion between stars, sky background, instrumental noise, and DCR. For the latter, we discuss a method to reduce the noise further to 66 micro-arcsecond for observations within 1 hr of the meridian and 45 deg of zenith. We discuss optimal lists of target stars taken from the latest Gliese & Jahreiss catalog of nearby stars with the largest potential astrometric signals, declination limits for both telescope accessibility and reduced DCR, and galactic latitude limits for a sufficiant number of reference stars. Two samples are described from which one can perform statistically significant searches for gas giant planets around nearby stars. One sample contains 100 "solar class" stars with an average stellar mass of 0.82 solar mass; the other maximizes the number of stars, 574, by searching mainly low-mass M stars. We perform Monte Carlo simulations of the statistical significance of the expected results by using measured and estimated noise quantities. We show the semimajor axis parameter spaces that are searched for each star and how an increase in the length of the observing program expands these spaces. The search over semimajor axis parameter space relates to the theory of gas giant planet formation.

Pravdo, Steven H.↗

RCEMIP-II: mock-Walker simulations as phase II of the radiative–convective equilibrium model intercomparison project

Abstract. The radiative–convective equilibrium (RCE) model intercomparison project (RCEMIP) leveraged the simplicity of RCE to focus attention on moist convective processes and their interactions with radiation and circulation across a wide range of model types including cloud-resolving models (CRMs), general circulation models (GCMs), single-column models, global cloud-resolving models, and large-eddy simulations. While several robust results emerged across the spectrum of models that participated in the first phase of RCEMIP (RCEMIP-I), two points that stand out are (1) the strikingly large diversity in simulated climate states and (2) the strong imprint of convective self-aggregation on the climate state. However, the lack of consensus in the structure of self-aggregation and its response to warming is a barrier to understanding. Gaining a deeper understanding of convective aggregation and tropical climate will require reducing the degrees of freedom with which convection can vary. Therefore, we propose phase II of RCEMIP (RCEMIP-II) that utilizes a prescribed sinusoidal sea surface temperature (SST) pattern to provide a constraint on the structure of convection and move one critical step up the model hierarchy. This so-called “mock-Walker” configuration generates features that resemble observed tropical circulations. The specification of the mock-Walker protocol for RCEMIP-II is described, along with example results from one CRM and one GCM. RCEMIP-II will consist of five required simulations: three simulations with the same three mean SSTs as in RCEMIP-I but with an SST gradient and two additional simulations at one of the mean SSTs with different values of the SST gradients. We also test the sensitivity to the imposed SST gradient and the domain size. Under weak SST gradients, unforced self-aggregation emerges across the entire domain, similar to what was found in RCEMIP. As the SST gradient increases, the convective region narrows and is more confined to the warmest SSTs. At warmer mean SSTs and stronger SST gradients, low-frequency variability in the convective aggregation emerges, suggesting that simulations of at least 200 d may be needed to achieve robust equilibrium statistics in this configuration. Simulations with different domain sizes generally have similar mean statistics and convective structures, depending on the value of the SST gradient. The prescribed SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I, which enables comparison between the two; however, we also welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.

Geology↗

Atmospheric turbulence simulation techniques with application to flight analysis

Statistical modeling of atmospheric turbulence is discussed. The statistical properties of atmospheric turbulence, in particular the probability distribution, the spectra, and the coherence are reviewed. Different atmospheric turbulence simulation models are investigated, and appropriate statistical analyses are carried out to verify their validity. The models for simulation are incorporated into a computer model of aircraft flight dynamics. Statistical results of computer simulated landings for an aircraft having characteristics of a DC-8 are reported for the different turbulence simulation techniques. The significance of various degrees of sophistication in the turbulence simulation techniques on the landing performance of the aircraft is discussed.

Wang, S. T.↗

Simple Building Calculator

Whole building energy modeling is a powerful tool for analyzing energy use in buildings. For large buildings benefits of modeling services can easily be quantified due to the short payback period of the implemented measures. For small buildings, the upfront costs of developing robust energy models often deter building owners from investing in energy modeling services. To address this shortcoming, a ”Simple Building Calculator” (SBC) was developed. SBC uses pre-simulated results of a range of common measures over a wide range of efficiency inputs. It combines whole building simulation results with statistical modeling techniques to predict energy impact of measures very quickly. In this paper we present the modeling methodology used to develop the data supporting SBC.

Building Energy Simulation, commercial building, E↗

Soil carbon storage capacity of drylands under altered fire regimes

The determinants of fire-driven changes in soil organic carbon (SOC) across broad environmental gradients remains unclear, especially in global drylands. Here we combined datasets and field sampling of fire-manipulation experiments to evaluate where and why fire changes SOC and compared our statistical model to simulations from ecosystem models. Drier ecosystems experienced larger relative changes in SOC than humid ecosystems—in some cases exceeding losses from plant biomass pools—primarily explained by high fire-driven declines in tree biomass inputs in dry ecosystems. Many ecosystem models underestimated the SOC changes in drier ecosystems. Upscaling our statistical model predicted that soils in savannah–grassland regions may have gained 0.64 PgC due to net-declines in burned area over the past approximately two decades. Consequently, ongoing declines in fire frequencies have probably created an extensive carbon sink in the soils of global drylands that may have been underestimated by ecosystem models.

54 ENVIRONMENTAL SCIENCES↗

First Muon Neutrino on Argon Cross Section Measurements With an Off-Axis Beam

While over ninety years have passed since neutrinos were first predicted in 1930 by Wolfgang Pauli, much about them is still unknown [20]. How many types of neutrinos there are, how precisely neutrinos change flavor, and how neutrinos interact with various types of matter are largely unanswered questions that have implications for fields beyond particle physics, such as cosmology and astrophysics. The final question can be answered partially by differential cross section measurements, two of which constitute the primary result of this thesis. The NeUtrinos at the Main Injector (NuMI) Beam at Fermilab in Batavia, IL provides a source of muon neutrinos and muon antineutrinos to the MicroBooNE experiment, an 85 ton liquid argon time projection chamber (LArTPC). The technology allows for high precision particle tracking and energy reconstruction that can be used to make cross section measurements. Necessary work for making these measurements for which the author of this thesis made an essential contribution is presented in detail in References [21] [22] [11] [23]. The measurement presented in this thesis is made with 2.187e20 Protons on Target (POT), the standard metric for an amount of beam exposure, of data when the beam was in forward horn current mode. Single-differential cross section measurements in terms of muon candidate kinetic energy and direction are made, and their statistical fits with simulation are 5.8 (12 d.o.f.) and 11.4 (10 d.o.f.), respectively. This indicates that the muon kinematics for muon neutrino and muon antineutrino charged-current (CC) events primarily in the energy range [0-3] GeV can be trusted for oscillation experiments, including DUNE, because of the agreement to within the associated systematic and statistical uncertainties between the measurement and the model in every bin of both single-differential cross sections. This result adds to the corresponding result for the Booster Neutrino Beam (BNB) because it covers a completely different energy region than that measurement, a double-differential cross section in terms of muon momentum and direction [14]. This result contributes to reducing the uncertainty in near-far event rate comparisons in long baseline neutrino experiments and also adds to existing knowledge of the neutrino-nucleus interaction by describing the outgoing muon kinematics in interactions of muon neutrinos and muon antineutrinos. This result is the first-ever measurement of differential cross sections of muon neutrinos and muon antineutrinos in liquid argon made with an off-axis beam.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Direct Numerical Simulation of a Weakly Stratified Turbulent Wake

Direct numerical simulation (DNS) is used to investigate a time-dependent turbulent wake evolving in a stably stratified background. A large initial Froude number is chosen to allow the wake to become fully turbulent and axisymmetric before stratification affects the spreading rate of the mean defect. The uncertainty introduced by the finite sample size associated with gathering statistics from a simulation of a time-dependent flow is reduced, compared to earlier simulations of this flow. The DNS reveals the buoyancy-induced changes to the turbulence structure, as well as to the mean-defect history and the terms in the mean-momentum and turbulence-kinetic-energy budgets, that characterize the various states of this flow - namely the three-dimensional (essentially unstratified), non-equilibrium (or 'wake-collapse') and quasi-two-dimensional (or 'two-component') regimes observed elsewhere for wakes embedded in both weakly and strongly stratified backgrounds. The wake-collapse regime is not accompanied by transfer (or 'reconversion') of the potential energy of the turbulence to the kinetic energy of the turbulence, implying that this is not an essential feature of stratified-wake dynamics. The dependence upon Reynolds number of the duration of the wake-collapse period is demonstrated, and the effect of the details of the initial/near-field conditions of the wake on its subsequent development is examined.

Redford, J. A.↗

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗

Orbital flight simulation utility software unit specifications

The HP PASCAL source code contained in pages 6 through 104 was developed for the Mission Planning and Analysis Division (MPAD) and takes the place of detailed flow charts defining the specifications for a Utility Software Unit designed to support orbital flight simulators such as MANHANDLE and GREAS (General Research and Engineering Analysis Simulator). Besides providing basic input/output, mathematical, vector, matrix, quaternion, and statistical routines for such simulators, one of the primary functions of the Utility Software Unit is to isolate all system-dependent code in one well-defined compartment, thereby facilitating transportation of the simulations from one computer to another. Directives to the PASCAL compilers of the HP-9000 Series 200 PASCAL 3.0 operating system and the HP-9000 Series 500 HP-UX 5.0 operations systems are also provided.

Wilson, S. W.↗