Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 793 records · Page 44

Densification, microstructure, and mechanical properties of Mo–30W alloys fabricated from conditioned powders

Refractory alloys, such as molybdenum-based systems, are attracting growing interest for applications in extreme environments, such as in the nuclear and aerospace industries. Recent advances in sintering technologies, coupled with mechanical alloying, have enabled the tailored design of these alloys by leveraging powder characteristics to control final microstructures and mechanical properties. In this study, Mo-30W alloys were fabricated using electric field-assisted sintering (EFAS) from ball-milled powders with and without hydrogen treatment to investigate the influence of surface oxides on material properties and sintering behavior. The results revealed that samples processed from as-ball-milled powder contained a high density of oxides within the microstructure, whereas oxide presence was significantly reduced in samples fabricated from hydrogen-treated powders. Interestingly, the two powder types led to opposite trends in grain size distribution: samples from untreated powders exhibited grain refinement from sample periphery to the center, while samples from hydrogen-treated powders showed grain coarsening toward the center. This behavior is attributed to temperature gradients present during sintering due to electrical percolation pathway differences during Joule heating. The powder surface oxides may have influenced the temperature distribution and grain evolution. Microhardness profiles measured along both axial and thickness directions were consistent with the grain size distribution. Furthermore, oxide films on powder surfaces have delayed densification by hindering particle necking and atomic diffusion during sintering.

36 - MATERIALS SCIENCE↗

ATMOSPHERIC HEAT REDISTRIBUTION ON HOT JUPITERS

Infrared light curves of transiting hot Jupiters present a trend in which the atmospheres of the hottest planets are less efficient at redistributing the stellar energy absorbed on their daysides—and thus have a larger day-night temperature contrast—than colder planets. To this day, no predictive atmospheric model has been published that identifies which dynamical mechanisms determine the atmospheric heat redistribution efficiency on tidally locked exoplanets. Here we present a shallow-water model of the atmospheric dynamics on synchronously rotating planets that explains why heat redistribution efficiency drops as stellar insolation rises. Our model shows that planets with weak friction and weak irradiation exhibit a banded zonal flow with minimal day-night temperature differences, while models with strong irradiation and/or strong friction exhibit a day-night flow pattern with order-unity fractional day-night temperature differences. To interpret the model, we develop a scaling theory which shows that the timescale for gravity waves to propagate horizontally over planetary scales, τ{sub wave}, plays a dominant role in controlling the transition from small to large temperature contrasts. This implies that heat redistribution is governed by a wave-like process, similar to the one responsible for the weak temperature gradients in the Earth's tropics. When atmospheric drag can be neglected, the transition from small to large day-night temperature contrasts occurs when τ{sub wave}∼√(τ{sub rad}/Ω), where τ{sub rad} is the radiative relaxation time and Ω is the planetary rotation frequency. Alternatively, this transition criterion can be expressed as τ{sub rad} ∼ τ{sub vert}, where τ{sub vert} is the timescale for a fluid parcel to move vertically over the difference in day-night thickness. These results subsume the more widely used timescale comparison for estimating heat redistribution efficiency between τ{sub rad} and the horizontal day-night advection timescale, τ{sub adv}. Only because τ{sub adv} ∼ τ{sub vert} for hot Jupiters does the commonly assumed timescale comparison between τ{sub rad} and τ{sub adv} yield approximately correct predictions for the heat redistribution efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗

Broad band X-ray Telescope observations of the hot interstellar media in NGC 1399 and NGC 4472

We present our analysis and interpretation of spatially resolved X-ray spectroscopy of the elliptical galaxies NGC 1399 and NGC 4472 obtained with the Broad Band X-Ray Telescope. The X-ray emission from both galaxies is dominated by thermal emission from a hot interstellar medium. The temperature of the plasma in NGC 1399 is 1.0-1.2 keV with a mild positive temperature gradient; for NGC 4472 the average temperature is 0.7-1.0 keV. The temperature of NGC 1399, when combined with other X-ray and optical data, implies that, if the gas is gravitationally confined, about 70 percent of the mass inside a radius of about 8.6 arcmin (38 kpc) is nonluminous. The derived metallicities and metallicity gradients are consistent with optically determined stellar metallicities, and allow surprisingly small upper limits to be placed on Type Ia supernovae rates in both galaxies. The implications of this for the gasdynamical evolution of the two galaxies are discussed. Limits on the oxygen-to-iron ratio and the line-of-sight column density are derived. There is a possible detection of cold gas in the line of sight to NGC 1399 that exceeds the Galactic value. Constraints on the contribution from a harder spectral component are placed, and likely candidates for the origin of this secondary component are evaluated.

Serlemitsos, P. J.↗

ENSO Dynamics in the E3SM-1-0, CESM2, and GFDL-CM4 Climate Models

We report that this study examines historical simulations of ENSO in the E3SM-1-0, CESM2, and GFDL-CM4 climate models, provided by three leading U.S. modeling centers as part of the Coupled Model Intercomparison Project phase 6 (CMIP6). These new models have made substantial progress in simulating ENSO’s key features, including: amplitude; timescale; spatial patterns; phase-locking; spring persistence barrier; and recharge oscillator dynamics. However, some important features of ENSO are still a challenge to simulate. In the central and eastern equatorial Pacific, the models’ weaker-than-observed subsurface zonal current anomalies and zonal temperature gradient anomalies serve to weaken the nonlinear zonal advection of subsurface temperatures, leading to insufficient warm/cold asymmetry of ENSO’s sea surface temperature anomalies (SSTA). In the western equatorial Pacific, the models’ excessive simulated zonal SST gradients amplify their zonal temperature advection, causing their SSTA to extend farther west than observed. The models underestimate both ENSO’s positive dynamic feedbacks (due to insufficient zonal wind stress responses to SSTA) and its thermodynamic damping (due to insufficient convective cloud shading of eastern Pacific SSTA during warm events); compensation between these biases leads to realistic linear growth rates for ENSO, but for somewhat unrealistic reasons. The models also exhibit stronger-than-observed feedbacks onto eastern equatorial Pacific SSTAs from thermocline depth anomalies, which accelerates the transitions between events and shortens the simulated ENSO period relative to observations. Implications for diagnosing and simulating ENSO in climate models are discussed.

54 ENVIRONMENTAL SCIENCES↗

Effect of q-nonextensive parameter and saturation time on electron density steepening in electron-positron-ion plasmas

The effect of q-nonextensive parameter and saturation time on the electron density steepening in electron-positron-ion plasmas is studied by particle in cell method. Phase space diagrams show that the size of the holes, and consequently, the number of trapped particles strongly depends on the q-parameter and saturation time. Furthermore, the mechanism of the instability and exchange of energy between electron-positron and electric field is explained by the profiles of the energy density. Moreover, it is found that the q-parameter, saturation time, and electron and positron velocities affect the nonlinear evolution of the electron density which leads to the steepening of its structure. The q-nonextensive parameter or degree of nonextensivity is the relation between temperature gradient and potential energy of the system. Therefore, the deviation of q-parameter from unity indicates the degree of inhomogeneity of temperature or deviation from equilibrium. Finally, using the kinetic theory, a generalized q-dispersion relation is presented for electron-positron-ion plasma systems. It is found that the simulation results in the linear regime are in good agreement with the growth rate results obtained by the kinetic theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deformation analysis of SiC-SiC channel box for BWR applications

Silicon carbide fiber-reinforced silicon carbide matrix (SiC-SiC) composites are being considered as components in light water reactor cores to improve accident tolerance, including channel boxes and fuel cladding. In the nuclear reactor environment, core components like a channel box will be exposed to neutron and other radiation damage and temperature gradients. To ensure reliable and safe operation of a SiC-SiC channel box, it is important to assess its deformation behavior under in-reactor conditions including the expected neutron flux and temperature distributions. In particular, this work has evaluated the effect of non-uniform dimensional changes caused by spatially varying neutron flux and temperatures on the deformation behavior of the channel box over the course of one year. These analyses have been performed using the fuel performance modeling code BISON and the commercial finite element analysis code Abaqus, based on fast flux and temperature boundary conditions that have been calculated using the neutronics and thermal-hydraulics codes Serpent and CTF, respectively. The dependence of dimensions and thermophysical properties on fast flux and temperature has been incorporated into the material models. These initial results indicate significant bowing of the channel box with a lateral displacement greater than 6.5 mm. The channel box bowing behavior is time dependent and driven by the temperature dependence of the SiC irradiation-induced swelling and the neutron flux/fluence gradients. The bowing behavior gradually recovers during the course of the operating cycle as the swelling of the SiC-SiC material saturates. Furthermore, the bending relaxation due to temperature gradients does not fully recover and residual bending remains after the swelling saturates in the entire channel box.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Unconventional ballooning structures for toroidal drift waves

With strong gradients in the pedestal of high confinement mode (H-mode) fusion plasmas, gyrokinetic simulations are carried out for the trapped electron and ion temperature gradient modes. A broad class of unconventional mode structures is found to localize at arbitrary poloidal positions or with multiple peaks. It is found that these unconventional ballooning structures are associated with different eigen states for the most unstable mode. At weak gradient (low confinement mode or L-mode), the most unstable mode is usually in the ground eigen state, which corresponds to a conventional ballooning mode structure peaking in the outboard mid-plane of tokamaks. However, at strong gradient (H-mode), the most unstable mode is usually not the ground eigen state and the ballooning mode structure becomes unconventional. This result implies that the pedestal of H-mode could have better confinement than L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and numerical validation of an ultrasonic visualization system for nuclear reactor application

This paper deals with the design of ultrasonic imaging systems for the next generation of nuclear reactors cooled by liquid metal. Indeed, a generation IV research nuclear reactor is being developed by the Belgian Nuclear Research Center (SCK-CEN) in the frame of the MYRRHA project (for Multipurpose hYbrid Research Reactor for High-tech Applications). This sub-critical/critical reactor is cooled by Lead-Bismuth Eutectic (LBE). The opacity of this liquid metal requires the development of an ultrasonic visualization system for internal inspection and object detection. But due to the peculiar conditions met in the core of the reactor, velocity as well as temperature gradients are expected and are likely to affect directly the ultrasonic propagation. The objective of this work is to validate the ultrasonic imaging strategy by tests performed in a dedicated test rig and by numerical simulations using a ray-tracing method. The experimental investigations have been performed on a specific water facility reproducing conditions similar to those encountered in the core of the MYRRHA reactor. These conditions include the propagation over large distance, and the presence of temperature and velocity gradients. In the MYRRHA reactor application, the distance to be travelled by the acoustic waves of the visualization system is about 5 m, including a reflection. The acoustic absorption, the scattering losses, the beam divergence and the transmitted energy during reflection have been determined as a function of the travelled distance. The experimental values are compared with the literature for validation. The presence of temperature and velocity gradients in the core of the reactor is due to the coolant circulation. These gradients are about 5 K over 0.1 m and 1 m/s over 0.2 m, respectively, and are reproduced in the facility for investigating their influence on the propagation of ultrasounds. The experimental data are used for improving and validating a ray-tracing algorithm. This algorithm is then applied to the full scale application for validating the ultrasonic visualization system for nuclear reactor application. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Full-f gyrokinetic simulation over a confinement time

A long time ion temperature gradient driven turbulence simulation over a confinement time is performed using the full-f gyrokinetic Eulerian code GT5D. The convergence of steady temperature and rotation profiles is examined, and it is shown that the profile relaxation can be significantly accelerated when the simulation is initialized with linearly unstable temperature profiles. In the steady state, the temperature profile and the ion heat diffusivity are self-consistently determined by the power balance condition, while the intrinsic rotation profile is sustained by complicated momentum transport processes without momentum input. The steady turbulent momentum transport is characterized by bursty non-diffusive fluxes, and the resulting turbulent residual stress is consistent with the profile shear stress theory [Y. Camenen et al., “Consequences of profile shearing on toroidal momentum transport,” Nucl. Fusion 51, 073039 (2011)] in which the residual stress depends not only on the profile shear and the radial electric field shear but also on the radial electric field itself. Based on the toroidal angular momentum conservation, it is found that in the steady null momentum transport state, the turbulent residual stress is cancelled by the neoclassical counterpart, which is greatly enhanced in the presence of turbulent fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Oscillatory/Chaotic Thermocapillary Flow Induced by Radiant Heating

There is a continuing need to understand the fluid physics occurring under low gravity conditions in processes such as crystal growth, materials processing, and the movement of bubbles or droplets. The fluid flow in such situations is often caused by a gradient in interfacial tension. If a temperature gradient is created due to a heat source, the resulting flow is called thermocapillary flow, a special case of Marangoni Convection. In this study, an experimental investigation was conducted using silicone oil in cylindrical containers with a laser heat source at the free surface. It was desired to determine the conditions under which steady, axisymmetrical thermocapillary flow becomes unstable and oscillatory three-dimensional flow states develop. The critical Marangoni number for each observed oscillatory state was measured as a function of the container aspect ratio and the dynamic Bond number, a measure of buoyant force versus ii thermocapillary force. Various oscillatory modes were observed during three- dimensional convection, and chaotic flow was reached in one test condition. The critical Marangoni numbers are compared with those measured in previous studies, and the power spectra and phase trajectories of the instantaneous surface temperature distributions are used to characterize the routes of transitions to the chaotic flow state. Results show that only superharmonic modes appear in the routes to chaos while infinite number of subharmonic modes occur in flow transitions for pure Rayleigh convection.

DeWitt, Kenneth J.↗

Multi-temperature representation of electron velocity distribution functions. I. Fits to numerical results

Electron energy distribution functions are expressed as a sum of 6-12 Maxwellians or a sum of 3, but each multiplied by a finite series of generalized Laguerre polynomials. We fitted several distribution functions obtained from the finite difference Fokker-Planck code 'FPI'[Matte and Virmont, Phys. Rev. Lett. 49, 1936 (1982)] to these forms, by matching the moments, and showed that they can represent very well the coexistence of hot and cold populations, with a temperature ratio as high as 1000. This was performed for two types of problems: (1) the collisional relaxation of a minority hot component in a uniform plasma and (2) electron heat flow down steep temperature gradients, from a hot to a much colder plasma. We find that the multi-Maxwellian representation is particularly good if we accept complex temperatures and coefficients, and it is always better than the representation with generalized Laguerre polynomials for an equal number of moments. For the electron heat flow problem, the method was modified to also fit the first order anisotropy f{sub 1}(x,v,t), again with excellent results. We conclude that this multi-Maxwellian representation can provide a viable alternative to the finite difference speed or energy grid in kinetic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of injection and confinement improvement through impurity induced profile modifications on the Wendelstein 7-X stellarator

Pulsed injections of boron carbide granules into Wendelstein 7-X stellarator (W7-X) plasmas transiently increase the plasma stored energy and core ion temperatures above the reference W7-X experimental programs by up to 30%. In a series of 4 MW electron cyclotron resonance heating experiments, the PPPL Probe Mounted Powder Injector provided 50 ms bursts of 100 μm granules every 350 ms at estimated quantities ranging from approximately 1 mg/pulse to over 30 mg/pulse. For each injection, the stored energy was observed to initially drop and the radiated power transiently increased, while the radial electron density profile rose at the edge as material was assimilated. Once the injected boron carbide was fully absorbed, the density rise transitioned to the core while the stored energy increased above the previous baseline level by an amount linearly correlated with the injection quantity. During the injection, the ion temperature gradient steepened with peak core ion temperatures observed to increase from a nominal 1.7 keV to over 2.6 keV for the largest injection amounts. Enhanced performance is accompanied by a reversal of the radial electric field at ρ < 0.3, indicating that the core transport has switched to the ion root. These observations are suggestive of a change in transport and provide further evidence that externally induced profile modifications provide a possible path to enhanced W7-X performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Metrology for femtosecond pulsed x-ray heating in diamond anvil cell experiments at the European XFEL: Revisiting the iron phase diagram up to 150 GPa

The development of pulsed intense x-ray sources, such as free electron laser, offers new avenues for high pressure experiments. Here, we study the feasibility and metrology of x-ray heating in diamond anvil cells at the European x-ray free electron laser. This method enables one to volumetrically heat the sample while inhibiting chemical migration and probing the crystallographic structure of the sample throughout the heating with a high repetition rate. We focus our study on iron, whose phase diagram is well established up to 100 GPa, to explore the possibilities and limitations of this technique. We volumetrically heat iron samples at starting pressures ranging from 10 to 138 GPa, using the x-ray beam pulsed at 4.5 MHz in a serial pump-and-probe experimental design. Experimental challenges arise from temperature gradients within the sample, changes in temperature at the 100 ns timescale, the difficulty of direct temperature estimates, the effect of thermal pressure, and the presence of metastable crystallites due to rapid cycles of heating and cooling. Hence, we develop a multi-crystal-like data processing method that allows us to account for sample heterogeneity in probed conditions. We then calibrate our measurements using known physical properties of iron under pressure. Thermal pressure in our experiments increases from 4% of the isochoric prediction at 10 GPa to 23% at 138 GPa, and we show that our data are in agreement with most previous observations of iron in this pressure range. The method can now be implemented at higher pressures and temperatures and on materials with unknown phase diagrams.

Materials science↗

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microwave-assisted, performance-advantaged electrification of propane dehydrogenation

Nonoxidative propane dehydrogenation (PDH) produces on-site propylene for value-added chemicals. While commercial, its modest selectivity and catalyst deactivation hamper the process efficiency and limit operation to lower temperatures. We demonstrate PDH in a microwave (MW)–heated reactor over PtSn/SiO 2 catalyst pellets loaded in a SiC monolith acting as MW susceptor and a heat distributor while ensuring comparable conditions with conventional reactors. Time-on-stream experiments show active and stable operation at 500°C without hydrogen addition. Upon increasing temperature or feed partial pressure at high space velocity, catalysts under MWs show resistance in coking and sintering, high activity, and selectivity, starkly contrasting conventional reactors whose catalyst undergoes deactivation. Mechanistic differences in coke formation are exposed. Gas-solid temperature gradients are computationally investigated, and nanoscale temperature inhomogeneities are proposed to rationalize the different performances of the heating modes. The approach highlights the great potential of electrification of endothermic catalytic reactions.

42 ENGINEERING↗

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗

Reduced African Easterly Wave Activity with Quadrupled CO 2 in the Superparameterized CESM

African easterly wave (AEW) activity is examined in quadrupled CO 2 experiments with the superparameterized CESM (SP-CESM). The variance of 2–10-day filtered precipitation increases with warming over the West African monsoon region, suggesting increased AEW activity. The perturbation enstrophy budget is used to investigate the dynamic signature of AEW activity. The northern wave track becomes more active associated with enhanced baroclinicity, consistent with previous studies. The southern track exhibits a surprising reduction of wave activity associated with less frequent occurrence of weak waves and a slight increase in the occurrence of strong waves. These changes are connected to changes in the profile of vortex stretching and tilting that can be understood as interconnected consequences of increased static stability from the lapse rate response, weak temperature gradient balance, and the fixed anvil temperature hypothesis.

58 GEOSCIENCES↗

Processing eutectics in space

The examination of the Ni-Mo-Al system is discussed. Experiments were performed to examine the role of thermal gradient magnitude at the solidifying interface on the geometrical regularity of the microstructure. The difficulty in measuring the thermal gradient in the melt using superheat temperatures as high as 1800 C resulted in the examination of three liquid temperature gradient ranges. Examination of the microstructure from directionally solidified ingots using transverse sections indicated that the highest gradient runs produced microstructure of no greater regularity than runs made at lower gradients. For the Ni-Al-Mo monovariant eutectic system, the magnitude of the thermal gradient at the solidifying interface is not a controlling parameter with respect to the geometrical regularity of the microstructure, provided it is sufficient to satisfy the gradient/rate ratio needed to produce a planar interface growth.

Douglas, F. C.↗