Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermal diffusivity”

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 523 records · Page 29

A Hotspot’s Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics

Shockwave interactions with a material’s microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion. Furthermore, the origin of the PE hotspot is traced to large intramolecular strains, storing energy in modes readily available for chemical decomposition.

36 MATERIALS SCIENCE↗

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗

Dynamic substrate reactions during room temperature heavy ion irradiation of CoCrCuFeNi high entropy alloy thin films

Abstract High entropy alloys (HEAs) are promising materials for various applications including nuclear reactor environments. Thus, understanding their behavior under irradiation and exposure to different environments is important. Here, two sets of near-equiatomic CoCrCuFeNi thin films grown on either SiO 2 /Si or Si substrates were irradiated at room temperature with 11.5 MeV Au ions, providing similar behavior to exposure to inert versus corrosion environments. The film grown on SiO 2 had relatively minimal change up to peak damage levels above 500 dpa, while the film grown on Si began intermixing at the substrate–film interface at peak doses of 0.1 dpa before transforming into a multi-silicide film at higher doses, all at room temperature with minimal thermal diffusion. The primary mechanism is radiation-enhanced diffusion via the inverse Kirkendall and solute drag effects. The results highlight how composition and environmental exposure affect the stability of HEAs under radiation and give insights into controlling these behaviors.

36 MATERIALS SCIENCE↗

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Internal gravity waves in massive stars: II. Frequency analysis across stellar mass

Stars that are over 1.6 solar masses are generally known to possess convective cores and radiative envelopes, which allows for the propagation of outwardly travelling internal gravity waves (IGWs). Here, we study the generation and propagation of IGWs in such stars using two-dimensional, fully non-linear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics. Compared to previous similar works, this study utilises radius-dependent thermal diffusivity profiles for five different stellar masses at the middle of the main sequence: 3 M ⊙ , 5 M ⊙ , 7 M ⊙ , 10 M ⊙ , and 13 M ⊙ . From the simulations, we find that the surface perturbations are larger for higher masses, but no noticeable trends are observed for the frequency slopes with different stellar masses. The slopes are also similar to the results from previous works. We compared our simulation results with stellar photometric data from a recent survey and we found that for frequency intervals above 8 μHz, there is a good agreement between the temperature frequency slopes from the simulations and the surface brightness variations of these observed stars. This indicates that the brightness variations are caused by core-generated IGWs.

79 ASTRONOMY AND ASTROPHYSICS↗

Extending the validation of multi-mode model for anomalous transport to high beta poloidal tokamak scenario in DIII-D

The Multi-Mode Model (MMM7.1) for anomalous transport is tested in predictive modeling of temperature profiles of a high beta poloidal DIII-D discharge. This new H-mode plasma regime, with high beta poloidal and high bootstrap currents, has been studied in DIII-D tokamak discharges. The role of instabilities that can drive the anomalous transport described by MMM7.1 is investigated. The temperature profiles for a high beta poloidal DIII-D discharge are computed using the NCLASS model for the neoclassical transport and the Weiland and Electron Temperature Gradient (ETG) components of the MMM7.1 model for the anomalous transport. Here, the neoclassical transport is found to be the main contributor to the ion thermal transport in the plasma core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of thermal annealing on the redistribution of alkali metals in Cu(In,Ga)Se 2 solar cells on glass substrate

The precise control of alkali-metal concentrations in Cu(In,Ga)Se 2 (CIGS) solar cells via post deposition treatment (PDT) has recently attracted attention. When PDT is performed at an elevated temperature, an accompanying annealing effect is expected. Here, we investigate how thermal annealing affects the redistribution of alkali metals in CIGS solar cells on glass substrates and the properties of the solar cells. In addition, we investigate the origin of non-homogeneous alkali-metal depth profiles that are typical of CIGS grown using a three-stage process. In particular, we use secondary-ion mass spectrometry measurements of the ion concentration as a function of distance from the CIGS surface to investigate the impact of thermal annealing on the distribution of alkali metals (Na, Ka, and Rb) and constituent elements (Ga and In) in the CIGS absorbers. We find that the depth profiles of the alkali metals strongly reflect the density of sites that tend to accommodate alkali metals, i.e., vacancies. Annealing at elevated temperature caused a redistribution of the alkali metals. The thermal-diffusion kinetics of alkali metals depends strongly on the species involved. We introduced low flux potassium fluoride (KF) to study a side effect of KF-PDT, i.e., Na removal from CIGS, separately from its predominant effects such as surface modification. When sufficient amounts of Na are supplied from the soda lime glass via annealing at an elevated temperature, the negative effect was not apparent. Conversely, when the Na supply was not sufficient, it caused a deterioration of the photovoltaic properties.

14 SOLAR ENERGY↗

Narrow band photoacoustic lamb wave generation for nondestructive testing using candle soot nanoparticle patches

The generation of ultrasonic surface waves with a photoacoustic-laser-source has become useful for the noncontact nondestructive testing and evaluation (NDT&E) of materials and structures. In this work, a hybrid ultrasound based NDT&E method is proposed based on the photoacoustic-laser-source as a noncontact Lamb wave generator by incorporating a line-arrayed patterned candle soot nanoparticle-polydimethylsiloxane (CSNPs-PDMS) patch as the signal amplifier and with a narrow bandwidth. The CSNP-PDMS composite has been investigated as the functional patch for its laser energy absorption efficiency, fast thermal diffusion, and large thermoelastic expansion capabilities. The signal amplitude (in mW) from the CSNP-PDMS patch exhibits 2.3 times higher amplitude than the no patch condition and a narrower bandwidth than other conditions. Furthermore, improvement in the sensitivity is also achieved through the selection of the aluminum nitride sensing system. The overall combination of the Lamb wave generation and receiver-sensing system in this study is found to be very promising for a broad range of noncontact NDT&E applications.

Acoustic transducers↗

Calculating RF current condensation with consistent ray-tracing and island heating

Finally, by exploiting the nonlinear amplification of the power deposition of RF waves, current condensation promises new pathways to the stabilization of magnetic islands. We present a numerical analysis of current condensation, coupling a geometrical optics treatment of wave propagation and damping to a thermal diffusion equation solver in the island. Taking into account the island geometry and relativistic damping, previous analytical theory can be made more precise and specific scenarios can be realistically predicted. With this more precise description, bifurcations and associated hysteresis effects could be obtained in an ITER-like scenario at realistic parameter values. Moreover, it is shown that dynamically varying the RF wave launching angles can lead to hysteresis and help to avoid the nonlinear shadowing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Linear ion-scale microstability analysis of high and low-collisionality NSTX discharges and NSTX-U projections

Linear gyrokinetic simulations were conducted to investigate ion-gyroradius-scale micro-instability predictions for high-beta NSTX discharges and NSTX-U projections that span over an order of magnitude variation in collisionality. A complex mix of microtearing modes and hybrid trapped electron modes/kinetic ballooning modes (TEM/KBM) is predicted for all experimental or projected conditions. Ion temperature gradient (ITG) instabilities are typically stable in the NSTX discharges investigated, consistent with the observed neoclassical ion thermal transport. ITG thresholds inferred from the simulations are typically much higher than the experimental NSTX gradients, as well as the projected gradients in the NSTX-U scenario, which assumed ion temperatures limited by neoclassical transport only. The analysis suggests ITG instabilities are unlikely to contribute significant anomalous thermal losses in high-beta, lower collisionality NSTX-U scenarios. On the other hand, the NSTX experimental profiles and NSTX-U projections are predicted to be very close to the predicted onset of unstable KBM at most radii investigated. The proximity of the various discharges to the KBM instability threshold implies it may play an important role in setting profile shapes and limiting global energy confinement. Here, it remains to be understood and predicted how KBM contributes to multi-channel transport (thermal and particle transport, for both ions and electrons) in a way that is consistent with experimental inferences.

36 MATERIALS SCIENCE↗

Role of edge-localized neoclassical tearing modes in quiescent H-mode plasmas in the DIII-D tokamak

The edge-harmonic oscillations (EHOs) in standard quiescent H-mode (QH-mode) plasmas in DIII-D are consistent with edge-localized neoclassical tearing modes (NTMs) based on nonlinear two-fluid MHD simulations. Using kinetic equilibria constrained by edge profile measurements, the MHD simulations show that the n = 1 NTM and its harmonics can be destabilized at the pedestal top of QH-mode plasma by the edge bootstrap current. The simulations further show that the unstable NTMs can saturate either at small (<2% ψ N ) or large (>4% ψ N ) island width depending on the magnitude of the edge bootstrap current, where ψ N is the normalized radius in poloidal flux. The onset of the EHO also results in a prompt decrease in the pedestal width and height, consistent with simulation results for the onset of the NTM at the top of the QH-mode pedestal. This suggests that the avoidance of edge-localized modes (ELMs) in QH-mode can be attributed to the enhanced local transport induced by the NTM that is sufficient to prevent the expansion of the pedestal to an unstable width, analogous to the mechanism explored for ELM suppression by resonant magnetic perturbations. Nonlinear MHD simulations scanning the E × B frequency and the ratio of parallel and perpendicular thermal diffusivity (χ | /χ ⊥ ) at the pedestal top show that edge-localized NTMs are destabilized for conditions of high E × B frequency, high pedestal temperature, and low pedestal density, qualitatively consistent with experimental conditions required for observing the EHO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Use of intrinsic hysteresis for the active control of internal transport barriers in magnetically confined fusion plasmas

In magnetically confined fusion devices, control of internal transport barriers (ITBs) is important both to enhance and suppress the turbulent transport to improve confinement control. Barrier control should allow for the improvement of confinement to aid in achieving the needed fusion criteria while also permitting the degradation of confinement to control profiles and clean the device by moving out the impurities accumulated near the core. In this work, we present a novel control scenario that takes advantage of the hysteresis intrinsic to transport barriers to easily cycle through enhanced and degraded confinement regimes. The control scenario is illustrated using a five-field simplified transport model for an ITB using typical parameters of a neutral beam injection-heated DIII-D tokamak discharge. Pellets and ion cyclotron resonance frequency power are used as control knobs for this active control scenario. These knobs adequately modify at will the local gradients and, therefore, the growth rates and shearing rates, allowing for easy and efficient control of the barrier by taking advantage of the barrier hysteresis. The result is a control cycle that could be operated with a relatively small amount of power in high performance regimes which, nowadays, typically require large power to control. It may also have advantages to avoid, or at least ameliorate, the appearance of magnetohydrodynamic instabilities in the barrier region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

A Diethylene Glycol Condensation Particle Counter for Rapid Sizing of sub-3 nm Atmospheric Clusters

AbstractThis paper describes the modification of a laminar flow, thermally diffusive universal-fluid condensation particle counter (standard operation: 50% detection efficiency at 5 nm) to rapidly measure the size distribution of sub 3 nm aerosol. Sub 3 nm detection was achieved by using diethylene glycol as the working fluid, which enabled high instrument super-saturations while minimizing homogenous nucleation of the working fluid; a detection efficiency of 50% was achieved at 1.6 nm with laboratory-generated ammonium sulfate aerosol. Rapid aerosol sizing beneath 3 nm was achieved by inverting the measured grown droplet size distribution (1 second sampling) to recover the sampled aerosol size distribution. The developed inversion algorithm utilizes analytical kernel functions determined from the instrument response to pseudo-monodisperse ammonium sulfate aerosol from 1.5 nm to 20 nm, generated by a high resolution DMA and a nano DMA. The inversion algorithm was tested numerically with assumed, idealized aerosol size distributions consistent with observed new particle formation events, yielding reasonable agreement between inverted and assumed aerosol size distributions below 3 nm. Furthermore, this technique provides a measure of the aerosol size assuming an aerosol composition identical to that of the aerosol used to generate the experimentally determined kernel function.

54 ENVIRONMENTAL SCIENCES↗

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↗