Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Simultaneous measurement of specific heat and thermal conductivity in pulsed magnetic fields

Here we report an experimental setup for simultaneously measuring specific heat and thermal conductivity in feedback-controlled pulsed magnetic fields of 50 ms duration at cryogenic temperatures. A stabilized magnetic field pulse obtained by the feedback control, which dramatically improves the thermal stability of the setup and sample, is used in combination with the flash method to obtain absolute values of thermal properties up to 37.2 T in the 22–16 K temperature range. We describe the experimental setup and demonstrate the performance of the present method with measurements on single-crystal samples of the geometrically frustrated quantum spin-dimer system SrCu 2 (BO 3 ) 2 . Our proof-of-principle results show excellent agreement with data taken using a standard steady-state method, confirming the validity and convenience of the present approach.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Field reversal in low pressure, unmagnetized radio frequency capacitively coupled argon plasma discharges

In general, the radio frequency (rf) electric field within a sheath points toward the metal electrode in low pressure, unmagnetized rf electropositive capacitively coupled plasma (CCP) glow discharges. This is due to the large ratio of electron to ion mobility and the formation of an ion sheath. In this work, we studied, using fully kinetic particle-in-cell simulations, a reversed electric field induced by the strong secondary electron emission during the phase of sheath collapse in a high-voltage rf-driven low pressure CCP glow discharge. We explored the transition behavior of the formation of field reversal as a function of driving voltage amplitude and found that field reversal starts to form at around 750 V, for a discharge with an electrode spacing of 4 cm at 10 mTorr argon pressure driven at 13.56 MHz. Accordingly, the energy distribution function of electrons incident on the electrode shows peaks from around 3 to 10 eV while varying the driving voltage from 150 to 2000 V, showing potentially beneficial effects in plasma material processing where relatively directional electrons are preferred to solely thermal diffusion electrons.

Physics↗
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 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 OSTI2026

Statistical inference of anomalous thermal transport with uncertainty quantification for interpretive 2D SOL models

The critical task of inferring anomalous cross-field transport coefficients is addressed in simulations of boundary plasmas with fluid models. A workflow for parameter inference in the UEDGE fluid code is developed using Bayesian optimization with parallelized sampling and integrated uncertainty quantification. In this workflow, transport coefficients are inferred by maximizing their posterior probability distribution, which is generally multidimensional and non-Gaussian. Uncertainty quantification is integrated throughout the optimization within the Bayesian framework that combines diagnostic uncertainties and model limitations. As a concrete example, we infer the anomalous electron thermal diffusivity $\chi_\perp$ from an interpretive 2D model describing electron heat transport in the conduction-limited region with radiative power loss. The workflow is first benchmarked against synthetic data and then tested on H-, L-, and I-mode discharges to match their midplane temperature and divertor heat flux profiles. We demonstrate that the workflow efficiently infers diffusivity and its associated uncertainty, generating 2D profiles that match 1D measurements. Future efforts will focus on incorporating more complicated fluid models and analyzing transport coefficients inferred from a large database of experimental results.

Bayesian optimization↗
DOE OSTI2026

Routes to high-performance operation in Wendelstein 7-X: turbulence suppression with shaping of the density profile

Steep density gradients generally lead to improved plasma performance in the neoclassically optimized stellarator W7-X. This is evident in the global energy confinement time as well as in the ion temperature and can be explained by a strong reduction of the ion temperature gradient turbulence. Such conditions can experimentally be realized by several methods: injection of cryogenic hydrogen pellets, appropriate combination of neutral beam and electron cyclotron resonance heating (ECRH) and, in some cases, with low power of ECRH after preconditioning of the first wall. The duration of the improved phases is determined by the ability to sustain the steep density gradient, by technical limitations of the involved systems and, eventually, by the plasma stability. This paper gives an overview of relevant experimental results and presents example discharges where the improved confinement conditions could be extended to multiple seconds: up to 4 s using neutral beam injection and from 14 to 40 s with steady state pellet injection. In these plasmas the turbulent thermal diffusivity is reduced by a factor of 3 to 4 in a broad radial range, which allows high ion temperatures of up to 3 keV at the densities of about 1.5 • 10 20 m −3 .

high performance plasma↗
DOE OSTI2020

Magnetic shear effect on plasma transport at T e / T i ~ 1 through electron cyclotron heating in DIII-D plasmas

The effect of magnetic shear on plasma transport for an electron to ion temperature ratio ( T e / T i ) near unity has been explored in DIII-D utilizing electron cyclotron heating (ECH). Previous reports showed that significant confinement degradation occurred at T e / T i ~ 1 in positive shear (PS) plasmas in DIII-D, whereas reduced confinement degradation was observed in negative central shear (NCS) plasmas. In this study, plasma transport in weak magnetic shear (WS) plasmas with ECH is investigated and compared with that in NCS and PS plasmas. Here the magnetic shears ($\hat s$) are $\hat s$ > 0.5, ~0 and <-0.1 in the core region (ρ~ 0.3–0.4) of PS, WS and NCS plasmas, respectively, and flat or negative inside ρ~ 0.4 in the WS and NCS plasmas. Weak magnetic shear is found to be effective in minimizing degradation of ion thermal confinement as T e / T i increases through ECH application, and an improved confinement factor of H 98y2 ~ 1.2 is maintained, similar to NCS plasmas. At T e / T i ~ 1, the ion thermal diffusivity around an internal transport barrier decreases when changing the magnetic shear from positive to weak or negative shear. Also, reduced local particle and momentum transport was indicated by steeper density and toroidal rotation profiles in the weak and negative shear regimes. Linear gyrokinetic simulations predict little change in growth rates of low-k turbulence with ECH application in the WS and NCS plasmas, which is consistent with the transport and profile analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗
DOE OSTI2024

Synergistic effects of W–Ta alloys exposed to various ELM-like fusion events

Abstract Recently tungsten alloys have shown capabilities in suppressing morphological changes when exposed to harsh fusion conditions. Morphological changes due to He + and D + ions at high temperatures (1000–2000 K) can lead to easier erosion and sputtering of tungten (W) particles towards quenching the fusion plasma. This study investigated W–3Ta alloy as potential plasma facing materials (PFMs) exposed to different edge-localized mode (ELM)-like energies and ion loadings. Performance of W–3Ta has been tested as a function of ELM-like energies (0.6–1.5 GW m −2 ) and varying He + and D + ion concentrations (100% He + , 50% He + :50% D + ,10% He + :90% D + (similar ELM-like reactor condition), and 100% D + ions). Early-stage fuzz formation and surface pores were observed during 100% He + ion irradiation and mixed ion loadings but not observed during 100% D + ion loadings. Pores decreased in quantity and increased in size with increasing transient heat loading. Grain boundaries acted as trapping sites of helium and deuterium atoms during dual beam ion irradiations, as observed through increased pore size along the grain boundary. Qualitative erosion estimates were measured using witness plates. Increased erosion was observed with increasing D + ion concentration, which is attributed to D supersaturation in the near surface region affecting thermal and mechanical properties. More erosion was observed for W–3Ta samples compared to pure W, evident through large (>10 µ m) cracks on the surface, and an increase in W metal x-ray photoelectron spectroscopy peaks under similar loading conditions. The slight enhancement in erosion is attributed to the slightly decreased thermal properties, ∼10% thermal diffusivity drop, of W–3Ta at high temperatures. The damage threshold was impacted by ion loading, where distinct cracking was observed for 100% He + at a lower heat loading. This matrix study using reactor-like loading supports the need for extensive studies on future W alloy PFMs before selecting one as an alternative.

Physics↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗
DOE OSTI2023

A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations

ABSTRACT We perform two-dimensional (2D) numerical simulations of core convection for zero-age main-sequence stars covering a mass range from 3 to 20 M⊙. The simulations are performed with the fully compressible time-implicit code music. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 M⊙ models. The simulations cover hundreds to thousands of convective turnover time-scales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance (dov) scales with the stellar luminosity (L) and the convective core radius (rconv). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$, which is implemented in a one-dimensional stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses ≳10 M⊙. Our 2D simulations show the formation of a nearly adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent three-dimensional simulations of convection. The most luminous models show a growth in size with time of the nearly adiabatic layer. This growth seems to slow down as the upper edge of the nearly adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive time-scale in the overshooting layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dynamic evolution of ultracold electron beams from nanostructured cathodes

Dynamic Coulomb expansion of dense, nonuniform particle beamlets can lead to causticlike formations in the charge density profile of the resultant beam. Such features only occur for cold beams where thermal diffusion of the particles is minimized. These phenomena have been observed experimentally for select cases with ion beams but not yet for electron beams as the electron dynamics occur over far faster timescales. We simulate the interaction of electron bunches while varying the initial transverse temperature and density profiles to determine the thresholds that characterize this pattern formation. Furthermore, we extend the simulations to geometric arrangements that are more pertinent to the nanotip cathode application and consider the effects of asymmetries and the impact of a low-density halo on the overall process. The simulations are conducted using a previously documented high-fidelity algorithm for collisional particle dynamics. These studies identify the temperature and density thresholds required for caustic formation and provide guidance for optimizing the design of nanoscale structured cathodes. We also demonstrate that the correlation between the final density patterns and the underlying beam parameters enables a virtual diagnostic tool for measuring ultralow emittance beams solely based on a transverse charge density profile.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Nuclear fission reaction simulations in compact stars

Type-Ia supernovas are powerful stellar explosions that provide important distance indicators in cosmology. Recently, we proposed a new Type-Ia supernova mechanism that involves a nuclear fission chain-reaction in an isolated white dwarf [Phys. Rev. Lett. 126, 131101 (2021)]. Here, in this work, we perform novel reaction network simulations of the actinide-rich first solids in a cooling white dwarf. The network includes neutron-capture and fission reactions on a range of U and Th isotopes with various possible values for 235 U enrichment. We find, for modest 235 U enrichments, neutron capture on 238 U and 232 Th can breed additional fissile nuclei so that a significant fraction of all U and Th nuclei may fission during the chain reaction. Finally, we compute the energy release from the fission chain reaction for various uranium enrichments; a novel result that is a necessary input for thermal diffusion simulations of carbon ignition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Convection-Induced Compositional Patterning at Grain Boundaries in Irradiated Alloys

We consider the stability of precipitates formed at grain boundaries (GBs) by radiation-induced segregation in dilute alloys subjected to irradiation. Here, the effects of grain size and misorientation of symmetric-tilt GBs are quantified using phase field modeling. A novel regime is identified where, at long times, GBs are decorated by precipitate patterns that resist coarsening. Maps of the chemical Péclet number indicate that arrested coarsening takes place when solute advection dominates over thermal diffusion right up to the precipitate/matrix interface, preventing interfacial local equilibrium and overriding capillary effects. This contrasts with liquid-solid mixtures where convection always accelerates coarsening.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Challenges to Reliable Production Nitrogen Doping of Nb for SRF Accelerating Cavities

Over the last several years, alloying of the surface layer of niobium SRF cavities has been demonstrated to beneficially lower the superconducting RF surface resistance. Nitrogen, titanium, and oxygen have all been demonstrated as effective alloying agents, occupying interstitial sites in the niobium lattice within the RF penetration depth and even deeper, when allowed to thermally diffuse into the surface at appropriate temperatures. The use of nitrogen for this function has been often termed ’nitrogen doping’ and is being applied in the LCLS-II and LCLS-II HE projects. We report characterization studies of the distribution of nitrogen into the exposed niobium surface and how such distribution is affected by the quality of the vacuum furnace environment in which the doping takes place, and the complexity of nitride crystal growth on different grain orientations of surface niobium. Using state-of-the-art quantification methods by dynamic secondary ion mass spectrometry (SIMS) depth profiling in niobium, we identify several phenomena involving furnace-sourced contamination. We also highlight a potential issue with N₂ flow constraints from the flange ’caps’ used during heat treatments.

Reece, Charles↗
DOE OSTI2023

SIMS CHARACTERIZATION OF NITROGEN DOPING OF LCLS-II-HE PRODUCTION CAVITIES

The thermal diffusion of nitrogen into the surface of niobium has been shown to yield superior low-loss SRF performance. An effective solution was identified and promptly employed in the production of cryomodules for LCLS-II. With added experience and R&D, a modified process was chosen for use in the upgrade for LCLS-II-HE. Largely motivated by this circumstance, supporting research has significantly refined the technique for making calibrated secondary ion mass spectrometry (SIMS) measurements of the N concentration depth profiles produced by production processes. Standardized reference samples were included with four HE production cavities in their N-doping furnace runs. We report the calibrated dynamic SIMS depth profiles of N, C, and O for these samples, together with the cryogenic acceptance test performance of the associated cavities. Interpretation and comparison with similar samples acquired in other furnaces highlights the importance of intentional process quality control of furnace conditions.

Reece, Charles↗