Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusers”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Diffusion of chromium and aluminum in Ni-20Cr and TDNiCr /Ni-20Cr-2ThO2/.

Diffusion coefficients have been measured for Cr51 in fine- and coarse-grained TDNiCr (Ni-20Cr-2ThO2) and in fine-grained Ni-20Cr in the temperature range from 1038 to 1200 C. Selective diffusivities have also been determined for specimens of these alloys which were aluminized to give an initial surface concentration of 5.8 wt % Al. Finally, diffusion coefficients for interdiffusion of aluminum in TDNiCr and Ni-20Cr have been obtained from electron probe microanalysis of the aluminized specimens. For a given grain size and temperature there is no difference in diffusivities for chromium diffusion in TDNiCr or Ni-20Cr. Diffusion coefficients increase with decreasing grain size for both alloys. Comparison of aluminum diffusion data obtained from electron microprobe profiles with radiotracer chromium diffusivities suggests that aluminum diffuses approximately three times faster than chromium in TDNiCr and Ni-20Cr.

Seltzer, M. S.↗
NASA NTRS1973

Directional diffusion coefficients of solar protons inside and outside the bow shock.

The directional diffusion coefficients of low-energy (greater than or equal to 0.3 MeV) solar protons inside and outside the bow shock are examined during the solar flare event of Jan. 24, 1969. The data are derived from simultaneous observations obtained by Explorer 33 inside the magnetosheath and by Explorer 35 in the interplanetary medium. Although the gross properties of the spin-averaged intensities on a diffusion-type plot appear to be the same in both media, the directional intensities show significant variations. It is shown that directional intensities of low-energy protons can be described reasonably well by anisotropic diffusion with an associated diffusion coefficient. Directional diffusion coefficients are found to differ by a factor of as much as three among different directions in space, and from the spin-averaged diffusion coefficient. This suggests that anisotropic diffusion does indeed take place and that so called 'isotropic' diffusion coefficients derived in the past from spin-averaged intensities may actually be directional diffusion coefficients in cases where substantial anisotropies (greater than 50%) exist.

Verzariu, P.↗
NASA NTRS1978

Horizontal eddy diffusivities in the northern hemispheric stratosphere and lower mesosphere

Meteorological rocket soundings, representative of low, middle, and high altitudes, are employed with the use of the statistical theory of diffusion, to determine the zonal and meridional component of eddy diffusivity between 30 and 55 km as a function of season, latitude, and altitude. A comparison is also made between annualry-averaged eddy diffusivities above and below 30 km. It is shown that the zonal component of eddy diffusivity is approximately three to five times as large as the meridional component, in most cases. Both components of eddy diffusivity varied greatly with season, latitude, and altitude. Highest eddy diffusivities occurred in the vicinity of the winter westerly jet located near the tropopause and in the stratosphere at 35 km, 40 deg N. Annually, a minimum of eddy diffusivity occurred near 20 km. In most cases, the zonal component of the eddy diffusivity was approximately three to five times as large as the meridional component. The meridional cross-sections of the zonal and meridional components of the eddy diffusivity in the northern hemisphere are presented.

Kao, S. K.↗
NASA NTRS1978

Effect of concentration dependence of the diffusion coefficient on homogenization kinetics in multiphase binary alloy systems

Diffusion calculations were performed to establish the conditions under which concentration dependence of the diffusion coefficient was important in single, two, and three phase binary alloy systems. Finite-difference solutions were obtained for each type of system using diffusion coefficient variations typical of those observed in real alloy systems. Solutions were also obtained using average diffusion coefficients determined by taking a logarithmic average of each diffusion coefficient variation considered. The constant diffusion coefficient solutions were used as reference in assessing diffusion coefficient variation effects. Calculations were performed for planar, cylindrical, and spherical geometries in order to compare the effect of diffusion coefficient variations with the effect of interface geometries. In most of the cases considered, the diffusion coefficient of the major-alloy phase was the key parameter that controlled the kinetics of interdiffusion.

Tenney, D. R.↗
NASA NTRS1979

Experiments on tandem diffusers with boundary-layer suction applied in between

Experiments were performed on conical diffusers of various configurations with the same, but rather unusually large, 16:1 area ratio. Because available performance data on diffusers fall short of very large area ratio configurations, an unconventional design, consisting of two diffusers following each other in tandem, was proposed. Both diffusers had the same area ratio of 4:1, but had different taper angles. While for the first diffuser (called leading) the angle remained constant, for the second (called follower), the taper angle was stepped up to higher values. Boundary layer control, by way of suction, was applied between the diffusers, and a single slot suction ring was inserted between them. The leading diffuser had an enclosed nominal divergence angle 2 theta = 5 degrees, while the follower diffusers had either 10, 20, 30, or 40 degrees, respectively, giving 4 combinations. The experiments were performed at four different Reynolds numbers with various suction rates. The rates indicate a general improvement in the performance of all diffusers with boundary layer suction. It appears that the improvement of the pressure recovery depends on both the Reynolds number and the suction rate, and the largest increase, 0.075, was found at the lowest R sub e when the follower divergence was 2 theta = 40 degrees.

Barna, P. S.↗
NASA NTRS1979

Characteristics of acoustic-gravity waves in a diffusively separated atmosphere

The effects of wave-induced diffusion on the characteristics of plane acoustic gravity waves in the earth's thermosphere are investigated by means of a two fluid model atmosphere in diffusive equilibrium. At wave periods greater than the shortest characteristic diffusion time for the minor gas, diffusion induced by the wave is considered to act to eliminate amplitude and phase differences between the wave-induced density fluctuations of individual gases and restore the perturbed fluid to diffusive equilibrium. Vertical diffusion is most important at large scale sizes, but horizontal diffusion dominates for horizontal wavelengths of less than several hundred kilometers. As a result of wave-induced diffusion, AE-C satellite measurements of neutral density fluctuations of thermospheric constituents are only compatible with relatively small scale internal gravity waves. Velocity and temperature differences for acoustic-gravity waves are greatest at periods near the mean collision time and diffusion time and lead to the dissipation of wave energy. The predictions of the theory are consistent with a high latitude, low altitude source for the medium scale waves and a more localized random source for the waves observed by AE-C.

Del Genio, A. D.↗
NASA NTRS1984

Ambipolar diffusion in the middle atmosphere

In the middle atmosphere above 60 km, the electron concentration increases with altitude, reaching values of 10 to the 10th per cu m in the daytime ionospheric E region near 100 km. The electrons are more mobile than the ions and diffuse more rapidly through the neutral atmosphere. The electron diffusion polarizes the medium, causing an electric field to develop that acts to retard the electron diffusion and enhance the conduction current of ions. A global zonally averaged numerical model of atmosheric electricity from the ground to 100 km is used to examine the effect of ambipolar diffusion and the earth's geomagnetic field on the currents and fields in the middle atmosphere. The results show that above about 65 km, ambipolar diffusion generates local electric fields and conduction currents that balance electron diffusion currents. The electric fields and conduction currents are a few orders of magnitude larger than the vertical fields and currents calculated from the downward mapping of the ionospheric potential without taking electron diffusion into account. Ambipolar diffusion does not alter the total current flowing in the global circuit. It is a local effect where enhanced conduction currens flow to balance the electron diffusion current.

Tzur, I.↗
NASA NTRS1991

Diffusion of a multi-species component and its role in oxygen and water transport in silicates

The diffusion of a multispecies component is complicated by the different diffusion coefficient of each species and the interconversion reactions among the species. A diffusion equation is derived that incorporates the diffusive fluxes of all species contributing to the component's concentration. The effect of speciation on diffusion is investigated experimentally by measuring concentration profiles of all species developed during diffusion experiments. Data on water diffusion in rhyolitic glasses indicate that H2O molecules predominate over OH groups as the diffusing species at very low to high water concentrations. A simple theoretical relationship is drawn between the effective total oxygen diffusion coefficient and the total water concentration of silicates at low water content.

Zhang, Youxue↗
NASA NTRS1999

Investigation of Diffusion Flame Tip Thermodiffusive and Hydrodynamic Instability under Microgravity Conditions

We employ the opposed flow flame-spread configuration in order to examine flame-front instability of diffusion flames near cold, solid boundaries. The thermo-diffusive and hydrodynamic instabilities can transform an initially planar flame front into an irregularly curved, corrugated, possibly fragmented front. Under ordinary 1-g conditions, the buoyancy-induced flow masks the thermo-diffusive and hydrodynamic instabilities and produces planar flames. Such stable spreading flames have been observed for decades in laboratory experiments. Experiments in zero gravity are necessary to produce unstable flame fronts. The thermo-diffusive/hydrodynamic microgravity instability appears in diffusion flames such as, for example: the candle flame oscillations observed by Dietrich et al.; smolder instabilities on a recent Space Shuttle flight. Drs. T. Kashiwagi and S. Olson have attributed the latter to a lowered oxygen transport rate to the hot, reactive surface. Consider a burning surface near the flame extinction limit. The flow, or stretch, induced by the diffusion flame is weak, hence buoyancy plays a small role, thereby enabling previously secondary mechanisms, such as differential thermo-diffusion, to become the most important mechanisms. The flame leading edge becomes unstable; and diffusion flame breakup, oscillation, and rejoining all occur at a measurable frequency of approximately O(1 Hz). This project has only begun in January of this year, 1999. To date, there have been no flight experiments on flame spread instabilities. However, we have made numerous experiments in the NASA 2.2 and 5 second drop towers on flame spread over very thin cellulosic fuels. We have been very fortunate through a combination of factors, to be explained, to obtain some interesting, perhaps even compelling, results on diffusion flame instability in the presence of heat losses to cold surfaces.

Wichman, I. S.↗
NASA NTRS1998

Temperature Dependence of Diffusivities in Liquid Elements (LMD)

This research was to advance the understanding of diffusion mechanisms in liquid metals and alloys through accurate diffusivity measurements over a wide range of temperatures, including the proximity of the materials melting points. Specifically, it was driven towards developing a methodology (and subsequent flight hardware) to enable several diffusion coefficient measurements (i.e., at several different temperatures) to be performed using a single sample. The Liquid Metal Diffusion (LMD) was funded as a Flight Definition Project in February 1993 in response to NRA 91-OSSA-20 (Microgravity Science and Applications Division). The Science Concept Review for LAID was held during April 1994. In January 1995 we were informed that we had failed this review and the project was change to ground-based activities only. A new proposal was submitted for the next NRA addressing the panels concerns. As part of NASA's Risk Mitigation program, a scaled-down version of the hardware was funded in July of 1995 for a flight opportunity utilizing experiment on the Microgravity Isolation Mount. This experiment was to determine the self-diffusivity of indium at 185 C. The LMD was transferred to the Mir Space Station in STS-81 and returned on STS-84 (January - May 1997). Three, out of five, self-diffusion data sets were returned. A description of this experiment/hardware is included below. This summary is only intended to give the reader an overview of the results obtained for the tasks outlined in the original proposal. Research that was not published is explained in more detail. At the end of this report is a list of refereed publications and invited talks that were given as a result of this work. The reader is directed to these for further details. Attachment: Real-time diffusivity measurements in liquids at several temperatures with one sample, On the insensitivity of liquid diffusivity measurements to deviations from 1D transport, and Numerical simulations of the convective contamination of diffusivity measurements in liquids.

Banish, R. Michael↗
NASA NTRS2002

Flow/Soot-Formation Interactions in Nonbuoyant Laminar Diffusion Flames

This is the final report of a research program considering interactions between flow and soot properties within laminar diffusion flames. Laminar diffusion flames were considered because they provide model flame systems that are far more tractable for theoretical and experimental studies than more practical turbulent diffusion flames. In particular, understanding the transport and chemical reaction processes of laminar flames is a necessary precursor to understanding these processes in practical turbulent flames and many aspects of laminar diffusion flames have direct relevance to turbulent diffusion flames through application of the widely recognized laminar flamelet concept of turbulent diffusion flames. The investigation was divided into three phases, considering the shapes of nonbuoyant round laminar jet diffusion flames in still air, the shapes of nonbuoyant round laminar jet diffusion flames in coflowing air, and the hydrodynamic suppression of soot formation in laminar diffusion flames.

Dai, Z.↗
NASA NTRS2007

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm) (Fig. 1, left). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated (Fig. 1, right). A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗
NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗
NASA NTRS2006

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated. A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗
DOE OSTI2019

A Prediction of Uranium(VI) Diffusion in Montmorillonite at Various Chemical Solution Conditions

Most nuclear waste disposal options include compacted bentonite, predominately consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and a potential contaminant of water resources. Uranium sorption onto clay and its slow diffusive transport away from waste canisters are expected to limit U(VI) mobility, and possibly also control nuclear fuel degradation rates based on mass action considerations. A prediction of U(VI) diffusion in these systems is complicated by two main factors: (1) the complex mineralogical structure of montmorillonite clay, leading to two types of clay porosities and surface sites, and (2) the dependence of U(VI) sorption on U(VI) solution speciation and the specific chemical system conditions. For instance, depending on pH and ionic strength, U(VI) may sorb onto montmorillonite due to weak ion exchange reactions or the formation of stable surface complexes. Furthermore, in contrast to conventional porous media, U(VI) diffusion in montmorillonite is not only limited by contaminant sorption, but also by a (potential) lack of access to the full clay porosity for some U(VI) solution species. For instance, at the high degrees of clay compaction expected in future barrier systems, the partial or full exclusion of U(VI) from negatively-charged clay interlayer spaces may substantially decrease U(VI) fluxes. In this modeling study, we specifically investigated how the charge of U(VI) solution species affects U(VI) sorption and diffusion behavior. For this purpose, we first simplified complex chemical solution speciation diagrams to charge speciation diagrams. Then, we coupled an existing U(VI) surface complexation model (Tournassat et al., Geochim. Cosmochim. Acta, 2018) with an analytical solution of the diffusion equation to simulate U(VI) mobility in lab-scale through-diffusion experiments for various chemical conditions. Our results suggest that there are two important parameters that govern U(VI) diffusion in these systems. First, U(VI) surface complexation on clay edge surfaces leads to contaminant retardation. On the field scale, this will determine how early a plume may cross site boundaries and pose risks to drinking water resources. Second, the degree of access of U(VI) solution species to diffuse layers in clay interlayer spaces largely controls U(VI) diffusive fluxes. In the field, this parameter will govern the level of U(VI) concentrations that can be expected in contaminated waters.

Pistorino, Jonathan C.↗
DOE OSTI2018

Neoclassical Diffusion of Radiation-Belt Electrons Across Very Low L -shells

In the presence of drift-shell splitting intrinsic to the IGRF magnetic field model, pitch-angle scattering from Coulomb collisions experienced by radiation-belt electrons in the upper atmosphere and ionosphere produces extra radial diffusion, a form of neoclassical diffusion. The strength of the neoclassical radial diffusion at L < 1.2 exceeds that expected there from radial-diffusion mechanisms traditionally considered, and decreases with increasing L-shell. In this study we construct a numerical model for this coupled (radial and pitch-angle) collisional diffusion process and apply it to simulate raw count-rate data observed aboard the Gemini spacecraft for several years after the 1962 Starfish nuclear detonation. The data show apparent lifetimes 10-100 times as long as would have been expected from collisional pitch-angle diffusion and Coulomb drag alone. Our model reproduces apparent lifetimes for >0.5-MeV electrons in the region 1.14 < L < 1.26 to within a factor of two (comparable to the uncertainty quoted for the observations). We conclude that neoclassical radial diffusion (resulting from drift-shell splitting intrinsic to IGRF's azimuthal asymmetries) mitigates the decay expected from collisional pitch-angle diffusion and inelastic energy loss alone and thus contributes importantly to the long apparent lifetimes observed at these low L-shells.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Probing the Impact of Vacancy Diffusion on Void Dynamics at the Lithium Metal–Solid Electrolyte Interface

Lithium (Li) metal-based solid-state batteries (SSBs) are considered promising candidates for next-generation energy storage due to their superior energy density and enhanced safety compared to conventional Li-ion systems. However, their practical application is limited by challenges such as void formation at the Li-solid electrolyte (SE) interface, which disrupts ion transport and accelerates interfacial degradation. This work investigates how the coupled effects of electro-dissolution kinetics and surface diffusion at the Li metal surface govern the evolution of interfacial morphology during stripping. This work examines the influence of three distinct surface diffusion modes, which are terrace diffusion, step diffusion, and interlayer diffusion, on maintaining interfacial stability. In addition, how the dominant surface diffusion mechanism can overcome the contact loss due to high reaction kinetics is explored. Furthermore, the roughness of the Li metal anode surface is quantified, and the influence of different diffusion mechanisms on the evolution of the dynamic solid–solid interface is examined. The critical role of temperature in enhancing Li surface diffusivity and expanding the regime of stable contact is highlighted. By identifying distinct regimes of interface stability, this study analyzes how non-uniform electrochemical dynamics dictate void morphology evolution and interfacial contact. These insights offer guiding principles for engineering robust Li–SE interfaces in SSBs.

lithium metal anode↗
DOE OSTI2014

Finite-difference schemes for anisotropic diffusion

In fusion plasmas diffusion tensors are extremely anisotropic due to the high temperature and large magnetic field strength. This causes diffusion, heat conduction, and viscous momentum loss, to effectively be aligned with the magnetic field lines. This alignment leads to different values for the respective diffusive coefficients in the magnetic field direction and in the perpendicular direction, to the extent that heat diffusion coefficients can be up to 10{sup 12} times larger in the parallel direction than in the perpendicular direction. This anisotropy puts stringent requirements on the numerical methods used to approximate the MHD-equations since any misalignment of the grid may cause the perpendicular diffusion to be polluted by the numerical error in approximating the parallel diffusion. Currently the common approach is to apply magnetic field-aligned coordinates, an approach that automatically takes care of the directionality of the diffusive coefficients. This approach runs into problems at x-points and at points where there is magnetic re-connection, since this causes local non-alignment. It is therefore useful to consider numerical schemes that are tolerant to the misalignment of the grid with the magnetic field lines, both to improve existing methods and to help open the possibility of applying regular non-aligned grids. To investigate this, in this paper several discretization schemes are developed and applied to the anisotropic heat diffusion equation on a non-aligned grid.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗