Search NASA⌕ Search

SEARCH · Search NASA

Results for “temperature gradient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

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.↗

Millimeter-wave experiments for cometary space missions

Predicted brightness temperatures, computed by means of radiative transfer techniques adapted from the modeling of terrestrial ice and snow fields, are given for cometary nucleus models consisting of homogeneous layers of water ice and refractory grain mixtures presented as functions of wavelength. The computed millimeter-wave spectra are sensitive to the values of such physically significant nucleus parameters as (1) crust thickness, (2) subsurface temperature gradient, and (3) sublimating surface boundary temperature. Although antenna beam dilution is a major obstacle for ground-based molecular spectral line radio observations of comets, a millimeter-wave radiometer in the vicinity of the comet would be immune to this effect and able to make observations of several candidate parent molecules in the gas phase.

Hobbs, R. W.↗

Evolution of density in solar system ices

Pores present in ices in the solar system do not remain unchanged. In isothermal conditions they shrink, while in a thermal gradient they migrate towards the higher temperature and escape so that the ice densifies. This motion has been investigated for pure H2O- and CO2-ices in a very simple one-dimensional model assuming uniform thermal conductivity and temperature gradient. The results indicate that the densification of H2O-ice is so slow that it could be significant only for icy satellites having an internal heat source. On the other hand, CO2-ice densifies orders of magnitude faster and the effect should be important for the CO2 component of cometary nuclei. No effect is expected for icy planetary rings.

Smoluchowski, R.↗

Development of BEM for ceramic composites

It is evident that for proper micromechanical analysis of ceramic composites, one needs to use a numerical method that is capable of idealizing the individual fibers or individual bundles of fibers embedded within a three-dimensional ceramic matrix. The analysis must be able to account for high stress or temperature gradients from diffusion of stress or temperature from the fiber to the ceramic matrix and allow for interaction between the fibers through the ceramic matrix. The analysis must be sophisticated enough to deal with the failure of fibers described by a series of increasingly sophisticated constitutive models. Finally, the analysis must deal with micromechanical modeling of the composite under nonlinear thermal and dynamic loading. This report details progress made towards the development of a boundary element code designed for the micromechanical studies of an advanced ceramic composite. Additional effort has been made in generalizing the implementation to allow the program to be applicable to real problems in the aerospace industry.

Henry, D. P.↗

Large-Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulations↗

Large Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulation↗

Building the Next Generation of Regional Geothermal Potential Maps: Examples from the Great Basin Region, Western USA

Over the past decade research has greatly advanced our understanding of structurally controlled geothermal systems and the potential for “deep stratigraphic reservoirs” in the Great Basin region, western USA. For example, studies have shown that structurally controlled geothermal systems are associated with specific structural settings, such as terminations of major normal faults, accommodation zones, pull-aparts in strike-slip faults, displacement transfer zones, or step-overs in range-front faults. However, we know that not every structure (e.g., termination of a major normal fault or step-over in a range-front fault) hosts an active geothermal system. Recent research has also shown that regional strain rate model as determined from GPS geodesy, the composition of that strain (extensional versus shear or compressive strain), recency of Quaternary faulting, and slip and dilation tendency of faults correlate positively with the distribution of both high and low-temperature, structurally controlled geothermal systems within the Great Basin. Combining these strain data and fault slip history data sets with existing temperature gradient data, heat flow data, spring temperatures, fluid geothermometry, and geophysical data provides an effective platform for evaluating the geothermal potential of respective structural plays. Recent studies focused on deep stratigraphic reservoirs indicate the best conditions for economically viable development are associated with thick sedimentary basins with low-thermal-conductivity sedimentary cap rocks overlying permeable host rocks (e.g. carbonates) and co-located with areas of high heat flow. Some of the factors that contribute to the potential for deep stratigraphic reservoirs and conventional structurally controlled geothermal systems apply to both types of resources, and in some places the potential for these two types of resources overlap in space, suggesting the potential for hybrid reservoirs (structurally controlled deep stratigraphic reservoirs). Effective and complete evaluation of geothermal potential in the Great Basin should include parameters that cover both types of resources. In this study, two areas with a combined area of 26,000 km 2 were evaluated for geothermal resource potential, equating to <4% of the total 680,000 km 2 size of the Great Basin region. Within these two study areas, 73 potential structural targets and 8 potential deep sedimentary target areas were identified. Four of the 73 potential structural target areas are associated with known thermal anomalies. All 73 structural target areas were ranked according to multiple data sets.

15 GEOTHERMAL ENERGY↗