Engineering Papers⌕ Search

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.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Intra-Specific Variability in Plant Hydraulic Parameters Inferred From Model Inversion of Sap Flux Data

Understanding plant hydraulic regulation is critical for predicting plant and ecosystem responses to projected increases in drought stress. Plant hydraulic regulation is controlled by observable, diverse plant hydraulic traits that can vary as much across individuals of the same species as they do across different species. Direct measurements of plant hydraulic traits from a range of ecosystems remain limited in comparison to other, more readily measured traits (e.g., specific leaf area). Furthermore, plant hydraulic trait measurements, often made at leaf or branch levels, are not easily scaled to whole-plant values that are typically used to predict plant and ecosystem fluxes. In this study, multiple whole-plant hydraulic parameters are inferred from observations of plant water use (i.e., sap flow), soil properties, and meteorological data. We use a Markov Chain Monte Carlo model inversion approach to obtain the best estimates and uncertainty of plant hydraulic parameters that capture whole-plant effective embolism resistance and stomatal sensitivity to decreasing plant water potential. We then use the inferred values in the model to estimate whole-tree water use and isohydricity. This approach reliably infers whole-plant parameter values with enough specificity to resolve inter- and intra-specific differences, and thus supplements time- and labor-intensive direct measurements of traits.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Solar Neutrino Detection with a Pixelated Liquid-Argon Time Projection Chamber

This thesis presents a study of low-energy solar neutrino detection using large-scale LArTPCs, focussing on novel pixelated readout technologies. Solar neutrinos offer a unique probe of fundamental neutrino properties and solar physics, but their detection in the MeV range is challenged by backgrounds. We investigate two complementary technologies: SoLAr, which integrates LArPix-based pixelated charge collection with Silicon Photomultipliers (SiPMs) in a hybrid anode design for simultaneous charge and light detection; and Q-Pix, a triggerless pixelated readout architecture based on charge integrate-reset circuits with local clocks, where Reset Time Differences encode ionisation waveforms via time-to-charge conversion. Two SoLAr prototypes were developed and operated, demonstrating VUV-sensitive SiPM performance in liquid argon and accurate charge-light signal matching with a charge detection threshold of $\sim 100 \mathrm{keV}$. We also implement a complete simulation and reconstruction framework, incorporating realistic detector geometry, electron transport, readout response, and detailed signal and background models, including intrinsic argon and radon progeny, as well as site-specific $\gamma$-ray and neutron fluxes. For Q-Pix, we demonstrate that with a pixel size of $4\times 4$ mm$^2$ and a reset threshold of 1 fC ($\sim 0.1475$ MeV), full-scale operation produces data volumes below 1 PB per 10 ktonne-year. For SoLAr, assuming a shielded DUNE-like detector and 100 ktonne-year exposure, we project uncertainties of $0.90\times 10^{-5}$ eV$^2$ on $\Delta m^2_{21}$ and 0.033 on $\sin^2\theta_{12}$, improving to $0.46\times 10^{-5}$ eV$^2$ and 0.025 with 400 kilotonne-year. At this higher exposure, we also obtain a day–night flux asymmetry at the level of $( -5.6 \pm 3.6 ) \%$. Combining Monte Carlo modelling, hardware validation, and advanced reconstruction techniques, this work establishes a path toward next-generation ktonne-scale LArTPCs as observatories for precision solar neutrino physics.

Ruiz Ferreira, Guilherme [Manchester U.] (ORCID:00↗
DOE OSTI2022

Study of dijet events with large rapidity separation in proton-proton collisions at $\sqrt{s}$ = 2.76 TeV

The cross sections for inclusive and Mueller-Navelet dijet production are measured as a function of the rapidity separation between the jets in proton-proton collisions at $\sqrt{s}$ = 2.76 TeV for jets with transverse momentum p T > 35 GeV and rapidity |y| < 4.7. Various dijet production cross section ratios are also measured. A veto on additional jets with p T > 20 GeV is introduced to improve the sensitivity to the effects of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution. The measurement is compared with the predictions of various Monte Carlo models based on leading-order and next-to-leading-order calculations including the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi leading-logarithm (LL) parton shower as well as the LL BFKL resummation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Measurements of Lund subjet multiplicities in 13 TeV proton-proton collisions with the ATLAS detector

This Letter presents a differential cross-section measurement of Lund subjet multiplicities, suitable for testing current and future parton shower Monte Carlo algorithms. This measurement is made in dijet events in 140 fb -1 of $\sqrt{s}$ =13 TeV proton–proton collision data collected with the ATLAS detector at CERN's Large Hadron Collider. The data are unfolded to account for acceptance and detector-related effects, and are then compared with several Monte Carlo models and to recent resummed analytical calculations. The experimental precision achieved in the measurement allows tests of higher-order effects in QCD predictions. Most predictions fail to accurately describe the measured data, particularly at large values of jet transverse momentum accessible at the Large Hadron Collider, indicating the measurement's utility as an input to future parton shower developments and other studies probing fundamental properties of QCD and the production of hadronic final states up to the TeV-scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Unraveling the Early-Stage Ordering of Krypton Solid Bubbles in Molybdenum: A Multimodal Study

Self-organization of defects such as fission gas bubbles in materials can lead to high inventory capacity for fission gas storage and help mitigate swelling caused by fission gases in nuclear fuel materials under radiation in nuclear reactors. Here, we report the physical mechanism of self-organization of krypton (Kr) gas bubbles in molybdenum (Mo) under ion implantation. The ion fluence and temperature-dependent formation of Kr solid bubble superlattice (SBS) in Mo were investigated by using both synchrotron-based small-angle X-ray scattering (SAXS) and transmission electron microscopy (TEM). Early stage self-organization of gas bubbles is observed at a fluence of 2.5 × 10 16 ions/cm 2 at temperatures of 300–400 °C. The bubble lattice constant increases with increasing implantation temperature from 300 to 400 °C. Both experiments and atomic kinetic Monte Carlo modeling indicate that the Kr solid bubbles are weakly ordered in comparison to previously studied helium (He) gas bubble superlattice (GBS) while the lattice constant are relatively smaller for Kr SBS compared to that of He GBS. The irradiation conditions suggest that spinodal decomposition, which is a form of phase separation, probably precedes gas bubble ordering in Mo. Altogether, our work sheds light on the formation mechanism of noble gas superlattice toward the development of radiation-tolerant materials which are important for the design of advanced nuclear reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles Dissociation Pathways of BCl 3 on the Si(100)-2 × 1 Surface

BCl 3 is a promising acceptor precursor for atomic-precision δ-doping of silicon, as it has been observed to rapidly dissociate into boron doped into the silicon surface and surface chlorine, which can be removed upon annealing. The chemical pathway and the resulting kinetics, through which BCl 3 adsorbs and dissociates on silicon, however, have only been partially explained. Here, in this work, we use density functional theory to expand the dissociation reactions of BCl 3 to include reactions that take place across multiple silicon dimer rows and reactions which end in a bare B atom either at the surface, substituted for a surface silicon, or in a subsurface position. We further simulate the resulting scanning tunneling microscopy images for each of these BCl x dissociation fragments, demonstrating that they often display distinct features that may allow for relatively confident experimental identification. Finally, we input the full dissociation pathway for BCl 3 into a kinetic Monte Carlo model, which simulates realistic reaction pathways as a function of environmental conditions, such as the pressure and temperature of dosing. We find that BCl 2 is broadly dominant at low temperatures, while high temperatures and ample space on the silicon surface for dissociation encourage the formation of bridging BCl fragments and B substitutions on the surface. This work provides the chemical mechanisms for understanding atomic-precision doping of Si with B, enabling a number of relevant quantum applications, such as bipolar nanoelectronics, acceptor-based qubits, and superconducting Si.

Campbell, Quinn T. [Sandia National Laboratories (↗
DOE OSTI2023

Unexpected Thermomechanical Behavior of Off-Stoichiometry Epoxy/Amine Materials

Recent studies on off-stoichiometric thermosets reveal unique viscoelastic behavior derived from increased free volume and physical interactions between chain ends. To understand structural characteristics arising from cure and its effect on properties, we developed a Monte Carlo model based on step-growth polymerization. Our model accurately predicted structure–property trends for a two-component system of EPON 828 (EPON) and ethylenediamine. A second epoxy monomer, D.E.R. 732 (DER), was investigated to modulate T g . Binary mixtures of EPON and DER in off-stoichiometric, amine-rich formulations resulted in nonlinear evolution of thermomechanical properties with respect to initial formulation stoichiometry. Modifying our model with kinetic parameters allowing for differential epoxide/amine reaction kinetics only partially accounted for trends in T g , suggesting that spatiotemporal contributions-not captured by our model-were significant determinants of material properties compared to polymer architecture for three-component systems. Here these findings underpin the importance of spatial awareness in modeling to inform the development of dynamic thermosets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Overcoming the Entropy Penalty of Direct Air Capture for Efficient Gigatonne Removal of Carbon Dioxide

Atmospheric carbon poses an existential threat to civilization via global climate change. Hundreds of gigatonnes of carbon dioxide must be removed from earth’s atmosphere in the next three decades, necessitating a low-cost, energy-efficient process to extract low concentrations of carbon dioxide for conversion to a stable material permanently stored for thousands of years. In this work, the challenge of removing gigatonnes of CO 2 is described via the scale of effort and the thermodynamics of collecting and reducing this diffuse chemical, the accumulation of which imparts a substantial entropy penalty on any atmospheric carbon capture process. The methods of CO 2 reduction combined with upstream direct air capture (DAC) including absorption, membrane separation, and adsorption are compared with biomass torrefaction and permanent burial (BTB). A Monte Carlo model assesses the mass, energy, and economics of the full process of biomass torrefaction from biomass collection and transport to stable carbon burial to determine that 95% of scenarios could remove carbon for less than $200 per CO 2 -tonne-equivalent. Torrefied carbon is further discussed for its long-term stability and availability at the scale required to substantially mitigate the threat of climate change.

biomass↗
DOE OSTI2024

Amorphous Oxyhalide Matters for Achieving Lithium Superionic Conduction

The recently surged halide-based solid electrolytes (SEs) are great candidates for high-performance all-solid-state batteries (ASSBs), due to their decent ionic conductivity, wide electrochemical stability window, and good compatibility with high-voltage oxide cathodes. In contrast to the crystalline phases in halide SEs, amorphous components are rarely understood but play an important role in Li-ion conduction. Here, we reveal that the presence of amorphous component is common in halide-based SEs that are prepared via mechanochemical method. The fast Li-ion migration is found to be associated with the local chemistry of the amorphous proportion. Taking Zr-based halide SEs as an example, the amorphization process can be regulated by incorporating O, resulting in the formation of corner-sharing Zr–O/Cl polyhedrons. This structural configuration has been confirmed through X-ray absorption spectroscopy, pair distribution function analyses, and Reverse Monte Carlo modeling. The unique structure significantly reduces the energy barriers for Li-ion transport. As a result, an enhanced ionic conductivity of (1.35 ± 0.07) × 10 –3 S cm –1 at 25 °C can be achieved for amorphous Li 3 ZrCl 4 O 1.5 . In addition to the improved ionic conductivity, amorphization of Zr-based halide SEs via incorporation of O leads to good mechanical deformability and promising electrochemical performance. In conclusion, these findings provide deep insights into the rational design of desirable halide SEs for high-performance ASSBs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Chemical Framework to Design Linear-like Relaxors toward Capacitive Energy Storage

ABO 3 -type perovskite relaxor ferroelectrics (RFEs) have emerged as the preferred option for dielectric capacitive energy storage. However, the compositional design of RFEs with high energy density and efficiency poses significant challenges owing to the vast compositional space and the absence of general rules. Here, we present an atomic-level chemical framework that captures inherent characteristics in terms of radius and ferroelectric activity of ions. By categorizing A/B-site ions as host framework, rattling, ferroelectrically active, and blocking ions and assembling these four types of ions with specific criteria, linear-like relaxors with weak locally correlated and highly extendable unit-cell polarization vectors can be constructed. As example, we demonstrate two new compositions of Bi 0.5 K 0.5 TiO 3- based and BaTiO 3- based relaxors, showing extremely high recoverable energy densities of 17.3 and 12.1 J cm –3 , respectively, both with a high efficiency of about 90%. Further, the role of different types of ions in forming heterogeneous polar structures is identified through element-specific local structure analysis using neutron total scattering combined with reverse Monte Carlo modeling. Further, our work not only opens up new avenues toward rational compositional design of high energy storage performance lead-free RFEs but also sheds light on atomic-level manipulation of functional properties in compositionally complex ferroelectrics.

25 ENERGY STORAGE↗
DOE OSTI2023

Local inversion-symmetry breaking in a bismuthate high-Tc superconductor

Abstract The doped perovskite BaBiO 3 exhibits a maximum superconducting transition temperature ( T c ) of 34 K and was the first high- T c oxide to be discovered, yet pivotal questions regarding the nature of both the metallic and superconducting states remain unresolved. Although it is generally thought that superconductivity in the bismuthates is of the conventional s -wave type, the pairing mechanism is still debated, with strong electron-phonon coupling and bismuth valence or bond disproportionation possibly playing a role. Here we use diffuse x-ray scattering and Monte Carlo modeling to study the local structure of Ba 1-x K x BiO 3 across its insulator-metal boundary. We find no evidence for either long- or short-range disproportionation, which resolves a major conundrum, as disproportionation and the related polaronic effects are likely not relevant for the metallic and superconducting states. Instead, we uncover nanoscale structural correlations that break inversion symmetry, with far-reaching implications for the electronic physics. This unexpected finding furthermore establishes that the bismuthates belong to the broader classes of materials with hidden spin-orbit coupling and a tendency towards inversion-breaking displacements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Neutron spectroscopy of plutonium using a handheld detection system

The ability to distinguish multiple forms of plutonium from one another, such as oxide and metal, is paramount in areas of nuclear nonproliferation and international safeguards. In its metal form, plutonium can be readily used in a nuclear weapon, while oxide forms are associated with nuclear reactor fuel. Oxide-based plutonium forms emit neutrons with an energy spectrum that is significantly different from the fission neutrons that are emitted from plutonium metal. Organic scintillation detectors output pulses that are proportional to the neutron energy deposited, and therefore present a means of distinguishing these plutonium forms based on their energy spectra. In this work, metal and oxide forms of plutonium were measured using a handheld detection system based on an organic glass scintillator. Monte Carlo modeling of these experiments was performed to provide insight into the origin of the features in the observed light output spectra. Through analysis of multiple regions of these spectra, in a matter of minutes we were able to unambiguously discriminate oxide and metal plutonium forms from one another and from a plutonium-beryllium neutron source, which was considered for comparison because these sources are commonly used in industrial applications. The ability to discriminate weapons-usable material from nuclear reactor fuel has applications in nuclear treaty verification and safeguards.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Structural regulation-induced Li-electron disentanglement for stabilized oxygen redox of Li-excess disordered rock-salt cathode materials

Since the discovery of its electrochemical activity, Li-excess disordered rock-salt (DRX) cathode material has received worldwide attention as it sets up a new way to exploit oxygen redox beyond the conventional layered structure with late-3d transition metals. However, the intricate structure-function relationship in the disordered lattice of the DRX material fogs the researcher's lens on the underlying redox mechanisms. Here, in this study, we employ a synergistic approach combining neutron total scattering with reverse Monte Carlo modeling and density functional theory calculations to unravel the landscape of oxygen redox reactions in DRX. Redox activities are evaluated in diverse oxygen clusters (OLi x TM 6-x ) and the spatial distribution of these clusters in the model DRX structure (Li 1.16 Ti 0.37 Ni 0.37 Nb 0.1 O 2 and Li 1.2 Ti 0.35 Ni 0.35 Nb 0.1 O 1.8 F 0.2 ) is explicitly determined. The results unveil that by regulating the short-range ordering between cations, fluorine atoms can effectively decouple the location of Li extraction and electron depletion. Such disentanglement between the Li reservoir and electron reservoir in the DRX lattice could play a pivotal role in protecting the oxidized oxygen and preserving the lattice framework during cycling. Through a tentatively designed non-fluorinated DRX oxide realizing similar Li-electron decoupling, an obvious enhancement of the cycling capability can be achieved without compromising the capacity release.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Intermolecular correlations of liquid and glassy CS 2 studied by synchrotron radiation x-ray diffraction

How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Combinatorial deposition of Au–Bi alloys via high-rate magnetron sputtering

Gold–bismuth (Au–Bi) alloy films are promising candidate materials for inertial confinement fusion (ICF) hohlraums due to their high laser-to-x-ray conversion efficiency, particularly compared with Au, Ta–Au, and Bi hohlraums. However, the fabrication of uniform and dense Au–Bi alloy films remains a challenge. Here, we use a combination of Monte-Carlo modeling and experiments to demonstrate that the microstructure and properties of Au–Bi alloy films can be greatly improved when direct-current magnetron sputtering with high deposition rates of ⩾5 μm h -1 is used. Resultant films are ~90% of their maximum theoretical densities and have low O content of <1 at.%. Films with Bi content above 40 at. % exhibit high electrical resistivity >100μΩ cm, making them suitable for both magnetized and non-magnetized ICF schemes.

Materials science↗
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 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↗