SEARCH · Engineering Papers

Results for “Direct simulation 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 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Model-free Rayleigh weight from x-ray Thomson scattering measurements

X-ray Thomson scattering (XRTS) has emerged as a powerful tool for the diagnostics of matter under extreme conditions. In principle, it gives one access to important system parameters such as the temperature, density, and ionization state, but the interpretation of the measured XRTS intensity usually relies on theoretical models and approximations. In this context, a key property is given by the Rayleigh weight that describes the electronic localization around the ions. Here, we show that it is possible to extract the Rayleigh weight directly from the experimental data without the need for any model calculations or simulations. As a practical application, we consider an experimental measurement of strongly compressed Be at the National Ignition Facility [Döppner et al., Nature 618, 270–275 (2023)]. We demonstrate that experimental results for the Rayleigh weight open up new avenues for the interpretation of XRTS experiments by matching the measurement with ab initio simulations such as density functional theory or path integral Monte Carlo. Interestingly, this new procedure leads to significantly lower density compared to previously used chemical models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Energy Gap Closure of Crystalline Molecular Hydrogen with Pressure

Here, we study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by coupled electron ion Monte Carlo methods. Nuclear zero point effects cause a large reduction in the gap (~2 eV). Depending on the structure, the fundamental indirect gap closes between 380 and 530 GPa for ideal crystals and 330–380 GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to ~450–500 GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1973

Decision-directed detector for overlapping PCM/NRZ signals.

A decision-directed (DD) technique for the detection of overlapping PCM/NRZ signals in the presence of white Gaussian noise is investigated. The performance of the DD detector is represented by probability of error Pe versus input signal-to-noise ratio (SNR). To examine how much improvement in performance can be achieved with this technique, Pe's with and without DD feedback are evaluated in parallel. Further, analytical results are compared with those found by Monte Carlo simulations. The results are in good agreement.

Wang, C. D.↗

NASA NTRS1981

Monte Carlo simulation of charged particles upstream of the earth's bow shock

A simple model of the earth's bow shock is presented and compared with ISEE-1 observations of diffuse ions upstream of the shock. Ion behavior is modeled in a test particle simulation using Monte Carlo techniques for a plane parallel shock at which cold solar wind ions are heated sufficiently to allow diffusion back upstream. Fitting of the model to the observations indicates that the diffusion coefficient must increase faster than linearly with velocity, and provides an upper limit to the acceleration efficiency defined as the proportion of incoming solar wind ions accelerated and observed upstream. Results suggest that astrophysical shocks are capable of accelerating ions directly from the thermal plasma without the need for a distinct preenergizing injection stage.

Ellison, D. C.↗

NASA NTRS1982

Monte Carlo studies of ocean wind vector measurements by SCATT: Objective criteria and maximum likelihood estimates for removal of aliases, and effects of cell size on accuracy of vector winds

The scatterometer on the National Oceanic Satellite System (NOSS) is studied by means of Monte Carlo techniques so as to determine the effect of two additional antennas for alias (or ambiguity) removal by means of an objective criteria technique and a normalized maximum likelihood estimator. Cells nominally 10 km by 10 km, 10 km by 50 km, and 50 km by 50 km are simulated for winds of 4, 8, 12 and 24 m/s and incidence angles of 29, 39, 47, and 53.5 deg for 15 deg changes in direction. The normalized maximum likelihood estimate (MLE) is correct a large part of the time, but the objective criterion technique is recommended as a reserve, and more quickly computed, procedure. Both methods for alias removal depend on the differences in the present model function at upwind and downwind. For 10 km by 10 km cells, it is found that the MLE method introduces a correlation between wind speed errors and aspect angle (wind direction) errors that can be as high as 0.8 or 0.9 and that the wind direction errors are unacceptably large, compared to those obtained for the SASS for similar assumptions.

Pierson, W. J.↗

NASA NTRS2007

Extinction-to-Backscatter Ratios of Lofted Aerosol Layers Observed During the First Three Months of CALIPSO Measurements

Case studies from the first three months of the Cloud and Aerosol Lidar and Infrared Pathfinder Spaceborne Observations (CALIPSO) measurements of lofted aerosol layers are analyzed using transmittance [Young, 1995] and two-wavelength algorithms [Vaughan et al., 2004] to determine the aerosol extinction-to-backscatter ratios at 532 and 1064 nm. The transmittance method requires clear air below the layer so that the transmittance through the layer can be determined. Suitable scenes are selected from the browse images and clear air below features is identified by low 532 nm backscatter signal and confirmed by low depolarization and color ratios. The transmittance and two-wavelength techniques are applied to a number of lofted layers and the extinction-to-backscatter ratios are compared with values obtained from the CALIPSO aerosol models [Omar et al., 2004]. The results obtained from these studies are used to adjust the aerosol models and develop observations based extinction-to-backscatter ratio look-up tables and phase functions. Values obtained by these techniques are compared to Sa determinations using other independent methods with a goal of developing probability distribution functions of aerosol type-specific extinction to backscatter ratios. In particular, the results are compared to values determined directly by the High Spectral Resolution Lidar (HSRL) during the CALIPSO CloudSat Validation Experiments (CCVEX) and Sa determined by the application of the two-wavelength lidar Constrained Ratio Aerosol Model-fit (CRAM) retrieval approach [Cattrall et al., 2005; Reagan et al., 2004] to the HSRL data. The results are also compared to values derived using the empirical relationship between the multiple-scattering fraction and the linear depolarization ratio by using Monte Carlo simulations of water clouds [Hu et al., 2006].

Omar, Ali H.↗

NASA NTRS2010

Absorbed Dose and Dose Equivalent Calculations for Modeling Effective Dose

While in orbit, Astronauts are exposed to a much higher dose of ionizing radiation than when on the ground. It is important to model how shielding designs on spacecraft reduce radiation effective dose pre-flight, and determine whether or not a danger to humans is presented. However, in order to calculate effective dose, dose equivalent calculations are needed. Dose equivalent takes into account an absorbed dose of radiation and the biological effectiveness of ionizing radiation. This is important in preventing long-term, stochastic radiation effects in humans spending time in space. Monte carlo simulations run with the particle transport code FLUKA, give absorbed and equivalent dose data for relevant shielding. The shielding geometry used in the dose calculations is a layered slab design, consisting of aluminum, polyethylene, and water. Water is used to simulate the soft tissues that compose the human body. The results obtained will provide information on how the shielding performs with many thicknesses of each material in the slab. This allows them to be directly applicable to modern spacecraft shielding geometries.

Welton, Andrew↗

DOE OSTI2021

Proof that Combining the Forced-collision and DXTRAN Monte Carlo Variance-reduction Techniques is Fair

This report provides mathematical proof that the MCNP DXTRAN (also known in other codes as forced-flight) and forced-collision variance-reduction techniques do not bias the expected value of Monte Carlo simulation estimates when combined. Proof is also provided that the techniques are unbiased when used independently. To prove that the techniques are unbiased, this report derives the first History Score Moment Equation (HSME) for non-multiplying media employing only forced collisions, only DXTRAN, and then both the variance-reduction techniques combined, as defined next. The HSMEs are found by forming the History Score Probability Density Functions (HSPDFs) and then taking the first score moment. Following its derivation, each HSME with a variance-reduction technique employed is reduced to the HSME for an analog simulation. Because the HSME represents the expected contribution to estimators in the simulation, reducing the HSME with variance reduction to the analog HSME shows that the simulation is unbiased despite the variance-reduction technique considered in the HSME. Analysis regarding higher-score moments of each technique is the subject of prior work and is not addressed herein. Throughout this work, the phase space p is defined to be the particle position x, direction-of-flight unit vector $\hat{Ω}$, energy E, and statistical weight w, $$p ≡ (x; \hat{Ω}; E; w).$$ A reduced phase-space excluding the statistical weight of the particle, $$r ≡ (x; \hat{Ω}; E),$$ is also used.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

The cosmic Mach number - Direct comparisons of observations and models

A robust method is used to derive the cosmic Mach number, M, from existing peculiar velocity data sets. For the IR Tully-Fisher distances of spirals in the Local Supercluster of Aaronson et al. (1979, 1982), M of 1.03 and a characteristic distance from the Local Group of 1464 km/s are found. For the elliptical galaxy sample of Faber et al. (1987, 1988), M of 0.57, with a characteristic distance of 2572 km/s is found. These results are compared with Monte Carlo simulations of the observational realizations drawn from numerical simulations of the universe based on various scenarios, including the standard cold dark matter (CDM) scenario. The effect of velocity bias on the derived Mach number is found to be small. Only 5 percent of the CDM simulations have M as large as or larger than that derived for the Aaronson et al. sample.

Strauss, Michael A.↗

DOE OSTI2025

High-throughput multimodal exploration of a nanocrystalline Cu-Ag library

Sputter-deposited, nanocrystalline Cu-Ag thin films produced across a broad compositional and deposition-parameter space were evaluated to unravel the process-structure-property relationships important for creating hard, conductive electrical contacts and coatings. Combinatorial deposition involving pulsed direct current magnetron sputtering of elemental targets enabled swift examination of nearly the full range of alloy compositions and a relevant portion of deposition atomistics. Several high-throughput characterization modalities were employed to evaluate the chemistry, structure, and properties of the films. The resultant hardness, modulus, film density, crystal texture, and resistivity were analyzed in terms of key deposition characteristics (incident atom kinetic energy and incidence angle) predicted by binary-collision, kinematic Monte Carlo simulations. The study revealed improved hardness, parabolic resistivity dependence on composition, and compositional and process dependencies of film tarnishing. The results are discussed in the context of variations in microstructure and film density. Transmission electron microscopy and X-ray diffraction demonstrate several forms of compositional variation including solute segregation to grain boundaries as well as periodic, intragranular compositional modulations. Annealing of a Cu-rich alloy film exhibiting grain boundary segregation showed that this as-deposited, compositional variation is not stable above 100 °C. Finally, the Cu-Ag system is shown to have potential for hard, conductive, tarnish-resistant and room temperature-stable nanocrystalline thin films across the composition space.

36 MATERIALS SCIENCE↗

NASA NTRS1971

Study of automatic flare and decrab guidance and control system for the space shuttle

The automatic landing of an unpowered space shuttle vehicle was studied along with the design of the guidance and flight controls, selection of the ground reference component of the landing system (signals in space), and evaluation of the performance of the overall system. The study was carried out by means of an all-digital simulation which incorporated wind shears, gusts and error models for the hardware components. The performance of the landing system operating in winds and turbulence and with hardware errors was determined by Monte Carlo techniques. It was found that the interactions of turbulence and vehicle characteristics were the limiting factors on system performance. The touchdown errors due to the measurement errors of the ground reference system were less significant. The guidance designed was predictive in the longitudinal direction and fixed path in the lateral direction.

Cockayne, W.↗

NASA NTRS2014

Solar Wind Implantation into Lunar Regolith: Hydrogen Retention in a Surface with Defects

Solar wind protons are implanted directly into the top 100 nm of the lunar near-surface region, but can either quickly diffuse out of the surface or be retained, depending upon surface temperature and the activation energy, U, associated with the implantation site. In this work, we explore the distribution of activation energies upon implantation and the associated hydrogen-retention times; this for comparison with recent observation of OH on the lunar surface. We apply a Monte Carlo approach: for simulated solar wind protons at a given local time, we assume a distribution of U values with a central peak, U(sub c) and width, U(sub w), and derive the fraction retained for long periods in the near-surface. We find that surfaces characterized by a distribution with predominantly large values of U (greater than 1 eV) like that expected at defect sites will retain implanted H (to likely form OH). Surfaces with the distribution predominantly at small values of U (less than 0.2 eV) will quickly diffuse away implanted H. However, surfaces with a large portion of activation energies between 0.3 eV less than U less than 0.9 eV will tend to be H-retentive in cool conditions but transform into H-emissive surfaces when warmed (as when the surface rotates into local noon). These mid-range activation energies give rise to a diurnal effect with diffusive loss of H at noontime.

regolith↗

NASA NTRS1985

Shock acceleration of diffuse ions at the earth's bow shock Acceleration efficiency and A/Z enhancement

Observations of particle spectra, intensity, and enhancement of alpha particles over protons at diffuse ion events at the quasi-parallel earth bow shock are compared to a Monte Carlo simulation of diffusive shock acceleration. The simulation includes the back reaction of accelerated particles on the shock structure, particle escape at an upstream free escape boundary, and a low energy per nucleon threshold for thermal leakage of downstream, shock-heated particles into the upstream region. The simulation assumes that the same scattering operator that gives rise to shock acceleration can also describe a viscous shock governed by hydrodynamic turbulence. This implies that accelerated ions can be drawn directly from the thermal solar wind with no separate superthermal seed population. Good agreement between the simulation and observations made during nearly radial magnetic field configurations lends support to thermal leakage of downstream, shock-heated ions as the mode of injection for diffusion ion events.

Ellison, D. C.↗

DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2011

Calcium in Mercury's Exosphere: Modeling MESSENGER Data

Mercury is surrounded by a surface-bounded exosphere comprised of atomic species including hydrogen, sodium, potassium, calcium, magnesium, and likely oxygen. Because it is collisionless. the exosphere's composition represents a balance of the active source and loss processes. The Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface. Space ENvironment. GEochemistry. and Ranging (MESSENGER) spacecraft has made high spatial-resolution observations of sodium, calcium, and magnesium near Mercury's surface and in the extended, anti-sunward direction. The most striking feature of these data has been the substantial differences in the spatial distribution of each species, Our modeling demonstrates that these differences cannot be due to post-ejection dynamics such as differences in photo-ionization rate and radiation pressure. but instead point to differences in the source mechanisms and regions on the surface from which each is ejected. The observations of calcium have revealed a strong dawn/dusk asymmetry. with the abundance over the dawn hemisphere significantly greater than over the dusk. To understand this asymmetry, we use a Monte Carlo model of Mercury's exosphere that we developed to track the motions of exospheric neutrals under the influence of gravity and radiation pressure. Ca atoms can be ejected directly from the surface or produced in a molecular exosphere (e.g., one consisting of CaO). Particles are removed from the system if they stick to the surface or escape from the model region of interest (within 15 Mercury radii). Photoionization reduces the final weighting given to each particle when simulating the Ca radiance. Preliminary results suggest a high temperature ( I-2x 10(exp 4) K) source of atomic Ca concentrated over the dawn hemisphere. The high temperature is consistent with the dissociation of CaO in a near-surface exosphere with scale height <= 100 km, which imparts 2 eV to the freshly produced Ca atom. This source region and energy are consistent with data from the three MESSENGER flybys; whether this holds true for the data obtained in orbit is under investigation.

Burger, Matthew H.↗

DOE OSTI2021

Dynamic life-cycle carbon analysis for fast pyrolysis biofuel produced from pine residues: implications of carbon temporal effects

Abstract Background Woody biomass has been considered as a promising feedstock for biofuel production via thermochemical conversion technologies such as fast pyrolysis. Extensive Life Cycle Assessment studies have been completed to evaluate the carbon intensity of woody biomass-derived biofuels via fast pyrolysis. However, most studies assumed that woody biomass such as forest residues is a carbon–neutral feedstock like annual crops, despite a distinctive timeframe it takes to grow woody biomass. Besides, few studies have investigated the impacts of forest dynamics and the temporal effects of carbon on the overall carbon intensity of woody-derived biofuels. This study addressed such gaps by developing a life-cycle carbon analysis framework integrating dynamic modeling for forest and biorefinery systems with a time-based discounted Global Warming Potential (GWP) method developed in this work. The framework analyzed dynamic carbon and energy flows of a supply chain for biofuel production from pine residues via fast pyrolysis. Results The mean carbon intensity of biofuel given by Monte Carlo simulation across three pine growth cases ranges from 40.8–41.2 g CO 2 e MJ −1 (static method) to 51.0–65.2 g CO 2 e MJ −1 (using the time-based discounted GWP method) when combusting biochar for energy recovery. If biochar is utilized as soil amendment, the carbon intensity reduces to 19.0–19.7 g CO 2 e MJ −1 (static method) and 29.6–43.4 g CO 2 e MJ −1 in the time-based method. Forest growth and yields (controlled by forest management strategies) show more significant impacts on biofuel carbon intensity when the temporal effect of carbon is taken into consideration. Variation in forest operations and management (e.g., energy consumption of thinning and harvesting), on the other hand, has little impact on the biofuel carbon intensity. Conclusions The carbon temporal effect, particularly the time lag of carbon sequestration during pine growth, has direct impacts on the carbon intensity of biofuels produced from pine residues from a stand-level pine growth and management point of view. The carbon implications are also significantly impacted by the assumptions of biochar end-of-life cases and forest management strategies.

09 BIOMASS FUELS↗

NASA NTRS1988

Power spectral density analysis of wind-shear turbulence for related flight simulations

Meteorological phenomena known as microbursts can produce abrupt changes in wind direction and/or speed over a very short distance in the atmosphere. These changes in flow characteristics have been labelled wind shear. Because of its adverse effects on aerodynamic lift, wind shear poses its most immediate threat to flight operations at low altitudes. The number of recent commercial aircraft accidents attributed to wind shear has necessitated a better understanding of how energy is transferred to an aircraft from wind-shear turbulence. Isotropic turbulence here serves as the basis of comparison for the anisotropic turbulence which exists in the low-altitude wind shear. The related question of how isotropic turbulence scales in a wind shear is addressed from the perspective of power spectral density (psd). The role of the psd in related Monte Carlo simulations is also considered.

Laituri, Tony R.↗

NASA NTRS1999

Diffuse Reflection of Laser Light From Clouds

Laser light reflected from an aqueous suspension of particles or "cloud" with known thickness and particle size distribution defines the "cloud radiative Green's function", G. G is sensitive to cloud thickness, allowing retrieval of that important quantity. We describe a laboratory simulation of G, useful in design of an offbeam Lidar instrument for remote sensing of cloud thickness. Clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The size distribution extends from 0.5 microns to 25 microns, roughly lognormal, similar to real clouds. Density of suspended spheres is adjusted so photon mean-free-path is about 10 cm, 1000 times smaller than in real clouds. The light source is a Nd:YAG laser at 530 nm. Detectors are flux and photon-counting PMTs, with a glass probe for precise positioning. A Labview 5 VI controls position and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider,and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns. Step size is selectable. Far from the beam, the rate of exponential increase in the beam direction scales as expected from diffusion theory, linearly with cloud thickness, and inversely as the square root of the reduced optical thickness, independent of particle size. Nearer the beam the signal increases faster than exponential and depends on particle size. Results verify 3D Monte Carlo simulations that demonstrate detectability of remotely sensed offbeam returns, without filters at night, with narrow bandpass filter in day.

Cahalan, R. F.↗