Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Atomic masses with machine learning for the astrophysical r process

The astrophysical r process plays a vital role in the production of heavy elements. Modeling of the r process is sensitive to masses and further requires knowledge of masses beyond current experimental reach. Therefore, simulations of the r process offer a unique test bed for predicting mass extrapolations. We take a Machine-Learning (ML) approach to model the masses across the entire chart of nuclides. For the first time, we simulate r-process nucleosynthesis with a physics-based ML mass model. We compare simulated abundances to solar data in order to evaluate the model's performance far from stability. The resulting r-process abundances up to thorium and uranium qualitatively match those of the observed solar system abundance pattern, with the characteristic peaks well positioned. We propagate the mass uncertainties obtained from the ML model to r-process abundance yields to estimate an uncertainty band associated with our approach. The size of the uncertainty band is approximately one order of magnitude which aligns with the uncertainty reported using alternative techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗
DOE OSTI2022

Reconciling Observed and Predicted Tropical Rainforest OH Concentrations

We present OH observations made in Amazonas, Brazil during the Green Ocean Amazon campaign (GoAmazon2014/5) from February to March of 2014. The average diurnal variation of OH peaked with a midday (10:00–15:00) average of 1.0 × 10 6 (±0.6 × 10 6 ) molecules cm -3 . This was substantially lower than previously reported in other tropical forest photochemical environments (2–5 × 10 6 molecules cm -3 ) while the simulated OH reactivity was lower. The observational data set was used to constrain a box model to examine how well current photochemical reaction mechanisms can simulate observed OH. We used one near-explicit mechanism (MCM v3.3.1) and four condensed mechanisms (i.e., RACM2, MOZART-T1, CB05, CB6r2) to simulate OH. A total of 14 days of analysis shows that all five chemical mechanisms were able to explain the measured OH within instrumental uncertainty of 40% during the campaign in the Amazonian rainforest environment. Future studies are required using more reliable NO x and VOC measurements to further investigate discrepancies in our understanding of the radical chemistry in the tropical rainforest.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Solar cell contacts: quantifying the impact of interfacial layers on selectivity, recombination, charge transfer, and V oc

Interfacial layers (IFLs) are ubiquitous in solar cells, but their precise impact on carrier transfer rates and the relation of these rates to performance metrics and the concepts of selectivity and recombination is lacking. We report the use of a well-defined interdigitated back-contact (IBC) silicon solar cell to determine the precise role of such IFLs. We characterize the action of the common IFL spiro-OMeTAD by making it a third contact to the IBC cell. This architecture creates three solar cells within a single structure that, with numerical simulation, provide the exchange current densities (i.e., charge transfer rates) for electrons ( J 0n ) and holes ( J 0p ) and the quasi Fermi-level splitting in the absorber, which measures the balance of generation and recombination. Further, we describe the relation of V oc to contact recombination, the asymmetry in electron/hole collection rates at a single contact (contact selectivity), and the asymmetry in collection rates of the same carrier at separate contacts (carrier selectivity). Relative to bare gold, neat spiro-OMeTAD reduces J 0n and J 0p (their geometric mean ( J 0n J 0p ) 0.5 decreases by 10 4 ), decreasing contact recombination. Addition of the common dopant Li-TFSI and air increase J 0p / J 0n by 10 6 with little effect on ( J 0n J 0p ) 0.5 , increasing contact hole selectivity. The significant increase in V oc observed by introducing spiro-OMeTAD/Li-TFSI IFLs into the cells studied, however, is due to an increase in the carrier selectivity rather than the contact selectivity or recombination of the spiro-OMeTAD-modified contact. Operando measurements further show voltage-dependent changes in the J 0 s, demonstrating that spiro-OMeTAD contributes to hysteresis. Broader context : Solar cells operate by photogenerating excess charge carriers in an absorber material and, in competition with recombination, asymmetrically extracting them at so-called carrier-selective contacts, one that ideally collects only electrons and the other, only holes. Particularly in emerging photovoltaics such as perovskites, thin layers of organic semiconductors or related materials are introduced between the absorber and contact to improve power conversion efficiency. In terms of interfacial charge transfer, a prevailing view is that such interfacial layers improve performance by helping block the collection of the undesired carrier, considered a form of recombination. We use a novel platform to study the simultaneous impact of spiro-OMeTAD, a common interfacial layer, on the collection of electrons, the collection of holes, and on the recombination of electrons and holes. We quantitatively demonstrate that spiro-OMeTAD layers indeed passivate the contact toward recombination, but that larger improvements in the open-circuit voltage, a key cell metric, can come not only from this but also from increasing the asymmetry of the collection of electrons in the system. Further, operando measurements show transient changes in the properties of spiro-OMeTAD which suggest that it contributes to hysteresis phenomena commonly observed in perovskite and other solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Understanding and predicting damage and failure at grain boundaries in BCC Ta

Understanding the effect of grain boundaries (GBs) on the deformation and spall behavior is critical to designing materials with tailored failure responses under dynamic loading. This understanding is hampered by the lack of in situ imaging capability with the optimum spatial and temporal resolution during dynamic experiments, as well as by the scarcity of a systematic data set that correlates boundary structure to failure, especially in BCC metals. To fill in this gap in the current understanding, molecular dynamics simulations are performed on a set of 74 bi-crystals in Ta with a [110] symmetric tilt axis. Our results show a correlation between GB misorientation angle and spall strength and also highlight the importance of GB structure itself in determining the spall strength. Specifically, we find a direct correlation between the ability of the GB to plasticity deform through slip/twinning and its spall strength. Additionally, a change in the deformation mechanism from dislocation-meditated to twinning-dominated plasticity is observed as a function of misorientation angles, which results in lowered spall strengths for high-angle GBs.

Chen, J. (ORCID:000000018633229X)↗
DOE OSTI2020

Perspectives on the generation of electron beams from plasma-based accelerators and their near and long term applications

This article first gives the authors' perspectives on how the field of plasma-based acceleration (PBA) developed and how the current experiments, theory, and simulations are motivated by long term applications of PBA to a future linear collider and an x-ray free electron laser. Here, we then focus on some early applications that will likely emerge from PBA research such as electron beam radiotherapy, directional but incoherent x-ray beams for science and technology, near single cycle continuously tunable infrared pulses for spectroscopy, and non-perturbative quantum electrodynamics enabled by PBA electron beams. In our opinion, these near term applications could be developed within the next decade with a concerted effort by the community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulation of an inductively coupled plasma with a two-dimensional Darwin particle-in-cell code

A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Classification of dynamical Lie algebras generated by spin interactions on undirected graphs

Dynamical Lie algebras (DLAs) are a versatile tool for various topics that span from the expressibility-trainability of variational quantum algorithms (VQAs), to simulation of many body Hamiltonians. Quantum gates and most of the Hamiltonians of interest consist of local interactions; therefore, the analysis of all possible DLAs generated by 1- and 2-local operators is crucial for quantum simulation and VQAs on current hardware. Previously in [R. Wiersema et al ., npj Quantum Inf. 10 , 110 (2024)], we analyzed the DLAs on linear, circular and all-to-all topologies, and obtained results about their dimensions and algebraic structure. Here, in this work, we extend our analysis into any possible hardware topology and provide a classification of all DLAs generated by Pauli strings on any undirected interaction graph. Our results indicate that the DLAs depend solely on whether the connectivity or interaction graph is bipartite or not. In addition, we find that the non-trivial polynomially scaling DLAs appear only on 1D line or circle topologies, and all other DLAs have dimensions scaling exponentially with the system size. Together with the current VQA literature, our results imply that either the majority of VQAs are non-trainable, or we are yet to understand the role of DLAs on the trainability of VQAs.

Algebraic structures↗
DOE OSTI2020

A Bayesian analysis of nuclear deformation properties with Skyrme energy functionals

In spite of numerous scientific and practical applications, there is still no comprehensive theoretical description of the nuclear fission process based solely on protons, neutrons and their interactions. Furthermore, the most advanced simulations of fission are currently carried out within nuclear density functional theory (DFT). In spite of being fully quantum-mechanical and rooted in the theory of nuclear forces, DFT still depends on a dozen or so parameters characterizing the energy functional. Calibrating these parameters on experimental data results in uncertainties that must be quantified for applications. This task is very challenging because of the high computational cost of DFT calculations for fission. In this paper, we use Gaussian processes to build emulators of DFT models in order to quantify and propagate statistical uncertainties of theoretical predictions for a range of nuclear deformations relevant to describing the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

3D numerical investigation of a free-burning argon arc with metal electrodes using a novel sheath coupling procedure

In this study, a novel coupling procedure describing the complicated plasma-electrode interaction process has been developed and applied into the 3D finite volume simulation of a direct current tungsten inert gas (TIG) welding system. This is achieved by making the space charge layer (sheath) incorporated into the computation domain and interact with both bulk plasma and cathode through the effective electrical conductivity. Both chemical and thermal nonequilibrium phenomena as well as the self-induced magnetic fields have been taken into consideration by the model to ensure a realistic numerical description of a non-thermal arc. The applicability of this coupling procedure is further improved by calculating the real electric potential, which is capable of accounting for the effects of the complicated drift and diffusion processes. Numerical results of both 100 and 200A discharge currents are presented, field reversal is obtained at near-anode regions in both cases, which is followed by the negative anode sheath potential drop. The region of the strongest electron overpopulation appears at the intersection of plasma fringes and electrode surface. Finally, numerical results of plasma temperature and voltage show good agreement with experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Atacama Cosmology Telescope: DR6 power spectrum foreground model and validation

We discuss the model of astrophysical emission at millimeter wavelengths used to characterize foregrounds in the multi-frequency power spectra of the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6), expanding on Louis et al. (2025) (2503.14452). We detail several tests to validate the capability of the DR6 parametric foreground model to describe current observations and complex simulations, and show that cosmological parameter constraints are robust against model extensions and variations. We demonstrate consistency of the model with pre-DR6 ACT data and observations from Planck and the South Pole Telescope. We evaluate the implications of using different foreground templates and extending the model with new components and/or free parameters. In all scenarios, the DR6 ΛCDM and ΛCDM+N eff cosmological parameters shift by less than 0.5σ relative to the baseline constraints. Some foreground parameters shift more; we estimate their systematic uncertainties associated with modeling choices. From our constraint on the kinematic Sunyaev-Zel'dovich power, we obtain a conservative limit on the duration of reionization of Δz rei < 4.4, assuming a reionization midpoint consistent with optical depth measurements and a minimal low-redshift contribution, with varying assumptions for this component leading to tighter limits. Finally, we analyze realistic non-Gaussian, correlated microwave sky simulations containing Galactic and extragalactic foreground fields, built independently of the DR6 parametric foreground model. Processing these simulations through the DR6 power spectrum and likelihood pipeline, we recover the input cosmological parameters of the underlying cosmic microwave background field, a new demonstration for small-scale CMB analysis. These tests validate the robustness of the ACT DR6 foreground model and cosmological parameter constraints.

CMBR experiments↗
DOE OSTI2023

Magnetic island formation and rotation braking induced by low-Z impurity penetration in an EAST plasma

Abstract Recent observations of the successive formations of the 4 / 1 , 3 / 1 , and 2 / 1 magnetic islands as well as the subsequent braking of the 2 / 1 mode during a low- Z impurity penetration process in EAST experiments are well reproduced in our 3 D resistive MHD simulations. The enhanced parallel current perturbation induced by impurity radiation predominately contributes to the tearing mode growth, and the 2 / 1 island rotation is mainly damped by the impurity accumulation as results of the influence from high n modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

What if GW190425 did not produce a black hole promptly?

ABSTRACT It is widely believed that the binary neutron star merger GW190425 produced a black hole promptly upon merger. Motivated by the potential association with the fast radio burst FRB 20190425A, which took place 2.5 h after the merger, we revisit the question of the outcome of GW190425 by means of numerical relativity simulations. We show that current laboratory and astrophysical constraints on the equation of state of dense matter do not rule out the formation of a long-lived remnant. However, the formation of a stable remnant would have produced a bright kilonova, in tension with upper limits by ZTF at the location and time of FRB 20190425A. Moreover, the ejecta would have been optically thick to radio emission for days to months, preventing a putative FRB from propagating out. The predicted dispersion measure is also several orders of magnitude larger than that observed for FRB 20190425A. Our results indicate that FRB 20190425A and GW190425 are not associated. However, we cannot completely rule out the formation of a long-lived remnant, due to the incomplete coverage of the relevant sky regions. More observations of GW190425-like events, including potential upper limit, have the potential to constrain nuclear physics. To this aim, it is important that follow-up observational campaigns of gravitational wave events are informed by the properties of the source, such as their chirp mass, and we urge the LIGO-Virgo-KAGRA collaboration to promptly release them publicly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Finite-temperature tensor network study of the Hubbard model on an infinite square lattice

Here, the Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network—an infinite projected entangled pair state—evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With U⁡(1)×U⁡(1) symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators—measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit—provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Design and validation of a fission-product source for radioactive beams

A neutron-generator-based source of fission products (nuCARIBU) is being developed to replace the existing source at Argonne National Laboratory. Neutrons generated by reaction are thermalized in a compact polyethylene [(moderator and impinge on a thin 235 U foil to produce fission products. Experiments have been conducted to verify the nuCARIBU concept and simulation results using a 4 MeV, 0.16 µ⁢A proton beam from the ATLAS accelerator irradiating a LiF neutron production target and two different simplified versions of the nuCARIBU (C 2 ⁢H 4 ) n moderator with thicknesses of 3.81 cm and 1.91 cm separating it from the enriched 235 U fission foil. Here, the emitted γ-ray peaks of 1383.93 keV ( 92 Sr), 358.0 keV ( 104 Tc), 565.992 keV ( 134 Te), and 258.411 keV ( 138 Xe) obtained from fission products are used to determine the neutron-induced fission rates in the fission foil. The measured fission rates are in good agreement with those simulated using the PHITS Monte Carlo code with user-defined neutron production cross sections based on well-known cross section data. These experimental results, consistent with the simulations, show that the current 252 Cf source (CARIBU) can be replaced by the neutron induced fission source (nuCARIBU).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Distribution of blackouts in the power grid and the Motter and Lai model

Carreras, Dobson, and colleagues have studied empirical data on the sizes of the blackouts in real grids and modeled them with computer simulations using the direct current approximation. They have found that the resulting blackout sizes are distributed as a power law and suggested that this is because the grids are driven to the self-organized critical state. In contrast, more recent studies found that the distribution of cascades is bimodal resulting in either a very small blackout or a very large blackout, engulfing a finite fraction of the system. Here we reconcile the two approaches and investigate how the distribution of the blackouts changes with model parameters, including the tolerance criteria and the dynamic rules of failure of the overloaded lines during the cascade. Finally, we study the same problem for the Motter and Lai model and find similar results, suggesting that the physical laws of flow on the network are not as important as network topology, overload conditions, and dynamic rules of failure.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Two-component dynamics and the liquid-like to gas-like crossover in supercritical water

Molecular-scale dynamics in sub- to super-critical water is studied with inelastic X-ray scattering and molecular dynamics simulations. The obtained longitudinal current correlation spectra can be decomposed into two main components: a low-frequency (LF), gas-like component and a high-frequency (HF) component arising from the O--O stretching mode between hydrogen-bonded molecules, reminiscent of the longitudinal acoustic mode in ambient water. With increasing temperature, the hydrogen-bond network diminishes and the spectral weight shifts from HF to LF, leading to a transition from liquid-like to gas-like dynamics with rapid changes around the Widom line.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗