Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forcing fields”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Entropy-based feature selection for capturing impacts in Earth system models with abrupt forcing

This paper presents the development of a new entropy-based feature selection method for identifying and quantifying impacts. Here, impacts are defined as statistically significant differences in spatio-temporal fields when comparing datasets with and without an external forcing in an Earth system model. Temporal feature selection is performed by first computing the cross-fuzzy entropy to quantify similarity of patterns between two datasets and then applying changepoint detection to identify regions of statistically constant entropy. The method is used to capture temperate north surface cooling from a 9-member simulation ensemble of the Mt. Pinatubo volcanic eruption, which injected 10 Tg of SO 2 into the stratosphere. The results estimate a mean difference decrease in near surface air temperature of -0.560 K with a 99% confidence interval between -0.864 K and -0.257 K between April and November of 1992, one year following the eruption. A sensitivity analysis with decreasing SO 2 injection revealed that the impact is statistically significant at 5 Tg but not at 3 Tg. Using identified features, a dependency graph model based on a 9-day lag had significantly fewer nodes than a graph based on monthly means. Furthermore, this demonstrates our method’s ability to perform dimension reduction while still uncovering source-to-impact pathways.

Changepoint detection↗
DOE OSTI2024

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Implementation and Results of Extended-MHD Post-processing of xRAGE Simulations

Magnetic fields are known to affect the evolution of conducting fluids via the J× B force, where J is the current density and B is the magnetic field. Typically, J× B in high energy density plasmas (HEDP) systems is considered small compared to the plasma pressure gradients. However, many experiments lie in a part of the parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate the so-called Biermann magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental observations of multiple modes of filamentary structures in the magnetized dusty plasma experiment (MDPX) device

Due to the small charge-to-mass ratio of dust particles, it is often necessary to use large magnetic fields of B ≥ 1 T, in order to observe the influence of magnetic forces in laboratory dusty plasmas. However, when experiments are performed at high magnetic fields in capacitively coupled, radio frequency discharges used for these dusty plasma experiments, the plasma is often observed to form filamentary structures between the electrodes that are aligned to the external magnetic field which disrupt the uniformity of the plasma and adversely impact some of our dusty plasma experiments. Recent experiments performed in the Magnetized Dusty Plasma Experiment device seek to identify and characterize these filamentary structures. This paper discusses the morphology of several distinct filamentary modes that are formed in low temperature argon plasmas. Three different spatial modes of filaments that appear to be strongly dependent on the magnetic field and neutral pressure are identified. There is strong evidence that each spatial mode has a threshold condition that is dependent on the ion Hall parameter. The criteria for the formation of the filaments are shown to be consistent with predictions of recent numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Pinning Characteristics of Zr and Hf- Added REBCO Coated Conductors Made by Advanced MOCVD in Low-to-High Magnetic Fields

BaMO 3 (M=Zr, Hf) pinning centers introduced in REBa 2 Cu 3 O 7-x (REBCO and RE = rare earth) coated conductors yield superior performance in high magnetic fields. We present the critical current density J c over a temperature range of 4.2-77 K and magnetic fields of 0-14 T (B || c-axis), of 5-15 mol.% Zr- and Hf-added REBCO 4+ μm thick film tapes fabricated by advanced metal organic chemical vapor deposition (A-MOCVD). The morphology of self-assembled BMO nanorods aligned along c-axis is found to be dependent on the Zr/Hf content. We also observe a correlation between the density of RE 2 O 3 in-plane nano-precipitates and the continuity and concentration of BMO nanorods in films of different Ba content. It is found that the (Ba+M)/Cu (M=Zr, Hf) content in REBa 2 Cu 3 O 7-x affects the shape of pinning force density curves over a wide magnetic field range and temperatures below 20 K. Another remarkable observation is the similarity in the critical current properties of Hf and Zr -added REBCO films as function of (Ba+M)/Cu content at intermediate to high range of (Ba+M)/Cu content. As a result, a quantitative analysis of pinning efficiency and correlation with the microstructure of Zr and Hf- added REBCO coated conductors is discussed.

2G high-temperature superconductor (2G-HTS)↗
DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗
DOE OSTI2023

Magnetohydrodynamics in free surface liquid metal flow relevant to plasma-facing components

Abstract While flowing Liquid Metal (LM) Plasma-Facing Components (PFCs) represent a potentially transformative technology to enable long-pulse operation with high-power exhaust for fusion reactors, Magnetohydrodynamic (MHD) drag in the conducting LM will reduce the flow speed. Experiments have been completed in the linear open-channel LMX-U device [Hvasta et al 2018 Nucl. Fusion 58 01602] for validation of MHD drag calculations with either insulating or conducting walls, with codes similar to those used to design flowing LM PFCs for a Fusion Nuclear Science Facility [Kessel et al 2019 Fusion Sci. Technol . 75 886]. We observe that the average channel flow speed decreased with the use of conducting walls and the strength of the applied transverse magnetic field. The MHD drag from the retarding Lorentz force resulted in an increase of the LM depth in the channel that ‘piled up’ near the inlet, but not the outlet. As reproduced by OpenFOAM and ANSYS CFX calculations, the magnitude and characteristics of the pileup in the flow direction increased with the applied traverse magnetic field by up to 120%, as compared to the case without an applied magnetic field, corresponding to an average velocity reduction of ∼45%. Particle tracking measurements confirmed a predicted shear in the flow speed, with the surface velocity increasing by 300%, despite the 45% drop in the average bulk speed. The MHD effect makes the bulk flow laminarized but keeps surface waves aligned along the magnetic field lines due to the anisotropy of MHD drag. The 3D fringe field and high surface velocity generate ripples around the outlet region. It was also confirmed that the MHD drag strongly depends on the conductivity of the channel walls, magnetic field, and volumetric flow rate, in agreement with the simulations and a developed analytical model. These validated models are now available to begin to determine the conditions under which the ideal LM channel design of a constant flow speed and fluid depth could be attained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

Plasma instabilities dominate radioactive transients magnetic fields: the self-confinement of leptons in Type Ia and core-collapse supernovae, and kilonovae

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta ^+$ and $\beta ^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

instabilities↗
DOE OSTI2026

Multiscale Characterization of Electrode-Induced Degradation in Perovskite Solar Cells

The stability of metal-halide-perovskite (MHP) solar cells must be understood and improved for the commercial viability of MHP technologies. Here, we apply multiscale characterization methods to study degradation modes, specifically electrode corrosion, for p-i-n MHP partial device stacks and full devices that are stored in the dark under an inert atmosphere. Our multiscale characterization approaches include full-device electro-optical performance using current-voltage (JV) curves and spatial imaging with electroluminescence (EL) and photoluminescence (PL). We further correlate interface properties using cross-sectional Kelvin probe force microscopy, which maps the nanoscale electric field properties, and electron microscopy, which demonstrates structural and chemical features. Devices stored as a full device stack degrade primarily by metal (Ag) electrode diffusion into the absorber, with formation of AgI byproducts and Ag accumulation near the indium tin oxide (ITO) contact. This causes decomposition of the perovskite absorber domains, loss of the potential drop at the electron transport layer (ETL)/perovskite interface near the metal contact, and increased equivalent resistance at the perovskite/hole transport layer (HTL) interface near the ITO contact. The devices stored without metal show a different degradation pathway dominated by corrosion of the ITO, creating voids at the ITO electrode surface with diffusion of In and Sn into the absorber. We conclude that metal electrode-induced degradation is the most severe degradation pathway under dark storage, but that ITO corrosion and absorber instability must also be mitigated. We further demonstrate mitigation of these degradation pathways by changes to the device stack, including a SnO x blocking layer at the ETL side and replacing ITO with FTO at the HTL side. These results provide a useful demonstration of specific dark degradation pathways at each electrode interface, as well as a unique multiscale example that links degradation of chemical, structural, and electrical interface properties to the full-device electro-optical characteristics.

14 SOLAR ENERGY↗
DOE OSTI2022

Element-specific first order reversal curves measured by magnetic transmission x-ray microscopy

The first-order reversal curve (FORC) method is a macroscopic measurement technique that can be used to extract quantitative and microscopic properties of hysteretic systems. Using magnetic transmission x-ray microscopy (MTXM), local element-specific FORC measurements are performed on a 20 nm thick film of CoTb. The FORCs measured with microscopy reveal a step-by-step domain evolution under the magnetic field cycling protocol and provide a direct visualization of the mechanistic interpretation of FORC diagrams. They are compared with magnetometry FORCs and show good quantitative agreement. Furthermore, the high spatial resolution and element-specific sensitivity of MTXM provide new capabilities to measure FORCs in small regions or specific phases within multicomponent systems, including buried layers in heterostructures. The ability to perform FORCs on very small features is demonstrated with the MTXM-FORC measurement of a rectangular microstructure with vortex-like Landau structures. This work demonstrates the confluence of two uniquely powerful techniques to achieve quantitative insight into nanoscale magnetic behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Equation of State of Neutron-Rich Matter at Fourth Order of Chiral Effective Field Theory and the Radius of a Medium-Mass Neutron Star

We report neutron star predictions based on our most recent equations of state. These are derived from chiral effective field theory, which allows for a systematic development of nuclear forces, order by order. We utilize high-quality two-nucleon interactions and include all three-nucleon forces up to fourth order in the chiral expansion. Our ab initio predictions are restricted to the domain of applicability of chiral effective field theory. However, stellar matter in the interior of neutron stars can be up to several times denser than normal nuclear matter at saturation, and its composition is essentially unknown. Following established practices, we extend our microscopic predictions to higher densities matching piecewise polytropes. The radius of the average-size neutron star, about 1.4 solar masses, is sensitive to the pressure at normal densities, and thus it is suitable to constrain ab initio theories of the equation of state. For this reason, we focus on the radius of medium-mass stars. We compare our results with other theoretical predictions and recent constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Electron beam irradiation effects on bulk metals: a comparative study of polycrystalline versus single-crystalline structures

This study investigates the effects of electron beam (e-beam) irradiation on the mechanical and structural properties of eight bulk metallic samples, comprising both polycrystalline (PC) and single-crystalline (SC) forms of Ni, Cr, V, and Ti. These metals were evaluated as potential candidates for beam exit windows in high-power (MW-class) particle accelerators. The primary objective is to identify metals capable of withstanding the conditions of high-power/MW-class e-beam accelerators and serve effectively as exit windows. Selection criteria were based on each metal’s intrinsic properties, power dissipation capability, and irradiation-induced changes in mechanical behavior, including hardness, elastic modulus, and defect density. Comprehensive characterization was conducted using field-emission scanning electron microscopy (FESEM), atomic force microscopy (AFM), scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM–EDS), X-ray diffraction (XRD), high-resolution transmission electron microscopy (HRTEM), selected area electron diffraction (SAED), and nanoindentation, performed both before and after exposure to a ~¿66 kGy dose from a 10 MeV e-beam accelerator. Results revealed that e-beam irradiation induced hardening in PC Ni, whereas PC Ti, commonly used in beam exit windows, exhibited softening. The observed softening in PC Ti is attributed to grain coarsening, elongation, and the formation of twins and twin boundaries, in contrast to the smaller, compressed grains in the pristine (Pr) PC Ti samples, consistent with the Hall–Petch relationship. The stresses due to twinning are small and insignificant in influencing the overall hardening of the PC Ti irradiated sample when compared to the stresses due to the dislocation density. Conversely, SC Ti samples exhibited irradiation-induced hardening. The SC Ti irradiated samples developed additional irradiation-induced modifications in crystallographic texture of (100), (101), (110), (200), (112), (004), and (211) as evidenced from the XRD results, which could probably explain the hardening effect that is caused by irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

MHD flow in liquid metal blankets: Major design issues, MHD guidelines and numerical analysis

The design of breeding blankets represents the major challenge for fusion reactor engineering because of performance requirements and severe operating conditions in terms of heat load and neutron flux. Liquid metal alloys such as lead-lithium, due to their lithium content, can be used to breed tritium, one of the plasma fuel components, and owing to their high thermal conductivity, they may serve as coolants. On the other hand, there are technical issues related to the fact that the liquid metals are electrically conducting and interact with the plasma-confining magnetic field. Induced electric currents and generated electromagnetic forces affect velocity and pressure distribution in the blankets. Magnetohydrodynamic (MHD) flows for fusion applications have been often investigated in simplified geometries, such as pipes, ducts, bends, with focus on their fundamental features. These analyses are essential, since results remain valid as background for the development of blanket designs, even when a concept is dismissed. However, the conceptual study of fusion blankets requires to take into account the global multiple effects, that arise when the full system is considered. Progress made in fusion-related MHD research results from combined numerical and experimental activities. In this paper we review and summarize features of 2D and 3D MHD flows that are typical in liquid metal blankets, together with available correlations for MHD pressure losses. This knowledge can provide simple design MHD guidelines that support a preliminary estimate of MHD effects in a blanket concept, in terms of pressure drop and flow distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

Electric-Field-Induced Assembly of an Ionic Liquid–Water Interphase Enables Efficient Heavy Metal Electrosorption

Controlling ion desolvation, transport, and charge transfer at the electrode-electrolyte interface (EEI) is critical to enabling the rational design of efficient and selective separation of targeted heavy metals and decontamination of industrial waste water. One challenge is to sufficiently resolve and interrogate the intermediate transformation steps between solvated metal cations and bare cations for electroreduction at the EEI and establish pathways to modulate these steps to achieve efficient energy transfer for targeted reactive separations. Herein, we obtained a predictive understanding of the effect of adding a hydrophobic ionic liquid, 1-ethyl-3-methylimidazolium chloride (EMIMCl), to an aqueous electrolyte on modulating the desolvation and electrosorption of Pb 2+ cations using a combination of experimental and theoretical techniques. Here, we reveal formation of a compact interphase layer consisting of EMIMCl-Pb complex clusters under an applied electric field using operando electrochemical Raman spectroscopy, atomic force microscopy, and electrochemical impedance spectroscopy measurements combined with classical molecular dynamics simulations. The application of a larger negative potential is shown to result in formation of a well oriented layer with the positive imidazolium ring of the EMIMCl lying on the electrode and the long hydrophobic alkyl chain extending into the bulk electrolyte. This oriented layer is demonstrated to facilitate desolvation of incoming solvated Pb 2+ cations and decrease the charge transfer resistance for Pb electrodeposition, which has important implications for selective removal of Pb from contaminated mixtures. Overall, our studies open up new opportunities to modulate ion desolvation using hydrophobic ionic liquids in aqueous electrolytes for efficient heavy metal separation.

36 MATERIALS SCIENCE↗