SEARCH · Engineering Papers

Results for “variational Monte Carlo”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Three-dimensional rarefied flow simulations for the Aeroassist Flight Experiment vehicle

Hypersonic rarefied flow about the Aeroassist Flight Experiment vehicle has been investigated using a three-dimensional direct simulation Monte Carlo method. Calculations are performed for the transitional flows encountered during the vehicle's atmospheric entry for altitudes of 110 and 100 km with an entry velocity of 9.9 km/s. The simulations are performed using a five-species reacting gas model that account for rotational and vibrational internal energies. The solutions indicate that dissociation is important at altitudes of 110 km and below. Results are presented for surface pressures, convective heating, flowfield structure, and aerodynamic coefficient variations with altitude.

Celenligil, M. Cevdet↗

NASA NTRS2008

Uncertainty and Sensitivity Analyses of a Two-Parameter Impedance Prediction Model

This paper presents comparisons of predicted impedance uncertainty limits derived from Monte-Carlo-type simulations with a Two-Parameter (TP) impedance prediction model and measured impedance uncertainty limits based on multiple tests acquired in NASA Langley test rigs. These predicted and measured impedance uncertainty limits are used to evaluate the effects of simultaneous randomization of each input parameter for the impedance prediction and measurement processes. A sensitivity analysis is then used to further evaluate the TP prediction model by varying its input parameters on an individual basis. The variation imposed on the input parameters is based on measurements conducted with multiple tests in the NASA Langley normal incidence and grazing incidence impedance tubes; thus, the input parameters are assigned uncertainties commensurate with those of the measured data. These same measured data are used with the NASA Langley impedance measurement (eduction) processes to determine the corresponding measured impedance uncertainty limits, such that the predicted and measured impedance uncertainty limits (95% confidence intervals) can be compared. The measured reactance 95% confidence intervals encompass the corresponding predicted reactance confidence intervals over the frequency range of interest. The same is true for the confidence intervals of the measured and predicted resistance at near-resonance frequencies, but the predicted resistance confidence intervals are lower than the measured resistance confidence intervals (no overlap) at frequencies away from resonance. A sensitivity analysis indicates the discharge coefficient uncertainty is the major contributor to uncertainty in the predicted impedances for the perforate-over-honeycomb liner used in this study. This insight regarding the relative importance of each input parameter will be used to guide the design of experiments with test rigs currently being brought on-line at NASA Langley.

Jones, M. G.↗

DOE OSTI2025

Observation of Cosmic-Ray Anisotropy in the Southern Hemisphere with 12 yr of Data Collected by the IceCube Neutrino Observatory

We analyzed the 7.92 × 10 11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗

DOE OSTI2023

VarIabiLity seLection of AstrophysIcal sources iN PTF (VILLAIN): I. Structure function fits to 71 million objects

Light-curve variability is well-suited to characterizing objects in surveys with high cadence and a long baseline. This is especially relevant in view of the large datasets to be produced by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We aim to determine variability parameters for objects in the Palomar Transient Factory (PTF) and explore differences between quasars (QSOs), stars, and galaxies. We relate variability and colour information in preparation for future surveys. We fit joint likelihoods to structure functions (SFs) of 71 million PTF light curves with a Markov chain Monte Carlo method. For each object, we assume a power-law SF and extract two parameters: the amplitude on timescales of one year, A, and a power-law index, γ. With these parameters and colours in the optical (Pan-STARRS1) and mid-infrared (WISE), we identify regions of parameter space dominated by different types of spectroscopically confirmed objects from SDSS. Candidate QSOs, stars, and galaxies are selected to show their parameter distributions. QSOs show high-amplitude variations in the R band, and the highest γ values. Galaxies have a broader range of amplitudes and their variability shows relatively little dependency on timescale. With variability and colours, we achieve a photometric selection purity of 99.3% for QSOs. Even though hard cuts in monochromatic variability alone are not as effective as seven-band magnitude cuts, variability is useful in characterizing object subclasses. Through variability, we also find QSOs that were erroneously classified as stars in the SDSS. We discuss perspectives and computational solutions in view of the upcoming LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Beam-induced backgrounds measured in the ATLAS detector during local gas injection into the LHC beam vacuum

Inelastic beam-gas collisions at the Large Hadron Collider (LHC), within a few hundred metres of the ATLAS experiment, are known to give the dominant contribution to beam backgrounds. These are monitored by ATLAS with a dedicated Beam Conditions Monitor (BCM) and with the rate of fake jets in the calorimeters. These two methods are complementary since the BCM probes backgrounds just around the beam pipe while fake jets are observed at radii of up to several metres. In order to quantify the correlation between the residual gas density in the LHC beam vacuum and the experimental backgrounds recorded by ATLAS, several dedicated tests were performed during LHC Run 2. Local pressure bumps, with a gas density several orders of magnitude higher than during normal operation, were introduced at different locations. The changes of beam-related backgrounds, seen in ATLAS, are correlated with the local pressure variation. In addition the rates of beam-gas events are estimated from the pressure measurements and pressure bump profiles obtained from calculations. Using these rates, the efficiency of the ATLAS beam background monitors to detect beam-gas events is derived as a function of distance from the interaction point. These efficiencies and characteristic distributions of fake jets from the beam backgrounds are found to be in good agreement with results of beam-gas simulations performed with theFluka Monte Carlo programme.

43 PARTICLE ACCELERATORS↗

NASA NTRS1993

DSMC and continuum analyses of low-density nozzle flow

Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chan-Hong↗

NASA NTRSDate not supplied

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗

NASA NTRS1994

Possible consequences of absence of Jupiters in planetary systems

The formation of the gas giant planets Jupiter and Saturn probably required the growth of massive approximately 15 Earth-mass cores on a time scale shorter than the approximately 10(exp 7) time scale for removal of nebular gas. Relatively minor variations in nebular parameters could preclude the growth of full-size gas giants even in systems in which the terrestrial planet region is similar to our own. Systems containing 'failed Jupiters,' resembling Uranus and Neptune in their failure to capture much nebular gas, would be expected to contain more densely populated cometary source regions. They will also eject a smaller number of comets into interstellar space. If systems of this kind were the norm, observation of hyperbolic comets would be unexpected. Monte Carlo calculations of the orbital evolution of region of such systems (the Kuiper belt) indicate that throughout Earth history the cometary impact flux in their terrestrial planet regions would be approximately 1000 times greater than in our Solar System. It may be speculated that this could frustrate the evolution of organisms that observe and seek to understand their planetary system. For this reason our observation of these planets in our Solar System may tell us nothing about the probability of similar gas giants occurring in other planetary systems. This situation can be corrected by observation of an unbiased sample of planetary systems.

Wetherill, George W.↗

DOE OSTI2022

Influence of temporal variations in plasma conditions on the electric potential near self-organized dust chains

Self-organization of dust grains into stable filamentary dust structures (or “chains”) largely depends on dynamic interactions between individual charged dust grains and complex electric potential arising from the distribution of charges within a local plasma environment. Recent studies have shown that the positive column of the gas discharge plasma in the Plasmakristall-4 (PK-4) experiment at the International Space Station supports the presence of fast-moving ionization waves, which lead to variations of plasma parameters by up to an order of magnitude from the average background values. The highly variable environment resulting from ionization waves may have interesting implications for the dynamics and self-organization of dust particles, particularly concerning the formation and stability of dust chains. Here, we investigate the electric potential surrounding dust chains in the PK-4 experiment by employing a molecular dynamics model of the dust and ions with boundary conditions supplied by a particle-in-cell with Monte Carlo collision simulation of the ionization waves. The model is used to examine the effects of the plasma conditions within different regions of the ionization wave and compare the resulting dust structure to that obtained by employing the time-averaged plasma conditions. Here, the comparison between simulated dust chains and experimental data from the PK-4 experiment shows that the time-averaged plasma conditions do not accurately reproduce observed results for dust behavior, indicating that more careful treatment of plasma conditions in the presence of ionization waves is required. It is further shown that commonly used analytic forms of the electric potential do not accurately describe the electric potential near charged dust grains under these plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1988

Direct simulation of three-dimensional flow about the AFE vehicle at high altitudes

Three-dimensional hypersonic rarefied flow about the Aeroassist Flight Experiment (AFE) vehicle was studied using the direct simulation Monte Carlo (DSMC) technique. Results are presented for the transitional flow regime encountered between 120 and 200 km altitudes with a reentry velocity of 9.92 km/s. In the simulations, a five-species reacting real-gas model that accounts for internal energies (rotational and vibrational) is used. The results indicate that the transitional effects are significant even at an altitude of 200 km and influence the overall vehicle aerodynamics. For the cases considered, the aerodynamic coefficients, surface pressures, convective heating, and flow field structure variations with rarefaction effects are presented.

Celenligil, M. Cevdet↗

NASA NTRS1989

Direct simulation of three-dimensional flow about the AFE vehicle at high altitudes

Celenligil, M. Cevdet↗

NASA NTRS1994

Monte Carlo prediction of radiative heat transfer in inhomogeneous, anisotropic, nongray media

A Monte Carlo solution technique has been formulated to predict the radiative heat transfer in three-dimensional, inhomogeneous participating media which exhibit spectrally dependent emission and absorption and anisotropic scattering. Details of the technique and selected numerical sensitivities are discussed. The technique was applied to a problem involving a medium composed of a gas mixture of carbon dioxide and nitrogen and suspended carbon particles. A homogeneous medium was modeled to examine the effect of total pressure and carbon-particle concentration on radiative heat transfer. Variation in total pressure, over the range studied, had minimal effect on the amount of heat radiated to the enclosure walls and on the radiative-flux distribution within the medium. Increases in the carbon particle concentration produced significantly higher heat fluxes at the boundaries and altered the radiative flux distribution. The technique was then applied to an inhomogeneous medium to examine effects of specific temperature and carbon particle concentration distributions on radiative heat transfer. For the inhomogeneous conditions examined, the largest radiative flux divergence occurs near the center of the medium and the regions near some enclosure walls act as energy sinks.

HEAT FLUXES↗

NASA NTRS2004

Angle-of-Attack-Modulated Terminal Point Control for Neptune Aerocapture

An aerocapture guidance algorithm based on a calculus of variations approach is developed, using angle of attack as the primary control variable. Bank angle is used as a secondary control to alleviate angle of attack extremes and to control inclination. The guidance equations are derived in detail. The controller has very small onboard computational requirements and is robust to atmospheric and aerodynamic dispersions. The algorithm is applied to aerocapture at Neptune. Three versions of the controller are considered with varying angle of attack authority. The three versions of the controller are evaluated using Monte Carlo simulations with expected dispersions.

Queen, Eric M.↗

DOE OSTI2025

On the Interplay of Electronic and Lattice Screening on Exciton Binding in Two-Dimensional Lead Halide Perovskites

We use path integral Monte Carlo to study the energetics of excitons in layered, hybrid organic–inorganic perovskites in order to elucidate the relative contributions of dielectric confinement and electron–phonon coupling. While the dielectric mismatch between polar perovskite layers and nonpolar ligand layers significantly increases the exciton binding energy relative to their three-dimensional bulk crystal counterparts, formation of exciton polarons attenuates this effect. The contribution from polaron formation is found to be a nonmonotonic function of the lead halide layer thickness, which is clarified by a general variational theory. Accounting for both of these effects provides a description of exciton binding energies in good agreement with experimental measurements. By studying isolated layers and stacked layered crystals of various thicknesses, with ligands of varying polarity, we provide a systematic understanding of the excitonic behavior of this class of materials and how to engineer their photophysics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1972

Degree and plane of polarization of multiple scattered light. 2: Earth's atmosphere with aerosols

The degree of polarization, as well as the direction of the plane of polarization, were calculated by a Monte Carlo method for the reflected and transmitted photons from the earth's atmosphere. The solar photons were observed during multiple collisions with aerosols and the Rayleigh scattering centers in the atmosphere. The aerosol number density, as well as the ratio of aerosol to Rayleigh scattering, varies with height. The proportion of aerosol to Rayleigh scattering was appropriately chosen at each wavelength 0.4 microns and 0.7 microns; ozone absorption was included where appropriate. Three different aerosol number densities were used to study the effects of aerosol variations. Results are given for a solar zenith angle of 81.37 deg and a surface albedo of zero. The polarization of the reflected and transmitted photons was found to be sensitive to the amount of aerosols in the atmosphere at certain angles of observation.

Plass, G. N.↗

NASA NTRS1979

Monte Carlo simulation of lunar megaregolith and implications

A realistic Monte Carlo model closely simulating the evolution of the lunar megaregolith over a large area of 67 million sq. km of the front surface of the moon is presented. Craters larger than 100 km in diameter observed over the entire surface of the moon and those less than 100 km lying in the referenced area are included in the simulation. A total of 21,664 craters are processed. The model predicts the average thickness of the megaregolith to be about 1.9-2.0 km. Curves for the variation of the regolith thickness across the simulated area are given and show that about 50% of the area is covered with regolith less than 1 km thick. The model produces crater structures similar to the ones observed in the lunar highlands, it partially supports the layering theory for crater structures that the variations in strength and density of target materials may be responsible for the observed differences in the morphologies of lunar craters, and rules out the possibility that all craters when formed are bowl-shaped with a fixed depth/diameter ratio characteristic of small craters.

Aggarwal, H. R.↗

DOE OSTI2026

Role of electron correlation on the adenine dimer interaction for non-equilibrium geometries: a benchmark Quantum Monte Carlo study

The accurate description of non-covalent interactions is critical for understanding the structure, dynamics, and eventual function of biomolecules. The adenine dimer serves as a benchmark system for computational methods due to its role in nucleic acid structures and its rich conformational landscape. In this study, we employ benchmark diffusion quantum Monte Carlo (DMC) methods to investigate the relative energies and role of electron correlation on a set of adenine dimer conformations generated via a search of the potential energy landscape using the global optimizer algorithm. Relative DMC energies are compared against a wide range of density functional theory (DFT) approximation results. We find that although most of the DFT functionals perform well for low-energy structures, their accuracy varies significantly for higher-energy conformations, including stacked and T-shaped structures. A large fraction of the variation is due to the treatment of the van der Waals interaction. BLYP, B3LYP, and PBE0 significantly improve with added D4 dispersion, while the recent r2SCAN-D4 and ωB97M-V functionals show the least scatter and closest agreement with the DMC. These findings highlight the delicate nature of these interactions in biomolecular systems and provide guidance for simulations of their structure and dynamics and for the development of machine learned interatomic potentials.

Washburn, Laurel [ORNL] (ORCID:0000000324179335)↗