Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Design and simulation of a muon detector to characterize geological overburden

This study presents the design, construction, and simulation of a mobile muon detector tailored for geological overburden characterization. The detector employs plastic scintillator paddles with silicon photomultipliers (SiPMs) and a QuarkNet data acquisition system, offering a portable solution suitable for remote field deployment. The simulator’s modular aluminum frame allows for adjustable geometry and directional sensitivity, while its battery system supports over a week of autonomous operation. Preliminary experimental tests confirmed that its muon flux measurements were consistent with theoretical expectations. A comprehensive simulation framework using Geant4 and CORSIKA was developed to model detector response and overburden effects. Analytical and Monte Carlo methods were used to assess quadrant resolution and infer muon directionality. This work lays the foundation for future overburden mapping and supports the development of reconstruction algorithms for geological applications.

72 - PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Fixed Source Sensitivity Calculations for Inertial Confinement Fusion Applications

A numerical code library was developed for the radiation transport code MCNP6.3 to calculate generalized response sensitivity coefficients for fixed source neutron transport problems with applications to inertial confinement fusion (ICF) experiments. The new MCNP6.3 dependency is used to generate a novel time convolution response that represents a neutron time-of-flight (nToF) signal. The traditional suite of macroscopic cross-section sensitivities and constrained fixed source probability distribution sensitivities are available for both the standard and the new response tallies in this library. However, novel sensitivity coefficients for the constrained hyperparameters of analytic fixed source probability distributions are emphasized in this work for their connection to ICF neutron transport models. Particularly, advanced Monte Carlo methods are developed for calculating the sensitivity of a nToF signal to perturbations in an ICF plasma’s ion temperature and burn history as well as perturbations in the target liner mass density and the shape parameters of the nToF detector’s impulse response function. Together, these capabilities form an advanced suite of computational tools that can be used to analyze and extract information from any ICF experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Wang–Landau approach to the simulation of water clusters

The Wang–Landau Monte Carlo method [Wang FG, Landau DP. Efficient, multiple-range random walk algorithm to calculate the density of states. Phys Rev Lett. 2001;86:2050-2053; Determining the density of states for classical statistical models: a random walk algorithm to produce a flat histogram. Phys Rev E. 2001;64:056101] has been established as an effective computational approach for studying complex systems. In this study, we review the properties of 6 popular empirical water potentials, including a recently proposed explicit 3-body potential, obtained by applying the Wang–Landau method to the simulation of water clusters. Finally, comparing to numerous classical and first principle calculations, we show that several classical water models can provide a consistent description of the structural and thermodynamical properties of water clusters partially in agreement of ab initio results, although they were designed to reproduce macroscopic properties of water at ambient conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Comparison of Stack Measurement Data from R&D Facilities to Regulatory Criteria. A Case Study from PNNL

Chemical emissions from research and development (R&D) activities are difficult to estimate because of the large number of chemicals used and the potential for continual changes in processes. In this case study, stack measurements taken from R&D facilities at Pacific Northwest National Laboratory (PNNL) were examined, including extreme worst-case emissions estimates and alternate analyses using a Monte Carlo method that takes into account the full distribution of sampling results. The results from these analyses were then compared to emissions estimated from chemical inventories. Results showed that downwind ambient air concentrations calculated from the stack measurement data were below acceptable source impact levels (ASILs) for almost all compounds, even under extreme worst-case analyses. However, for compounds with averaging periods of a year, the unrealistic but simplifying extreme worst-case analysis often resulted in exceedances of lower level regulatory criteria used to determine modeling requirements or to define trivial releases. Compounds with 24-hour averaging periods were nearly all several orders of magnitude below all, including the trivial release, criteria. The alternate analysis supplied a more realistic basis of comparison and an ability to explore effects under different operational modes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Improved estimates of the Milky Way’s stellar mass and star formation rate from hierarchical Bayesian meta-analysis

We present improved estimates of several global properties of the Milky Way, including its current star formation rate (SFR), the stellar mass contained in its disk and bulge+bar components, as well as its total stellar mass. We do so by combining previous measurements from the literature using a hierarchical Bayesian (HB) statistical method that allows us to account for the possibility that any value may be incorrect or have underestimated errors. We show that this method is robust to a wide variety of assumptions about the nature of problems in individual measurements or error estimates. Ultimately, our analysis yields an SFR for the Galaxy of M-dot {sub ⋆}=1.65±0.19 M{sub ⊙} yr{sup −1}, assuming a Kroupa initial mass function (IMF). By combining HB methods with Monte Carlo simulations that incorporate the latest estimates of the Galactocentric radius of the Sun, R{sub 0}, the exponential scale length of the disk, L{sub d}, and the local surface density of stellar mass, Σ{sub ⋆}(R{sub 0}), we show that the mass of the Galactic bulge+bar is M{sub ⋆}{sup B}=0.91±0.07×10{sup 10} M{sub ⊙}, the disk mass is M{sub ⋆}{sup D}=5.17±1.11×10{sup 10} M{sub ⊙}, and their combination yields a total stellar mass of M{sub ⋆}=6.08±1.14×10{sup 10} M{sub ⊙} (assuming a Kroupa IMF and an exponential disk profile). This analysis is based upon a new compilation of literature bulge mass estimates, normalized to common assumptions about the stellar IMF and Galactic disk properties, presented herein. We additionally find a bulge-to-total mass ratio for the Milky Way of B/T=0.150{sub −0.019}{sup +0.028} and a specific SFR of M-dot {sub ⋆}/M{sub ⋆}=2.71±0.59×10{sup −11} yr{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dipolar focusing in laser-assisted positronium formation

Abstract We consider positronium formation in collisions of positrons with excited hydrogen atoms H(n) in an infrared laser field theoretically. This process is assisted by the dipolar focusing effect: a positron moving in a superposition of a laser field and the dipolar field can approach the atomic target even if its trajectory starts with a very large impact parameter, leading to a significant enhancement of the Ps formation cross section. The classical trajectory Monte Carlo method, which is justified for n ⩾ 3 , allows efficient calculation of this enhancement. A similar effect can occur in collisions of positrons with other atoms in excited states, which can lead to improvements in the efficiency of positronium formation.

Optics↗
DOE OSTI2019

Fluctuations and phase transitions of uniaxial and biaxial liquid crystals using a theoretically informed Monte Carlo and a Landau free energy density

In this work, we explore fluctuations during phase transitions of uniaxial and biaxial liquid crystals using a phenomenological free energy functional. We rely on a continuum-level description of the liquid crystal ordering with a tensorial parameter and a temperature dependent Landau polynomial expansion of the tensor's invariants. The free energy functional, over a three-dimensional periodic domain, is integrated with a Gaussian quadrature and minimized with a theoretically informed Monte Carlo method. We reconstruct analytical phase diagrams, following Landau and Doi 's notations, to verify that the free energy relaxation reaches the global minimum. Importantly, our relaxation method is able to follow the thermodynamic behavior provided by other non-phenomenological approaches; we predict the first order character of the isotropic-nematic transition, and we identify the uniaxial-biaxial transition as second order. Finally, we use a finite-size scaling method, using the nematic susceptibility, to calculate the transition temperatures for 4-Cyano-4'-pentylbiphenyl (5CB) and N-(4-methoxybenzylidene)-4-butylaniline (MBBA). Our results show good agreement with experimental values, thereby validating our minimization method. Our approach is an alternative towards the relaxation of temperature dependent continuum-level free energy functionals, in any geometry, and can incorporate complicated elastic and surface energy densities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Monte-Carlo fluid approaches to detached plasmas in non-axisymmetric divertor configurations

Fluid transport modeling in three-dimensional boundaries of toroidal confinement devices is reviewed with the emphasis on a Monte-Carlo approach to simulate detached plasmas. The loss of axisymmetry in such configurations presents a major challenge for numerical implementation of the standard fluid model widely applied to fusion experimental devices. A large-scale effort has been made to address this problem under complementary aspects including different magnetic topologies and numerical techniques. Here in this paper, we give a brief review of the different strategies pioneered and the challenges involved. A more detailed description is provided for the Monte-Carlo code—EMC3-Eirene, where the physics model and the basic idea behind the applied Monte-Carlo method are presented. The focus is put on its applications to detachment studies for stellarators and tokamaks. Here, major achievements and difficulties encountered are described. Model limitations and further development plans are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Trigonometric Parallaxes of Two T Dwarfs With Keck and ShaneAO Astrometry

We present trigonometric parallax and proper motion measurements for two T-type brown dwarfs. We derive our measurements from infrared laser guide star adaptive optics observations spanning five years from the ShaneAO/SHARCS and NIRC2/medium-cam instruments on the Shane and Keck telescopes, respectively. To improve our astrometric precision, we measure and apply a distortion correction to our fields for both instruments. We also transform the Keck and ShaneAO astrometric reference frames onto the ICRS using five-parameter parallax and proper motion solutions for background reference stars from Gaia DR2. Fitting for parallax and proper motion, we measure parallaxes of 73.5 ± 9.2 mas and 70.1 ± 6.7 mas for WISEJ19010703+47181688 (WISE1901) and WISEJ21543294+59421370 (WISE2154), respectively. We utilize Monte Carlo methods to estimate the error in our sparse field methods, taking into account overfitting and differential atmospheric refraction. Comparing to previous measurements in the literature, all of our parallax and proper motion values fall within 2σ of the published measurements, and 4 of 6 measurements are within 1σ. These data are among the first parallax measurements of these T dwarfs and serve as precise measurements for calibrating stellar formation models. These two objects are the first results of an ongoing survey of T dwarfs with Keck/NIRC2 and the Shane Adaptive Optics system at Lick Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Re-entrance in nuclei: competitive phenomena

Using the shell-model Monte Carlo method, we investigate how temperature and rotation affect pairing properties for nuclei in the fp - gds region. The re-entrance of pairing correlations with temperature is predicted at high rotational frequencies. It manifests through an anomalous behavior of the specific heat and level density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

State-selective charge transfer cross sections for light ion impact of atomic hydrogen

Owing to the utility of diagnosing plasma properties such as impurity concentration and spatial distribution, and plasma temperature and rotation, by detection of photon emission following capture of electrons from atomic hydrogen to excited states of multiply charged ions, new calculations of state-selective charge transfer involving light ions have been carried out using the atomic orbital close-coupling and the classical trajectory Monte Carlo methods. By comparing these with results of other approaches applicable in a lower impact energy regime, and by benchmarking them using key experimental data, knowledge of the cross sections can be made available across the range parameters needed by fusion plasma diagnostics.

Schultz, D. R.↗
DOE OSTI2022

Inference of molecular divertor density from filtered camera analysis of molecularly induced Balmer line emission during detachment in JET L-mode plasmas

Here a previously presented Monte Carlo method for estimating local plasma conditions in 2D based on intensity ratios of deuterium Balmer D α , D γ and D ε lines was amended to consider also the D α and D γ emission contributions arising from molecular processes. The obtained estimates were used to infer the molecular divertor density with the help of the molecular databases of EIRENE. The method was benchmarked against EDGE2D-EIRENE simulations and observed to reproduce the molecularly induced emission fractions and the molecular divertor densities primarily within 25% of the references. Experimental analysis of a JET L-mode density scan suggested molecularly induced D α and D γ contributions of up to 60–70% and 20%, respectively, during the process of detachment. The independent estimates of the molecular divertor density inferred from the obtained molecularly induced D α and D γ intensities agree within uncertainties with each other. Both estimates show the molecular density increasing up to approximately 1.0–2.0 x 10 20 m -3 at the outer strike point in deep detachment with its ratio to the local electron density agreeing with EDGE2D-EIRENE predictions within the scatter of the experimental data.

interferometry↗
DOE OSTI2022

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2022

SDSS-IV MaStar: [α/Fe] for the MaNGA Stellar Library from synthetic model spectra

ABSTRACT We calculate the α-enhancement ratio [α/Fe] for the Mapping Nearby Galaxies at APO (MaNGA) Stellar Library (MaStar) while also fitting for the fundamental atmospheric parameters effective temperature, surface gravity, and metallicity – Teff, log g, [Fe/H]. This approach builds upon a previous catalogue of stellar parameters, whereby only the fundamental atmospheric parameters are fit with solar-scaled models. Here, we use the same Markov Chain Monte Carlo method with the additional free parameter [α/Fe]. Using the full spectral fitting code pPXF, we are able to fit multiple lines sensitive to [α/Fe] for a more robust measurement. Quality flags based on the convergence of the sampler, errors in [α/Fe] and a cut in the χ2 of the model fit are used to clean the final catalogue, returning 17 214 spectra and values in the range of −0.25 < [α/Fe] < 0.48. Comparing our calculated [α/Fe] with literature values reveals a degeneracy in cool stars with log g ≥ ∼4; this comparison is then used to create an alternative and calibrated parameter set. We also plot the final catalogue in an [Fe/H] versus [α/Fe] diagram and recover the expected result of increasing [α/Fe] with decreasing [Fe/H] for Milky Way disc-halo stars. We apply our method to a subsample of spectra of uniform resolution and higher signal to noise that finds that our results are independent of this higher signal to noise. In the context of stellar population models, we are able to cover a parameter space for the creation of intermediate to old age models at solar-scaled [α/Fe], high [Fe/H] and enhanced [α/Fe], low [Fe/H].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The information content of projected galaxy fields

ABSTRACT The power spectrum of the non-linearly evolved large-scale mass distribution recovers only a minority of the information available on the mass fluctuation amplitude. We investigate the recovery of this information in 2D ‘slabs’ of the mass distribution averaged over ≈100 h−1 Mpc along the line of sight, as might be obtained from photometric redshift surveys. We demonstrate a Hamiltonian Monte Carlo method to reconstruct the non-Gaussian mass distribution in slabs, under the assumption that the projected field is a point-transformed Gaussian random field, Poisson-sampled by galaxies. When applied to the Quijote N-body suite at z = 0.5 and at a transverse resolution of 2 h−1 Mpc, the method recovers ∼30 times more information than the 2D power spectrum in the well-sampled limit, recovering the Gaussian limit on information. At a more realistic galaxy sampling density of 0.01 h3 Mpc−3, shot noise reduces the information gain to a factor of 5 improvement over the power spectrum at resolutions of 4 h−1 Mpc or smaller.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Chemical order-disorder nanodomains in Fe 3 Pt bulk alloy

Chemical ordering is a common phenomenon and highly correlated with the properties of solid materials. By means of the redistribution of atoms and chemical bonds, it invokes an effective lattice adjustment and tailors corresponding physical properties. To date, however, directly probing the 3D interfacial interactions of chemical ordering remains a big challenge. In this work, we deciphered the interlaced distribution of nanosized domains with chemical order/disorder in Fe 3 Pt bulk alloy. HAADF-STEM images evidence the existence of such nanodomains. The reverse Monte Carlo method with the X-ray pair distribution function data reveal the 3D distribution of local structures and the tensile effect in the disordered domains at the single-atomic level. The chemical bonding around the domain boundary changes the bonding feature in the disordered side and reduces the local magnetic moment of Fe atoms. This results in a suppressed negative thermal expansion and extended temperature range in Fe 3 Pt bulk alloy with nanodomains. Our study demonstrates a local revelation for the chemical order/disorder nanodomains in bulk alloy. The understanding gained from atomic short-range interactions within the domain boundaries provides useful insights with regard to designing new functional compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Variations of the physical parameters of the blazar Mrk 421 based on analysis of the spectral energy distributions

Abstract We report on the variations of the physical parameters of the jet observed in the blazar Mrk 421, and discuss the origin of X-ray flares in the jet, based on analysis of several spectral energy distributions (SEDs). The SEDs are modeled using the one-zone synchrotron self-Compton model, its parameters determined using a Markov chain Monte Carlo method. The lack of data at TeV energies means many of the parameters cannot be uniquely determined and are correlated. These are studied in detail. We find that the optimal solution can be uniquely determined only when we apply a constraint to one of four parameters: the magnetic field (B), the Doppler factor, the size of the emitting region, and the normalization factor of the electron energy distribution. We used 31 sets of SEDs from 2009 to 2014 with optical–UV data observed with UVOT/Swift and the Kanata telescope, X-ray data with XRT/Swift, and γ-ray data with the Fermi Large Area Telescope. The result of our SED analysis suggests that, in the X-ray faint state, the emission occurs in a relatively small area (∼1016 cm) with a relatively strong magnetic field (B ∼ 10−1 G). The X-ray bright state shows a tendency opposite to that of the faint state, that is, a large emitting area (∼1018 cm), probably downstream of the jet, and a weak magnetic field (B ∼ 10−3 G). The high X-ray flux was due to an increase in the maximum energy of electrons. On the other hand, the presence of two kinds of emitting areas implies that the one-zone model is unsuitable for reproducing at least part of the observed SEDs.

Yamada, Yurika↗
DOE OSTI2022

Negative string tension of a higher-charge Schwinger model via digital quantum simulation

We study some properties of generalized global symmetry for the charge- q Schwinger model in the Hamiltonian formalism, which is the (1 + 1)D quantum electrodynamics with a charge- q Dirac fermion. This model has the |$\mathbb {Z}_q\, 1$| -form symmetry, which is a remnant of the electric |$U(1)\, 1$| -form symmetry in the pure Maxwell theory. It is known that, if we put the theory on closed space, then the Hilbert space is decomposed into q distinct sectors, called universes, and some states with higher energy density do not decay to the ground state due to the selection rule of the 1-form symmetry. Even with open boundaries, we can observe the stability of such states by seeing a negative string tension behavior, meaning that opposite charges repel each other. In order to see negative string tensions, the vacuum angle θ has to be large enough and the standard path-integral Monte Carlo method suffers from the sign problem. We develop a method based on the adiabatic state preparation to see this feature with digital quantum simulation and confirm it using a classical simulator of quantum devices. In particular, we measure the local energy density and see how it jumps between the inside and outside of the insertion of the probe charges. We explicitly see that the energy density inside is lower than that outside. This is a clear signature of the negative string tension.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗