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 109 records · Page 6

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↗

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↗

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↗

High-pressure annealing driven nanocrystal formation in Zr 50 Cu 40 Al 10 metallic glass and strength increase

Pressure-induced structural changes in metallic glasses have been of great interest as they are expected to open new ways to synthesize novel materials with unexpected properties. Here, we investigated the effect of simultaneous high-pressure and high-temperature treatment on the structure and properties of a Zr 50 Cu 40 Al 10 metallic glass by in situ X-ray structure measurement and property analysis of the final material. We found the unusual formation of Cu-rich nanocrystals at high pressure and temperature, accompanied by significant strength and hardness enhancement. Based on reverse Monte Carlo modeling and molecular dynamics simulations, the structure of the metallic glass changed to a densely packed, chemically uniform configuration with high short-range and medium-range ordering at high pressure and temperature. These results show that high-pressure annealing processes provide a new way to improve and control properties without changing their composition.

36 MATERIALS SCIENCE↗

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↗

Laser-produced plasmas as drivers of laboratory collisionless quasi-parallel shocks

The creation of a repeatable collisionless quasi-parallel shock in the laboratory would provide a valuable platform for experimental studies of space and astrophysical shocks. However, conducting such an experiment presents substantial challenges. Scaling the results of hybrid simulations of quasi-parallel shock formation to the laboratory highlights the experimentally demanding combination of dense, fast, and magnetized background and driver plasmas required. One possible driver for such experiments is high-energy laser-produced plasmas (LPPs). Preliminary experiments at the University of California, Los Angeles, have explored LPPs as drivers of quasi-parallel shocks by combining the Phoenix Laser Laboratory [Niemann et al., J. Instrum. 7, P03010 (2012)] with a large plasma device [Gekelman et al., Rev. Sci. Instrum. 87, 025105 (2016)]. Beam instabilities and waves characteristic of the early stages of shock formation are observed, but spatial dispersion of the laser-produced plasma prematurely terminates the process. This result is illustrated by experimental measurements and Monte Carlo calculations of LPP density dispersion. The experimentally validated Monte Carlo model is then applied to evaluate several possible approaches to mitigating LPP dispersion in future experiments.

Heuer, P. V. (ORCID:0000000150506606)↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Beam-energy dependence of identified two-particle angular correlations in $\sqrt{S_{NN}} = 7.7 – 200$ GeV Au + Au collisions

The two-particle angular correlation functions, $R_2$ , of pions, kaons, and protons in Au + Au collisions at $\sqrt{S_{NN}}$ = 7.7 , 11.5 , 14.5 , 19.6 , 27 , 39 , 62.4 , and 200 GeV were measured by the STAR experiment at the BNL Relativistic Heavy Ion Collider. These correlations were measured for both like-sign and unlike-sign charge combinations and versus the centrality. The correlations of pions and kaons show the expected near-side (i.e., at small relative angles) peak resulting from short-range mechanisms. The amplitudes of these short-range correlations decrease with increasing beam energy. However, the proton correlation functions exhibit strong anticorrelations in the near-side region. This behavior is observed for the first time in an A + A collision system. The observed anticorrelation is $p_T$ independent and decreases with increasing beam energy and centrality. The experimental results are also compared to the Monte Carlo models UrQMD, Hijing, and AMPT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

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↗

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↗