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 1,639 records · Page 91

A Blueprint for Demonstrating Quantum Supremacy with Superconducting Qubits

Long coherence times and high fidelity control recently achieved in scalable superconducting circuits paved the way for the growing number of experimental studies of many-qubit quantum coherent phenomena in these devices. Albeit full implementation of quantum error correction and fault tolerant quantum computation remains a challenge the near term pre-error correction devices could allow new fundamental experiments despite inevitable accumulation of errors. One such open question foundational for quantum computing is achieving the so called quantum supremacy, an experimental demonstration of a computational task that takes polynomial time on the quantum computer whereas the best classical algorithm would require exponential time and/or resources. It is possible to formulate such a task for a quantum computer consisting of less than a 100 qubits. The computational task we consider is to provide approximate samples from a non-trivial quantum distribution. This is a generalization for the case of superconducting circuits of ideas behind boson sampling protocol for quantum optics introduced by Arkhipov and Aaronson. In this presentation we discuss a proof-of-principle demonstration of such a sampling task on a 9-qubit chain of superconducting gmon qubits developed by Google. We discuss theoretical analysis of the driven evolution of the device resulting in output approximating samples from a uniform distribution in the Hilbert space, a quantum chaotic state. We analyze quantum chaotic characteristics of the output of the circuit and the time required to generate a sufficiently complex quantum distribution. We demonstrate that the classical simulation of the sampling output requires exponential resources by connecting the task of calculating the output amplitudes to the sign problem of the Quantum Monte Carlo method. We also discuss the detailed theoretical modeling required to achieve high fidelity control and calibration of the multi-qubit unitary evolution in the device. We use a novel cross-entropy statistical metric as a figure of merit to verify the output and calibrate the device controls. Finally, we demonstrate the statistics of the wave function amplitudes generated on the 9-gmon chain and verify the quantum chaotic nature of the generated quantum distribution. This verifies the implementation of the quantum supremacy protocol.

Kechedzhi, Kostyantyn↗

An Evaluation of Teleconnections Over the United States in an Ensemble of AMIP Simulations with the MERRA-2 Configuration of the GEOS Atmospheric Model

The atmospheric general circulation model that is used in NASA's Modern Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2) is evaluated with respect to the relationship between large-scale teleconnection patterns and daily temperature and precipitation over the United States (US) using a ten-member ensemble of simulations, referred to as M2AMIP. A focus is placed on four teleconnection patterns that are known to influence weather and climate in the US: El Nino Southern Oscillation, the Pacific Decadal Oscillation, the North Atlantic Oscillation, and the Pacific-North American Pattern. The monthly and seasonal indices associated with the patterns are correlated with daily temperature and precipitation statistics including: (i) monthly mean 2 m temperature and precipitation, (ii) the frequency of extreme temperature events at the 90th, 95th, and 99th percentiles, and (iii) the frequency and intensity of extreme precipitation events classified at the 90th, 95th, and 99th percentiles.Correlations obtained with M2AMIP data and thus the strength of teleconnections in the free-running model are evaluated through comparison against corresponding correlations computed from observations and from MERRA-2. Overall, the strongest teleconnections in all datasets occur during the winter and coincide with the largest agreement between the observations, MERRA-2, and M2AMIP. When M2AMIP does capture the correlation seen in observations, there is a tendency for the spatial extent to be exaggerated. The weakest agreement between the data sources, for all teleconnection patterns, is in the correlation with extreme precipitation; however there are discrepancies between the datasets in the number of days with at least 1 mm of precipitation: M2AMIP has too few days with precipitation in the Northwest and the Northern Great Plains and too many days in the Northeast. In JJA, M2AMIP has too few days with precipitation in the western two-thirds of the country and too many days with precipitation along the east coast.

JJA↗

Unifying Radiative Transfer Models in Computer Graphics and Remote Sensing, Part II: A Differentiable, Polarimetric Forward Model and Validation

The constellation of Earth-observing satellites continuously collects measurements of scattered radiance, which must be transformed into geophysical parameters in order to answer fundamental scientific questions about the Earth. Retrieval of these parameters requires highly flexible, accurate, and fast forward and inverse radiative transfer models. Existing forward models used by the remote sensing community are typically accurate and fast, but sacrifice flexibility by assuming the atmosphere or ocean is composed of plane-parallel layers. Monte Carlo forward models can handle more complex scenarios such as 3D spatial heterogeneity, but are relatively slower. We propose looking to the computer graphics community for inspiration to improve the statistical efficiency of Monte Carlo forward models and explore new approaches to inverse models for remote sensing. In Part 2 of this work, we demonstrate that Monte Carlo forward models in computer graphics are capable of sufficient accuracy for remote sensing by extending Mitsuba 3, a forward and inverse modeling framework recently developed in the computer graphics community, to simulate simple atmosphere-ocean systems and show that our framework is capable of achieving error on par with codes currently used by the remote sensing community on benchmark results.

Radiative transfer↗

Effect of simulated lunar impact on the survival of bacterial spores.

In order to test the effect of impact on organisms, the survival of bacterial spores after being propelled at high velocity in Pyrex and plastic beads into crushed basalt was measured. The beads were fired into sterilized canisters by both a conventional powder and a light gas gun. Results indicate that at the minimum (2.4 km/sec) lunar capture velocity, the number of colony forming units (CFUs) decreased by five orders of magnitude, and at 5.5 km/sec, statistically a more probable capture velocity, no CFUs were found. The decrease in CFUs observed with increasing velocity indicates that the spores were most probably killed by the impact.

Whitfield, O.↗

Data simulation studies for field programs, satellite system tests, and regional weather prediction

In the present paper, the status of the Drexel/NCAR regional numerical weather forecast system is summarized to show the extent of our experience and our immediate requirements for data within the United States on scales smaller than rawinsonde spacing (approximately 300 km). Preliminary experiments are described, and data assimilation problems are examined. It is shown that data requirements, assimilation methods, and forecast sensitivity differ for regional and large-scale models. Forecast verification statistics suggest that middle-troposphere temperature is well known but relative humidity is not. In some cases, regional forecasts require much better humidity data than available from rawinsondes. Numerics, lateral boundary conditions, model physics, and computing power are adequate for near real-time experimental regional prediction with horizontal mesh sizes down to 35 km. Subsynoptic data and assimilation are the primary challenges facing regional NWP.

Kreitzberg, C. W.↗

Cosmology with the Roman Space Telescope - Multiprobe Strategies

We simulate the scientific performance of the Nancy Grace Roman Space Telescope High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6-yr HLS Reference survey is currently envisioned to image 2000 deg2 in multiple bands to a depthof∼26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multiband photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g. weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper, we explore multiprobe strategies that can be implemented, given the telescope’s instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the HLS survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the HLS reference survey alone can achieve a standard dark energy FoM of>300 when including all probes. This assumes no information from external data sets, we assume a flat universe however, and includes realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community-driven effort to simulate and optimize the science return of the Roman Space Telescope.

Cosmological parameters↗

The use of a simulator in a synthetic aperture radar data processing tradeoff study.

A computer simulation of a synthetic aperture radar system is described. The objective in developing the simulator was to make it possible to assess the imaging performance of various configurations against a variety of target types. The simulation involves targets, the radar system itself, and the processor; several analysis and output options are provided. Targets may be internally generated or come from external sources, or externally produced (e.g., real) radar signals can be used. A number of parameters of the radar system can be varied. The characteristics of the processor can be adjusted to a considerable extent. Statistical analyses may be performed. Several output options are available, including imagery.

Lybanon, M.↗

Adaptive statistical pattern classifiers for remotely sensed data

A technique for the adaptive estimation of nonstationary statistics necessary for Bayesian classification is developed. The basic approach to the adaptive estimation procedure consists of two steps: (1) an optimal stochastic approximation of the parameters of interest and (2) a projection of the parameters in time or position. A divergence criterion is developed to monitor algorithm performance. Comparative results of adaptive and nonadaptive classifier tests are presented for simulated four dimensional spectral scan data.

Gonzalez, R. C.↗

Numerical study of sink-flow boundary layers

Direct numerical simulations of sink-flow boundary layers, with acceleration parameters K between 1.5 x 10 to the -6th and 3.0 x 10 to the -6th, are presented. The three-dimensional, time-dependent Navier-Stokes equations are solved numerically, using a spectral method, with about one million degrees of freedom. The flow is assumed to be statistically steady, and self-similar. A multiple-scale approximation and periodic conditions are applied to the fluctuations. The turbulence is studied using instantaneous and statistical results. Good agreement with the experiments of Jones and Launder (1972) is observed. The two effects of the favorable pressure gradient are to extend the logarithmic layer, and to alter the energy balance of the turbulence near the edge of the boundary layer. At low Reynolds number the logarithmic layer is shortened and slightly displaced, but wall-layer streaks are present even at the lowest values of R(theta) for which turbulence can be sustained. Large quiescent patches appear in the flow. Relaminarization occurs at K = 3.0 x 10 to the -6th, corresponding to a Reynolds number R(theta) of about 330.

Spalart, Philippe R.↗

Classification of multifrequency multilook synthetic aperture radar data

A technique for segmentation of multifrequency multilook intensity synthetic aperture radar (SAR) data into regions of homogeneous and similar backscatter characteristics is presented. Two statistical models, one for the multifrequency multilook SAR intensities and the other for the distribution of the region labels, are combined to obtain the a posteriori probability distribution function of the region labels given the multifrequency speckled intensities. As the maximization of the posterior distribution is computationally intensive, a suboptimal technique for region labeling is proposed. Several examples using both simulated and real multifrequency multilook imagery are given to illustrate the performance of the algorithm.

Rignot, Eric↗

Comparison of current Shuttle and pre-Challenger flight suit reach capability during launch accelerations

The Challenger accident prompted the creation of a crew escape system which replaced the former Launch Entry Helmet (LEH) ensemble with the current Launch Entry Suit (LES). However, questions were raised regarding the impact of this change on crew reach capability. This study addressed the question of reach capability and its effects on realistic ground-based training for Space Shuttle missions. Eleven subjects performed reach sweeps in both the LEH and LES suits during 1 and 3 Gx acceleration trials in the Brooks AFB centrifuge. These reach sweeps were recorded on videotape and subsequently analyzed using a 3D motion analysis system. The ANOVA procedure of the Statistical Analysis System program was used to evaluate differences in forward and overhead reach. The results showed that the LES provided less reach capability than its predecessor, the LEH. This study also demonstrated that, since there was no substantial difference between 1 and 3 Gx reach sweeps in the LES, realistic Shuttle launch training may be accomplished in ground based simulators.

Bagian, James P.↗

LAT Perspectives in Detection of High Energy Cosmic Ray Electrons

The GLAST Large Area Telescope (LAT) science objectives and capabilities in the detection of high energy electrons in the energy range from 20 GeV to approx. 1 TeV are presented. LAT simulations are used to establish the event selections. It is found that maintaining the efficiency of electron detection at the level of 30% the residual hadron contamination does not exceed 2-3% of the electron flux. LAT should collect approx. ten million of electrons with the energy above 20 GeV for each year of observation. Precise spectral reconstruction with high statistics presents us with a unique opportunity to investigate several important problems such as studying galactic models of IC radiation, revealing the signatures of nearby sources such as high energy cutoff in the electron spectrum, testing the propagation model, and searching for KKDM particles decay through their contribution to the electron spectrum.

Moiseev, Alexander↗

Uncertainty in the response of sudden stratospheric warmings and stratosphere- troposphere coupling to quadrupled CO2 concentrations in CMIP6 models

Major sudden stratospheric warmings (SSWs), vortex formation and final breakdown dates are key highlight points of the stratospheric polar vortex. These phenomena are relevant for stratosphere-troposphere coupling, which explains the interest in understanding their future changes. However, up to now, there is not a clear consensus on which projected changes to the polar vortex are robust, particularly in the Northern Hemisphere, possibly due to short data record or relatively moderate CO2 forcing. The new simulations performed under the Coupled Model Intercomparison Project, Phase 6, together with the long daily data requirements of the DynVarMIP project in preindustrial and quadrupled CO2 (4xCO2 ) forcing simulations provide a new opportunity to revisit this topic by overcoming the limitations mentioned above. In this study, we analyze this new model output to document the change, if any, in the frequency of SSWs under 4xCO2 forcing. Our analysis reveals a large disagreement across the models as to the sign of this change, even though most models show a statistically significant change. The models, however, are in good agreement as to the impact of SSWs over the North Atlantic: there is no indication of a change under 4xCO2 forcing. Over the Pacific, however, the change is more uncertain. Finally, the models show robust changes to the seasonal cycle in the stratosphere. Specifically, we find a longer duration of the stratospheric polar vortex, and thus a longer season of stratosphere-troposphere coupling.

sudden stratospheric warmings (SSWs)↗

Comparison of Clouds and Cloud Feedback between AMIP5 and AMIP6

We examine the changes in clouds and cloud feedback between Phase 5 (AMIP5) and Phase 6 (AMIP6) of the Atmospheric Model Intercomparison Project. Each model is perturbed by uniformly increasing the sea surface temperature by 4 K. The simulated cloud fraction, the perturbed states and cloud radiative kernels are used to derive cloud feedback in the shortwave (SW), longwave (LW) and their sum (Net). Compared to AMIP5, the cloud fraction in AMIP6 increases by 9.1%, while the perturbation leads to a 0.25% decrease. The Net cloud feedback at the top of the atmosphere (TOA) is almost double (174%). Statistical tests support that this change is mainly due to an increase in the surface SW cloud feedback caused by optically thick, middle and low clouds. The contribution of the atmospheric Net component (12%) stems from the increase in the atmospheric LW cloud feedback, likely to play a role in weakening (strengthening) the northward (southward) meridional atmospheric energy transport, while the opposite is true for the surface LW and Net cloud feedback in the meridional oceanic energy transport. The substantial increase in cloud feedback at the TOA primarily contributes to the higher climate sensitivity. The cloud feedback spread in AMIP6 is comparable to that in AMIP5.

Clouds↗

Uncertainty Quantification and Sensitivity Analysis in Process-Structure-Property Simulations for Laser Powder Bed Fusion Additive Manufacturing

Process variations and process-induced defects like porosity cause significant uncertainty in the microstructure and mechanical behavior of additively manufactured metals. Establishing process-structure-property (PSP) relationships and quantifying uncertainty using experiments alone is costly, especially for structural applications where mechanical allowables must be established for qualification and certification. This work presents a PSP simulation framework for laser powder bed fusion with a focus on uncertainty quantification through probabilistic calibration and multi-fidelity uncertainty propagation. Motivated by phenomenological input parameters related to grain nucleation and growth that are difficult to characterize, a global sensitivity analysis (GSA) is completed. Through GSA, the most important input parameters are identified based on their influence on the statistical distributions of microstructural metrics that influence mechanical behavior, including grain size, morphology, and crystallographic texture. The results provide insight on what experiments are necessary to quantify and control PSP uncertainties, particularly those associated with the more challenging input parameters.

additive manufacturing↗

Image statistics decoding for convolutional codes

It is a fact that adjacent pixels in a Voyager image are very similar in grey level. This fact can be used in conjunction with the Maximum-Likelihood Convolutional Decoder (MCD) to decrease the error rate when decoding a picture from Voyager. Implementing this idea would require no changes in the Voyager spacecraft and could be used as a backup to the current system without too much expenditure, so the feasibility of it and the possible gains for Voyager were investigated. Simulations have shown that the gain could be as much as 2 dB at certain error rates, and experiments with real data inspired new ideas on ways to get the most information possible out of the received symbol stream.

Pitt, G. H., III↗

On solar cosmic rays

The report summarizes the major results obtained over the last four years in the following areas: (1) transport in interplanetary space which includes solar particle observations combined with theoretical efforts to develop models which simulate the dynamics of the event; (2) transport in the solar corona and release of particles in the interplanetary medium; and (3) characteristics of source spectra and acceleration models. Modelling the propagation of solar cosmic rays through interplanetary space and establishing the magnitude and the radial and energy dependence of the diffusion coefficient K are considered. An important result due to statistical observation is that the coronal propagation must be independent of or weakly dependent on both rigidity and energy. Increasing evidence is found for prolonged injection at the sun.

Van Hollebeke, M. A. I.↗

Monte Carlo simulation of vibrational relaxation in nitrogen

Monte Carlo simulation of nonequilibrium vibrational relaxation of (rotationless) N2 using transition probabilities form an extended SSH theory is presented. For the range of temperatures considered, 4000-8000 K, the vibrational levels were found to be reasonably close to an equilibrium distribution at an average vibrational temperature based on the vibrational energy of the gas. As a result, they do not show any statistically significant evidence of the bottleneck observed in earlier studies of N2. Based on this finding, it appears that, for the temperature range considered, dissociation commences after all vibrational levels equilibrate at the translational temperature.

Olynick, David P.↗