Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DIFFUSION”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Understanding Large-Scale Dynamos in Unstratified Rotating Shear Flows

For this work, we combine simulations with new analyses that overcome previous pitfalls to explicate how nonhelical mean-field dynamos grow and saturate in unstratified, magnetorotationally driven turbulence. Shear of the mean radial magnetic field amplifies the azimuthal component. Radial fields are regenerated by velocity fluctuations that induce shear of radial magnetic fluctuations, followed by Lorentz and Coriolis forces that source a negative off-diagonal component in the turbulent diffusivity tensor. We present a simple schematic to illustrate this dynamo growth. A different part of the Lorentz force forms a third-order correlator in the mean electromotive force that saturates the dynamo.

accretion disk & black-hole plasma↗
DOE OSTI2025

Q3 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q3: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis. • Benchmark between two first wall heat flux mapping methods, identify importance of various heat flux sources and physics impact of using fully coupled CESOL vs post-analysis evaluation. 2. Generate medium fidelity parametrized CAD. • Generate parametrized CAD components for the CAT example case via either user-defined modules called within the geometry generation or by defeatured/parametrized CAD, including DCLL blanket matched to divertor boundary and magnets. Define materials, labels, and boundary conditions for passing the mesh to CFD tools. 3. Demonstrate multiphysics magnet analysis. • Demonstrate magnet analysis workflow called from the FREDA workflow, and 4. Demonstrate nuclear analysis. • Add model to OpenFOAM and/or other codes possibly including Diablo to account for tritium diffusion in solids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Analysis of static Wilson line correlators from lattice QCD at finite temperature with T -matrix approach

The thermodynamic T-matrix approach is used to study Wilson line correlators (WLCs) for a static quark-antiquark pair in the quark-gluon plasma (QGP). Selfconsistent results that incorporate constraints from the QGP equation of state can approximately reproduce WLCs computed in 2+1-flavor lattice-QCD (lQCD), provided the input potential exhibits less screening than in previous studies. Utilizing the updated potential to calculate pertinent heavylight T-matrices we evaluate thermal relaxation rates of heavy quarks in the QGP. We find a more pronounced temperature dependence for low-momentum quarks than in our previous results (with larger screening), which turns into a weaker temperature dependence of the (temperature-scaled) spatial diffusion coefficient, in fair agreement with the most recent lQCD data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Ohmic contacts to nitrogen-incorporated n-type ultrananocrystalline diamond film grown on intrinsic single crystal diamond substrate

Nitrogen-incorporated ultrananocrystalline diamond (n-UNCD) films offer tremendous potential for diverse electronic applications. However, the absence of a reliable Ohmic contact at room temperature limits their practical integration and broader applicability. Here, in this study, we investigate the room temperature specific contact resistivity (ρ c ) of Ti/Pt/Au metal stack deposited on n-UNCD films grown on an intrinsic single crystal diamond substrate using a microwave plasma chemical vapor deposition system. We employ a circular transfer length model (c-TLM) and find the room temperature ρ c to be ∼4.67 × 10 −5 Ω cm 2 , which is among the lowest reported value for n-UNCD films. High temperature vacuum annealing conducted at 700 and 800 °C results in an initial improvement, followed by a minor degradation in ρ c values, respectively. The electrical contacts remain highly Ohmic for all measurements. Furthermore, cross-sectional transmission electron microscopy analysis suggests formation of conductive titanium carbide layer with no significant inter metallic diffusion. Overall, the electrical contacts demonstrate robust thermal stability, both of which are critical for attaining high-performance nanocrystalline diamond-based electronic devices.

Low resistance contacts↗
DOE OSTI2025

A novel closed-form inversion of the convection–diffusion equation for rapid convection, diffusion, and source profile estimation

To simplify and routinize particle transport analysis in fusion devices, a novel closed form linear inversion of the 1-D convection diffusion equation to estimate diffusion and convection profiles D(r ⃗ ), v(r ⃗ ) and source distribution s(r ⃗ ), of a single species from measured data is derived and demonstrated on synthetic data. Profile estimates of D(r ⃗ ), v(r ⃗ ), s(r ⃗ ) and their uncertainties are given as a matrix expression constructed directly from the incoming density data of the transported species in space and time, as well as physics assumptions such as particle conservation and experimental geometry. The derived matrix expression can be applied to a pumped or non-pumped recycling species, or a non-recycling species that is effectively “pumped” by plasma-facing surfaces.

Hinson, Edward [ORNL] (ORCID:000000019713140X)↗
DOE OSTI2024

3D printed carbon fiber reinforced carbon as an energy efficient alternative to graphite for EFAS tooling

As Electric Field Assisted Sintering (EFAS) gains more industrial acceptance and use, it becomes more important to develop more efficient means to implement this technology. To this aim, 3D printed continuous carbon fiber reinforced carbon (CCC) was manufactured and fabricated into tooling for EFAS systems as an alternative to traditional graphite tooling. The impact of fiber orientation on the thermal and electrical properties of the CCC was characterized. Sample material was sintered in Tokai G535 graphite tooling, under common processing conditions and compared with CCC tooling. There was nearly 50 % energy savings compared to graphite while maintaining equivalent sample density and microstructure plus keeping ram temperatures 39 % cooler. This is due to spatial control of generated heat and thermal diffusivity within the molds, by means of fiber orientation anisotropy. Finite element modeling of the tooling design supported the experimental results as well as displays the effect of optimization of this 3D printed CCC material.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A snapshot of high-entropy alloy processing techniques and their effects on resulting mechanical properties

High-entropy alloys (HEAs) exhibit exceptional strength, corrosion resistance, and thermal stability, making them promising candidates for nuclear, aerospace, and other extreme applications. While most prior work has focused on compositional design, manufacturing techniques themselves can alter microstructure and mechanical properties as dramatically as alloy chemistry. This review compiles and compares the effects of processing routes—including arc melting, induction melting, mechanical alloying with spark plasma sintering, and additive manufacturing—on the structure and properties of HEAs. Quantitative comparisons highlight, for example, that SPS-processed alloys can achieve ∼20–45 % higher yield and tensile strength than arc-melted counterparts, while Bridgman solidification produces nearly single-crystal structures with elongation to failure exceeding 80 %. Additive manufacturing routes such as selective laser melting offer fine microstructures but also introduce anisotropy and porosity, leading to yield strengths spanning 100–600 MPa for the same composition. By synthesizing such results, this review provides actionable insight into how processing routes interact with HEA core effects (high entropy, lattice distortion, sluggish diffusion, and cocktail effect) to determine performance, thereby offering a practical guide for optimizing manufacturing strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chemically Amplified, Dry-Develop Poly(aldehyde) Photoresist

The catalytic decomposition of poly(phthalaldehyde) with a photoacid generator can be used as dry-develop photoresist, where the exposed film depolymerizes into small molecules to allow the development of features via controlled vaporization. Higher temperatures enabled shorter dry-development times, but also promoted faster photoacid diffusion that compromised pattern fidelity. Trihexylamine was used as a base quencher to counteract acid diffusion in a phthalaldehyde-propanal co-polymer photoresist. The propanal co-monomer in the polymer improves the vaporization rate because it has a higher vapor pressure than phthalaldehyde. Addition of the base quencher was found to improve the contrast, pattern fidelity, and ease-of-handling of the dry-develop resist in a direct-write UV lithography tool. The dry-development of 4 μ m features was achieved with no appreciable residue. For large area features, a spatially variable exposure method was used to direct the residue away from the exposed area. The gradient exposure method was used to produce 100 μ m features. Plasma etching after dry-development was also used to achieve residue-free dry-developed patterns. These results show the benefits of incorporating base additives into a dry-develop depolymerizable resist system and highlight the need for addressing residue formation.

Materials Science↗
DOE OSTI2025

Vertically Aligned Carbon Nanotubes Grown on Copper Foil as Electrodes for Electrochemical Double Layer Capacitors

This study reports a binder-free, catalyst-free method for fabricating vertically aligned carbon nanotubes (VACNTs) directly on copper (Cu) foil using plasma-enhanced chemical vapor deposition (PECVD) for electrochemical double-layer capacitor (EDLC) applications. This approach eliminates the need for catalyst layers, polymeric binders, or substrate pre-treatments, simplifying electrode design and enhancing electrical integration. The resulting VACNTs form a dense, uniform, and porous array with strong adhesion to the Cu substrate, minimizing contact resistance and improving conductivity. Electrochemical analysis shows gravimetric specific capacitance (C grav ) and areal specific capacitance (C areal ) of 8 F g −1 and 3.5 mF cm −2 at a scan rate of 5 mV/s, with low equivalent series resistance (3.70 Ω) and charge transfer resistance (0.48 Ω), enabling efficient electron transport and rapid ion diffusion. The electrode demonstrates excellent rate capability and retains 92% of its initial specific capacitance after 3000 charge–discharge cycles, indicating strong cycling stability. These results demonstrate the potential of directly grown VACNT-based electrodes for high-performance EDLCs, particularly in applications requiring rapid charge–discharge cycles and sustained energy delivery.

25 ENERGY STORAGE↗
DOE OSTI2026

Nonlinear causality of Israel-Stewart theory with diffusion

We present the first fully nonlinear causality constraints in D = 3 + 1 dimensions for Israel-Stewart theory in the presence of energy and number diffusion in the Eckart and Landau hydrodynamic frames, respectively. These constraints are algebraic inequalities that make no assumption on the underlying geometry of the spacetime or the equation of state. In order to highlight the distinct physical and structural behavior of the two hydrodynamic frames, we discuss the special ultrarelativistic ideal gas equation of state considered in earlier literature in D = 1 + 1 dimensions, and show that our general D = 3 + 1 constraints reduce to their results upon an appropriate choice of angles. For this equation of state in both D = 1 + 1 and D = 3 + 1 dimensions one can show that: (i) there exists a region allowed by nonlinear causality in which the baryon current transitions into a spacelike vector in the Landau frame, and (ii) an analogous argument shows that the solutions of the Eckart frame equations of motion never violate the dominant energy condition, assuming nonlinear causality holds. Furthermore, we then compare our results with those from linearized Israel-Stewart theory and show that the linear causality bounds fail to capture the new physical constraints on energy and number diffusion that are successfully obtained through our nonlinear causality approach.

Quark-gluon plasma↗
DOE OSTI2025

Ohm’s Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm’s law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets with a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Self-Consistent Relativistic Electron Scattering using the Sherlock Scattering Model for X-ray Diagnostics

We present on a new, self-consistent, arbitrary-temperature Romberg integration scheme for modeling electron scattering in materials in a LANL Lagrangian Shock Hydro (LSH) code. Electron beam-target interactions are fundamental to a wide range of scientific and technological applications. When high-energy electron beams hit their target, they may scatter, deposit energy, or ionize the source. These processes govern the behavior and outcomes in nanotechnology manufacturing, electron microscopy, and modern X-ray diagnostics. Simulating these interactions is essential for interpreting experimental results, predicting material responses, and designing efficient tools and experiments. At Los Alamos, this is done using a LSH code, which is a multi-dimension, multi-material, massively parallel, multi-physics code used to simulate applications from asteroid impacts to electron beam interactions. By effectively and efficiently modeling the way that electrons scatter from the beam we can bolster these simulations and more accurately predict experimental outcomes. The model currently implemented in the LSH of interest is based on work by Papp and does not self-consistently preserve momentum in the slightly relativistic regime; here we adopt a model proposed by Braams and Karney and implement a Romberg integration scheme to compute the diffusion tensor. In this paper we will provide background on the Braams-Karney diffusion tensor as well as the Romberg integration scheme we employed to numerically solve for it. We will show that our integration scheme is accurate in solving for the set of scalar potentials used to re-express the diffusion tensor in differential form, and in solving for the diffusion coefficients in the larger LSH code. By using this diffusion tensor rather than the existing Papp one, and numerically integrating it with a Romberg method, we produce much more accurate, self-consistent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Lattice B -field correlators for heavy quarks

We analyze the color-magnetic (or “ B ”) field two-point function that encodes the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The simulations are done on fine isotropic lattices in the quenched approximation at 1.5 T c , using a range of gradient flow times for noise suppression and operator renormalization. The continuum extrapolation is performed at fixed flow time followed by a second extrapolation to zero flow time. Perturbative calculations to next-to-leading order of this correlation function, matching gradient-flowed correlators to MS ¯ , are used to resolve nontrivial renormalization issues. We perform a spectral reconstruction based on perturbative model fits to estimate the coefficient κ B of the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The approach we present here yields high-precision data for the correlator with all renormalization issues incorporated at next-to-leading order and is also applicable for actions with dynamical fermions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

How does negative triangularity mitigate ITG turbulence and transport?

Improved confinement in negative triangularity (NT) experiments is attributed to reduced fluxes driven by micro-turbulence. The physical mechanism of why thermal confinement improves in NT relative to PT is unknown. This study employs gyrokinetic flux tube simulations using the GENE code with local Miller equilibrium to elucidate the physical mechanisms behind the beneficial effects of NT flux surface shapes. The focus is on collisionless ion temperature gradient (ITG) driven turbulence with adiabatic electrons. The kinetic profiles are held fixed across a scan of triangularity values, thus enabling comparisons on a level playing field. The reduced linear growth rates for NT is shown to be due to a reduced eigenmode averaged magnetic drift frequency and a wider, stronger negative local magnetic shear region about the outboard mid-plane. The nonlinear heat flux is lower for NT than that for PT, due to reduced radial correlation length and increased correlation time (τ c ) of fluctuations. These, in turn, are due to a comparatively higher level of self-generated zero-frequency E × B zonal shearing rate ω E in NT as compared to PT. Though the linear zonal potential residual is lower for NT, the nonlinearly generated E × B zonal shearing rate is higher for NT than for PT. This outcome is linked to the distinctive features of the radial wavenumber spectra of the zonal potential and the zonal shearing rate. The dimensionless parameter ω E τ c is suggested as a figure of merit. This is higher for NT than for PT. Thus, the reduced heat diffusivity for NT is linked to increased ω E τ c . Self-generated temperature corrugations (i.e. zonal temperature gradients) are much weaker than the background mean temperature gradient. Nevertheless, temperature corrugations are more pronounced in NT than in PT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

High-Throughput In-Line Deposition of Silicon Oxide for Polycrystalline Silicon Passivating Contacts

Polycrystalline silicon passivating contacts rely on an ultrathin (1–2 nm) silicon oxide layer to minimize recombination at the wafer/oxide interface and regulate dopant diffusion. Traditionally formed by thermal or chemical oxidation, this oxide is herein replaced by silicon oxide deposited via aerosol impact-driven assembly (AIDA), enabling high wafer-per-hour throughput and precise thickness control. In this study, we show that AIDA coatings conformally cover planar or textured substrates and achieve a SiO x /poly-Si(n) structure with an implied open-circuit voltage (iV oc = 726 mV) and contact saturation current density (J 0 = 8.8 fA/cm 2 ). Furthermore, annealing AIDA SiO x films at elevated temperatures desorbs hydroxyl groups while the stoichiometry transitions toward SiO 2 , improving passivation quality. Together, these results highlight AIDA’s potential for scalable, high-throughput manufacturing of advanced passivating contacts, offering a cost-effective alternative to conventional low-pressure chemical vapor deposition and plasma-enhanced chemical vapor deposition-based silicon and oxide processes.

TOPcon↗
DOE OSTI2025

Disentangling core and edge mechanisms of the density limit in DIII-D negative triangularity plasmas

The density limit is investigated in the DIII-D negative triangularity plasmas which lack a standard H-mode edge. We find the limit may not be a singular disruptive boundary but a multifaceted density saturation phenomenon governed by distinct core and edge transport mechanisms. Sustained, non-disruptive operation is achieved at densities up to 1.8 times the Greenwald limit ($n_\mathrm{G}$) until the termination of auxiliary heating. Systematic power scans show distinct power scalings for the core ($n_e \propto P_\mathrm{SOL}^{0.27\pm0.03}$) and edge ($n_e \propto P_\mathrm{SOL}^{0.42\pm0.04}$) density limits. The edge density saturation is triggered by the onset of a non-disruptive, high-field side radiation front and the associated cooling, which clamps the edge density below . In contrast, the core density continues to rise until it saturates, a state characterized by enhanced core turbulence. Core transport evolves from a diffusive to an intermittent, avalanche-like state, as indicated by heavy-tailed probability density functions (kurtosis ≈ 6), increased Hurst exponents, and a $1/f$-type power spectrum. These findings suggest that the density limit in the low-confinement regime is determined by a combination of edge radiative cooling and core turbulent transport. This distinction provides separate targets for control strategies aimed at extending the operational space of future fusion devices.

density limit↗