Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 847 records · Page 47

Tectonic plates, D (double prime) thermal structure, and the nature of mantle plumes

It is proposed that subducting tectonic plates can affect the nature of thermal mantle plumes by determining the temperature drop across a plume source layer. The temperature drop affects source layer stability and the morphology of plumes emitted from it. Numerical models are presented to demonstrate how introduction of platelike behavior in a convecting temperature dependent medium, driven by a combination of internal and basal heating, can increase the temperature drop across the lower boundary layer. The temperature drop increases dramatically following introduction of platelike behavior due to formation of a cold temperature inversion above the lower boundary layer. This thermal inversion, induced by deposition of upper boundary layer material to the system base, decays in time, but the temperature drop across the lower boundary layer always remains considerably higher than in models lacking platelike behavior. On the basis of model-inferred boundary layer temperature drops and previous studies of plume dynamics, we argue that generally accepted notions as to the nature of mantle plumes on Earth may hinge on the presence of plates. The implication for Mars and Venus, planets apparently lacking plate tectonics, is that mantle plumes of these planets may differ morphologically from those of Earth. A corollary model-based argument is that as a result of slab-induced thermal inversions above the core mantle boundary the lower most mantle may be subadiabatic, on average (in space and time), if major plate reorganization timescales are less than those acquired to diffuse newly deposited slab material.

Lenardic, A.↗

Optimal Inlet Shape Design of N2B Hybrid Wing Body Configuration

The N2B hybrid wing body aircraft was conceptually designed to meet environmental and performance goals for the N+2 generation transport set by the Subsonic Fixed Wing project of NASA Fundamental Aeronautics Program. In the present study, flow simulations are conducted around the N2B configuration by a Reynolds-averaged Navier-Stokes flow solver using unstructured meshes. Boundary conditions at engine fan face and nozzle exhaust planes are provided by the NPSS thermodynamic engine cycle model. The flow simulations reveal challenging design issues arising from boundary layer ingestion offset inlet and airframe-propulsion integration. Adjoint-based optimal designs are then conducted for the inlet shape to minimize the airframe drag force and flow distortion at fan faces. Design surfaces are parameterized by NURBS, and the cowl lip geometry is modified by a spring analogy approach. By the drag minimization design, flow separation on the cowl surfaces are almost removed, and shock wave strength got remarkably reduced. For the distortion minimization design, a circumferential distortion indicator DPCP(sub avg) is adopted as the design objective and diffuser bottom and side wall surfaces are perturbed for the design. The distortion minimization results in a 12.5 % reduction in the objective function.

Kim, Hyoungjin↗

Efficiency and reliability enhancements in propulsion flowfield modeling

The implementation of traditional CFD algorithms in practical propulsion related flowfields often leads to dramatic reductions in efficiency and/or robustness. The present research is directed at understanding the reasons for this deterioration and finding methods to circumvent it. Work to date has focussed on low Mach number regions, viscous dominated regions, and high grid aspect ratios. Time derivative preconditioning, improved definition of the local time stepping, and appropriate application of boundary conditions are employed to decrease the required time to obtain a solution, while maintaining accuracy. A number of cases having features typical of rocket engine flowfields are computed to demonstrate the improvement over conventional methods. These cases include laminar and turbulent high Reynolds number flat plate boundary layers, flow over a backward-facing step, a diffusion flame, and wall heat-flux calculations in a turbulent converging-diverging nozzle. Results from these cases show convergence that is virtually independent of the local Mach number and the grid aspect ratio, which translates to a convergence speed-up of up to several orders of magnitude over conventional algorithms. Current emphasis is in extending these results to three-dimensional flows with highly stretched grids.

Buelow, Philip E. O.↗

Conversion Kinetics and Ionic Conductivity in Na-β”-Alumina + YSZ (Naβ”AY) Sodium Solid Electrolyte via Vapor Phase Conversion Process

Sodium ion batteries have been receiving increasing attention and may see potential revival in the near future, particularly in large-scale grid energy storage coupling with wind and solar power generation, due to the abundant sodium resources, low cost, and sufficiently high energy density. Among the known sodium ion conductors, the Na-β”-alumina electrolyte remains highly attractive because of its high ionic conductivity. This study focuses on the vapor phase synthesis of a Na-β”-Alumina + YSZ (Naβ”AY) composite sodium electrolyte, which has higher mechanical strength and stability than conventional single phase β”-Alumina. The objectives are the measurement of conversion kinetics through a newly developed weight-gain based model and the determination of sodium ionic conductivity in the composite electrolyte. Starting samples contained ~70 vol% α-Alumina and ~30 vol% YSZ (3 mol% Y 2 O 3 stabilized Zirconia) with and without a thin alumina surface layer made by sintering in air at 1600 °C. The sintered samples were placed in a powder of Na-β”-alumina and heat-treated at 1250 °C for various periods. Sample dimensions and weight were measured as a function of heat treatment time. The conversion of α-Alumina in the α-Alumina + YSZ composite into Naβ”AY occurred by coupled diffusion of sodium ions through Na-β”-alumina and of oxygen ions through YSZ, effectively diffusing Na2O. From the analysis of the time dependence of sample mass and dimensions, the effective diffusion coefficient of Na 2 O through the sample, D eff , was estimated to be 1.74 x 10 -7 cm2 s -1 , and the effective interface transfer parameter, k eff , was estimated as 2.33 x 10 -6 cm s -1 . By depositing a thin alumina coating layer on top of the bulk composite, the chemical diffusion coefficient of oxygen through single phase Na-β”-alumina was estimated as 4.35 x 10 -10 cm 2 s -1 . An AC impedance measurement was performed on a fully converted Naβ”AY composite, and the conductivity of the composite electrolyte was 1.3 x 10 -1 S cm -1 at 300 °C and 1.6 x 10 -3 S cm -1 at 25 °C, indicating promising applications in solid state or molten salt batteries at low to intermediate temperatures.

36 MATERIALS SCIENCE↗

Progress toward the analysis of supersonic inlet flows

A zonal flow analysis procedure was developed to predict the flow through the supersonic diffuser of an axisymmetric mixed compression inlet at angle-of-attack. In this analysis, the inlet flow is divided into three types of regions, each with different dominant flow phenomena. These are the inviscid supersonic core, boundary layer, and shock/boundary layer interaction flows. An appropriate analysis was selected or developed for the three-dimensional flow in each type of region. Procedures were developed to interface these analyses for the overall inlet flow analysis. This analysis was applied to an inlet operating at M = 2.58 at several angle-of-attack conditions. Comparisons are presented between computed and measured flow properties for the inlet and for the component analysis flows. Extensions of the present procedure to include the terminal shock and subsonic diffuser flows are recommended. Desirable experiments for evaluation of the inlet analysis procedure or the component analyses and to support improved modeling or extension of the inlet analysis are defined and recommended.

Paynter, G. C.↗

Jupiter System Data Analysis Program: Mechanisms, Manifestation, and Implications of Cryomagmatism on Europa

The objectives of the work completed under NASA Grant NAG5-8898 were (i) to document and characterize the low-albedo diffuse surfaces associated with triple bands and lenticulae, (ii) to determine their mechanisms of formation, and (iii) to assess the implications of these features for the resurfacing (in space and time) of the Europa and the nature of the Europan interior. The approach involved a combination of processing and analysis of Solid State Imaging data returned by the Galileo spacecraft during the primary and extended mission phases, together with numerical modeling of the physical processes interpreted to the observed features. We have modeled the formation of Europan triple explosive venting of cryoclastic material from bands and lenticulae halos by two processes: (i) a liquid layer in the Europan interior, and (ii) lag deposit formation by the thermal influence of subsurface cryomagmatic intrusions. We favor the latter hypothesis for explaining these features, and further suggest that a liquid water or brine intrusion is required to provide sufficient lateral heating of surface ice to explain the 25 km size of the largest features. (Solid ice diapirs, even under the most favorable conditions, become thermally exhausted before they heat significant lateral distances). We argue that water circulating in open fractures, or repeated cryomagmatic 'diking' events would provide sufficient thermal input to produce the observed features. Thus our work argues for the existence of a liquid beneath Europa's surface. Our results might most easily be explained by the presence of a continuous liquid layer (the putative Europan ocean); this would concur with the findings of the Galileo magnetometer team. However, we cannot rule out the possibility that discrete liquid pockets provide injections of fluid closer to the surface.

Fagents, Sarah A.↗

NASA Tech Briefs, November 2003

Topics covered include: Computer Program Recognizes Patterns in Time-Series Data; Program for User-Friendly Management of Input and Output Data Sets; Noncoherent Tracking of a Source of a Data-Modulated Signal; Software for Acquiring Image Data for PIV; Detecting Edges in Images by Use of Fuzzy Reasoning; A Timer for Synchronous Digital Systems; Prototype Parts of a Digital Beam-Forming Wide-Band Receiver; High-Voltage Droplet Dispenser; Network Extender for MIL-STD-1553 Bus; MMIC HEMT Power Amplifier for 140 to 170 GHz; Piezoelectric Diffraction-Based Optical Switches; Numerical Modeling of Nanoelectronic Devices; Organizing Diverse, Distributed Project Information; Eigensolver for a Sparse, Large Hermitian Matrix; Modified Polar-Format Software for Processing SAR Data; e-Stars Template Builder; Software for Acoustic Rendering; Functionally Graded Nanophase Beryllium/Carbon Composites; Thin Thermal-Insulation Blankets for Very High Temperatures; Prolonging Microgravity on Parabolic Airplane Flights; Device for Locking a Control Knob; Cable-Dispensing Cart; Foam Sensor Structures Would be Self-Deployable and Survive Hard Landings; Real-Gas Effects on Binary Mixing Layers; Earth-Space Link Attenuation Estimation via Ground Radar Kdp; Wedge Heat-Flux Indicators for Flash Thermography; Measuring Diffusion of Liquids by Common-Path Interferometry; Zero-Shear, Low-Disturbance Optical Delay Line; Whispering-Gallery Mode-Locked Lasers; Spatial Light Modulators as Optical Crossbar Switches; Update on EMD and Hilbert-Spectra Analysis of Time Series; Quad-Tree Visual-Calculus Analysis of Satellite Coverage; Dyakonov-Perel Effect on Spin Dephasing in n-Type GaAs; Update on Area Production in Mixing of Supercritical Fluids; and Quasi-Sun-Pointing of Spacecraft Using Radiation Pressure.

Source record↗

Design of PGM-free cathodic catalyst layers for advanced PEM fuel cells

Here, the design of cathodic catalysts layer (CCL) consisted of Platinum Group Metal-free (PGM-free) electrocatalysts was done by catalyst coated membrane approach. Three different Fe-Mn-N-C compounds were synthesized with Fe:Mn ratio of 1:1, 2:1 and 2:1 with modified heat treatment profile. The catalysts were characterized by X-ray photoelectron spectroscopy, X-ray powder diffraction, pore and particle size distribution, zeta potential and transmission electron microscopy. Electrocatalysts were integrated into membrane electrode assembly and evaluated by electrochemical methods. Electrochemical impedance spectroscopy in combination with modeling were used for estimation of proton conductivity of CCL and its oxygen diffusivity. It was found that all CCLs possess extremely high proton conductivity, which was demonstrated for the first time for these types of PGM-free catalysts. The observed ORR mechanism was predominantly 4e- due to peroxide/radicals scavenging effect of Mn.

impedance↗

Stationary Oxygen Vacancy Construction toward a Superior-Performance Ultrahigh Nickel Single-Crystal Cathode

Oxygen vacancies exert a complex and profound influence on the layered cathodes, especially those with ultrahigh nickel content. They can facilitate lithium-ion transport and enhance electronic conductivity, while aggressive oxygen vacancy formation causes structural degradation and electrolyte decomposition. Herein, taking ultrahigh nickel single-crystal LiNi 0.92 Co 0.06 Mn 0.02 O 2 (SC-Ni92) as a model material, we propose a pinning strategy to harness the benefits of oxygen vacancies while mitigating their detrimental effects. Through a carefully controlled thermal process, both oxygen vacancies and pinning atoms are successfully introduced into the surface region. The resulting anchored oxygen vacancies, capitalizing on their inherent advantages, improve conductivity and lithium-ion diffusion. Simultaneously, the neighboring pinning atoms effectively increase the migration barrier and suppress the adverse effects of these vacancies, including electrolyte decomposition and structural degradation during long-term electrochemical cycling. Consequently, oxygen vacancy-anchored single-crystal LiNi 0.92 Co 0.06 Mn 0.02 O 2 (SC-Ni92-OV) demonstrates significantly improved high-voltage electrochemical performance, with 86.16% capacity retention after 200 cycles at 4.6 V and 1 C in a half-cell and 90.71% after 300 cycles at 4.5 V and 1 C in a full cell. Furthermore, this study not only provides valuable insights into the chemistry of oxygen vacancy but also introduces a viable strategy for leveraging oxygen vacancies to achieve stable high-voltage performance in ultrahigh nickel single-crystal cathodes.

defects in solids↗

Gravitational Influences on Flame Propagation Through Non-Uniform Premixed Gas Systems

Flame propagation through non-uniformly premixed gases occurs in several common combustion situations. As summarized in a previous conference paper, non-uniform premixed gas combustion has received scant attention compared to the more usual limiting cases of diffusion or uniformly premixed flames. It is the goal of this research to further our knowledge of layered combustion, in which a fuel concentration gradient exists normal to the direction of flame spread, in particular by focusing on the role that gravity plays. Gravity can affect flame propagation in at least three ways: through a hydrostatic pressure gradient, by altering the initial distribution of fuel vapor, and through buoyantly induced flows once ignition has occurred. An understanding of the phenomena involved is important to fire safety, especially aboard spacecraft since no microgravity data exist. The data obtained will also be useful to verify theoretical models of this problem, which are easier to implement if buoyancy is neglected.

Miller, Fletcher J.↗

Ionomer Optimization for Hydroxide-Exchange-Membrane Water Electrolyzers Operated with Distilled Water: A Modeling Study

The hydroxide-exchange-membrane water electrolyzer (HEMWE) is a promising means to store intermittent renewable energy in the form of hydrogen chemical energy. The hydroxide-exchange ionomer (HEI) in the gas-evolving electrodes and the hydroxide-exchange membrane (HEM) are key components of HEMWE. In this work, we simulate the cell and examine explicitly the impact of HEI and HEM properties with a focus on improving HEMWE performance when operated with distilled water ( i.e ., no supporting electrolyte). The tradeoff between the ionic conductivity gain and electrochemically active surface area (ECSA) loss is studied. For a constant catalyst loading, distributing more catalyst next to the HEM or making thinner but denser catalyst layer is beneficial for HEMWE performance. The results demonstrate that a higher water diffusion coefficient is desired for HEM to supply reactant water to the cathode. In contrast, a lower water diffusion coefficient is preferred for the cathode HEI to retain the water in the regions with high reaction rates. Overall, the findings provide important insights to optimizing HEI/HEM materials for improved HEMWE performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensing depths in frequency domain thermoreflectance

In this work, a method is developed to calculate the length into a sample to which a Frequency Domain Thermoreflectance (FDTR) measurement is sensitive. Sensing depth and sensing radius are defined as limiting cases for the spherically spreading FDTR measurement. A finite element model for FDTR measurements is developed in COMSOL multiphysics and used to calculate sensing depth and sensing radius for silicon and silicon dioxide samples for a variety of frequencies and laser spot sizes. The model is compared to experimental FDTR measurements. Design recommendations for sample thickness are made for experiments where semi-infinite sample depth is desirable. For measurements using a metal transducer layer, the recommended sample thickness is three thermal penetration depths, as calculated from the lowest measurement frequency.

36 MATERIALS SCIENCE↗

Transient heat transfer of superfluid He 4 in nonhomogeneous geometries: Second sound, rarefaction, and thermal layer

Transient heat transfer in superfluid 4 He (He II) is a complex process that involves the interplay of the unique counterflow heat-transfer mode, the emission of second-sound waves, and the creation of quantized vortices. Many past researches focused on homogeneous heat transfer of He II in a uniform channel driven by a planar heater. Here, we report our systematic study of He II transient heat transfer in nonhomogeneous geometries that are pertinent to emergent applications. By solving the He II two-fluid equations of motion coupled with Vinen's equation for vortex-line density, we examine and compare the characteristics of transient heat transfer from planar, cylindrical, and spherical heaters in He II. Our results show that as the heater turns on, an outgoing second-sound pulse emerges, within which the vortex-line density grows rapidly. These vortices attenuate the second sound and result in a heated He II layer in front of the heater, i.e., the thermal layer. In the planar case where the vortices are created throughout the space, the second-sound pulse is continuously attenuated, leading to a thick thermal layer that diffusely spreads following the heat pulse. On the contrary, in the cylindrical and the spherical heater cases, vortices are created mainly in a thin thermal layer near the heater surface. As the heat pulse ends, a rarefaction tail develops following the second-sound pulse, in which the temperature drops. This rarefaction tail can promptly suppress the thermal layer and take away the deposited thermal energy. The effects of the heater size, heat flux, pulse duration, and temperature on the thermal-layer dynamics are discussed. We also show how the peak heat flux for the onset of boiling in He II can be studied in our model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave-particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel and Engelmann (Phys. Fluids vol. 9, 1966, pp. 2377- 2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that pitch angle scatter them out of and into resonance. As a result, the resonant particle-wave interactions occur in the center of a narrow collisional boundary when the collision frequency nu is very small compared to the wave frequency omega. The diffusive nature of the pitch angle scattering combined with the wave-particle resonance condition enhances the collision frequency by (omega/nu)2/3 >>1, resulting in an effective resonant particle collision time of tau_int ~ (nu /omega)2/3 nu <<1/ nu . A rigorous collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel-Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Decay of the zincate concentration gradient at an alkaline zinc cathode after charging

The transport of the zincate ion to the alkaline zinc cathode was studied by observing the decay of the zincate concentration gradient at a horizontal zinc cathode after charging. This decay was found to approximate first order kinetics as expected from a proposed boundary layer model. The concentrations were calculated from polarization voltages. The decay half life was shown to be a linear function of the thickness of porous zinc deposit on the cathode indicating a very rapid transport of zincate through porous zinc metal. The rapid transport is attributed to an electrochemical mechanism. From the linear dependence of the half life on the thickness the boundary layer thickness was found to be about 0.010 cm when the cathode was at the bottom of the cell. No significant dependence of the boundary layer thickness on the viscosity of electrolyte was observed. The data also indicated a relatively sharp transition between the diffusion and convection transport regions. When the cathode was at the top of the cell, the boundary layer thickness was found to be roughly 0.080 cm. The diffusion of zincate ion through asbestos submerged in alkaline electrolyte was shown to be comparable with that predicted from the bulk diffusion coefficient of the zincate ion in alkali.

Kautz, H. E.↗

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗