SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Correlations of baryon and charge stopping in heavy ion collisions*

Baryon numbers are theorized to be carried by valence quarks in the standard QCD picture of the baryon structure. Another theory proposed an alternative baryon number carrier, a non-perturbative Y-shaped configuration of the gluon field, called the baryon junction in the 1970s. However, neither of these theories has been verified experimentally. Recently, searching for the baryon junction by investigating the correlation of net-charge and net-baryon yields at midrapidity in heavy-ion collisions has been suggested. Here, this paper presents studies of such correlations in collisions of various heavy ions from oxygen to uranium with the UrQMD Monte Carlo model. The UrQMD model implements valence quark transport as the primary means of charge and baryon stopping at midrapidity. Detailed studies are also conducted for isobaric $_{40}^{96}{\rm{Zr}}$ + $_{40}^{96}{\rm{Zr}}$ and $_{44}^{96}{\rm{Ru}}$ + $_{44}^{96}{\rm{Ru}}$ collisions. We found a universal trend of charge stopping with respect to baryon stopping and discovered that the charge stopping is always greater than the baryon stopping. This study provides a model baseline in valence quark transport for what is expected in net-charge and net-baryon yields at the midrapidity of relativistic heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Beam modulation and bump-on-tail effects on Alfvén eigenmode stability in DIII-D

Beam modulation effects on Alfvén eigenmode stability have been investigated in a recent DIII-D experiment and show that variations in neutral beam modulation period can have an impact on the beam driven Alfvén eigenmode spectrum and resultant fast ion transport despite similar time-averaged input power. The experiment was carried out during the current ramp phase of L-mode discharges heated with sub-Alfvénic 50–80 kV deuterium neutral beams that drive a variety of Alfvén eigenmodes unstable. The modulation period of two interleaved beams with different tangency radii was varied from shot to shot in order to modify the relative time dependent mix of the beam pitch angle distribution as well as the persistence of a bump-on-tail feature near the injection energy (a feature confirmed by imaging neutral particle analyzer measurements). As the beam modulation period is varied from 7 ms to 30 ms on/off (typical full energy slowing down time of τ slow ≈ 50 ms at mid-radius), toroidicity-induced Alfvén eigenmodes (TAEs) located in the outer periphery of the plasma become intermittent and coincident with the more tangential beam. Core mode activity changes from reversed shear Alfvén eigenmodes (RSAEs) to a mix of RSAE and beta-induced Alfvén eigenmodes. Discharges with 30 ms on/off period do not have a persistent bump-on-tail feature, have the lowest average mode amplitude and least fast ion transport. Furthermore, detailed analysis of an individual TAE using TRANSP kick modeling (Monte Carlo evolution of the distribution function with probabilistic 'kicks' by the AEs) and the resistive MHD code with kinetic fast ions, MEGA, find no strong role of energy gradient drive due to bump-on-tail features. Instead, the observed TAE modulation with interleaved beams is likely a pitch angle dependent result combined with slowing down of the tangential beam between pulses. For the conditions investigated, bump-on-tail contributions to TAE drive were found to be 5% or less of the total drive at any given time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Verification of global gyrokinetic simulation of low frequency mode excited by thermal plasma in spherical tokamak

A low frequency Alfven eigenmode in the ST40 spherical tokamak (ST) is found to be resonantly excited by thermal ions from global gyrokinetic GTC simulations. The analyses of phase-space resonances and wave-particle energy exchanges show that the Beta-induced Alfvén-Acoustic Eigenmode (BAAE) is primarily driven unstable by thermal ions, but is stabilized by energetic particles (EPs). When all kinetic effects are artificially suppressed, GTC simulations in the ideal magnetohydrodynamic (MHD) limit confirm the existence of the MHD eigenmode with a frequency in the BAAEs gap, as previously found by the ideal MHD code NOVA. The kinetic effects of thermal ions induce a frequency downshift compared to the stable BAAE in the MHD simulations. These results support the verification for global gyrokinetic simulations of the BAAE in the ST. However, the simulations do not explain the experimental observations of the BAAE in the ST40, possibly owing to the limitation of the collisionless gyrokinetic Monte Carlo model for the EPs.

Alfvén eigenmode↗

DOE OSTI2021

Measurement of the angular distribution of wavelength-shifted light emitted by TPB

We present measurements of the angular distribution of re-emitted light from tetraphenyl butadiene thin films when exposed to 128 nm light, which is the peak of the liquid Argon (LAr) scintillation spectrum, in vacuum. Films ranging from 250 nm to 5.5 μm in thickness are measured. All films were fabricated by evaporation deposition on ultraviolet transmitting (UVT) acrylic substrates. Preliminary comparisons of the angular distribution to that produced by a detailed Monte Carlo model are also presented. The current shortcomings of the model are discussed and future plans briefly outlined.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Above-threshold ionization at laser intensity greater than 10 20 W / cm 2

Here we present experimental observation of above-threshold ionization (ATI) electrons produced by ionization of the neon K shell in a laser field where intensity exceeds 10 20 W/cm 2 . An array of plastic scintillating calorimeter detectors was used to measure the high-energy electrons at four angles in the laser forward direction. Coarse energy resolution was obtained using aluminum filters of several thicknesses to block lower-energy electrons. A threshold intensity around 2×10 20 W/cm 2 is observed for production of energetic ATI electrons in the laser forward direction, with maximum electron energy exceeding 10 MeV. L-shell electrons with energies <1.4 MeV are scattered farther forward along the laser direction than expected. We present comparisons of the measured total electron energies to the predictions of Monte Carlo models employing the ADK-PPT ionization model and the Augst barrier suppression ionization model.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Flow-plane decorrelations in heavy-ion collisions with multiple-plane cumulants

The azimuthal correlations between local flow planes at different (pseudo)rapidities (η) may reveal important details of the initial nuclear matter density distributions in heavy-ion collisions. Extensive experimental measurements of a factorization ratio (r 2 ) and its derivative (F 2 ) have shown evidence of the longitudinal flow-plane decorrelation. However, nonflow effects also affect this observable and prevent a quantitative understanding of the phenomenon. In this paper, to distinguish decorrelation and nonflow effects, we propose a new cumulant observable, T 2 , which largely suppresses nonflow. Here, the technique sensitivity to different initial-state scenarios and nonflow effects are tested with a simple Monte Carlo model, and in the end, the method is applied to events simulated by a multiphase transport model (AMPT) for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. We also emphasize that a distinct decorrelation signal requires not only the right sign of an observable but also its proper dependence on the η window of the reference flow plane to be consistent with the pertinent decorrelation picture.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Exclusive J / ψ , ψ ( 2 s ) , and e + e - pair production in Au + Au ultraperipheral collisions at the BNL Relativistic Heavy Ion Collider

Measurements of exclusive J/ψ, ψ(2⁢s), and electron-positron (e + ⁢e - ) pair photoproduction in Au+Au ultraperipheral collisions are reported by the STAR experiment at $\sqrt{s_{NN}}$ = 200GeV. Here, we report several first measurements at the BNL Relativistic Heavy Ion Collider, which are (i) J/ψ photoproduction with large momentum transfer up to 2.2⁢(GeV/c) 2 , (ii) coherent J/ψ photoproduction associated with neutron emissions from nuclear breakup, (iii) the rapidity dependence of incoherent J/ψ photoproduction, (iv) the ψ⁡(2⁢s) photoproduction cross section at midrapidity, and (v) e + ⁢e - pair photoproduction up to high invariant mass of 6GeV/c 2 . For measurement (ii), the coherent J/ψ total cross section of γ + Au → J/ψ + Au as a function of the center-of-mass energy W γN has been obtained without photon energy ambiguities. The data are quantitatively compared with the Monte Carlo models STARlight, Sartre, BeAGLE, and theoretical calculations of gluon saturation with color glass condensate, nuclear shadowing with leading twist approximation, quantum electrodynamics, and the next-to-leading-order perturbative QCD. At the photon-nucleon center-of-mass energy of 25.0 GeV, the coherent and incoherent J/ψ cross sections of Au nuclei are found to be 71% ± 10% and 36% ± 7%, respectively, of that of free protons. These data provide an important experimental constraint for nuclear parton distribution functions and a unique opportunity to advance the understanding of the nuclear modification effect at the top RHIC energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

RMCProfile7 : reverse Monte Carlo for multiphase systems

This work introduces a completely rewritten version of the programRMCProfile(version 7), big-box, reverse Monte Carlo modelling software for analysis of total scattering data. The major new feature ofRMCProfile7is the ability to refine multiple phases simultaneously, which is relevant for many current research areas such as energy materials, catalysis and engineering. Other new features include improved support for molecular potentials and rigid-body refinements, as well as multiple different data sets. An empirical resolution correction and calculation of the pair distribution function as a back-Fourier transform are now also available.RMCProfile7is freely available for download at https://rmcprofile.ornl.gov/.

Chemistry↗

DOE OSTI2022

Land‐based climate solutions for the United States

Abstract Meeting end‐of‐century global warming targets requires aggressive action on multiple fronts. Recent reports note the futility of addressing mitigation goals without fully engaging the agricultural sector, yet no available assessments combine both nature‐based solutions (reforestation, grassland and wetland protection, and agricultural practice change) and cellulosic bioenergy for a single geographic region. Collectively, these solutions might offer a suite of climate, biodiversity, and other benefits greater than either alone. Nature‐based solutions are largely constrained by the duration of carbon accrual in soils and forest biomass; each of these carbon pools will eventually saturate. Bioenergy solutions can last indefinitely but carry significant environmental risk if carelessly deployed. We detail a simplified scenario for the United States that illustrates the benefits of combining approaches. We assign a portion of non‐forested former cropland to bioenergy sufficient to meet projected mid‐century transportation needs, with the remainder assigned to nature‐based solutions such as reforestation. Bottom‐up mitigation potentials for the aggregate contributions of crop, grazing, forest, and bioenergy lands are assessed by including in a Monte Carlo model conservative ranges for cost‐effective local mitigation capacities, together with ranges for (a) areal extents that avoid double counting and include realistic adoption rates and (b) the projected duration of different carbon sinks. The projected duration illustrates the net effect of eventually saturating soil carbon pools in the case of most strategies, and additionally saturating biomass carbon pools in the case of forest management. Results show a conservative end‐of‐century mitigation capacity of 110 (57–178) Gt CO 2 e for the U.S., ~50% higher than existing estimates that prioritize nature‐based or bioenergy solutions separately. Further research is needed to shrink uncertainties, but there is sufficient confidence in the general magnitude and direction of a combined approach to plan for deployment now.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Production of light-flavor hadrons in pp collisions at $\sqrt{s}~=~7\text { and }\sqrt{s} = 13 \, \text { TeV} $

The production of $\pi ^{\pm }$ , $\mathrm{K}^{\pm }$ , $\mathrm{K}^{0}_{S}$ , $\mathrm{K}^{*}(892)^{0}$ , $\mathrm{p}$ , $\phi (1020)$ , $\Lambda $ , $\Xi ^{-}$ , $\Omega ^{-}$ , and their antiparticles was measured in inelastic proton–proton (pp) collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV at midrapidity ( $|y|<0.5$ ) as a function of transverse momentum ( $p_{\mathrm{T}}$ ) using the ALICE detector at the CERN LHC. Furthermore, the single-particle $p_{\mathrm{T}}$ distributions of $\mathrm{K}^{0}_{S}$ , $\Lambda $ , and $\overline{\Lambda }$ in inelastic pp collisions at $\sqrt{s} = 7$ TeV are reported here for the first time. The $p_{\mathrm{T}}$ distributions are studied at midrapidity within the transverse momentum range $0\le p_{\mathrm{T}}\le 20$ GeV/ c , depending on the particle species. The $p_{\mathrm{T}}$ spectra, integrated yields, and particle yield ratios are discussed as a function of collision energy and compared with measurements at lower $\sqrt{s}$ and with results from various general-purpose QCD-inspired Monte Carlo models. A hardening of the spectra at high $p_{\mathrm{T}}$ with increasing collision energy is observed, which is similar for all particle species under study. The transverse mass and $x_{\mathrm{T}}\equiv 2p_{\mathrm{T}}/\sqrt{s}$ scaling properties of hadron production are also studied. As the collision energy increases from $\sqrt{s}$ = 7–13 TeV, the yields of non- and single-strange hadrons normalized to the pion yields remain approximately constant as a function of $\sqrt{s}$ , while ratios for multi-strange hadrons indicate enhancements. The $p_\mathrm{{T}}$ -differential cross sections of $\pi ^{\pm }$ , $\mathrm {K}^{\pm }$ and $\mathrm {p}$ ( $\overline{\mathrm{p}}$ ) are compared with next-to-leading order perturbative QCD calculations, which are found to overestimate the cross sections for $\pi ^{\pm }$ and $\mathrm{p}$ ( $\overline{\mathrm{p}}$ ) at high $p_\mathrm{{T}}$ .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Multiplicity dependence of (multi-)strange hadron production in proton-proton collisions at $\sqrt{s}$s = 13 TeV

The production rates and the transverse momentum distribution of strange hadrons at mid-rapidity ( | y | < 0.5 ) are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV as a function of the charged particle multiplicity, using the ALICE detector at the LHC. The production rates of K$^0_S$, Λ, Ξ , and Ω ncrease with the multiplicity faster than what is reported for inclusive charged particles. The increase is found to be more pronounced for hadrons with a larger strangeness content. Possible auto-correlations between the charged particles and the strange hadrons are evaluated by measuring the event-activity with charged particle multiplicity estimators covering different pseudorapidity regions. When comparing to lower energy results, the yields of strange hadrons are found to depend only on the mid-rapidity charged particle multiplicity. Several features of the data are reproduced qualitatively by general purpose QCD Monte Carlo models that take into account the effect of densely-packed QCD strings in high multiplicity collisions. However, none of the tested models reproduce the data quantitatively. This work corroborates and extends the ALICE findings on strangeness production in proton-proton collisions at 7 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Measurement of charged particle multiplicity distributions in DIS at HERA and its implication to entanglement entropy of partons

Charged particle multiplicity distributions in positron-proton deep inelastic scattering at a centre-of-mass energy √s = 319 GeV are measured. The data are collected with the H1 detector at HERA corresponding to an integrated luminosity of 136 pb -1 . Charged particle multiplicities are measured as a function of photon virtuality Q 2 , inelasticity y and pseudorapidity η in the laboratory and the hadronic centre-of-mass frames. Predictions from different Monte Carlo models are compared to the data. The first and second moments of the multiplicity distributions are determined and the KNO scaling behaviour is investigated. The multiplicity distributions as a function of Q 2 and the Bjorken variable x bj are converted to the hadron entropy S hadron , and predictions from a quantum entanglement model are tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Measurements of multiplicity fluctuations of identified hadrons in inelastic proton–proton interactions at the CERN Super Proton Synchrotron

Measurements of multiplicity fluctuations of identified hadrons produced in inelastic p+p interactions at 31, 40, 80, and 158 $\text {Ge}\text {V}/c$ beam momentum are presented. Three different measures of multiplicity fluctuations are used: the scaled variance $\omega $ and strongly intensive measures $\Sigma $ and $\Delta $. These fluctuation measures involve second and first moments of joint multiplicity distributions. Data analysis is preformed using the Identity method which corrects for incomplete particle identification. Strongly intensive quantities are calculated in order to allow for a direct comparison to corresponding results on nucleus–nucleus collisions. The results for different hadron types are shown as a function of collision energy. A comparison with predictions of string-resonance Monte-Carlo models: Epos, Smash and Venus, is also presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Jet-like correlations with respect to $\textrm{K}^{0}_\mathrm{{S}}$ and $\Lambda $ ($\overline{\Lambda }$) in pp and central Pb–Pb collisions at $\mathbf {\sqrt{s_\textrm{NN}}} = 5.02$ TeV

Two-particle correlations with $\textrm{K}^{0}_\mathrm{{S}}$, $\Lambda $ /$\overline{\Lambda }$, and charged hadrons as trigger particles in the transverse momentum range 8 < p T,trig < 16 GeV/c , and associated charged particles within 1 < p T,assoc < 8 GeV/c , are studied at midrapidity in pp and central Pb–Pb collisions at a centre-of-mass energy per nucleon–nucleon collision $\mathbf {\sqrt{s_\textrm{NN}}} = 5.02$ TeV with the ALICE detector at the LHC. After subtracting the contributions of the flow background, the per-trigger yields are extracted on both the near and away sides, and the ratio in Pb–Pb collisions with respect to pp collisions (I AA ) is computed. The per-trigger yield in Pb–Pb collisions on the away side is strongly suppressed to the level of I AA ≈0.6 for p T,assoc > 3 GeV/c as expected from strong in-medium energy loss, while an enhancement develops at low p T,assoc on both the near and away sides, reaching I AA ≈1.8 and 2.7 respectively. These findings are in good agreement with previous ALICE measurements from two-particle correlations triggered by neutral pions (π 0 –h) and charged hadrons (h–h) in Pb–Pb collisions at $\mathbf {\sqrt{s_\textrm{NN}}} = 2.76$ TeV. Moreover, the correlations with $\textrm{K}^{0}_\mathrm{{S}}$ mesons and $\Lambda $ /$\overline{\Lambda }$ baryons as trigger particles are compared to those of inclusive charged hadrons. The results are compared with the predictions of Monte Carlo models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Measurement of groomed event shape observables in deep-inelastic electron-proton scattering at HERA

The H1 Collaboration at HERA reports the first measurement of groomed event shape observables in deep inelastic electron-proton scattering (DIS) at $\sqrt{s} =319$ GeV, using data recorded between the years 2003 and 2007 with an integrated luminosity of 351 pb -1 . Event shapes provide incisive probes of perturbative and non-perturbative QCD. Grooming techniques have been used for jet measurements in hadronic collisions; this paper presents the first application of grooming to DIS data. The analysis is carried out in the Breit frame, utilizing the novel Centauro jet clustering algorithm that is designed for DIS event topologies. Events are required to have squared momentum-transfer $Q^2 > 150$ GeV 2 and inelasticity $0.2< y < 0.7$. We report measurements of the production cross section of groomed event 1-jettiness and groomed invariant mass for several choices of grooming parameter. Monte Carlo model calculations and analytic calculations based on Soft Collinear Effective Theory are compared to the measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Benchmarking of ENDF/B-VIII.0 cross sections for LWR ex-core transport applications

ENDF/B-VIII.0 cross section data is benchmarked on a simple Monte-Carlo model representative of the ex-core geometry of a PWR used for reactor vessel neutron fluence calculations. The benchmarking indicates that the ENDF/B-VIII.0 results show large differences as compared to previous evaluations, which have been extensively compared to measured data. This exercise suggests that ENDF/B-VIII.0 nuclear data is not suitable for ex-core LWR radiation transport applications. Evaluators are aware of the problems with the ENDF/B-VIII.0 evaluation of iron and work is underway to improve the evaluations for an upcoming release.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Joint Focal Mechanism Inversion Using Downhole and Surface Monitoring at the Decatur, Illinois, CO2 Injection Site

ABSTRACT The three-year CO2 injection period at the Illinois Basin - Decatur Project site (Decatur, Illinois, United States) produced a number of microseismic events distributed in very distinct spatiotemporal clusters with different orientations. Further characterization of the microseismicity encompasses the determination of the event source mechanisms. Initially, the microseismic monitoring network consisted solely of borehole sensors, but has been extended with surface sensors, thereby significantly improving the data coverage over the focal sphere. This article focuses on 23 events from the northernmost microseismic cluster (about 2 km from the injection point) and takes advantage of both, surface and downhole, recordings. The resulting strike-slip east–west-oriented focal planes are all consistent with the east–west orientation of the cluster in map view. The injection-related increase of pore pressure is far below the formation fracture pressure; however, small stress-field changes associated with the pore-pressure increase may reach as far as to the investigated cluster location. Monte Carlo modeling of the slip reactivation potential within this cluster showed that the observed maximum stress-field orientation of N068° is the optimum orientation for fault reactivation of the east–west-oriented cluster. Our results suggest that the east–west orientation of the investigated cluster is the main reason for its activation, even though the cluster is about 2 km away from the low-pressure injection point.

Geochemistry & Geophysics↗

DOE OSTI2020

Sodium Void Worth and Control Rod Worth Measurements in ZPPR 15

Sodium void worth and control rod worth were measured in Phases A, B, C and D of the ZPPR15 experimental program. Sodium voiding was simulated by replacing stainless steel cans containing sodium by empty stainless steel cans of the same sizes. The control rod worth measurements were measurements of simulated control rods (CRs) and simulated empty control rod positions (CRPs). A CRP consisted of a 2 x 2 array of drawers filled with sodium plates. A CR consisted of a 2 x 2 array of drawers containing alternating columns of boron carbide plates and sodium plates in the front 18 inches and just sodium plates in the back 18 inches. The ZPPR staff published the measured sodium void worths and control rod worths and the statistical uncertainties in the measured values. In reality, the published uncertainty for each of these measurements is only one component of the total uncertainty. There are additional contributions from uncertainties related to measurement technique, plate compositions and reproducibility. However, the dominant uncertainty for the sodium void worth and control rod worth measurements is the uncertainty in βeff which is required to convert from the natural measurement units, cents, to the calculation units, pcm or Δ/k 1 k 2 . A full uncertainty analysis was performed for the sodium void worth and control rod worth measurements. The significant uncertainties were quantified, and a total uncertainty was derived for each sodium void worth or control rod worth measurement. The ZPPR-15 experimental records were used to create detailed as-built Monte Carlo models of the sodium void worth and control rod worth configurations. These models were used to compute the sodium void worths and control rod worths. (C/E – 1) and the related uncertainty in (C/E – 1) were derived for the measured worths.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗