Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explainable”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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

Design nanoporous metal thin films via solid state interfacial dealloying

Thin-film solid-state interfacial dealloying (thin-film SSID) is an emerging technique to design nanoarchitecture thin films. The resulting controllable 3D bicontinuous nanostructure is promising for a range of applications including catalysis, sensing, and energy storage. Using a multiscale microscopy approach, we combine X-ray and electron nano-tomography to demonstrate that besides dense bicontinuous nanocomposites, thin-film SSID can create a very fine (5–15 nm) nanoporous structure. Not only is such a fine feature among one of the finest fabrications by metal-agent dealloying, but a multilayer thin-film design enables creating nanoporous films on a wider range of substrates for functional applications. Through multimodal synchrotron diffraction and spectroscopy analysis with which the materials’ chemical and structural evolution in this novel approach is characterized in details, we further deduce that the contribution of change in entropy should be considered to explain the phase evolution in metal-agent dealloying, in addition to the commonly used enthalpy term in prior studies. Here, the discussion is an important step leading towards better explaining the underlying design principles for controllable 3D nanoarchitecture, as well as exploring a wider range of elemental and substrate selections for new applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural evolution during solution-based synthesis of Li 7 P 3 S 11 solid electrolyte by synchrotron X-ray total scattering

Solution-based synthesis has emerged as a promising approach for scalable and efficient synthesis of solid electrolytes (SEs). However, the ionic conductivity of solution-synthesized SEs is much lower than that of SEs synthesized by solid state approaches. The structural evolution during the synthesis of SEs is crucial to explaining the decreased ionic conductivities but has not been well understood. Here, in this work, we use synchrotron X-ray total scattering combined with pair distribution function (PDF) analysis to understand the structural evolution during annealing of the solution-processed precursor for the synthesis of Li 7 P 3 S 11 SEs. Combining Bragg diffraction and diffuse scattering, X-ray total scattering enables detection of structural information including phase transition and crystallinity. The results show that the formation of Li 7 P 3 S 11 glass-ceramic starts at 230 °C. Increasing the annealing temperature helps to improve the crystallinity, but the degree of crystallization is still lower than 70.3% even after annealing at 290 °C when thermal decomposition of Li 7 P 3 S 11 to Li 4 P 2 S 6 has occurred. PDF analysis also revealed the existence of acetonitrile used in the solution-based synthesis in the SEs after annealing at high temperatures. The results rationalize the synthesis–property relationship, explain why annealing at 290 °C provides the highest ionic conductivity at the studied temperatures due to an optimal balance for achieving a high crystallinity while suppressing severe decomposition of Li 7 P 3 S 11 , and provide novel insights to develop approaches for synthesizing superior Li 7 P 3 S 11 SEs on a large scale.

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

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Wang, Xiaoyu↗
DOE OSTI2023

On the Low Risk of SSR in Type III Wind Turbines Operating in Grid-Forming Control: Preprint

We have shown in a previous work that the risk of sub-synchronous resonance (SSR) between a wind power plant with Type III wind turbines and series-compensated transmission lines is low when the wind turbines in the plant are operated in grid-forming mode, instead of the standard grid-following mode. The fundamental mechanism behind the improved damping characteristics is explained in this paper by modeling the positive and negative sequence impedances of Type III wind turbines for GFM operation mode. It is discovered that the GFM control naturally acts against the negative resistance of Type III wind turbines that results from an interaction between the proportional gain of the rotor-side converter current controller and negative slip. The developed sequence impedance models and improved damping behavior are verified using PSCAD simulations of a 2.5-MW Type III GFM wind turbine. The modeling predictions are also supported by experimentally measuring the sequence impedance response of a 2.5-MW Type III wind turbine during operation in GFL and GFM modes.

grid-following turbine↗
DOE OSTI2021

Accuracy of environmental tracers and consequences for determining the Type Ia supernova magnitude step

Type Ia Supernovae (SNe Ia) are standardizable candles that allow us to measure the recent expansion rate of the Universe. Due to uncertainties in progenitor physics, potential astrophysical dependencies may bias cosmological measurements if not properly accounted for. The dependency of the intrinsic luminosity of SNe Ia with their host-galaxy environment is often used to standardize SNe Ia luminosity and is commonly parameterized as a step function. This functional form implicitly assumes two-populations of SNe Ia. In the literature, multiple environmental indicators have been considered, finding different, sometimes incompatible, step function amplitudes. We compare these indicators in the context of a two-populations model, based on their ability to distinguish the two populations. We show that local Hα-based specific star formation rate (lsSFR) and global stellar mass are better tracers than, for instance, host galaxy morphology. We show that tracer accuracy can explain the discrepancy between the observed SNe Ia step amplitudes found in the literature. Using lsSFR or global mass to identify the two populations can explain all other observations, though lsSFR is favoured. As lsSFR is strongly connected to age, our results favour a prompt and delayed population model. In any case, there exists two populations that differ in standardized magnitude by at least 0.121 ± 0.010 mag.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Characterization and dynamics of the peculiar stream Jhelum: A tentative role for the Sagittarius dwarf galaxy

Stellar streams are a promising tool to study the Milky Way’s dark matter subhalo population, as interactions with subhalos are expected to leave visible imprints in the streams in the form of substructure. However, there may be other causes of substructure. Here we studied the kinematics and the unusual morphology of the stellar stream Jhelum. Using a combination of ground-based photometry and Gaia EDR3 astrometry, we characterized the morphology of Jhelum. We combined this new data with radial velocities from the literature to perform orbit integrations of the stream in static Galactic potentials. We also carried out N-body simulations in the presence of the Sagittarius dwarf galaxy. The new data reveal a previously unreported tertiary component in the stream, as well as several gaps and a kink-like feature in its narrow component. We find that for a range of realistic Galactic potentials, no single orbit is able to reproduce Jhelum’s radial velocity data entirely. A generic property of the orbital solutions is that they share a similar orbital plane to Sagittarius and this leads to repeated encounters with the stream. Using N-body simulations that include a massive Sagittarius, we explored its effect on Jhelum, and we show that these encounters are able to qualitatively reproduce the narrow and broad components in Jhelum, as well as create a tertiary component in some cases. We also find evidence that such encounters can result in an apparent increase in the velocity dispersion of the stream by a factor up to four due to overlapping narrow and broad components. Our findings suggest that the Jhelum stream is even more complex than once thought; however, its morphology and kinematics can tentatively be explained via the interactions with Sagittarius. In this scenario, the formation of Jhelum’s narrow and broad components occurs naturally, yet some of the smaller gap-like features remain to be explained.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The persistent shadow of the supermassive black hole of M87. II. Model comparisons and theoretical interpretations

The Event Horizon Telescope (EHT) observation of M87∗ in 2018 has revealed a ring with a diameter that is consistent with the 2017 observation. The brightest part of the ring is shifted to the southwest from the southeast. In this paper, we provide theoretical interpretations for the multi-epoch EHT observations for M87∗ by comparing a new general relativistic magnetohydrodynamics model image library with the EHT observations for M87∗ in both 2017 and 2018. The model images include aligned and tilted accretion with parameterized thermal and nonthermal synchrotron emission properties. The 2018 observation again shows that the spin vector of the M87∗ supermassive black hole is pointed away from Earth. A shift of the brightest part of the ring during the multi-epoch observations can naturally be explained by the turbulent nature of black hole accretion, which is supported by the fact that the more turbulent retrograde models can explain the multi-epoch observations better than the prograde models. The EHT data are inconsistent with the tilted models in our model image library. Assuming that the black hole spin axis and its large-scale jet direction are roughly aligned, we expect the brightest part of the ring to be most commonly observed 90 deg clockwise from the forward jet. This prediction can be statistically tested through future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗
DOE OSTI2020

Direct Oxidation of Primary Alcohols to Carboxylic Acids

Oxidation of primary alcohols to carboxylic acids is a fundamental transformation in organic chemistry, yet despite its simplicity, extensive use, and relationship to pH, it remains a subject of active research for synthetic organic chemists. Since 2013, a great number of new methods have emerged that utilize transition-metal compounds as catalysts for acceptorless dehydrogenation of alcohols to carboxylates. The interest in this reaction is explained by its atom economy, which is in accord with the principles of sustainability and green chemistry. Furthermore, the methods for the direct synthesis of carboxylic acids from alcohols is ripe for a modern survey, which we provide in this review.

acceptorless dehydrogenation↗
DOE OSTI2020

Axial magnetic flux amplification in Hall-magnetohydrodynamic simulations of externally magnetized z-pinches

It is shown that externally magnetized gas-puff z-pinches under compression can exhibit robust axial flux amplification under suitable axial boundary conditions. This effect relies upon the Hall term in the generalized Ohm's law to generate azimuthal currents in the presence of a driving axial electric field. Under dynamic compression, the total current tends to flow in a mostly force-free boundary layer, separating the predominately azimuthal field outside the layer and the predominately axial field inside the layer. The effect only appears to occur if the axial boundaries allow for outflow or absorption of electromagnetic energy. The effect is mitigated by imposing either periodic or conducting axial boundary conditions. A semi-analytic equilibrium analysis agrees with steady-state solutions of the time-dependent electron-magnetohydrodynamic equations and provides an estimate of the scaling of the boundary layer as well as suggesting a scenario for the formation of the boundary layer. When operative, the effect can significantly impede plasma compression due to the increase in axial flux that diffuses into the pre-compressed plasma or through the presence of conductors that inhibit movement of the generated axial flux. Several facilities have noted unusual implosion behavior in z-pinch experiments with applied axial fields that does not appear to be explainable within the standard magnetohydrodynamic model. In conclusion, it is suggested that these experiments can be explained by the axial flux amplification and concentration effect

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Cryogenic operation of planar ultrananocrystalline diamond field emission source in SRF injector

Here, we report a demonstration of electron beam generation in a superconducting radio frequency (SRF) Teraelectronvolt Energy Superconducting Linear Accelerator 1.3 GHz gun equipped with a field emission cathode and operated at 2 K. The cathode is a submicrometer film of nitrogen-incorporated ultrananocrystalline diamond [(N)UNCD] deposited atop a Nb RRR300 cathode plug that is 4 mm in diameter. Measurements of the output current showed that it increased exponentially as a function of the gun gradient. Our results demonstrate a feasible path toward simplified fully cryogenic SRF injector technology. At a maximal gradient of 0.9 MV/m, a current density of 22 μA/cm2 was obtained [calculated as a 2.75 μA current over the total cathode surface area of π(2 mm)2]. One important finding is that the electron emitter made of (N)UNCD, a material long been known as a highly efficient field emission material, demonstrated a record low turn-on gradient of 0.6 MV/m. A hypothesis explaining this behavior is proposed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Analytical model of low and high ablation regimes in carbon arcs

Graphite ablation by an electric arc or a laser/solar flux is widely used for the synthesis of carbon nanomaterials. Previously, it was observed in multiple arc experiments that the ablation rate is a strong nonlinear function of the arc current and it drastically increases at some threshold current. As such, we developed an analytical model explaining this transition in the rate of ablation by an electric arc or a laser/solar flux. The model not only explains the observations but can also accurately predict the experimentally observed ablation rates. The model takes into account redeposition of carbon back to the ablated surface, which is the key process responsible for the observed effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Wave-driven torques to drive current and rotation

In the classic Landau damping initial value problem, where a planar electrostatic wave transfers energy and momentum to resonant electrons, a recoil reaction occurs in the nonresonant particles to ensure momentum conservation. To explain how net current can be driven in spite of this conservation, the literature often appeals to mechanisms that transfer this nonresonant recoil momentum to ions, which carry negligible current. However, this explanation does not allow the transport of net charge across magnetic field lines, precluding E X B rotation drive. In this study, we show that in steady state, this picture of current drive is incomplete. Using a simple Fresnel model of the plasma, we show that for lower hybrid waves, the electromagnetic energy flux (Poynting vector) and momentum flux (Maxwell stress tensor) associated with the evanescent vacuum wave become the Minkowski energy flux and momentum flux in the plasma and are ultimately transferred to resonant particles. Thus, the torque delivered to the resonant particles is ultimately supplied by the electromagnetic torque from the antenna, allowing the nonresonant recoil response to vanish and rotation to be driven. We present a warm fluid model that explains how this momentum conservation works out locally, via a Reynolds stress that does not appear in the one-dimensional initial value problem. This model is the simplest that can capture both the nonresonant recoil reaction in the initial-value problem, and the absence of a nonresonant recoil in the steady-state boundary value problem, thus forbidding rotation drive in the former while allowing it in the latter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Laboratory study of the failed torus mechanism in arched, line-tied, magnetic flux ropes

Coronal mass ejections (CMEs) are some of the most energetic and violent events in our solar system. The prediction and understanding of CMEs are of particular importance due to the impact that they can have on Earth-based satellite systems and, in extreme cases, ground-based electronics. CMEs often occur when long-lived magnetic flux ropes (MFRs) anchored to the solar surface destabilize and erupt away from the Sun. One potential cause for these eruptions is an ideal magnetohydrodynamic (MHD) instability, such as the kink or torus instability. Previous experiments on the magnetic reconnection experiment revealed a class of MFRs that were torus-unstable but kink-stable, which failed to erupt. These “failed-tori” went through a process similar to Taylor relaxation, where the toroidal current was redistributed before the eruption ultimately failed. Herein we have investigated this behavior through additional diagnostics that measure the current distribution at the foot points and the energy distribution before and after an event. These measurements indicate that ideal MHD effects are sufficient to explain the energy distribution changes during failed torus events. This excludes Taylor relaxation as a possible mechanism of current redistribution during an event. A new model that only requires non-ideal effects in a thin layer above the electrodes is presented to explain the observed phenomena. This work broadens our understanding of the stability of MFRs and the mechanism behind the failed torus through the improved prediction of the torus instability and through new diagnostics to measure the energy inventory and current profile at the foot points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Theoretical examination of nuclear spin diffusion in light-induced spin coherences in photosystem I

Microsecond coherence times are predicted in spin dynamics calculations of a donor–acceptor electron spin pair P + A 1A - created upon light excitation of photosystem I (PSI). The effects of nuclear spin diffusion (NSD) due to the anisotropic protein environment on the predicted coherence times T M are studied. Closely positioned pairs and triples of protons located 5–8 Å from the electron spin are shown to largely control T M . Additionally, knowledge of the crystal structure of PSI allows for spin dynamics calculations where specific cofactors and amino acid residues are removed or replaced and the identification of anisotropic environmental features controlling electronic decoherence. Finally, we show that NSD alone cannot explain the >3× shorter experimentally observed coherence times and suggest that methyl groups at key protein sites may explain this discrepancy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Thinking Bayesian for plasma physicists

Bayesian statistics offers a powerful technique for plasma physicists to infer knowledge from the heterogeneous data types encountered. To explain this power, a simple example, Gaussian Process Regression, and the application of Bayesian statistics to inverse problems are explained. The likelihood is the key distribution because it contains the data model, or theoretic predictions, of the desired quantities. By using prior knowledge, the distribution of the inferred quantities of interest based on the data given can be inferred. Because it is a distribution of inferred quantities given the data and not a single prediction, uncertainty quantification is a natural consequence of Bayesian statistics. The benefits of machine learning in developing surrogate models for solving inverse problems are discussed, as well as progress in quantitatively understanding the errors that such a model introduces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗