Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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

Mantle devolatilization and rheology in the framework of planetary evolution

We investigate the thermal history of an Earth-like planet with the help of a parameterized mantle convection model including the volatile exchange between mantle and surface reservoirs. The weakening of mantle silicates by dissolved volatiles is described by a functional relationship between creep rate and water fugacity. We use flow law parameters of diffusion creep in olivine under dry and wet conditions. The mantle degassing rate is considered as directly proportional to the seafloor spreading rate, which is also dependent on the mantle heat flow. To calculate the spreading rate, we assume that the heat flow under the mid-ocean ridges is double the average mantle heat flow. The rate of regassing also depends on the seafloor spreading rate as well as on other factors like the efficiency of volatile recycling through island arc volcanism. Both mechanisms (de- and regassing) are coupled self-consistently with the help of the parameterized convection model under implementation of a temperature and volatile-content-dependent mantle viscosity.

Franck, S.↗

Critical review of electrical conductivity measurements and charge distribution analysis of magnesium oxide

The electrical conductivity sigma of MgO single crystals shows a sharp increase at 500-800 C, in particular of sigma surface, generally attributed to surface contamination. Charge Distribution Analysis (CDA), a new technique providing information on fundamental properties that was previously unavailable, allows for the determination of surface charges, their sign and associated internal electric field. Data on 99.99% purity, arc-fusion grown MgO crystals show that mobile charge carriers start to appear in the bulk of the MgO crystals between 200 and 400 C when sigma (measured by conventional techniques) is in t he 10(exp -14) to 10(exp -16) /omega/cm range. Above 500 C, as sigma increases to 10(exp -6) to 10(exp -7)/omega/cm, more charges appear giving rise to a strong positive surface charge supported by a strong internal field. This indicates that charges are generated in the bulk and diffuse to the surface by an internally controlled process. On the basis of their positive sign they are identified as holes (defect electrons). Because of the low cation content of these very pure MgO crystals, theses holes cannnot be associated with transition metal impurties. Instead, they are associated with the O(2-) sublattice, e.g. consist of O(-) states or positive holes. This conclusion is supported by magnetic susceptibility data showing the appearance of 1000 +/- 500 ppm paramagnetic species between 200-500 C. The magnetic data are consistent with strongly coupled, diamagnetic O(-) pairs below 200-500 C, chemically equivalent to peroxy anions, O2(2-), and probably associated with cation vacancies in the MgO matrix. The formation of O2(2-) in arc-fusion grown MgO crystals is very unexpected because of the highly reducing growth conditions. Their presence implies an internal redox reaction involving dissolved 'water' by which OH(-) pairs convert to O2(2-) plus H2 molecules. This redox conversion is supported by mass spectroscopic measurements of the H2 release from highly OH(-)-doped, finely divided MgO and by wet-chemical analysis of its oxidant concentration.

Freund, Friedemann↗

Technology Alignment and Portfolio Prioritization (TAPP): Advanced Methods in Strategic Analysis, Technology Forecasting and Long Term Planning for Human Exploration and Operations, Advanced Exploration Systems and Advanced Concepts

The Advanced Concepts Office (ACO) at NASA, Marshall Space Flight Center is expanding its current technology assessment methodologies. ACO is developing a framework called TAPP that uses a variety of methods, such as association mining and rule learning from data mining, structure development using a Technological Innovation System (TIS), and social network modeling to measure structural relationships. The role of ACO is to 1) produce a broad spectrum of ideas and alternatives for a variety of NASA's missions, 2) determine mission architecture feasibility and appropriateness to NASA's strategic plans, and 3) define a project in enough detail to establish an initial baseline capable of meeting mission objectives ACO's role supports the decision­-making process associated with the maturation of concepts for traveling through, living in, and understanding space. ACO performs concept studies and technology assessments to determine the degree of alignment between mission objectives and new technologies. The first step in technology assessment is to identify the current technology maturity in terms of a technology readiness level (TRL). The second step is to determine the difficulty associated with advancing a technology from one state to the next state. NASA has used TRLs since 1970 and ACO formalized them in 1995. The DoD, ESA, Oil & Gas, and DoE have adopted TRLs as a means to assess technology maturity. However, "with the emergence of more complex systems and system of systems, it has been increasingly recognized that TRL assessments have limitations, especially when considering [the] integration of complex systems." When performing the second step in a technology assessment, NASA requires that an Advancement Degree of Difficulty (AD2) method be utilized. NASA has used and developed or used a variety of methods to perform this step: Expert Opinion or Delphi Approach, Value Engineering or Value Stream, Analytical Hierarchy Process (AHP), Technique for the Order of Prioritization by Similarity to Ideal Solution (TOPSIS), and other multi­‐criteria decision-making methods. These methods can be labor-intensive, often contain cognitive or parochial bias, and do not consider the competing prioritization between mission architectures. Strategic Decision-Making (SDM) processes cannot be properly understood unless the context of the technology is understood. This makes assessing technological change particularly challenging due to the relationships "between incumbent technology and the incumbent (innovation) system in relation to the emerging technology and the emerging innovation system." The central idea in technology dynamics is to consider all activities that contribute to the development, diffusion, and use of innovations as system functions. Bergek defines system functions within a TIS to address what is actually happening and has a direct influence on the ultimate performance of the system and technology development. ACO uses similar metrics and is expanding these metrics to account for the structure and context of the technology. At NASA technology and strategy is strongly interrelated. NASA's Strategic Space Technology Investment Plan (SSTIP) prioritizes those technologies essential to the pursuit of NASA's missions and national interests. The SSTIP is strongly coupled with NASA's Technology Roadmaps to provide investment guidance during the next four years, within a twenty-year horizon. This paper discusses the methods ACO is currently developing to better perform technology assessments while taking into consideration Strategic Alignment, Technology Forecasting, and Long Term Planning.

Funaro, Gregory V.↗

Investigating Global Ion and Neutral Atom Populations with IBEX and Voyager

The main objective of this project was to investigate pickup ion (PUI) production in the solar wind and heliosheath (the region between the termination shock and the heliopause) and compute the distributed energetic neutral atom fluxes throughout the helioshpere. The simulations were constrained by comparing the model output against observations from Ulysses, New Horizons, Voyager 1 and 2, and IBEX space probes. As evidenced by the number of peer reviewed journal publications resulting from the project (13 plus three submitted) and their citation rate (156 citations over three years), the project has made a lasting contribution to the field. The outcome is a significant improvement of our understanding of the pickup ion production and distribution in the distant heliosphere. The team has accomplished the entire set of tasks A-H set forth in the proposal. Namely, the transport modeling framework has been augmented with two populations of pickup ions (PUIs), the boundary conditions for the plasma and interstellar neutral hydrogen were verified against Ulysses and New Horizons PUI and an optimal set of velocity diffusion parameters established. The multi-component fluxes of PUIs were computed and isotropic velocity distributions generated for each cell in the computer simulation that covered the heliosphere from 1.5 AU to the heliopause. The distributions were carefully compared with in situ measurements at 3 AU (Ulysses), 12 AU (New Horizons), and 80-90 AU (Voyager 1 and 2) as well as those inferred from ENA fluxes measured by Cassini and IBEX (Wu et al., 2016). Some examples of modeldata comparison are shown in Figure 1. We have used coupled MHD-plasma and kinetic-neutral code to investigate the likely range of plasma and magnetic field parameters in the local interstellar medium (LISM), based on the assumption that the shape of the IBEX ribbon could be used to determine the orientation of the interstellar magnetic field. While the magnetic field is believed to be oriented toward the center of the ribbon, constraining its strength requires comparing the model-predicted angular diameter and circularity of the ribbon with the observations. The study, published in Heerikhuisen et al. (2014), found that the most likely range for the LISM magnetic field strength is between 0.2 and 0.3 nT, which is less than previously thought. Figure 2 shows the IBEX data (left) and compares it to the simulation with a 0.2 nT interstellar magnetic field (center) and a 0.4 nT (right).

Florinski, Vladimir↗

A geostationary imaging spectrometer TOMS instrument

One design for a geostationary Total Ozone Mapping Spectrometer (TOMS) with many desirable features is an imaging spectrometer. A preliminary study makes use of a 0.25 m Czerny-Turner spectrometer with which the Earth is imaged on a charge-coupled device (CCD) in dispersed light. The wavelength is determined by a movable grating which can be set arbitrarily by ground control. The signal integration time depends on wavelength but this system allows arbitrary timing by command. Special circumstances such as a requirement to track a low-lying sulfur dioxide cloud or a need to discriminate high level ozone from total ozone at midlatitudes could be obtained by adding a particular wavelength to the normally pre-programmed time sequence. The incident solar irradiance is measured by deploying a diffuser plate in the field of view. Individual detector elements correspond to scene elements in which the several wavelengths are serially sampled and the Earth radiance is compared to the incident sunlight. Thus the problem of uncorrelated drift of multiple detectors is removed.

Krueger, Arlin J.↗

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.

Kennedy, Christopher A.↗

Mechanisms of Ionospheric Mass Ejection

Ionospheric outflows are directly responsive to solar wind disturbances, particularly in the dayside auroral cusp or cleft regions. Inputs of both electromagnetic energy (Poynting flux) and kinetic energy (particle precipitation) are closely correlated with these outflows. We assess the importance of processes thought to drive ionospheric outflows. These begin with the diffuse effects of photoionization and thermal equilibrium of the ionospheric topside, enhancing Jeans' escape, with ambipolar diffusion and acceleration. Auroral outflows begin with dayside reconnexion and resultant field-aligned currents and driven convection. These produce plasmaspheric plumes, collisional heating and wave-particle interactions, centrifugal acceleration, and auroral acceleration by parallel electric fields, including enhanced ambipolar fields from electron heating by precipitation particles. Solar wind energy dissipation is concentrated by the geomagnetic field into auroral regions with an amplification factor of 10-100, enhancing heavy species plasma and gas escape from gravity, and providing more current carrying capacity. Internal plasmas thus enable electromagnetic driving via coupling to the plasma and neutral gas. We assess the importance of each of these processes in terms of local escape flux production as well as global outflow, and suggest methods for their implementation within multi-species global simulation codes. We conclude by assessing outstanding obstacles to this objective.

Moore, Thomas Earle↗

Estimating the Distribution of Colored Dissolved Organic Matter During the Southern Ocean Gas Exchange Experiment Using Four-Dimensional Variational Data Assimilation

We diagnosed the effect of various physical processes on the distribution of mixed-layer colored dissolved organic matter (CDOM) and a sulfur hexauoride (SF6) tracer during the Southern Ocean Gas Exchange Experiment (SO GasEx). The biochemical upper ocean state estimate uses in situ and satellite biochemical and physical data in the study region, including CDOM (absorption coefcient and spectral slope), SF6, hydrography, and sea level anomaly. Modules for photobleaching of CDOM and surface transport of SF6 were coupled with an ocean circulation model for this purpose. The observed spatial and temporal variations in CDOM were captured by the state estimate without including any new biological source term for CDOM, assuming it to be negligible over the 26 days of the state estimate. Thermocline entrainment and photobleaching acted to diminish the mixed-layer CDOM with time scales of 18 and 16 days, respectively. Lateral advection of CDOM played a dominant role and increased the mixed-layer CDOM with a time scale of 12 days, whereas lateral diffusion of CDOM was negligible. A Lagrangian view on the CDOM variability was demonstrated by using the SF6 as a weighting function to integrate the CDOM elds. This and similar data assimilation methods can be used to provide reasonable estimates of optical properties, and other physical parameters over the short-term duration of a research cruise, and help in the tracking of tracer releases in large-scale oceanographic experiments, and in oceanographic process studies.

Del Castillo, C. 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.↗

Radio to Gamma-Ray Emission from Shell-Type Supernova Remnants: Predictions from Non-Linear Shock Acceleration Models

Supernova remnants (SNRs) are widely believed to be the principal source of galactic cosmic rays, produced by diffusive shock acceleration in the environs of the remnant's expanding blast wave. Such energetic particles can produce gamma-rays and lower energy photons via interactions with the ambient plasma. The recently reported observation of TeV gamma-rays from SN1006 by the CANGAROO Collaboration, combined with the fact that several unidentified EGRET sources have been associated with known radio/optical/X-ray-emitting remnants, provides powerful motivation for studying gamma-ray emission from SNRs. In this paper, we present results from a Monte Carlo simulation of non-linear shock structure and acceleration coupled with photon emission in shell-like SNRs. These non-linearities are a by-product of the dynamical influence of the accelerated cosmic rays on the shocked plasma and result in distributions of cosmic rays which deviate from pure power-laws. Such deviations are crucial to acceleration efficiency considerations and impact photon intensities and spectral shapes at all energies, producing GeV/TeV intensity ratios that are quite different from test particle predictions.

Baring, Matthew G.↗

The Effect of Stochastic Acceleration on the Dynamics of Solar Energetic Particles in the Heliosphere

Solar Energetic Particles (SEPs) are crucial in space weather phenomena, impacting space-based technologies and human activities in space. Understanding the dynamics of SEPs, including their acceleration and transport mechanisms, is essential for improving predictive models and mitigating adverse effects. This study investigates the impact of stochastic acceleration (SA) on the dynamics of SEPs within the heliosphere. By integrating stochastic acceleration into a SEP transport model, we assess its impact relative to diffusive shock acceleration (DSA). The investigation utilizes the Adaptive Mesh Particle Simulator (AMPS) within the Space Weather Modeling Framework (SWMF) to simulate SEP dynamics, incorporating solar wind parameters and turbulence from the Alfven Wave Solar Model (AWSoM). To do this, we solve the focused transport equation to model the propagation of SEPs along magnetic field lines, extending these lines to 5 AU to incorporate the potential effects of pitch angle scattering beyond 1 AU. This aspect could notably impact the decay phase of a SEP event. Modeling the transport of SEPs is coupled with the simulation of solar wind dynamics, the interplanetary magnetic field, and the characteristics of Alfven wave turbulence. This presentation covers the modeling approach used in this research, the characterization of the impact of the stochastic acceleration, and the role of pitch angle scattering at various heliocentric distances on the dynamics of the SEP events' decay phase.

SEP↗

The treatment of convected vortices in compressible potential flow

A method is described for incorporating line vortices into the three dimensional compressible potential flow equation. A modified Biot-Savart law is used to compute a vortical velocity field, which is added to the gradient of the potential to form a total velocity. A rapidly converging approximate factorization (AFZ) scheme is then used to compute a potential such that the modified potential flow equation as well as the appropriate boundary conditions, based on total velocity, are satisfied. As part of a coupled iteration procedure, the positions of the line vortices are computed so that they convect with the total flow. The method is used to compute the field due to a single line vortex convecting past a wing. This represents an approximation of the effect of a canard or other lifting surface ahead of the wing, which sheds a tip vortex. It is seen that the flow field is substantially modified by the passage of the vortex. Unlike Euler equation schemes, which are also used to compute these flows, the solutions exhibit no numerical diffusion: The convected vortices retain their initial upstream width. Euler solutions, on the other hand, involve a vorticity which is numerically convected in an Eulerian frame and, unless extensive adaptive grid refinement is used they result in vortices with spread as they convect. Also, the potential flow method requires approximately two orders of magnitude less computing time and much less computer storage than the Euler methods.

Steinhoff, J.↗

Slewing and vibration control of the SCOLE

A discussion of Slewing and Vibration Control makes the following conclusions: (1) A 2-stage approach is feasible and promising for rapid slewing and precision pointing of SCOLE; (2) Not all bang-bang type of time-minimized slew maneuvers will excite large structural vibrations in SCOLE; and (3) Modal dashpots can be a concentrated high-power vibration control, as well as the usual diffuse (broadband, low-power (low-authority) control. The following recommendations are made: (1) Limit the magnitude of applied forces on reflector to either the 25-lb limit of vernier thrusters on the real Space Shuttle or the 150-lb level equivalent to the cold-gas jets of laboratory SCOLE; (2) to complete stage 2, add an integrated design of LQF/LTR (Linear-Quadratic-Gaussian/Loop-Transfer Recovery) and Modal Dashpots; and, (3) Validate the 2-stage approach using the SCOLE laboratory facility with a comprehensive sequence of integrated designs and experiments coupling nonlinear rigid-body motions with flexible-body dynamics.

Lin, Jiguan Gene↗

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗

Analysis of Thermal and Reaction Times for Hydrogen Reduction of Lunar Regolith

System analysis of oxygen production by hydrogen reduction of lunar regolith has shown the importance of the relative time scales for regolith heating and chemical reaction to overall performance. These values determine the sizing and power requirements of the system and also impact the number and operational phasing of reaction chambers. In this paper, a Nusselt number correlation analysis is performed to determine the heat transfer rates and regolith heat up times in a fluidized bed reactor heated by a central heating element (e.g., a resistively heated rod, or a solar concentrator heat pipe). A coupled chemical and transport model has also been developed for the chemical reduction of regolith by a continuous flow of hydrogen. The regolith conversion occurs on the surfaces of and within the regolith particles. Several important quantities are identified as a result of the above analyses. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the particle Reynolds number, the Archimedes number, and the time needed for hydrogen to diffuse into the pores of the regolith particles. The analysis is used to determine the heat up and reaction times and its application to NASA s oxygen production system modeling tool is noted.

Hegde, U.↗

Analysis of Thermal and Reaction Times for Hydrogen Reduction of Lunar Regolith

System analysis of oxygen production by hydrogen reduction of lunar regolith has shown the importance of the relative time scales for regolith heating and chemical reaction to overall performance. These values determine the sizing and power requirements of the system and also impact the number and operational phasing of reaction chambers. In this paper, a Nusselt number correlation analysis is performed to determine the heat transfer rates and regolith heat up times in a fluidized bed reactor heated by a central heating element (e.g., a resistively heated rod, or a solar concentrator heat pipe). A coupled chemical and transport model has also been developed for the chemical reduction of regolith by a continuous flow of hydrogen. The regolith conversion occurs on the surfaces of and within the regolith particles. Several important quantities are identified as a result of the above analyses. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the particle Reynolds number, the Archimedes number, and the time needed for hydrogen to diffuse into the pores of the regolith particles. The analysis is used to determine the heat up and reaction times and its application to NASA s oxygen production system modeling tool is noted.

Hegde, U.↗

Xenon Fractionation, Hydrogen Escape, and the Oxidation of the Earth

Xenon in Earth's atmosphere is severely mass fractionated and depleted compared to any plausible solar system source material, yet Kr is unfractionated. These observations seem to imply that Xe has escaped from Earth. But to date no process has been identified that can cause Xe, which is heavier than Kr, to escape while Kr does not. Vigorous hydrodynamic hydrogen escape can produce mass fractionation in heavy gases. The required hydrogen flux is very high but within the possible range permitted by solar Extreme Ultraviolet radiation (EUV, which here means radiation at wavelengths short enough to be absorbed efficiently by hydrogen) heating when Earth was on the order of 100 Myrs old or younger. However this model cannot explain why Xe escapes but Kr does not. Recently, what appears to be ancient atmospheric xenon has been recovered from several very ancient (3-3.5 Ga) terrestrial hydrothermal barites and cherts. What is eye-catching about this ancient Xe is that it is less fractionated that Xe in modern air. In other words, it appears that a process was active on Earth some 3 to 3.5 billion years ago that caused xenon to fractionate. By this time the Sun was no longer the EUV source that it used to. If xenon was being fractionated by escape — currently the only viable hypothesis — it had to be in the less unfamiliar context of Earth’s Archean atmosphere and under rather modest levels of EUV forcing. This requires a new model. Here we address the circumstances in which Xe, but not Kr, could escape from Earth as an ion. In a hydrodynamically escaping hydrogen wind the hydrogen is partially photo-ionized. The key concepts are that ions are much more strongly coupled to the escaping flow than are neutrals (so that a relatively modest flow of H and H+ to space could carry Xe+ along with it), and that xenon alone among the noble gases is more easily ionized than hydrogen. This sort of escape is possible if not prevented by a planetary magnetic field. The best prospects for Earth are therefore escape along the polar field lines, although a very weak or absent magnetic field would likely work as well. As applied to the Archean Earth the discussion will be constrained by diffusion-limited hydrogen escape. The extended history of hydrogen escape implicit in Xe escape in the Archean is consistent with suggestions that hydrogen escape from the anoxic Archean atmosphere was considerable, because biogenic methane is expected to have been rather abundant. Hydrogen escape plausibly played the key role in creating oxidizing condition at the surface of the Earth and setting the stage for the creation of an O2 atmosphere.

photo-ionized gases↗

Vertical transport and photochemistry in the terrestrial mesosphere and lower thermosphere /50-120 km/

A study is conducted of the coupled effects of kinetics, solar cycle flux variations, and vertical transport on the distribution of long-lived hydrogen-carbon-oxygen compounds in the terrestrial mesosphere and lower thermosphere, using a one-dimensional aeronomy model. The calculations account for the important chemical reactions and use rocket measurements of the solar flux at solar minimum and maximum. Photodissociation rates appropriate for the mesosphere are determined with a spherical shell atmosphere formalism. Detailed corrections for the O2 Schumann-Runge bands and the temperature dependence of the CO2 cross sections are used. An eddy diffusion profile is derived which is in agreement with the Aladdin 74 mass spectral measurements of atomic O, O2, CO2, and Ar in the lower thermosphere and observations of the O3 minimum at about 80 km.

Allen, M.↗