Search NASA⌕ Search

SEARCH · Search NASA

Results for “Temperature Gradient”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Bifurcation physics of magnetic islands and stochasticity explored by heat pulse propagation studies in toroidal plasmas

Bifurcation physics of the magnetic island was investigated using the heat pulse propagation technique produced by the modulation of electron cyclotron heating. There are two types of bifurcation phenomena observed in LHD and DIII-D. One is a bifurcation of the magnetic topology between nested and stochastic fields. The nested state is characterized by the bi-directional (inward and outward) propagation of the heat pulse with slow propagation speed. The stochastic state is characterized by the fast propagation of the heat pulse with electron temperature flattening. The other bifurcation is between magnetic island with larger thermal diffusivity and that with smaller thermal diffusivity. The damping of toroidal flow is observed at the O-point of the magnetic island both in helical plasmas and in tokamak plasmas during a mode locking phase with strong flow shears at the boundary of the magnetic island. Associated with the stochastization of the magnetic field, the abrupt damping of toroidal flow is observed in LHD. The toroidal flow shear shows a linear decay, while the ion temperature gradient shows an exponential decay. Lastly, this observation suggests that this flow damping is due to the change in the non-diffusive term of momentum transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

A Monte Carlo model for 3D grain evolution during welding

Welding is one of the most wide-spread processes used in metal joining. However, there are currently no open-source software implementations for the simulation of microstructural evolution during a weld pass. Here we describe a Potts Monte Carlo based model implemented in the SPPARKS kinetic Monte Carlo computational framework. The model simulates melting, solidification and solid-state microstructural evolution of material in the fusion and heat-affected zones of a weld. The model does not simulate thermal behavior, but rather utilizes user input parameters to specify weld pool and heat-affect zone properties. Weld pool shapes are specified by Bezier curves, which allow for the specification of a wide range of pool shapes. Pool shapes can range from narrow and deep to wide and shallow representing different fluid flow conditions within the pool. Surrounding temperature gradients are calculated with the aide of a closest point projection algorithm. Furthermore, the model also allows simulation of pulsed power welding through time-dependent variation of the weld pool size. Example simulation results and comparisons with laboratory weld observations demonstrate microstructural variation with weld speed, pool shape, and pulsed-power.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Matching time and spatial scales of rapid solidification: Dynamic TEM experiments coupled to CALPHAD-informed phase-field simulations

A combination of dynamic transmission electron microscopy (DTEM) experiments and CALPHAD-informed phase-field simulations was used to study rapid solidification in Cu–Ni thin-film alloys. Experiments—conducted in the DTEM—consisted of in situ laser melting and determination of the solidification kinetics by monitoring the solid–liquid interface and the overall microstructure evolution (time-resolved measurements) during the solidification process. Modelling of the Cu–Ni alloy microstructure evolution was based on a phase-field model that included realistic Gibbs energies and diffusion coefficients from the CALPHAD framework (thermodynamic and mobility databases). DTEM and post mortem experiments highlighted the formation of microsegregation-free columnar grains with interface velocities varying from ~0.1 to ~0.6 m s –1 . After an 'incubation' time, the velocity of the planar solid–liquid interface accelerated until solidification was complete. In addition, a decrease of the temperature gradient induced a decrease in the interface velocity. The modelling strategy permitted the simulation (in 1D and 2D) of the solidification process from the initially diffusion-controlled to the nearly partitionless regimes. Lastly, results of DTEM experiments and phase-field simulations (grain morphology, solute distribution, and solid–liquid interface velocity) were consistent at similar time (μs) and spatial scales (μm).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Solidification and grain formation in alloys: a 2D application of the grand-potential-based phase-field approach

Solidification is a significant step in the forming of crystalline structures during various manufacturing and material processing techniques. Solidification characteristics and the microstructures formed during the process dictate the properties and performance of the materials. Hence, understanding how the process conditions relate to various microstructure formations is paramount. In this work, a grand-potential-based multi-phase, multi-component, multi-order-parameter phase-field model is used to demonstrate the solidification of alloys in 2D. This model has several key advantages over other multi-phase models such as it decouples the bulk energy from the interfacial energy, removes the constraints for the phase concentration variable, and prevents spurious third-phase formation at the two phase interfaces. Here, the model is implemented in a finite-element-based phase-field modeling code. The role of various modeling parameters in governing the solidification rate and the shape of the solidified structure is evaluated. It is demonstrated that the process conditions such as temperature gradient, thermal diffusion, cooling rate, etc, influence the solidification characteristics by altering the level of undercooling. Furthermore, the capability of the model to capture directional solidification and polycrystalline structure formation exhibiting various grain shapes is illustrated. In both these cases, the process conditions have been related to the growth rate and associated shape of the dendritic structure. As a result, this work serves as a stepping stone towards resolving the larger problem of understanding the process–structure–property–performance correlation in solidified materials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Implications of advanced collision operators for gyrokinetic simulation

In this paper, we explore both the potential improvements and pitfalls that arise when using advanced collision models in gyrokinetic simulations of plasma microinstabilities. Comparisons are made between the simple-but-standard electron Lorentz operator and specific variations of the advanced Sugama operator. The Sugama operator describes multi-species collisions including energy diffusion, momentum and energy conservation terms, and is valid for arbitrary wavelength. We report scans over collision frequency for both low and high ${k}_{\theta }{\rho }_{s}$ modes, with relevance for multiscale simulations that couple ion and electron scale physics. The influence of the ion–ion collision terms—not retained in the electron Lorentz model—on the damping of zonal flows is also explored. Collision frequency scans for linear and nonlinear simulations of ion-temperature-gradient instabilities including impurity ions are presented. Finally, implications for modeling turbulence in the highly collisional edge are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Fusion applications for lithium: wall conditioning in magnetic confinement devices

The application of lithium to plasma facing components (PFCs) has been widely used to improve plasma performance in magnetic confinement devices. The ability of lithium to retain hydrogen isotopes has been exploited to lower wall recycling. Low recycling leads to high edge plasma temperatures, and decreased temperature gradients reduce the free energy source for microinstabilities that cause transport. A variety of techniques have been developed to introduce lithium into fusion research devices, and its beneficial effects have been observed in their plasmas. The mechanism for the retention of hydrogen isotopes has been elucidated by detailed surface science studies, and 'first principles' molecular dynamics simulations. The remaining challenges primarily pertain to long pulse operation. Our knowledge of erosion, transport, and redeposition needs to improve for the mixed material PFCs lithium creates. As a result, flowing liquid lithium systems also must be demonstrated for the replenishable lithium walls that fusion reactors require.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Validation of gyrokinetic simulations of a National Spherical Torus eXperiment H-mode plasma and comparisons with a high-k scattering synthetic diagnostic

A new extensive validation study performed for a modest-beta National Spherical Torus eXperiment (NSTX) neutral beam injection-heated H-mode discharge predicts that electron thermal transport can be entirely explained by short-wavelength electron-scale turbulence fluctuations driven by the electron temperature gradient mode (ETG), both in conditions of strong and weak ETG turbulence drive. Quantitative comparisons between high-k fluctuation measurements and simulations are performed via a novel synthetic high-k diagnostic. Computationally intensive electron-scale simulations featuring an unusually large domain (L r , L θ ) ~ (20, 20)ρ s are shown to be required for accurate deployment of the synthetic diagnostic. Ion thermal transport is shown to be close to neoclassical levels, consistent with stable ion-scale turbulence simulations conducted with the GYRO code. Electron-scale GYRO simulations are shown to match the thermal power-balance estimates from TRANSP. The frequency spectra characteristics of electron-scale turbulence (spectral peak and width) can be consistently reproduced by the synthetic spectra, but these reveal not to be a critical constraint on the simulation model. The shape of the high-k wavenumber spectrum and the fluctuation level ratio between the strong and weak ETG conditions can also be simultaneously matched by electron-scale simulations within sensitivity scans about the experimental profile values, and result to be great discriminators of the turbulence models analyzed. Together, electron thermal power comparisons and quantitative agreement of electron-scale turbulence spectra give strong evidence supporting electron-scale ETG fluctuations as the main mechanism driving anomalous electron thermal transport in the two outer-core conditions of the modest-beta NSTX H-mode analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Recent progress on neoclassical impurity transport in stellarators with implications for a stellarator reactor

Accumulation of impurities in the core of the plasma is a potential problem for any fusion reactor, and is thought to be especially problematic for stellarators. We use the collisional transport code SFINCS to calculate, in a reactor stellarator scenario, the effects of flux-surface electrostatic potential variation on tungsten transport. The potential variation is found to have a large effect on the transport, reversing the sign of the transport coefficients related to the radial electric field and the bulk ion temperature gradient, and also increasing the magnitude of the former by two orders of magnitude. We optimize the electrostatic potential variation to minimize impurity accumulation, resulting in a 20%–50% reduction of the inward impurity flux. Furthermore, we investigate whether the electrostatic potential variation caused by fast particles due to ion-cyclotron resonance heating can be used to reduce the impurity accumulation, but find the effect to be small in the cases considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Zonal shear layer collapse and the power scaling of the density limit: old L-H wine in new bottles

Edge shear layer collapse causes edge cooling and aggravates radiative effects. This paper details on the microscopic dynamics of the emergence of power (Q) scaling of density limit (DL) from the shear layer collapse transport bifurcation scenario. The analysis is based on a novel 4-field model, which evolves turbulence energy, zonal flow energy, temperature gradient and density, including the neoclassical screening of zonal flow response. Bifurcation analysis yields power scaling of critical density for shear layer collapse as $n_{crit}\sim Q^{1/3}$. The favorable Q scaling of the DL emerges from the fact that the shear layer strength increases with Q, thus preventing shear layer collapse. This in turn reduces particle transport and improves particle confinement. RMP induced ambient stochastic fields degrade the shear layer by inducing decoherence in the Reynolds stress. As a result the particle transport increases and particle confinement degrades. This leads to the emergence of unfavorable stochastic field intensity ($b_{st}^{2}$) scaling of the critical density as $n_{crit}\sim(1+b_{st}^{2})^{-5/3}$. All fields, including zonal flow shear, exhibit hysteresis when the power (Q) is ramped cyclically across the bifurcation point. The hysteresis is due to dynamical delay in bifurcation on account of critical slowing down. Furthermore, the dynamical hysteresis here is fundamentally different from the hysteresis associated with the existence of bi-stable states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Gyrokinetic simulations of core turbulence and thermal transport in the high-β P discharge on EAST

The properties of core turbulence and thermal transport are investigated for EAST high-β P (β P ~ 3.1) plasmas with dominant electron heating (Te/Ti > 1) via gyrokinetic simulation with the NLT code. Linear simulations identify that the electrostatic η e -driven trapped electron mode (η e -TEM) dominates in the core region (ρ< 0.7) and the ion temperature gradient (ITG) mode dominates in the out region (ρ≥ 0.7), consistent with the linear threshold analysis of micro-instabilities. Sensitivity analysis shows that the normalized electron density gradient (R/L ne ) and ITG (R/L Ti ) are two effective parameters to stabilize TEM instability. Nonlinear simulations are also carried out and compared with the experimental results, which show that the electron thermal internal transport barrier (ITB) in EAST high-β P plasma is determined by the TEM-induced turbulence. A higher zonal flow shearing rate is observed in the ITB region (0.2 < ρ < 0.34), which can regulate energy transport induced by TEM turbulence and facilitate the formation of e-ITB. A plausible positive feedback mechanism to mitigate turbulence transport and improve energy confinement via enhanced ion heating is proposed for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of hydrogen isotope species on ITG microturbulence in LHD

Abstract The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.

Physics↗
DOE OSTI2024

Gyrokinetic simulation of pedestal degradation correlated with enhanced magnetic turbulence in a DIII-D ELMy H-mode discharge

Gyrokinetic simulation of a dedicated pedestal density ramping-up discharge on DIII-D can reproduce the enhancement of magnetic turbulence in the pedestal, which is identified to be caused by micro-tearing modes (MTMs). Here, an increase of MTM amplitude results in higher electron thermal diffusivity, consistent with experimentally observed lower electron temperature gradient and degraded pedestal height. Gyrokinetic simulation identifies the major cause of MTM enhancement to be the increase of collisionality, which has a significant impact on the MTM intensity and is beyond the description of any (quasi-)linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electron ITB formation in EAST high poloidal beta plasmas under dominant electron heating

Abstract Plasma confinement and transport in tokamaks play a crucial role in the development of high poloidal beta steady-state operation scenarios. Therefore, it is very important to study the relevant mechanism of internal transport barriers (ITBs), which can help plasmas to obtain better confinement in order to achieve higher fusion gain. This paper mainly introduces the analysis of the characteristics of electron heat transport of discharges with ITB in high β P operation regime in Experimental Advanced Superconducting Tokamak (EAST). Based on the statistical analysis of stable discharges with β P > 1.5, it is found that there is an obvious bifurcation of the normalised electron temperature gradient (ETG) ( R / L T e ) in the range of β P = 2–2.2. Then the discharges of lower β P ( β P < 2, where the value of β P is below the bifurcation threshold) and of higher β P ( β P > 2.2, where the value of β P is above the threshold) were selected for analysis. The diagnostic data provided by Thomson scattering, x-ray crystal spectrometry and charge exchange recombination spectroscopy are used to provide reliable parameter profiles and then combined with the data of external magnetic probe measurements and the polarisation interferometer diagnosis system to fully reconstruct the balance. A relevant plasma current calculation model is used to calculate and analyse the current density profiles and power deposition, and then the transport analysis is carried out. Interestingly, in the higher β P discharges, it is found that the turbulence intensity provided by the CO 2 laser collective scattering system gradually decreases and the normalised ETG R / L T e gradually increases. It is also found that the electron heat transport coefficient decreases in the discharge with higher β P and the growth rate of the electron-scale turbulence calculated by transport gyro-Landau fluid (TGLF) is significantly reduced. Meanwhile, similar conclusions are also obtained in the discharges when the β P is further increased.

Wang, Z. H. (ORCID:0009000203431036)↗
DOE OSTI2024

Examining transport and integrated modeling predictive capabilities for negative-triangularity scenarios

This paper investigates the predictive capabilities of TGYRO and TGLF models in assessing the performance of negative triangularity (NT) plasmas compared to positive triangularity (PT) plasmas in fusion devices. TGYRO predicts kinetic profiles, while TGLF analyzes turbulent transport. The study reveals that TGYRO reasonably predicts NT profiles similar to PT, although it overpredicts the high-power scenarios where there is increased experimental MHD activity. TGLF analysis finds reduced linear growth rates in NT and altered flux spectra relative to PT. Additionally, the TGLF SAT0 saturation model is observed to predict high-k transport and a reduction of particle transport with the electron temperature gradient. These findings are further corroborated by core-pedestal modeling using the Stability Transport Equilibrium Pedestal workflow, showing stronger confinement improvements in NT, particularly at higher power densities for the SAT0 saturation model. Furthermore, the study underscores the importance of accurately capturing turbulence saturation mechanisms for NT in order to project its performance accurately in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental characterization of turbulence properties in negative triangularity DIII-D plasmas

Negative Triangularity (NT) plasmas have demonstrated robustly ELM-free high-performance operation and provide a unique testbed to study how plasma shaping affects turbulent transport. The Beam Emission Spectroscopy (BES) diagnostic provides localized 2D measurements of low-k ($k_{\theta} \rho_s < 1$) density fluctuations. In a sweep of upper triangularity at fixed power, H-mode access is suppressed and an NT-edge is observed. The turbulence amplitude ($\tilde{n}/n$) is shown to decrease by $\sim50\%$ for $\rho<0.9$ in the NT-edge phase as compared to the H-mode phase. Additionally, low-velocity edge modes below 70 kHz are suppressed by triangularity and the dominant mode propagating in the electron diamagnetic direction is seen to broaden in wavenumber space in the NT-edge phase. In a strong NT plasma ($\delta_{avg}\sim-0.5$), low amplitude modes ($\tilde{n}/n <0.5\%$) propagating in the ion-diamagnetic direction with radially-poloidally symmetric eddy structure are identified for $\rho \sim 0.65-0.83$ consistent with Ion Temperature Gradient (ITG) turbulence. Modes consistent with Trapped Electron Mode (TEM) turbulence are observed propagating in the electron-diamagnetic direction for $\rho\gtrsim0.83$ with poloidally extended eddy structure and reduced amplitudes observed at the separatrix. The turbulence properties presented in this paper help validate our understanding of NT turbulence and help explain the improved confinement and unique edge features of NT plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Operation above the Greenwald density limit in high performance DIII-D negative triangularity discharges

The density limit in strongly-shaped negative triangularity (NT) discharges is studied experimentally in the DIII-D tokamak. Record-high Greenwald fractions f G are obtained, using gas puff injection only, with values up to near 2, where f G is defined as the ratio of the line-averaged density over n G = I p /(π a 2 ), with I p [MA] the plasma current and a[m] the plasma minor radius. A clear higher operational limit with higher auxiliary power is also demonstrated, with the ohmic density limit about two times lower than with additional neutral beam injection heating. The evolution of the electron density, temperature and pressure profiles are analyzed as well. The core density can be up to twice the Greenwald density and keeps increasing, while the value at the separatrix remains essentially constant and slightly below n G . The edge temperature gradient collapses to near zero and NT plasmas are shown to be resilient to such profiles in terms of disruptivity. We also present the time evolution of the inverse electron pressure scale length with the value at the last closed flux surface (LCFS) decreasing below the value at the normalized radius 0.9 near the density limit, demonstrating the clear drop of confinement starting from the edge. This inverse scale length ‘collapse’ at the LCFS also defines well the characteristic behavior of the kinetic profiles approaching a density limit.

density limit↗
DOE OSTI2026

Theory of zonal flow growth and propagation in toroidal geometry

The toroidal geometry of tokamaks and stellarators is known to play a crucial role in the linear physics of zonal flows (ZFs), leading to e.g. the Rosenbluth–Hinton residual and geodesic acoustic modes. However, descriptions of the nonlinear ZF growth from a turbulent background typically resort to simplified models of the geometry. We present a generalised theory of the secondary instability to model the ZF growth from turbulent fluctuations in toroidal geometry, demonstrating that the radial magnetic drift substantially affects the nonlinear ZF dynamics. In particular, the toroidicity gives rise to a new branch of propagating ZFs, the toroidal secondary mode, which is nonlinearly supported by the turbulence. We present a theory of this mode and compare the theory against gyrokinetic simulations of the secondary mode. The connection with other secondary modes—the ion-temperature-gradient and Rogers–Dorland–Kotschenreuther secondary modes—is also examined.

Nies, Richard [Princeton Univ., NJ (United States)↗
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗