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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2019

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↗

DOE OSTI2013

Thermal transpiration through single walled carbon nanotubes and graphene channels

Thermal transpiration through carbon nanotubes (CNTs) and graphene channels is studied using molecular dynamics (MD) simulations. The system consists of two reservoirs connected by a CNT. It is observed that a flow is developed inside the CNT from the low temperature reservoir to the high temperature reservoir when the two reservoirs are maintained at different temperatures. The influence of channel size and temperature gradient on the mean velocity is analysed by varying the CNT diameter and the temperature of one of the reservoirs. Larger flow rate is observed in the smaller diameter CNTs showing an increase in the mean velocity with increase in the temperature gradient. For the flow developed inside the CNTs, slip boundaries occur and the slip length is calculated using the velocity profile. We examine the effect of fluid-wall interaction strength (ε{sub fw}), diffusivity (D), and viscosity of the fluid (μ) on the temperature induced fluid transport through the CNTs. Similar investigations are also carried out by replacing the CNT with a graphene channel. Results show that the mean velocity of the fluid atoms in the graphene channel is lower than that through the CNTs. This can be attributed to the higher degree of confinement observed in the CNTs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Global gyrokinetic stability of collisionless microtearing modes in large aspect ratio tokamaks

Linear full radius gyrokinetic calculations show the existence of unstable microtearing modes (MTMs) in purely collisionless, high temperature, large aspect ratio tokamak plasmas. The present study takes into account fully gyrokinetic highly passing ions and electrons. The global 2-D structures of the collisionless mode with full radius coupling of the poloidal modes is obtained and compared with another electromagnetic mode, namely, the Alfvén Ion Temperature Gradient (AITG) mode (or Kinetic Ballooning Mode, KBM) for the same equilibrium profile. Several important characteristics of the modes are brought out and compared, such as a clear signature in the symmetry properties of the two modes, the plasma–β dependence, and radial and poloidal length scales of the electrostatic and magnetic vector potential fluctuations. Extensive parameter scans for this collisionless microtearing mode reveal the scaling of the growth rate with β and the electron temperature gradient η{sub e}. Scans at different β values show an inverse relationship between the η{sub e} threshold and β, leading to a stability diagram, and implying that the mode might exist at moderate to strong temperature gradients for finite β plasmas in large aspect ratio tokamaks. In contrast to small aspect ratio tokamaks where the trapped electron magnetic drift resonance is found to be important, in large aspect ratio tokamaks, a strong destabilization due to the magnetic drift resonance of passing electrons is observed and is identified as a possible collisionless drive mechanism for the collisionless MTM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

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↗

DOE OSTI2017

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↗

DOE OSTI2013

Magnetic island evolution in the presence of ion-temperature gradient-driven turbulence

Turbulence is known to drive and sustain magnetic islands of width equal to multiples of the Larmor radius. The nature of the drive is studied here by means of numerical simulations of a fluid electrostatic model in 2D (single helicity) sheared-slab geometry. The electrostatic model eliminates the coalescence of short wavelength islands as a mechanism for sustaining longer wavelength islands. In quiescent islands, the polarization current, which depends on the propagation velocity of the island through the plasma, plays a critical role in determining the growth or decay of island chains. For turbulent islands, the unforced propagation velocity is significantly changed by strong zonal flow. The simulations show, however, that the turbulent fluctuations in the current density are much larger and faster than those in the zonal flow, and that they dominate the steady-state perturbed current density. In order to distinguish the roles of the zonal flow from the direct action of the fluctuations on the islands, a new diagnostic is implemented. This new diagnostic separates the effects of all the sources of parallel current. These are the curvature (which drives Pfirsch-Schlüter currents) and the divergences of the viscous and Reynolds stresses (the latter driving polarization currents). The new diagnostic also enables the contributions from short and long wavelengths to be separated for each term. It shows that in the absence of curvature, the drive is dominated by the contributions to the polarization current from the short wavelength fluctuations, while the long-wavelength fluctuations play a stabilizing role. In the presence of unfavorable curvature, by contrast, the effects of the short- and long-wavelength contributions of the polarization current reverse roles but nearly cancel, leaving the Pfirsch-Schlüter current as the dominant drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

Free energy balance in gyrokinetic turbulence

Free energy plays an important role in gyrokinetic theory, since it is known to be a nonlinear invariant. Its evolution equations are derived and analyzed for the case of ion temperature gradient driven turbulence, using the formalism adopted in the Gene code. In particular, the ion temperature gradient drive, the collisional dissipation as well as entropy/electrostatic energy transfer channels represented by linear curvature and parallel terms are analyzed in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

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↗

NASA NTRS1975

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

NASA NTRS2011

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

NASA NTRS1993

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

DOE OSTI2013

THERMAL STRUCTURE OF CORONAL LOOPS AS SEEN WITH NORIKURA CORONAGRAPH

The thermal structure of a coronal loop, both along and across the loop, is vital in determining the exact plasma heating mechanism. High-resolution spectroscopic observations of the off-limb corona were made using the 25 cm Norikura coronagraph, located at Norikura, Japan. Observations on a number of days were made simultaneously in four forbidden iron emission lines, namely, the [Fe XI] 7892 A line, the [Fe XIII] 10747 A and 10798 A lines, and the [Fe XIV] 5303 A line and on some days made only in the [Fe XI] 7892 A and [Fe X] 6374 A lines. Using temperature sensitive emission line ratios [Fe XIV] 5303 A/[Fe XIII] 10747 A and [Fe XI] 7892 A/[Fe X] 6374 A, we compute the electron temperatures along 18 different loop structures observed on different days. We find a significant negative temperature gradient in all of the structures observed in Fe XIV and Fe XIII and a positive temperature gradient in the structures observed in Fe XI and Fe X. Combining these results with the previous investigations by Singh and his collaborators, we infer that the loop tops, in general, appear hotter when observed in colder lines and colder when observed in relatively hotter lines as compared to their coronal foot points. We suggest that this contrasting trend observed in the temperature variation along the loop structures can be explained by a gradual interaction of different temperature plasma. The exact mechanism responsible for this interaction must be investigated further and has the potential to constrain loop heating models.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2020

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↗

DOE OSTI2013

Global gyrokinetic simulations of the H-mode tokamak edge pedestal

Global gyrokinetic simulations of DIII-D H-mode edge pedestal show two types of instabilities may exist approaching the onset of edge localized modes: an intermediate-n, high frequency mode which we identify as the “kinetic peeling ballooning mode (KPBM),” and a high-n, low frequency mode. Our previous study [W. Wan et al., Phys. Rev. Lett. 109, 185004 (2012)] has shown that when the safety factor profile is flattened around the steep pressure gradient region, the high-n mode is clearly kinetic ballooning mode and becomes the dominant instability. Otherwise, the KPBM dominates. Here, the properties of the two instabilities are studied by varying the density and temperature profiles. It is found that the KPBM is destabilized by density and ion temperature gradient, and the high-n mode is mostly destabilized by electron temperature gradient. Nonlinear simulations with the KPBM saturate at high levels. The equilibrium radial electric field (E{sub r}) reduces the transport. The effect of the parallel equilibrium current is found to be weak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

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↗

NASA NTRS1989

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

DOE OSTI2013

Physics basis of Multi-Mode anomalous transport module

The derivation of Multi-Mode anomalous transport module version 8.1 (MMM8.1) is presented. The MMM8.1 module is advanced, relative to MMM7.1, by the inclusion of peeling modes, dependence of turbulence correlation length on flow shear, electromagnetic effects in the toroidal momentum diffusivity, and the option to compute poloidal momentum diffusivity. The MMM8.1 model includes a model for ion temperature gradient, trapped electron, kinetic ballooning, peeling, collisionless and collision dominated magnetohydrodynamics modes as well as model for electron temperature gradient modes, and a model for drift resistive inertial ballooning modes. In the derivation of the MMM8.1 module, effects of collisions, fast ion and impurity dilution, non-circular flux surfaces, finite beta, and Shafranov shift are included. The MMM8.1 is used to compute thermal, particle, toroidal, and poloidal angular momentum transports. The fluid approach which underlies the derivation of MMM8.1 is expected to reliably predict, on an energy transport time scale, the evolution of temperature, density, and momentum profiles in plasma discharges for a wide range of plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

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↗