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 343 records · Page 19

High spatial resolution temperature profile measurements of solid-oxide fuel cells

Temperature gradients resulting from local electrochemical reactions, current distribution and geometry of gas flow channels in solid oxide fuel cells (SOFCs) create thermal stresses, localized thermophysical property gradients and uneven property evolution, contributing to SOFC degradation. This paper presents a new method to perform temperature measurements (up to 800°C) at high spatial resolutions to monitor the operation of SOFCs. Using femtosecond laser irradiation, distributed fiber sensors were hardened for high temperature environment applications. Distributed fiber sensors were embedded in interconnected plates using an additive manufacturing method to perform temperature measurements with 4-mm spatial resolution during the operation of a planar fuel cell. The measurement revealed the impact of various H 2 fuel concentrations and current loads have on temperature profiles of the SOFC tested. Temperature variation on the anode side was found to be less than 5°C, and 3°C on the cathode side. The measurements were compared to results from a multiphysics fuel cell performance model simulating similar conditions. These simulations predicted similar temperature gradients, indicating the experimental data obtained is reasonable. The model also predicts that the effect of the embedded sensor has on the local temperature will be minimal and that the gradient of temperature in the gas channels will be captured despite the separation between the sensor and the gas flow. Finally, the high spatial resolution data harnessed by these distributed fiber sensors provides experimental support for model-based design and optimization to improve the operational efficiency and longevity of solid oxide fuel cells and fuel cell assemblies.

25 ENERGY STORAGE↗

Understanding Irrigation Impacts on Low-Level Jets over the Great Plains

Low-level jets (LLJs) over the Great Plains are well-known for their importance in transporting moisture from the Gulf of Mexico and forming precipitation in the Central US. However, the impact of irrigation practice on LLJs is rarely studied. To understand the irrigation impact on the Great Plains LLJs, the Weather Research and Forecasting (WRF) model is employed with and without an irrigation scheme. The simulated wind field is evaluated against the observed wind from Midlatitude Continental Convective Clouds Experiment (MC3E). The results show that adding irrigation improves the simulated wind, especially at lower altitudes below 800 hPa. Irrigation reduces the LLJ frequency over the Great Plains, mainly by reducing the intense LLJ frequency. The reduced LLJ frequency is caused by the irrigation-induced weakening of meridional wind. The favorable conditions for LLJs have been identified by examining the mechanism associated with LLJs: a ridge-like pattern with weak and calm synoptic-scale conditions in the mid-levels and a strong west-east temperature gradient at lower levels that promote meridional winds. Irrigation-induced cooling leads to a baroclinic pattern with high pressure at lower altitudes over the irrigation region and low pressure in the mid- and upper levels in the downstream region, resulting in a northerly component that hinders the formation of LLJs. Low-level surface cooling also decreases the west-east temperature gradient and hence weakens the thermal wind forcing. The irrigation-induced mid-level low and weaker west-east temperature gradient are responsible for the reduced meridional winds and the reduced LLJ frequencies.

Yang, Zhao↗

Geometric origin of the intrinsic transverse spin transport in a canted-antiferromagnet/heavy-metal heterostructure

We theoretically study the conditions under which an intrinsic spin Nernst effect–a transverse spin current induced by an applied temperature gradient–can occur in a canted-antiferromagnet insulator, such as LaFeO 3 and other materials of the same family. The spin Nernst effect may provide a microscopic mechanism for an experimentally observed anomalous thermovoltage in LaFeO 3 /Pt heterostructures, where spin is transferred across the insulator/metal interface when a temperature gradient is applied to LaFeO 3 parallel to the interface. We find that LaFeO 3 exhibits an intrinsic spin Nernst effect when inversion symmetry is broken on the axes parallel to both the applied temperature gradient and the direction of spin transport, which can result in a spin injection across the insulator/metal interface. Furthermore, our paper provides a general derivation of a symmetry-breaking-induced spin Nernst effect, which may open a path to engineering a finite spin Nernst effect in systems where it would otherwise not arise.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Global and Zonal-Mean Hydrological Response to Early Eocene Warmth

Earth's hydrological cycle is expected to intensify in response to global warming, with a “wet-gets-wetter, dry-gets-drier” response anticipated over the ocean. Subtropical regions (∼15°–30°N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterized by wetter conditions. Here we use an integrated data-modeling approach to reconstruct global and zonal-mean rainfall patterns during the early Eocene (∼56–48 million years ago). The Deep-Time Model Intercomparison Project (DeepMIP) model ensemble indicates that the mid- (30°–60°N/S) and high-latitudes (>60°N/S) are characterized by a thermodynamically dominated hydrological response to warming and overall wetter conditions. The tropical band (0°–15°N/S) is also characterized by wetter conditions, with several DeepMIP models simulating narrowing of the Inter-Tropical Convergence Zone. However, the latter is not evident from the proxy data. The subtropics are characterized by negative precipitation-evaporation anomalies (i.e., drier conditions) in the DeepMIP models, but there is surprisingly large inter-model variability in mean annual precipitation (MAP). Intriguingly, we find that models with weaker meridional temperature gradients (e.g., CESM, GFDL) are characterized by a reduction in subtropical moisture divergence, leading to an increase in MAP. These model simulations agree more closely with our new proxy-derived precipitation reconstructions and other key climate metrics and imply that the early Eocene was characterized by reduced subtropical moisture divergence. If the meridional temperature gradient was even weaker than suggested by those DeepMIP models, circulation-induced changes may have outcompeted thermodynamic changes, leading to wetter subtropics. This highlights the importance of accurately reconstructing zonal temperature gradients when reconstructing past rainfall patterns.

Margot J. Cramwinckel↗

Characteristics of Bay of Bengal Monsoon Depressions in the 21st Century

The controversy surrounding the observed trends in the frequency of South Asian Monsoon Depressions (MDs)1-4, the association of MDs with the summer season extremes5,6, and a moisture driven future precipitation response in Global Climate Models (GCMs)7,8 have elevated the importance of understanding climate change impacts on the characteristics of MDs. GCMs exhibit a significant relationship between the long-term mean occurrences of MDs and the strength of the background upper (lower) tropospheric meridional (zonal) winds and atmospheric moisture in the core genesis region of MDs. Likewise, the strength of the meridional tropospheric temperature gradient (MTG) in GCMs determines the trajectories of MDs over land. While monsoon dynamics progressively weakens in the future period, atmospheric moisture exhibits a strong increase, limiting the impact of changes in background dynamics on the frequency of MDs. Moreover, the weakening of the future MTG is many times weaker than its inherent underestimation in the GCMs. Hence the 21st century increase in radiative forcing does not significantly reduce the frequency of MDs or impact their trajectories over land in the GCMs. Although more wet extremes are expected to occur during the MD lifecycle, non-depression day extremes dominate the overall increase in extreme precipitation, which may render the monsoon extremes less predictable in the future.

Rastogi, Deeksha↗

Modeling electron temperature profiles in the pedestal with simple formulas for ETG transport

Abstract This paper reports on the refinement (building on (Hatch D.R. et al 2022 Phys. Plasmas 29 062501)) and application of simple formulas for electron heat transport from electron temperature gradient (ETG) driven turbulence in the pedestal. The formulas are improved by (1) improving the parameterization for certain key parameters and (2) carefully accounting for the impact of geometry and shaping in the underlying gyrokinetic simulation database. Comparisons with nonlinear gyrokinetic simulations of ETG transport in the MAST pedestal demonstrate the model’s applicability to spherical tokamaks in addition to standard aspect ratio tokamaks. We identify bounds for model applicability: the model is accurate in the steep gradient region, where the ETG turbulence is largely slab-like, but accuracy decreases as the temperature gradient becomes weaker in the pedestal top. We use the formula to model the electron temperature profile in the pedestal for four experimental scenarios while extensively varying input parameters to represent uncertainties. In all cases, the predicted electron temperature profile exhibits extreme sensitivity to separatrix temperature and density, which has implications for core-edge integration. The model reproduces the electron temperature profile for high η e = L n e / L T e scenarios but not for low η e scenarios in which microtearing modes (MTMs) have been identified. We develop a proof-of-concept model for MTM transport and explore the relative roles of ETG and MTM in setting the electron temperature profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Combined Structural Analysis of Core and Image Log of TGH 76-31 South East of Mt Baker, Washington State

Despite active volcanism, few geothermal energy resources have been developed in the Cascades Range. Temperature Gradient Hole 76-31 was drilled to ~440 m measured depth to probe for zones where fractures provide fluid conduits that transport deep volcanic heat to shallow depths that could support baseload, carbon neutral electrical generation. These zones were predicted by a Play Fairway assessment (PFA) of resource potential along a zone 11 km west-southwest of the summit of Mount Baker Volcano. Rock core, temperature logs, and an acoustic image log were obtained. By comparison to outcrops, the core has been initially interpreted as the Chilliwack group comprised of partially metamorphosed basaltic to andesitic volcanics. Core mapping reveals complex, steeply dipping networks of fractures and brecciation along slickensided strike slip faults; clay alteration is common in many of these structures. The majority of fractures are thoroughly healed by layers of chlorite and calcite; preservation of pyrite indicates anoxic conditions. The majority of fracture porosity resides in very dense fractures a few centimeters or less in length. The image log provides good insights into attitude of fractures that fully transect the core, but generally underestimates fracture density. The combination of complex, non-planar fracture zones containing many short fractures and healing promote misinterpretation of natural fracture attitude and density in the image log. The healing and anoxic minerals are consistent with the conductive temperature gradient measured in the well below a shallow isothermal zone, although, several fractures are open or only partially healed and resulted in fluid entries into the well. Here, the equilibrated measured temperature gradient of 64ºC/km and calculated heat flow of 145 mW/m2 is more than twice the regional average, indicating local influence of the Mt Baker magmatic system at the Little Park Creek TGH site.

Cascades↗

Gyrokinetic GENE simulations of DIII-D near-edge L-mode plasmas

We present gyrokinetic simulations with the GENE code addressing the near-edge region of an L-mode plasma in the DIII-D tokamak. At radial position ρ = 0.80, simulations with the ion temperature gradient (ITG) increased by 40% above the nominal value give electron and ion heat fluxes that are in simultaneous agreement with the experiment. This gradient increase is consistent with the combined statistical and systematic uncertainty σ of the charge exchange recombination spectroscopy measurements at the 1.6σ level. Multiscale simulations are carried out with a realistic mass ratio and geometry for the first time in the near-edge. These multiscale simulations suggest that the highly unstable ion temperature gradient (ITG) modes of the flux-matched ion-scale simulations suppress electron-scale transport, such that ion-scale simulations are sufficient at this location. At radial position ρ = 0.90, nonlinear simulations show a hybrid state of ITG and trapped electron modes, which was not expected from linear simulations. The nonlinear simulations reproduce the total experimental heat flux with the inclusion of E × B shear effects and an increase in the electron temperature gradient by ∼23%. This gradient increase is compatible with the combined statistical and systematic uncertainty of the Thomson scattering data at the 1.3σ level. These results are consistent with previous findings that gyrokinetic simulations are able to reproduce the experimental heat fluxes by varying input parameters close to their experimental uncertainties, pushing the validation frontier closer to the edge region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation and Life Cycle of Solar Spicules

The physical mechanism for the creation of solar spicules is proposed with three stages of their life cycle. It is assumed that at stage I the density hump is formed locally in the x-y plane in the lower chromosphere in the presence of temperature gradients of electrons and ions along the z-axis (the vertical direction). In this region, the density structure of quasi-neutral (n i ≃ n e = n) plasma after taking birth is accelerated in the vertical direction owing to the thermal force F th ∝ ∇n(x, y, t) × (∇T e + ∇T i ). The exact time-dependent analytical solution of two-fluid plasma equations is presented assuming that density is maximum at the center of the density structure and decays away from it gradually. The 2D density structure is created as a step function H(t) in time at the bottom of the chromosphere, and consequently, the vertical plasma velocity turns out to be the ramp function of time R(t) = tH(t), whereas the source term S(x, y, t) for the density follows the delta function δ(t) form. The upward acceleration a = a(x, y) $\hat{z}$ produced in this density structure is greater than the downward constant solar acceleration −g ⊙ in the chromosphere. In the transition region, the temperature gradients are steeper; therefore, the upward acceleration increases in magnitude g ⊙ ≪ a and the density hump spends less time there. This is stage II of its life cycle. In stage III, the density structure enters into the corona, where the gradients of temperatures vanish and the structure decelerates to zero velocity under the action of the solar gravitational force.

79 ASTRONOMY AND ASTROPHYSICS↗

Partnership for Edge Physics Simulation

A major goal of our participation in the Edge Physics Simulation project has been to contribute to the understanding of the self-organization of tokamak turbulence fluctuations resulting in the formation of a staircase structure in the ion temperature. A second important goal is to demonstrate how small scale turbulence in plasmas self-organizes with dynamically driven quasi-stationary flow shear. These goals have been accomplished through the analyses of the statistical properties of XGC1 flux driven Gyrokinetic electrostatic ion temperature gradient (ITG) turbulence simulation data in which neutrals are included. The ITG turbulence data, and in particular fluctuation data, were obtained from a massively parallel flux-driven gyrokinetic full-f particle-in-cell simulation of a DIII-D like equilibrium. Below some the findings are summarized. It was observed that the emergence of staircase structure is related to the variations in the normalized temperature gradient length (R/LT) and the poloidal flow shear. Average turbulence intensity is found to be large in the vicinity of minima in R/LTi, where ITG growth is expected to be lower. The distributions of the occurrences of potential fluctuation are found to be Gaussian away from the staircase-step locations, but they are found to be non-Gaussian in the vicinity of staircase-step locations. The results of analytically derived expressions for the distribution of the occurrences of turbulence intensity and intensity flux were compared with the corresponding quantities computed using XGC1 simulation data and good agreement is found. The derived expressions predicts inward and outward propagation of turbulence intensity flux in an intermittent fashion. The outward propagation of turbulence intensity flux occurs at staircase-step locations and is related to the change in poloidal flow velocity shear and to the change in the ion temperature gradient. The standard deviation, skewness and kurtosis for turbulence quantities were computed and found to be large in the vicinity of the staircase-step structures. Large values of skewness and kurtosis can be explained by a temporary opening and closing of the structure which allows turbulence intensity events to propagate. The staircase patterns may reduce the ion heat transport and a manipulation of these patterns may be used to optimize heat transport in tokamaks. An additional objective of the research in support of the Edge Physics Simulation initiative has been to improve the understanding of scrape-off layer thermal transport. In planning experiments and designing future tokamaks, it is important to understand the physical effects that contribute to divertor heat-load fluxes. The research accomplished will contribute to developing new models for the scrape-off layer region. The XGC0 code was used to compute the heat fluxes and the heat-load width in the outer divertor plates of C-Mod and DIII-D tokamaks. It was observed that the width of the XGC0 neoclassical heat-load was approximately inversely proportional to the total plasma current. Anomalous transport in the H-mode pedestal region of five Alcator C-Mod discharges, representing a collisionality scan, was analyzed. The understanding of anomalous transport in the pedestal region is important for the development of a comprehensive model for the H-mode pedestal slope. It was found that the electron thermal anomalous diffusivities at the pedestal top increase with the electron collisionality. This dependence can point to the DRIBM as the modes that drive the anomalous transport in the plasma edge of highly collisional discharges. The effects of plasma shaping on the H-mode pedestal structure was also investigated. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities were discussed. For the discharges with higher elongation, it was found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continued to evolve to higher pedestal pressures and bootstrap currents until the peeling ballooning stability conditions were satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes were triggered at earlier times. The plasma elongation was found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency was large, and the H-mode pedestal evolves rapidly. It was found that the temperature of neutrals in the scrape-off-layer region can affect the dynamics of the H-mode pedestal buildup. However, the final pedestal profiles were nearly independent of the neutral temperature. The elongation and triangularity affected the pedestal widths of plasma density and electron temperature profiles differently. This study illustrated a new mechanism for controlling the pedestal bootstrap current and the pedestal stability.

97 MATHEMATICS AND COMPUTING↗

Washington Geothermal Play Fairway Analysis Heat, Permeability, and Fracture Model Data

This submission contains raster and vector data for the entire state of Washington, with specific emphasis on the three geothermal play fairway sites: Mount St. Helens seismic zone (MSHSZ), Wind River valley (WRV), and Mount Baker (MB). Data are provided for 3 major geothermal models: heat, permeability, and fluid-filled fractures, and an additional infrastructure model. Both of the permeability and fluid-filled-fracture models are produced at 200 m and at 2 km depths; the heat model is only produced at the 200 m depth. Values are provided for both model favorability and model confidence. A combined model at 200m and 2 km depths is provided for favorability, confidence, and exploration risk. Raster data are provided in GeoTiff format and have a statewide coverage. Cell size is 104.355 ft; file type is unsigned 8-bit integer (0-255); 0 represents no favorability or confidence; 255 represents maximum favorability or confidence. The NAD83(HARN)/Washington South (ftUS) projection is used (EPSG:2927). Vector data are provided in shapefile or comma-delimited text file formats. Geographic coordinates, where provided, are in WGS84. A readme file accompanies each folder and provides an overview and description of the enclosed data. The heat model combines 5 intermediate raster layers (which are included in the download package): temperature gradient wells, young volcanic vents, hot springs, young intrusive volcanic rocks, and geothermometry. The permeability model combines 8 intermediate raster layers: density of mapped faults, 2D dilation tendency of mapped faults, 2D slip tendency of mapped faults, seismicity, 3D dilation tendency, 3D slip tendency, 3D maximum coulomb shear stress, and 3D slip gradients. The fluid-filled fracture model combines up to 4 intermediate rasters: resistivity from magneto-telluric 3D inversions, seismicity, Vp/Vs anomalies from passive seismic tomography, and Vs anomalies from ambient-noise tomography. A statewide infrastructure model is also provided that formalizes land-use constraints and restrictions relevant for geothermal prospecting and development. This model combines 10 intermediate rasters: areas off limits to drilling, existing or proposed geothermal leases, DNR-owned land, land-use restrictions along the Columbia River Gorge, areas inundated by water, availability of potential process water, proximity to existing roads, proximity to transmission lines, distance from urban areas, and snow-related elevation restrictions. Supporting vector data for the development of each raster layer is provided. For details on the areas of interest and modeling process please see the 'WA_State_Play_Fairway_Phase_2_Technical_Report' in the download package.

15 GEOTHERMAL ENERGY↗

Solid-propellant rocket motor ballistic performance variation analyses

Results are presented of research aimed at improving the assessment of off-nominal internal ballistic performance including tailoff and thrust imbalance of two large solid-rocket motors (SRMs) firing in parallel. Previous analyses using the Monte Carlo technique were refined to permit evaluation of the effects of radial and circumferential propellant temperature gradients. Sample evaluations of the effect of the temperature gradients are presented. A separate theoretical investigation of the effect of strain rate on the burning rate of propellant indicates that the thermoelastic coupling may cause substantial variations in burning rate during highly transient operating conditions. The Monte Carlo approach was also modified to permit the effects on performance of variation in the characteristics between lots of propellants and other materials to be evaluated. This permits the variabilities for the total SRM population to be determined. A sample case shows, however, that the effect of these between-lot variations on thrust imbalances within pairs of SRMs is minor in compariosn to the effect of the within-lot variations. The revised Monte Carlo and design analysis computer programs along with instructions including format requirements for preparation of input data and illustrative examples are presented.

Sforzini, R. H.↗

Sources of Variation in Creep Testing

Creep rupture is an important material characteristic for the design of rocket engines. It was observed during the characterization of GRCop-84 that the complete data set had nearly 4 orders of magnitude of scatter. This scatter likely confounded attempts to determine how creep performance was influenced by manufacturing. It was unclear if this variation was from the testing, the material, or both. Sources of variation were examined by conducting tests on identically processed specimens at the same specified stresses and temperatures. Significant differences existed between the five constant-load creep frames. The specimen temperature was higher than the desired temperature by as much as 43 C. It was also observed that the temperature gradient was up to 44 C. Improved specimen temperature control minimized temperature variations. The data from additional tests demonstrated that the results from all five frames were comparable. The variation decreased to 1/2 order of magnitude from 2 orders of magnitude for the baseline data set. Independent determination of creep rates in a reference load frame closely matched the creep rates determined after the modifications. Testing in helium tended to decrease the sample temperature gradient, but helium was not a significant improvement over vacuum.

Loewenthal, William S.↗

Protein crystal growth in low gravity

This Final Technical Report for NASA Grant NAG8-774 covers the period from April 27, 1989 through December 31, 1992. It covers five main topics: fluid flow studies, the influence of growth conditions on the morphology of isocitrate lyase crystals, control of nucleation, the growth of lysozyme by the temperature gradient method and graphoepitaxy of protein crystals. The section on fluid flow discusses the limits of detectability in the Schlieren imaging of fluid flows around protein crystals. The isocitrate lyase study compares crystals grown terrestrially under a variety of conditions with those grown in space. The controlling factor governing the morphology of the crystals is the supersaturation. The lack of flow in the interface between the drop and the atmosphere in microgravity causes protein precipitation in the boundary layer and a lowering of the supersaturation in the drop. This lowered supersaturation leads to improved crystal morphology. Preliminary experiments with lysozyme indicated that localized temperature gradients could be used to nucleate crystals in a controlled manner. An apparatus (thermonucleator) was designed to study the controlled nucleation of protein crystals. This apparatus has been used to nucleate crystals of materials with both normal (ice-water, Rochelle salt and lysozyme) and retrograde (horse serum albumin and alpha chymotrypsinogen A) solubility. These studies have lead to the design of an new apparatus that small and more compatible with use in microgravity. Lysozyme crystals were grown by transporting nutrient from a source (lysozyme powder) to the crystal in a temperature gradient. The influence of path length and cross section on the growth rate was demonstrated. This technique can be combined with the thermonucleator to control both nucleation and growth. Graphoepitaxy utilizes a patterned substrate to orient growing crystals. In this study, silicon substrates with 10 micron grooves were used to grow crystals of catalase, lysozyme and canavalin. In all cases, the crystals grew oriented to the substrate. The supersaturation needed for nucleation and growth was lower on the patterned substrates. In some cases, isolated, large crystals were grown.

Feigelson, Robert S.↗

Midlatitude Cloud Systems

In contrast to the tropics and subtropics, the middle latitudes are characterised by large meridional temperature gradients, created as a consequence of differential radiative heating between high and low latitudes. These meridional temperature gradients often concentrate in relatively narrow baroclinic zones that become unstable to wave-like perturbations called baroclinic eddies, or more commonly baroclinic storms or midlatitude cyclones. Baroclinic storms constitute the primary source of poleward energy transport at midlatitudes, which is accomplished through contrasting transports of warm air masses poleward (warm fronts) and cold air masses equatorward (cold fronts). Baroclinic storms also flux momentum into midlatitude regions, driving a region of enhanced westerly winds from the surface to the upper troposphere called the eddy-driven jet stream.

Midlatitude clouds↗

Late-time Cooling of Neutron Star Transients and the Physics of the Inner Crust

An accretion outburst onto a neutron star transient heats the neutron star’s crust out of thermal equilibrium with the core. After the outburst, the crust thermally relaxes toward equilibrium with the neutron star core, and the surface thermal emission powers the quiescent X-ray light curve. Crust cooling models predict that thermal equilibrium of the crust will be established ≈1000 days into quiescence. Recent observations of the cooling neutron star transient MXB 1659-29, however, suggest that the crust did not reach thermal equilibrium with the core on the predicted timescale and continued to cool after ≈2500 days into quiescence. Because the quiescent light curve reveals successively deeper layers of the crust, the observed late-time cooling of MXB 1659-29 depends on the thermal transport in the inner crust. In particular, the observed late-time cooling is consistent with a low thermal conductivity layer near the depth predicted for nuclear pasta that maintains a temperature gradient between the neutron star’s inner crust and core for thousands of days into quiescence. As a result, the temperature near the crust–core boundary remains above the critical temperature for neutron superfluidity, and a layer of normal neutrons forms in the inner crust. We find that the late-time cooling of MXB 1659-29 is consistent with heat release from a normal neutron layer near the crust–core boundary with a long thermal time. We also investigate the effect of inner crust physics on the predicted cooling curves of the accreting transient KS 1731-260 and the magnetar SGR 1627-41.

79 ASTRONOMY AND ASTROPHYSICS↗

Temperature profiles in high gradient furnaces

Accurate temperature measurement of the furnace environment is very important in both the science and technology of crystal growth as well as many other materials processing operations. A high degree of both accuracy and precision is acutely needed in the directional solidification of compound semiconductors in which the temperature profiles control the freezing isotherm which, in turn, affects the composition of the growth with a concomitant feedback perturbation on the temperature profile. Directional solidification requires a furnace configuration that will transport heat through the sample being grown. A common growth procedure is the Bridgman Stockbarger technique which basically consists of a hot zone and a cold zone separated by an insulator. In a normal growth procedure the material, contained in an ampoule, is melted in the hot zone and is then moved relative to the furnace toward the cold zone and solidification occurs in the insulated region. Since the primary path of heat between the hot and cold zones is through the sample, both axial and radial temperature gradients exist in the region of the growth interface. There is a need to know the temperature profile of the growth furnace with the crystal that is to be grown as the thermal load. However it is usually not feasible to insert thermocouples inside an ampoule and thermocouples attached to the outside wall of the ampoule have both a thermal and a mechanical contact problem as well as a view angle problem. The objective is to present a technique of calibrating a furnace with a thermal load that closely matches the sample to be grown and to describe procedures that circumvent both the thermal and mechanical contact problems.

Fripp, A. L.↗

Plasma transport in an Eulerian AMR code

A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions to flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unit of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗