Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion Coefficients”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
NASA NTRS2002

Seasonal to Decadal-Scale Variability in Satellite Ocean Color and Sea Surface Temperature for the California Current System

Support for this project was used to develop satellite ocean color and temperature indices (SOCTI) for the California Current System (CCS) using the historic record of CZCS West Coast Time Series (WCTS), OCTS, WiFS and AVHRR SST. The ocean color satellite data have been evaluated in relation to CalCOFI data sets for chlorophyll (CZCS) and ocean spectral reflectance and chlorophyll OCTS and SeaWiFS. New algorithms for the three missions have been implemented based on in-water algorithm data sets, or in the case of CZCS, by comparing retrieved pigments with ship-based observations. New algorithms for absorption coefficients, diffuse attenuation coefficients and primary production have also been evaluated. Satellite retrievals are being evaluated based on our large data set of pigments and optics from CalCOFI.

Mitchell, B. Greg↗
NASA NTRS2002

Many-body Effects in a Laterally Inhomogeneous Semiconductor Quantum Well

Many body effects on conduction and diffusion of electrons and holes in a semiconductor quantum well are studied using a microscopic theory. The roles played by the screened Hartree-Fock (SHE) terms and the scattering terms are examined. It is found that the electron and hole conductivities depend only on the scattering terms, while the two-component electron-hole diffusion coefficients depend on both the SHE part and the scattering part. We show that, in the limit of the ambipolax diffusion approximation, however, the diffusion coefficients for carrier density and temperature are independent of electron-hole scattering. In particular, we found that the SHE terms lead to a reduction of density-diffusion coefficients and an increase in temperature-diffusion coefficients. Such a reduction or increase is explained in terms of a density-and temperature dependent energy landscape created by the bandgap renormalization.

Ning, Cun-Zheng↗
NASA NTRS1995

Absorption of water and lubricating oils into porous nylon

Oil and water absorption from air into sintered porous nylon can be described by infiltration into the pores of the material. This process can be modeled by a diffusion-like mechanism. For water absorption, we find a formal diffusion coefficient of 1.5 x 10(exp -4)sq cm/min when the nylon is initially dry. The diffusion coefficient is 4 x 10(exp -6)sq cm/min when the nylon is oil-impregnated prior to air exposure. In a 52% RH atmosphere, dry nylon absorbs 3% w/w water, and oil-impregnated nylon absorbs 0.6% w/w water. For oil absorption there are three steps: (1) surface absorption and infiltration into (2) larger and (3) smaller pores. Surface absorption is too fast to be measured in these experiments. The diffusion coefficient for the second step is 6 x 10(exp -4)sq cm/min for SRG-60 oil into dry nylon and 4 x 10(exp -4)sq cm/min for air-equilibrated nylon. The diffusion coefficient for the third step is about 1 x 10(exp -6)sq cm/min for both cases. The total amount of oil absorbed is 31% w/w. The interaction between water and nylon is not as strong as that between water and cotton-phenolic: oil can replace water, and only a small amount of water can enter previously oil-impregnated nylon.

Bertrand, P. A.↗
NASA NTRS2001

A Hydrodynamic Theory for Spatially Inhomogeneous Semiconductor Lasers: Numerical Results - 2

We present numerical results of the diffusion coefficients (DCs) in the coupled diffusion model derived in the preceding paper for a semiconductor quantum well. These include self and mutual DCs in the general two-component case, as well as density- and temperature-related DCs under the single-component approximation. The results are analyzed from the viewpoint of free Fermi gas theory with many-body effects incorporated. We discuss in detail the dependence of these DCs on densities and temperatures in order to identify different roles played by the free carrier contributions including carrier statistics and carrier-LO phonon scattering, and many-body corrections including bandgap renormalization and electron-hole (e-h) scattering. In the general two-component case, it is found that the self- and mutual- diffusion coefficients are determined mainly by the free carrier contributions, but with significant many-body corrections near the critical density. Carrier-LO phonon scattering is dominant at low density, but e-h scattering becomes important in determining their density dependence above the critical electron density. In the single-component case, it is found that many-body effects suppress the density coefficients but enhance the temperature coefficients. The modification is of the order of 10% and reaches a maximum of over 20% for the density coefficients. Overall, temperature elevation enhances the diffusive capability or DCs of carriers linearly, and such an enhancement grows with density. Finally, the complete dataset of various DCs as functions of carrier densities and temperatures provides necessary ingredients for future applications of the model to various spatially inhomogeneous optoelectronic devices.

Li, Jianzhong↗
NASA NTRS1985

Middle atmosphere modeling

Breaking gravity waves generate and maintain a background level of turbulence which is capable of producing substantial cooling and/or heating in the upper mesosphere and lower thermosphere. The net thermodynamic effect of breaking gravity waves is critically dependent on the eddy Prandt number (P sub t) applicable to mesospheric turbulence. When P sub t is approximately 1, the calculations of the heat budget for the mesopause region imply that the globally averaged eddy or turbulent diffusion coefficient cannot exceed .000001 sq cm/s. This upper limit on turbulant diffusion applies to both potential temperature transport and chemically inert tracer transport when radiative damping is neglible. For chemically active species larger diffusion coefficients are permitted, because the effective eddy diffusion coefficient is increased by an additive term L/2 gamma (sup 2), where L is the chemical loss rate and gamma is the vertical wavenumber. For P sub t is approximately 4 to 6, the turbulent diffusion of momentum (D sub M) is sufficiently greater than the turbulent diffusion of heat (D sub H) that the conversion of gravity wave energy to heat with high efficiency nearly balances the divergence of the downward eddy heat flux in the wave breaking zone. Therefore the heat budget of the mesopause region would no longer provide a powerful and useful constraint on D sub H. If P sub t exceeds 6 with high efficiency for energy conversion to heat, gravity waves would heat the mesosphere throughout the wave breaking region.

Strobel, D. F.↗
NASA NTRS2000

Photochemistry of Saturn's Atmosphere: Hydrocarbon Chemistry and Comparisons with ISO Observations - 1

To investigate the details of hydrocarbon photochemistry on Saturn, we have developed a one-dimensional diurnally averaged model that couples hydrocarbon and oxygen photochemistry, molecular and eddy diffusion, radiative transfer, and condensation. The model results are compared with observations from the Infrared Space Observatory (ISO) to place tighter constraints on molecular abundances, to better define Saturn's eddy diffusion coefficient profile, and to identify important chemical schemes that control the abundances of the observable hydrocarbons in Saturn's upper atmosphere. From the ISO observations, we determine that the column 12 densities of CH3, CH3C2H, and C4H2 above 10 mbar are 4 (sup +2) (sub -1.5) x 10 (exp 13) cm (sup -2), (1.1 plus or minus 0.3) x 10 (exp 15) cm (exp -2), and (1.2 plus or minus 0.3) x 10 (exp 14) cm (sup -2), respectively. The observed ISO emission features also indicate C2H2 mixing ratios of 1.2 (sup +0.9) (sub -0.6) x 10 (exp -6) at 0.3 mbar and (2.7 plus or minus 0.8) x 10 (exp -7) at 1.4 mbar, and a C2H6 mixing ratio of (9 plus or minus 2.5) x 10 (exp -6) at 0.5 mbar. Upper limits are provided for C2H4, CH2CCH2, C3H8, and C6H2 sensitivity of the model results to variations in the eddy diffusion coefficient profile, the solar flux, the CH4 photolysis branching ratios, the atomic hydrogen influx, and key reaction rates are discussed in detail. We find that C4H2 and CH3C2H are particularly good tracers of important chemical processes and physical conditions in Saturn's upper atmosphere, and C2H6 is a good tracer of the eddy diffusion coefficient in Saturn's lower stratosphere. The eddy diffusion coefficient must be smaller than approximately 3 x 10 (exp 4) sq cm s (sup -1) at pressures greater than 1 mbar in order to reproduce the C2H6 abundance inferred from ISO observations. The eddy diffusion coefficients in the upper stratosphere could be constrained by observations of CH3 radicals if the low-temperature chemistry of CH3 were better understood. We also discuss the implications of our modeling for aerosol formation in Saturn's lower stratosphere-diacetylene, butane, and water condense between approximately 1 and 300 mbar in our model and will dominate stratospheric haze formation at nonauroral latitudes. Our photochemical models will be useful for planning observational sequences and for analyzing data from the upcoming Cassini mission.

Moses, Julianne I.↗
NASA NTRS2023

Analyses of ∼0.05-2 Mev Ions Associated With the 2022 February 16 Esp Event Observed by Parker Solar Probe

We present analyses of 0.05–2 MeV ions from the 2022 February 16 energetic storm particle event observed by Parker Solar Probe's (PSP) IS⊙IS/EPI-Lo instrument at 0.35 au from the Sun. This event was characterized by an enhancement in ion fluxes from a quiet background, increasing gradually with time with a nearly flat spectrum, rising sharply near the arrival of the coronal mass ejection (CME) driven shock, becoming nearly a power-law spectrum, then decaying exponentially afterwards, with a rate that was independent of energy. From the observed fluxes, we determine diffusion coefficients, finding that far upstream of the shock the diffusion coefficients are nearly independent of energy with a value of 10 20 cm 2 /s. Near the shock, the diffusion coefficients are more than one order of magnitude smaller and increase nearly linearly with energy. We also determine the source of energetic particles, by comparing ratios of the intensities at the shock to estimates of the quiet time intensity to predictions from diffusive shock acceleration theory. We conclude that the source of energetic ions is mostly the solar wind for this event. We also present potential interpretations of the near-exponential decay of the intensity behind the shock. One is that we suggest the shock was over-expanding when it crossed PSP and the energetic particle intensity decreased behind the shock to fill the expanding volume. Over-expanding CMEs could well be more common closer to the Sun, and this is an example of such a case.

Joe Giacalone↗
NASA NTRS2023

Analyses of ∼0.05–2 MeV Ions Associated with the 2022 February 16 Energetic Storm Particle Event Observed By Parker Solar Probe

We present analyses of 0.05–2 MeV ions from the 2022 February 16 energetic storm particle event observed by Parker Solar Probeʼs (PSP) ISeIS/EPI-Lo instrument at 0.35 au from the Sun. This event was characterized by an enhancement in ion fluxes from a quiet background, increasing gradually with time with a nearly flat spectrum, rising sharply near the arrival of the coronal mass ejection (CME)–driven shock, becoming nearly a power-law spectrum, then decaying exponentially afterward, with a rate that was independent of energy. From the observed fluxes, we determine diffusion coefficients, finding that far upstream of the shock the diffusion coefficients are nearly independent of energy, with a value of 10 20 cm 2 s −1 . Near the shock, the diffusion coefficients are more than 1 order of magnitude smaller and increase nearly linearly with energy. We also determine the source of energetic particles, by comparing ratios of the intensities at the shock to estimates of the quiet-time intensity to predictions from diffusive shock acceleration theory. We conclude that the source of energetic ions is mostly the solar wind for this event. We also present potential interpretations of the near-exponential decay of the intensity behind the shock. One possibility we suggest is that the shock was over expanding when it crossed PSP and the energetic particle intensity decreased behind the shock to fill the expanding volume. Over expanding CMEs could well be more common closer to the Sun, and this is an example of such a case.

Interplanetary particle acceleration↗
NASA NTRS2002

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗
NASA NTRS1986

Oxidation of commercial purity titanium

The oxidation kinetics of commercial purity Ti-A55 exposed to laboratory air in the 593-760 C temperature range were continuously monitored by thermogravimetric analysis. The oxide thickness was measured by microscopy, and the substrate contamination was estimated from microhardness measurements. The microhardness depth profiles were converted to oxygen composition profiles using calibration depth. The oxygen diffusion coefficient in alpha-Ti appears to be approximately concentration-independent in the 1-10 at. pct oxygen range. Diffusion coefficient for oxygen in TiO2 has been estimated as a function of temperature and is found to be about 50 times the value in alpha-Ti. The metallographically prepared cross sections of the oxidized specimens revealed a 'moving boundary' in the substrate, parallel to the oxide-metal interface. This boundary was associated with a specific oxygen level of 5.0 + or - 0.5 at. pct. It occurred at a distance from the oxide-metal interface which was correlatable with temperature and time of exposure. The diffusion coefficient corresponding to the composition of this moving boundary is in excellent agreement with the effective diffusion coefficient for the substrate contamination.

Unnam, J.↗
NASA NTRS1978

Partially averaged field approach to cosmic ray diffusion

A new nonlinear technique is used to derive the kinetic equation for particles interacting with turbulent fluctuations. Difficulties associated with quasi-linear theory are avoided. The new method evaluates the effects of the fluctuations along particle orbits which themselves include the effects of a statistically averaged subset of the possible configuration of the turbulence. As an illustration, the pitch-angle diffusion coefficient is calculated for particles interacting with 'slab model' magnetic turbulence, i.e., magnetic fluctuations linearly polarized transverse to a mean magnetic field. The pitch-angle diffusion coefficient is determined in the vicinity of 90-deg pitch angles where quasi-linear theory breaks down. The spatial diffusion coefficient parallel to a mean magnetic field is evaluated by use of the calculated pitch-angle diffusion coefficient. It is suggested that the partially averaged field method is not limited to small amplitude fluctuating fields and hence is not a perturbation theory.

Jones, F. C.↗
NASA NTRS1988

Venus' superrotation, mixing length theory and eddy diffusion - A parametric study

The concept of the Hadley mechanism is adopted to describe the axisymmetric circulation of the Venus atmosphere. It is shown that, for the atmosphere of a slowly rotating planet such as Venus, a form of the nonliner 'closure' (self-consistent solution) of the fluid dynamics system which constrains the magnitude of the eddy diffusion coefficients can be postulated. A nonlinear one-layer spectral model of the zonally symmetric circulation was then used to establish the relationship between the heat source, the meridional circulation, and the eddy diffusion coefficients, yielding large zonal velocities. Computer experiments indicated that proportional changes in the heat source and eddy diffusion coefficients do not significantly change the zonal velocities. It was also found that, for large eddy diffusion coefficients, the meridional velocity is virtually constant; below a threshold in the diffusion rate, the meridional velocity decreases; and, for large eddy diffusion and small heating rates, the zonal velocities decrease with decreasing planetary rotation rates.

Mayr, H. G.↗
NASA NTRS1992

Mean free paths of energetic particles at very large heliodistances (Pioneer 11 at 20 AU)

The parallel mean free path and the diffusion coefficient parallel to the magnetic field line are derived from magnetic field data at 20 AU to characterize heliospheric modulation and energetic-particle/magnetic-field interaction. The computational method of Moussas et al. (1975, 1982) is employed, and the values of the parallel mean free path are shown to be significantly larger than the values estimated in studies of up to 6 AU. The distance dependence of the parallel diffusion mean free path is found to follow a power law, and the diffusion coefficient dependence upon energy is determined by a constant mean free path and the velocity of the particle. The contribution of the diffusion coefficient perpendicular to the magnetic field is expected to dominate the radial diffusion coefficient of cosmic rays, although the contribution of the diffusion parallel to the field is important with respect to the small-scale structure of intensity gradients.

Moussas, X.↗
NASA NTRS2019

Improvements to an Explicit Algebraic Stress Model for Turbulent Jet Mixing Predictions

Modifications to key coefficients in a k − E based explicit algebraic stress model (EASM) are examined with the objective of improving the prediction of turbulent jet flows. The pressure strain coefficient, C2 and the turbulent diffusion coefficients, σk and σE were investigated. For a series of benchmark subsonic jets at heated and unheated conditions, lowering C2 from the default value of 0.36 to 0.10 resulted in a significant improvement in the jet mixing, when compared to experimental data. Changing σk and σE from default values of 1.00 and 1.4489, respectively, to 0.50 and 0.7244, respectively, improved the initial mixing rate, while reducing the farfield mixing rate and the peak turbulent kinetic energy along the centerline. A high-speed mixing layer was also investigated for performance of baseline and modified EASM coefficients, with similar results as for the jet cases. A flat plate boundary layer was briefly examined to determine the effects of changing the coefficients on the turbulent skin friction coefficient. The change to the pressure strain coefficient, C2 = 0.10 is recommended for future EASM calculation of jets flow; however, it is also recommended that the diffusion coefficients remain at their default values.

Turbulent Mixing↗
NASA NTRS1994

Long-term modulation of galactic cosmic radiation and its model for space exploration

As the human exploration of space has received new attention in the United States, studies find that exposure to space radiation could adversely impact the mission design. Galactic Cosmic Radiation (GCR), with its very wide range of charges and energies, is particularly important for a mission to Mars, because it imposes a stiff mass penalty for spacecraft shielding. Dose equivalent versus shielding thickness calculations, show a rapid initial drop on exposure with thickness, but an asymptotic behavior at a higher shielding thickness. Uncertainties in the radiobiology are largely unknown. For a fixed radiation risk, this leads to large uncertainties in shielding thickness for small uncertainties in estimated dose. In this paper we investigate the application of steady-state, spherically-symmetric diffusion-convection theory of solar modulation to individual measurements of differential energy spectra from 1954 to 1989 in order to estimate the diffusion coefficient, kappa(r,t), as a function of time. We have correlated the diffusion coefficient to the Climax neutron monitor rates and show that, if the diffusion coefficient can be separated into independent functions of space and time: kappa(r,t)=K(t) k(sub 0) beta Pk(sub 1) (r), where beta is the particle velocity and P the rigidity, then (i) The time dependent quantity 1/K(t), which is proportional to the deceleration potential, phi(r,t), is linearly related to the Climax neutron monitor counting rate. (ii) The coefficients obtained from hydrogen or helium intensity measurements are the same. (iii) There are different correlation functions for odd and even solar cycles. (iv) The correlation function for the Climax neutron monitor counting rate for given time, t, can be used to estimate mean deceleration parameter phi(t) to within +/- 15% with 90% confidence. We have shown that k(r,t) determined from hydrogen and/or helium data, can be used to fit the oxygen and iron differential energy spectra with a root mean square error of about +/- 10%, and essentially independent of the particle charge or energy. We have also examined the ion chamber and C-14 measurements which allow the analysis to be extended from the year 1906 to 1990. Using this model we have defined reference GCR spectra at solar minimum and solar maximum. These can be used for space exploration studies and provide a quantitative estimate of the error in dose due to changes in GCR intensities.

Badhwan, G. D.↗
NASA NTRS2019

EMIC Wave-Driven Bounce Resonance Scattering of Energetic Electrons in the Inner Magnetosphere

While electromagnetic ion cyclotron (EMIC) waves have been long studied as a scattering mechanism for ultrarelativistic (megaelectron volt) electrons via cyclotron‐resonant interactions, these waves are also of the right frequency to resonate with the bounce motion of lower‐energy (approximately tens to hundreds of kiloelectron volts) electrons. Here we investigate the effectiveness of this bounce resonance interaction to better determine the effects of EMIC waves on subrelativistic electron populations in Earth's inner magnetosphere. Using wave and plasma parameters directly measured by the Van Allen Probes, we estimate bounce resonance diffusion coefficients for four different events, illustrative of wave and plasma parameters to be encountered in the inner magnetosphere. The range of electron energies and pitch angles affected is examined to better assess the realistic effects of EMIC‐driven bounce resonance on energetic electron populations based on actual, locally observed event‐based parameters. Significant local diffusion coefficients (~ > 10(exp −6) s(exp −1)) for 50‐ to 100‐keV electrons are achieved for both H+ band wave events as well as He+ band, with diffusion coefficients peaking for near‐90° pitch angles but remaining elevated for intermediate ones as well. Diffusion coefficients for higher‐energy 200‐keV electrons are typically multiple orders of magnitude lower (ranging from 10(exp −11) to 10(exp −6) s(exp −1)) and often peak at lower pitch angles (~20–30°). These results suggest that both H+ and He+ band EMIC waves can play a role in shaping lower‐energy electron dynamics via bounce‐resonant interactions, in addition to their role in relativistic electron loss via cyclotron resonance.

Blum, L. W.↗
NASA NTRS1987

Oxygen-18 tracer study of the passive thermal oxidation of silicon

This work focuses on the thermal oxidation of silicon near 1273 K using the double-tracer oxidation method. The results confirm that oxidation occurs by the transport of electrically neutral nonnetwork oxygen through the interstitial space of the vitreous silica (v-SiO2) scale. Simultaneously, self- (or isotopic-) diffusion occurs in the network, resulting in characteristic isotopic fraction distributions near the gas-scale interface. The self-diffusion coefficients calculated from these profiles agree with those reported for tracer diffusion in v-SiO2, and the diffusion coefficient calculated from the scale growth is consistent with reported O2 permeation data. An important parameter that describes the double-oxidation behavior is a ratio relating the scale thickness grown during the second oxidation, the network self-diffusion coefficient for oxygen, and the time of the second oxidation.

Cawley, J. D.↗