SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Towards the ionizing radiation induced bond dissociation mechanism in oxygen, water, guanine and DNA fragmentation: a density functional theory simulation

The radiation-induced damages in bio-molecules are ubiquitous processes in radiotherapy and radio-biology, and critical to space projects. In this study, we present a precise quantification of the fragmentation mechanisms of deoxyribonucleic acid (DNA) and the molecules surrounding DNA such as oxygen and water under non-equilibrium conditions using the first-principle calculations based on density functional theory (DFT). Our results reveal the structural stability of DNA bases and backbone that withstand up to a combined threshold of charge and hydrogen abstraction owing to simultaneously direct and indirect ionization processes. We show the hydrogen contents of the molecules significantly control the stability in the presence of radiation. This study provides comprehensive information on the impact of the direct and indirect induced bond dissociations and DNA damage and introduces a systematic methodology for fine-tuning the input parameters necessary for the large-scale Monte Carlo simulations of radio-biological responses and mitigation of detrimental effects of ionizing radiation.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Controlling the Self-Assembly of DNA Origami Octahedra via Manipulation of Inter-Vertex Interactions

Recent studies have demonstrated novel strategies for the organization of nanomaterials into three-dimensional (3D) ordered arrays with prescribed lattice symmetries using DNA-based self-assembly strategies. In one approach, the nanomaterial is sequestered into DNA origami frames or “material voxels” and then coordinated into ordered arrays based on the voxel geometry and the corresponding directional interactions based on its valency. While the lattice symmetry is defined by the valency of the bonds, a larger-scale morphological development is affected by assembly processes and differences in energies of anisotropic bonds. To facilely model this assembly process, we investigate the self-assembly behavior of hard particles with six interacting vertices via theory and Monte Carlo simulations and exploration of corresponding experimental systems. We demonstrate that assemblies with different 3D crystalline morphologies, but the same lattice symmetry can be formed depending on the relative strength of vertex-to-vertex interactions in orthogonal directions. We observed three distinct assembly morphologies for such systems: cube-like, sheet-like, and cylinder-like. A simple analytical theory inspired by well-established ideas in the areas of protein crystallization, based on calculating the second virial coefficient of patchy hard spheres, captures the simulation results and thus represents a straightforward means of modeling this self-assembly process. To complement the theory and simulations, experimental studies were performed to investigate the assembly of octahedral DNA origami frames with varying binding energies at their vertices. In conclusion, x-ray scattering confirms the robustness of the formed nanoscale lattices for different binding energies, while both optical and electron microscopy imaging validated the theoretical predictions on the dependence of the distinct morphologies of assembled state on the interaction strengths in the three orthogonal directions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Partial deconfinement at strong coupling on the lattice

We provide evidence for partial deconfinement — the deconfinement of a SU(M) subgroup of the SU(N) gauge group — by using lattice Monte Carlo simulations. We take matrix models as concrete examples. By appropriately fixing the gauge, we observe that the M × M submatrices deconfine. This gives direct evidence for partial deconfinement at strong coupling. We discuss the applications to QCD and holography.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Nanofabrication and Demonstration of a Direct-Write Microevaporator

Direct-write vapor deposition is a new technique that would enable one-step 3D maskless nanofabrication on a variety of substrates. A novel silicon chip-based microevaporator is developed that allows evaporant to exit through 2000–300 nm nozzles while held at distances comparable to the nozzle diameter from the substrate by a three-axis nanopositioning stage in vacuum. This results in a localized deposition on the substrate, which may be scanned relative to the substrate to produce direct-write patterns. The performance of the microevaporator is tested by creating localized depositions of various materials and the line-writing potential is demonstrated. The relationship between linewidth and source-to-substrate distance is investigated by the application of Knudsen's cosine law and Monte-Carlo simulations, and then utilized to approximate the source-to-substrate distance from performed depositions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Observation of Cosmic-Ray Anisotropy in the Southern Hemisphere with 12 yr of Data Collected by the IceCube Neutrino Observatory

We analyzed the 7.92 × 10 11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

First measurement of prompt and non-prompt D* + vector meson spin alignment in pp collisions at $\sqrt{s}$ = 13 TeV

This letter reports the first measurement of spin alignment, with respect to the helicity axis, for D* + vector mesons and their charge conjugates from charm-quark hadronisation (prompt) and from beauty meson decays (non-prompt) in hadron collisions. The measurements were performed at midrapidity (|y| < 0.8) as a function of transverse momentum (p T ) in proton–proton (pp) collisions collected by ALICE at the centre-of-mass energy $\sqrt{s}$ = 13 TeV. The diagonal spin density matrix element $ρ$00 of D* + mesons was measured from the angular distribution of the D* + → D 0 (→ K - π + )π + decay products, in the D* + rest frame, with respect to the D* + momentum direction in the pp centre of mass frame. The $ρ$00 value for prompt D* + mesons is consistent with 1/3, which implies no spin alignment. However, for non prompt D* + mesons an evidence of $ρ$00 larger than 1/3 is found. The measured value of the spin density element is $ρ$00 = 0.455±0.022(stat.)±0.035(syst.) in the 5 < p T < 20 GeV/c interval, which is consistent with a Pythia 8 Monte Carlo simulation coupled with the EvtGen package, which implements the helicity conservation in the decay of D* + meson from beauty mesons. In non-central heavy-ion collisions, the spin of the D* + mesons may be globally aligned with the direction of the initial angular momentum and magnetic field. Based on the results for pp collisions reported in this letter it is shown that alignment of non-prompt D* + mesons due to the helicity conservation coupled to the collective anisotropic expansion may mimic the signal of global spin alignment in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Atomic-scale probing of short-range order and its impact on electrochemical properties in cation-disordered oxide cathodes

Abstract Chemical short-range-order has been widely noticed to dictate the electrochemical properties of Li-excess cation-disordered rocksalt oxides, a class of cathode based on earth abundant elements for next-generation high-energy-density batteries. Existence of short-range-order is normally evidenced by a diffused intensity pattern in reciprocal space, however, derivation of local atomic arrangements of short-range-order in real space is hardly possible. Here, by a combination of aberration-corrected scanning transmission electron microscopy, electron diffraction, and cluster-expansion Monte Carlo simulations, we reveal the short-range-order is a convolution of three basic types: tetrahedron, octahedron, and cube. We discover that short-range-order directly correlates with Li percolation channels, which correspondingly affects Li transport behavior. We further demonstrate that short-range-order can be effectively manipulated by anion doping or post-synthesis thermal treatment, creating new avenues for tailoring the electrochemical properties. Our results provide fundamental insights for decoding the complex relationship between local chemical ordering and properties of crystalline compounds.

25 ENERGY STORAGE↗

DOE OSTI2024

Simulated structure and thermodynamics of decagonal Al-Co-Cu quasicrystals

Atomic structures of Al-Co-Cu decagonal quasicrystals (dQCs) are investigated using empirical oscillating pair potentials (EOPP) in molecular dynamic (MD) simulations that we enhance by Monte Carlo (MC) swapping of chemical species and replica exchange. Predicted structures exhibit planar decagonal tiling patterns and are periodic along the perpendicular direction. We then recalculate the energies of promising structures using first-principles density functional theory (DFT), along with energies of competing phases. We find that our τ -inflated sequence of QC approximants (QCAs) are energetically unstable at low temperature by at least 3 meV/atom. Extending our study to finite temperatures by calculating harmonic vibrational entropy, as well as anharmonic contributions that include chemical species swaps and tile flips, our results suggest that the quasicrystal phase is entropically stabilized at temperatures in the range 600-800 K and above. It decomposes into ordinary (though complex) crystal phases at low temperatures, including a partially disordered B2-type phase. We discuss the influence of density and composition on QC phase stability; we compare the structural differences between Co-rich and Cu-rich quasicrystals; and we analyze the role of entropy in stabilizing the quasicrystal, concluding with a discussion of the possible existence of “high entropy” quasicrystals. Published by the American Physical Society 2024

Huang, Yang (ORCID:0009000045917347)↗

DOE OSTI2023

Improved thermonuclear rate of 42 Ti( p , γ ) 43 V and its astrophysical implication in the rp process

Context. Accurate 42 Ti(p,γ) 43 V reaction rates are crucial for understanding the nucleosynthesis path of the rapid capture process (rp process) that occurs in X-ray bursts. Aims. We aim to improve the thermonuclear rates of 42 Ti(p,γ) 43 V based on more complete resonance information and a more accurate direct component, together with the recently released nuclear masses data. We also explore the impact of the newly obtained rates on the rp process. Methods. We reevaluated the reaction rate of 42 Ti(p,γ) 43 V by the sum of the isolated resonance contribution instead of the Hauser-Feshbach statistical model. We used a Monte Carlo method to derive the associated uncertainties of new rates. The nucleosynthesis simulations were performed via the NuGrid post-processing code ppn. Results. The new rates differ from previous estimations due to the use of a series of updated resonance parameters and a direct S factor. Compared with the previous results from the Hauser-Feshbach statistical model, which assumes compound nucleus 43 V with a sufficiently high-level density in the energy region of astrophysical interest, large differences exist over the entire temperature region of rp-process interest, up to two orders of magnitude. We consistently calculated the photodisintegration rate using our new nuclear masses via the detailed balance principle, and found the discrepancies among the different reverse rates are much larger than those for the forward rate, up to ten orders of magnitude at the temperature of 10 8 K. Using a trajectory with a peak temperature of 1.95×10 9 K, we performed the rp-process nucleosynthesis simulations to investigate the impact of the new rates. Our calculations show that the adoption of the new forward and reverse rates result in abundance variations for Sc and Ca of 128% and 49%, respectively, compared to the variations for the statistical model rates. On the other hand, the overall abundance pattern is not significantly affected. The results of using new rates also confirm that the rp-process path does not bypass the isotope 43 V. Conclusions. Our study found that the Hauser-Feshbach statistical model is inappropriate to the reaction rate evaluation for 42 Ti(p,γ) 43 V. The adoption of the new rates confirms that the reaction path of 42 Ti(p,γ) 43 V(p,γ) 44 Cr(β + ) 44 V is a key branch of the rp process in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Monte Carlo calculations of the extraction of scintillation light from cryogenic $N$-type GaAs

We report the high scintillation luminosity of n-type GaAs at 10 °K is surprising because (1) with a refractive index of 3.5, escape is inhibited by total internal reflection and (2) narrow-beam experiments at 90 °K report infrared absorption coefficients of several per cm. This paper presents Monte Carlo calculations showing that the high luminosity at 10 °K can be explained if (1) the narrow-beam absorption is almost all optical scattering and (2) the absolute absorption coefficient is below 0.1 per cm. Sixteen surface reflector configurations are simulated for a range of internal scattering and absolute absorption coefficients, and these can guide the design of cryogenic scintillating GaAs targets for the direct detection of dark matter. The discussion section presents a possible infrared scattering mechanism based on the metallic nature of n-type GaAs. A supplement file describes (1) the Monte Carlo program steps in detail and (2) how narrow-beam and integrating sphere experiments can measure the cryogenic optical scattering and absolute absorption coefficients.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

High average-flux laser-driven neutron source

Laser-driven neutron generation is an attractive alternative to more established methods for compact, short-pulse-duration neutron sources with applications in medical science, material science and imaging. Despite extensive investigation of various techniques, achieving performance comparable to nuclear reactors or conventional accelerators remains challenging. In this work, we generate a stable, high-repetition-rate laser-driven neutron source reaching a record average flux of 7.8 × 10 7 n/sr/s, improving on other existing laser-based sources by more than one order of magnitude. Our approach is based on a two-step process where electrons are accelerated to relativistic energies via laser wakefield acceleration (LWFA), and subsequently generate neutrons through Bremsstrahlung emission followed by photonuclear reactions in a tungsten converter. Experimental results, supported by Monte Carlo simulations, show a neutron flux of 3.0 × 10 7 n/cm 2 /s near the target, on par with some compact accelerator-based neutron sources. Additionally, a direct comparison with the target-normal sheath acceleration (TNSA) pitcher-catcher scheme, performed on the same laser system, reveals a significantly higher total neutron yield of 3.9 × 10 8 neutrons per shot, outperforming the TNSA scheme by several orders of magnitude. These findings represent a significant advancement towards the development of practical laser-driven neutron sources and highlight the advantages of LWFA-based neutron generation for future applications.

Vallières, Simon [Institut National de la Recherch↗

DOE OSTI2025

Life Cycle Cost Analysis of Prestressed Concrete Poles Subjected to Wind, Surges, and Waves

Prestressed concrete (PC) poles are becoming popular choices to support coastal power transmission systems. However, the existing literature does not offer a detailed analysis of the effectiveness of PC poles in terms of long-term vulnerabilities and the direct and indirect costs. This is due to (1) lack of fragility models for PC poles and (2) lack of probabilistic wind, storm surge, and wave models in coastal settings. In this study, we address these gaps through a series of Monte Carlo simulations to estimate fragility of PC poles as a function of age and hazard (wind, surges, and waves) intensity, and the development of a probabilistic hazard model based on 10,000 years of synthetic tropical cyclone data. The probabilistic hazard model is used in conjunction with high-resolution hydrodynamic models to generate realizations of coastal wind, storm surges, and waves for the Louisiana and Mississippi coasts. A comprehensive life cycle cost analysis for a service life of 70 years considering direct and indirect losses is conducted to compare the performance of a transmission line located in Pascagoula, Mississippi, when wood poles are replaced by PC poles. Results showed that aging has a minor effect on the reliability of PC poles, highlighting the advantages of replacing wood poles with PC poles, especially in coastal areas. In addition, PC poles are significantly more cost-effective compared with wood poles over their life cycle, leading to an estimated saving of $11.55 million (68.17% reduction). The results of this study provide key insight to inform decision-making processes to keep the coastal power grids resilient and cost-effective against future storm hazards.

natural disasters↗

DOE OSTI2020

Ethanol Conversion over La 0.7 Sr 0.3 MnO 3–x (100): Autocatalysis, Adjacent O-Vacancies, Disproportionation, and Dehydrogenation

The mechanism for catalytic conversion of ethanol over La 0.7 Sr 0.3 MnO 3– x (100) surface to acetaldehyde and ethene was investigated. Pre-exposure temperature-programmed reaction (PE-TPR) experiments were performed in which ethanol was introduced to oxidized or reduced surfaces followed by heating. In particular, sequential PE-TPR experiments were conducted to incrementally and gradually reduce the surface. The products and their ratios were investigated as a function of surface reduction. The data show that acetaldehyde and ethene production is catalyzed with hydrogen abstraction and oxygen abstraction reactions occurring by intermediates in vacancies at various temperatures >400 K. Adsorption of acetaldehyde followed by a temperature-programmed reaction does not produce ethene, indicating that acetaldehyde is not an intermediate to ethene and that the hydrogen and oxygen abstraction from ethanol to ethene are decoupled steps. Further evidence for this mechanistic nuance was obtained using isotopically labeled ethanol (CD 3 CH 2 OH), which produces CD 3 CHO and CD 2 CH 2 . Additionally, the ratio of aldehyde production to alkene production increases with reduction, suggesting that aldehyde is produced from a disproportionation reaction between ethoxy species in adjacent O-vacancies, while ethene is produced from a dehydrogenation reaction with ethoxy species in vacancies without requiring adjacent O-vacancies. Counterintuitively, this finding indicates that the more oxygenated product (aldehyde vs ethene) is favored with more vacancies and that the net alcohol conversion is autocatalytic. Density functional theory calculations were able to find the previously unknown disproportionation pathway between ethoxies in adjacent O-vacancies, and kinetic Monte Carlo simulations support this interpretation by reproducing experimental selectivities. The activation energies for these pathways are estimated as 132 ± 10 kJ/mol for the disproportionation reaction (when occurring between ethoxies in adjacent vacancies) and as 148 ± 11 kJ/mol for the direct dehydrogenation reaction of an ethoxy in a vacancy. Based on these results, a mechanism with operative pathways based on elementary steps in O-vacancies is reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Improving Detector Systematic Uncertainties Through Data-Driven Machine Learning

Detector simulation in liquid argon time projection chambers (LArTPCs) is a constant challenge. In particular the modeling of electrons response on wires is highly nontrivial. However, new machine learning techniques exist which can be leveraged to ameliorate these concerns. We present a novel methodology to attempt to learn from cosmic muon data in the ICARUS detector how reconstructed wire signals are influenced by features of the hits such as location, particle direction, angle relative to the wire plane, etc. A model can then be generated which can apply the learned mapping to Monte Carlo events to create a more data-like simulation sample. By creating such a sample we expect to reduce the systematic uncertainties at ICARUS due to our detector modeling.

Hausner, Harry [Fermilab] (ORCID:0000000188932280)↗

DOE OSTI2026

Improving Detector Systematic Uncertainties Through Data-Driven Machine Learning

Hausner, Harry [Fermilab] (ORCID:0000000188932280)↗