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

Microstructural and material property changes in severely deformed Eurofer-97

Severe plastic deformation changes the microstructure and properties of steels, which may be favourable for their use in structural components of nuclear reactors. In this study, high-pressure torsion (HPT) was used to refine the grain structure of Eurofer-97, a ferritic/martensitic steel. Electron microscopy and X-ray diffraction were used to characterise the microstructural changes. Following HPT at room temperature to a maximum shear strain of 230, the average grain size reduced by a factor of ~30, with a marked increase in high-angle grain boundaries. Dislocation density also increased by more than one order of magnitude. The thermal stability of the deformed material was investigated via in-situ annealing during synchrotron X-ray diffraction. This revealed substantial recovery between 450 K – 800 K. Irradiation with 20 MeV Fe-ions to ~0.1 dpa caused a 20% reduction in dislocation density compared to the as-deformed material. However, HPT deformation prior to irradiation only had a minor effect in mitigating the irradiation-induced reductions in thermal diffusivity and surface acoustic wave velocity of the material. Microstructural and material property changes are dominated by deformation compared to irradiation. In light of this, the benefits of using HPT to improve the irradiation resistance of Eurofer-97 are limited. These results provide a multi-faceted view of the changes in ferritic/martensitic steels due to severe plastic deformation, and how these changes can be used to alter material properties.

ion-irradiation↗
DOE OSTI2023

Face mask integrated with flexible and wearable manganite oxide respiration sensor

Face masks are key personal protective equipment for reducing exposure to viruses and other environmental hazards such as air pollution. Integrating flexible and wearable sensors into face masks can provide valuable insights into personal and public health. The advantages that a breath-monitoring face mask requires, including multi-functional sensing ability and continuous, long-term dynamic breathing process monitoring, have been underdeveloped to date. Here, we design an effective human breath monitoring face mask based on a flexible La 0.7 Sr 0.3 MnO 3 (LSMO)/Mica respiration sensor. The sensor’s capabilities and systematic measurements are investigated under two application scenes, namely clinical monitoring mode and daily monitoring mode, to monitor, recognise, and analyse different human breath status, i.e., cough, normal breath, and deep breath. This sensing system exhibits super-stability and multi-modal capabilities in continuous and long-time monitoring of the human breath. We determine that during monitoring human breath, thermal diffusion in LSMO is responsible for the change of resistance in flexible LSMO/Mica sensor. Both simulated and experimental results demonstrate good discernibility of the flexible LSMO/Mica sensor operating at different breath status. Our work opens a route for the design of novel flexible and wearable electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Avalanche gain and its effect on energy resolution in GEM-based detectors

Here, we present avalanche gain and associated resolution measurements recorded with a 4 He:CO 2 (70:30) gas mixture and pure SF 6 , a Negative Ion (NI) gas. SF 6 is of particular interest to the directional dark matter detection community, as its low thermal diffusion helps to retain recoil ionization track features over long drift lengths. With the aid of a general form of the reduced first Townsend coefficient (RFTC), multiple GEM-based detector data sets are used to study the high-gain behavior of the 4 He:CO 2 gas mixture. The high-gain data is well described purely in terms of the reduced electric field strength and the number of GEMs, and the robust relationship between the RFTC and the average, reduced, electric field strength across the GEMs is emphasized. The associated (pulse-height) resolution measurements are used to discuss the variance of the avalanche distribution and to describe and estimate the lower limits of energy resolution one should expect to measure using a simple relationship with the RFTC. In the end, a description of avalanche gain, its effect on energy resolution, and the contributing experimental parameters in GEM-based detectors is developed over a broad parameter space for further use.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Optimization of target lifetime for production of 14 MeV neutrons

Two methods are examined for extending the life of tritium targets for production of 14 MeV neutrons by the 3 H( 2 H,n) 4 He nuclear reaction. With thick film targets the neutron production rate decreases with time due to isotope exchange of tritium in the film with implanted deuterium. Here, the target life is maximized by operating the target at elevated temperature where the implanted deuterium mixes by thermal diffusion throughout the entire thickness of the film. The number of neutrons obtained from a target is then proportional to the initial tritium content of the film. A novel thin-film target design was also developed and tested. With these thin-film targets, the incident deuterium is implanted through the tritide into the underlying substrate material. A thin permeation barrier layer between the tritide film and substrate, reduces the rate of tritium loss from the tritide film. Good thin-film target performance was achieved using W and Fe for the barrier and substrate materials respectively. Thin-film targets were fabricated and tested and shown to produce similar number of neutrons as thick-film targets while using only a small fraction of the amount of tritium.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Pool boiling critical heat flux studies of accident tolerant fuel cladding materials

Pool boiling experiments have been performed at atmospheric pressure on potential accident tolerant fuel (ATF) cladding materials (Zirlo®, Cr coated Zirlo®, FeCrAl coated Zirlo®, and monolithic SiC) to evaluate the Critical Heat Flux (CHF) points. Corrosion resistant coatings of Cr and FeCrAl alloy were deposited on flat Zirlo® samples using cold spray technology. The as-prepared samples after uniform surface polishing were subjected to autoclave tests at 360 °C water and 18.6 MPa for 360 h to simulate prototypic corrosion of the ATF cladding materials in LWR normal operation. Surface characteristics potentially influencing CHF such as surface morphology, roughness, static contact angle, and surface chemistry were characterized using a suite of characterization methods including scanning electron microscopy, 3D optical profilometry, optical contact angle measurements, and x-ray photoelectron spectroscopy (XPS). Thermo-physical properties of the samples such as density, thermal conductivity, and heat capacity were measured by differential scanning calorimetry and laser flash thermal diffusivity measurement. FeCrAl coatings showed CHF values comparable to bare Zirlo® samples, while slightly lower CHF values were observed for Cr coated samples and bulk SiC flats (produced by chemical vapor deposition, CVD). As-deposited Cr coatings showed 67% higher CHF than the as-polished Cr coatings due to its higher surface roughness level. The autoclave testing for the prepared samples generally increased the CHF except for the surface polished Cr coatings. Similar CHF values of bare Zirlo®, FeCrAl coating, and CVD SiC samples were observed after autoclave test but there were negligible change CHF of Cr coatings. XPS studies indicated that the formation or deposition of a few monolayers of hydrophobic carbonaceous species on surface can affect CHF values. The trends in CHF data are discussed in terms of evolution of the surface characteristics of the ATF materials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Structure–property relations in graphitic pebbles for nuclear applications

This work presents an analytical approach for holistically characterizing graphitic matrix pebbles for nuclear applications whereby the macrostructure, microstructure, and thermophysical properties of pebbles are determined. A systematic sectioning method was applied to several pebbles to describe the regional properties of the samples. Intact matrix-only spheres and sections of spheres fabricated by Kairos Power were characterized via optical imaging, x-ray computed tomography, x-ray diffraction, and ellipsometry to determine 2D and 3D macrostructure and anisotropy. The thermophysical properties of these materials were determined via measurements of density, specific heat, thermal expansion, and thermal diffusivity. The results of this study indicate that the pebble fabrication methods and their resultant effect on microstructure have a nontrivial effect on thermophysical properties, confirming the importance of robust characterization of these components. A discussion of the characterization approach and its applicability to nuclear fuel development activities is also included.

42 ENGINEERING↗
DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

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

Mechanics↗
DOE OSTI2021

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↗
DOE OSTI2025

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↗
DOE OSTI2022

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↗
DOE OSTI2020

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↗
DOE OSTI2023

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↗
DOE OSTI2020

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↗
DOE OSTI2022

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↗
DOE OSTI2022

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↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2023

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↗