Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chemical initiation”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Enabling Multireference Calculations on Multimetallic Systems with Graphic Processing Units

Modeling multimetallic systems efficiently enables faster prediction of desirable chemical properties and the design of new materials. This work describes an initial implementation for performing multireference wave function method localized active-space self-consistent field (LASSCF) calculations through the use of multiple graphics processing units (GPUs) to accelerate time-to-solution. Density fitting is leveraged to reduce memory requirements, and we demonstrate the ability to fully utilize multi-GPU compute nodes. Performance improvements of 5–10x in total application runtime were observed in LASSCF calculations for multimetallic catalyst systems up to 1200 AOs and an active space of (22e,40o) using up to four NVIDIA A100 GPUs. Furthermore, written with performance portability in mind, a comparable performance is also observed in early runs on the Aurora exascale system using Intel Max Series GPUs.

Algorithms↗
DOE OSTI2026

Electrochemically Initiated Depolymerization of Poly(Methyl Methacrylate)

Efficient depolymerization of polymers with all-carbon backbones under mild conditions would be valuable in chemically recycling commodity plastics. Poly(methyl methacrylate) (PMMA) is a commodity thermoplastic that is currently depolymerized under temperatures in excess of 400 °C. Herein, we lower the temperatures needed to achieve depolymerization of PMMA by performing radical generation and depropagation with orthogonal stimuli. This first demonstration of electrochemically initiated PMMA depolymerization relies on reduction of phthalimide esters that, upon subsequent decarboxylation, generate polymer-centered radicals. These radicals then spontaneously unzip the polymer back to its monomeric constituents at temperatures as low as 105 °C. We studied the mechanism and efficiency of this transformation as a function of phthalimide ester placement, incorporation density, and polymer molecular weight. We found that chain-end activation is effective for modest molecular weights but suffers diminished efficiency at higher degrees of polymerization. In contrast, pendent-group activation is more effective for depolymerizing higher molecular weight species. Integrating higher molar amounts of phthalimide ester pendants leads to more effective depolymerization, with >95% depolymerization in copolymers with 5 mol% phthalimide ester incorporation. We leveraged this understanding to create a custom electro-distillation apparatus that allowed us to simultaneously electrochemically depolymerize PMMA and directly distill methyl methacrylate in >22% yield, which could ultimately be repolymerized. These findings establish electrochemistry as a versatile and orthogonal stimulus for vinyl polymer depolymerization and provide a foundation for closed-loop electrochemical recycling of widely used plastics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

Space Weathering Experiments on Spacecraft Materials

A project to investigate space environment effects on specific materials with interest to remote sensing was initiated in 2016. The goal of the project is to better characterize changes in the optical properties of polymers and Mylar, specifically those found in multi-layered spacecraft insulation, due to electron bombardment. Previous analysis shows that chemical bonds break and potentially reform when exposed to high energy electrons. Among other properties these chemical changes altered the optical reflectance as documented in laboratory analysis. This paper presents results of the initial experiment results focused on the exposure of materials to various fluences of high energy electrons, used to simulate a portion of the geosynchronous space environment. The paper illustrates how the spectral reflectance changes as a function of time on orbit with respect to GEO environmental factors and investigates the survivability of the material after multiple electron doses. These results provide a baseline for analysis of aging effects on satellite systems used for remote sensing. They also provide preliminary analysis on what materials are most likely to encompass the high area-to-mass population of space debris in the geosynchronous environment. Lastly, the paper provides the results of the initial experimentation as a proof of concept for space aging on polymers and Mylar for conducting more experiments with a larger subset of spacecraft materials.

Cooper, R.↗
DOE OSTI2022

First-Principles Elucidation of Initial Dehydrogenation Pathways in Mg(BH 4 ) 2

Complex borohydrides such as Mg(BH 4 ) 2 offer one of highest capacities to chemically store hydrogen for onboard applications; however, it suffers greatly from kinetic constraints that prevent realization of full capacity and reversibility. Understanding these kinetic limitations solely from experiments is extremely challenging due to the unusual complexity of various competing elemental reaction steps involved during the de/rehydrogenation reaction. This work aims to map out the energetics associated with initial dehydrogenation of Mg(BH 4 ) 2 from first-principles simulations and to identify the preferred reaction pathways. Our calculations suggest the rate-limiting step during BH 4 – –B 3 H 8 – conversion is the formation of the B 2 H 7 – intermediate. We further emphasize and clarify that the B 3 H 8 – and H – intermediates, formed during initial Mg(BH 4 ) 2 decomposition, appear as molecular species that are embedded in the Mg–BH 4 –Mg matrix as evidenced in the nuclear magnetic resonance measurements and not as bulk MgH 2 and Mg(B 3 H 8 ) 2 as previously assumed in theoretical predictions of the thermodynamics.

08 HYDROGEN↗
DOE OSTI2026

Modelling the Sensitivity of Yukon River Biogeochemical Dynamics to Environmental and Chemical Drivers: Implications for Dissolved Organic Carbon

Riverine dissolved organic carbon (DOC) is a critical biogeochemical component that transmits information from Arctic soils to the Arctic Ocean, significantly influencing carbon dynamics in this unique ecosystem. As DOC travels downstream, it undergoes transformations that alter its composition and fate. The Yukon River serves as an effective testbed for modelling these dynamics, offering sufficient scale to capture key biogeochemical processes while having a simpler hydrology than other major Arctic rivers, as well as long-term DOC observational data for model validation. To investigate DOC transformations during transit in the Yukon River, we adapted our Arctic Riverine Organic Macromolecular Model by applying regional-specific parameterisations. Our model simulates the transport and transformation of 15 organic macromolecules, including CDOM (coloured dissolved organic matter), proteins, polysaccharides, lipids, lignin phenols, and humic substances. Initial DOC concentrations were derived from observed soil organic carbon stocks in the surrounding watershed, while chemical transformations and hydrological dynamics were modelled along the river's course. Sensitivity and uncertainty analyses were conducted using a Monte Carlo approach under two experimental setups. Results revealed that variability in DOC and CDOM concentrations at the river mouth were predominantly driven by initial DOC concentration (~70% of variability explained) and dilution at confluence points (~10%). The refractory fraction of DOC explained 21%–88% of the variability in 14 macromolecular concentrations and ranked in the top five sensitive parameters for all outputs when a uniform parameter distribution was assumed. However, when a more likely variability was applied to this parameter, its influence on DOC and CDOM decreased. Given that refractory DOC accounts for ~80% of total DOC in Arctic Rivers, this suggests that most DOC resists degradation and retains its chemical composition during transport to the coastal environment. River velocity, which determines residence time, explained 8%–47% of the variability in protein, polysaccharide, lipid, pigments, and lignin phenols at the river mouth. In contrast, chemical turnover times contributed only 1%–5% to output variability. Our findings underscore the need for improved land-specific headwater observations, including seasonal soil moisture and lateral transport dynamics that control the initial tributary-specific DOC inputs. With accelerated permafrost thaw and increasing river discharge, extending our model to other Arctic River systems and seasons will enhance understanding of Arctic riverine carbon fluxes and their contributions to the Arctic Ocean.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Chemical Analysis of the Ultrafaint Dwarf Galaxy Grus II. Signature of High-mass Stellar Nucleosynthesis

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultrafaint dwarf (UFD) galaxy Grus II. All stars exhibit a higher than expected [Mg/Ca] ratio compared to metal-poor stars in other UFD galaxies and in the Milky Way (MW) halo. Nucleosynthesis in high-mass ( 20 M ⊙) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (three) stars suggests the chemical enrichment of Grus II could have occurred through substantial high-mass stellar evolution, and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it cannot be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture (r-process) elements. The abundance pattern of the r-process elements in this star matches the scaled r-process pattern of the solar system and r-process enhanced stars in other dwarf galaxies and in the MW halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of r-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus II would increase the likelihood of any of these events occurring. The time delay between the α and r-process element enrichment of the galaxy favors a neutron star merger as the origin of the r-process elements in Grus II.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Some aspects of the ozone-depletion problem in the stratosphere

A preliminary investigation of O3 depletion by NOx in the exhaust of an SST wake is presented. Initial calculations were made using a model employing a two-step chemistry mechanism consisting of an NOx catalytic cycle and assuming total NOx and the ratio NO/NO2 to be conserved. The results indicate O3 depletion in the wake much less severe than reported in earlier investigations which used a much larger reaction rate constant. Subsequent investigation revealed, however, that even the current O3 depletion rates appear too large. Because of the apparent deficiencies in the two-step chemical model, a model is presented with a more realistic chemical mechanism without the requirement for NO conservation. Results indicate that total NOx and initial NO/NO2 ratio play an important role in determining O3 depletion. Furthermore, only in the case of rather extreme concentrations of water vapor or unrealistically large photolysis rates is there any significant conversion of NOx to HNO3.

Gupta, R. N.↗
DOE OSTI2020

Role of Intact Hydrogen-Bond Networks in Multiproton-Coupled Electron Transfer

The essential role of a well-defined hydrogen-bond network in achieving chemically reversible multiproton translocations triggered by one-electron electrochemical oxidation/reduction is investigated by using pyridylbenzimidazole–phenol models. The two molecular architectures designed for these studies differ with respect to the position of the N atom on the pyridyl ring. In one of the structures, a hydrogen-bond network extends uninterrupted across the molecule from the phenol to the pyridyl group. Experimental and theoretical evidence indicates that an overall chemically reversible two-protoncoupled electron-transfer process (E2PT) takes place upon electrochemical oxidation of the phenol. This E2PT process yields the pyridinium cation and is observed regardless of the cyclic voltammogram scan rate. In contrast, when the hydrogen-bond network is disrupted, as seen in the isomer, at high scan rates (~1000 mV s –1 ) a chemically reversible process is observed with an E 1/2 characteristic of a one-proton-coupled electron-transfer process (E1PT). At slow cyclic voltammetric scan rates (<1000 mV s –1 ) oxidation of the phenol results in an overall chemically irreversible two-proton-coupled electron-transfer process in which the second proton-transfer step yields the pyridinium cation detected by infrared spectroelectrochemistry. In this case, we postulate an initial intramolecular proton-coupled electron-transfer step yielding the E1PT product followed by a slow, likely intermolecular chemical step involving a second proton transfer to give the E2PT product. Insights into the electrochemical behavior of these systems are provided by theoretical calculations of the electrostatic potentials and electric fields at the site of the transferring protons for the forward and reverse processes. Furthermore, this work addresses a fundamental design principle for constructing molecular wires where protons are translocated over varied distances by a Grotthuss-type mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2005

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗
DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗
NASA NTRS2019

Detailed Infrared Study of Amorphous to Crystalline Propionitrile Ices Relevant to Observed Spectra of Titan's Stratospheric Ice Clouds

We have conducted a comprehensive study of propionitrile (C2H5CN) ice from the amorphous to crystalline phase in order to provide detailed information on this specific cyanide, which may potentially contribute to the chemical composition of the Haystack ice cloud observed in Titan’s stratosphere by the Cassini Composite InfraRed Spectrometer (CIRS). Infrared transmission spectra of thin films of pure propionitrile ices deposited at low temperature (30 – 160 K) were collected from 50 cm(exp -1) to 11700 cm(exp -1) (200 – 0.85 μm). The far-infrared spectral region was specifically targeted to compare with CIRS far-infrared limb spectra. The temperature and time evolution of C2H5CN ice was thoroughly investigated to better understand discrepancies reported in previously published laboratory studies on the crystalline phase of C2H5CN. Specifically, we observe peculiar temperature and time-driven ice phase transitions, revealed by significant spectral variations in the ice, which stabilizes once a complete crystalline phase is achieved. From these results, the crystalline phase of propionitrile ice was identified at deposition temperatures greater than or equal to 135 K and less than 140 K. Our findings corroborate previous studies that ruled out pure propionitrile ice as the sole chemical identity of Titan’s observed Haystack emission feature. As a result, we have initiated co-deposition experiments that incorporate mixtures of Titan relevant organics, many of which have corresponding vapors that are abundantly present in Titan's stratosphere. In this paper, we present the result of one example of a co-deposited trinary ice mixture containing 16% hydrogen cyanide (HCN), 23% C2H5CN, and 61% benzene (C6H6). Although this co-condensed ice mixture is the best fit thus far obtained to match the broad width of the Haystack, it is still not the appropriate chemical candidate. However, it reveals an intriguing result: the strong lattice mode of pure C2H5CN ice is drastically altered by the surrounding molecules as a result of mixing in a co-condensed phase. The laboratory results reported here on propionitrile ice may help to further constrain the chemical identification of Titan’s stratospheric Haystack ice cloud, as well as improve on the current state of knowledge of Titan's stratospheric ice cloud chemistry.

Delphine Nna-Mvondo↗
DOE OSTI2026

Summary of the Initial Post-Irradiation Characterization of HFIR-Irradiated Low-N and High-N HT-9 Steel

Reference cladding systems for sodium fast reactors are based on the historical steel, HT-9. HT-9 is a Fe12Cr ferritic/martensitic steel with additions of Mo, W, V, and other minor elements and demonstrates low irradiation swelling and adequate mechanical properties. Extensive irradiation literature exists on the use of HT-9 as cladding for metal fuel, primarily irradiation on the U-Zr/HT-9 system from the Experimental Breeder Reactor-II (EBR-II) and Fast Flux Test Facility (FFTF) sodium fast reactor, and as a structural material from experiments in the FFTF. The large amount of historical data makes the U-Zr/HT-9 system the reference fuel specification for many nuclear reactor vendors that seek to license modern sodium-cooled fast reactors in the United States. However, it is yet unclear how variations in impurity content within HT-9 fundamentally affect irradiation performance at various irradiation temperatures. Recent work suggests that impurity content may noticeably alter the production of helium through nuclear transmutation. For these reasons, High-Flux Isotope Reactor (HFIR) irradiation of HT-9 steels with known variations in the impurity content is particularly timely to generate data to enable more accurate refinement of the chemical specification for nuclear-grade HT-9 material. This report summarizes the initial transmission electron microscopy characterization of HFIR-irradiated HT-9 steels following mechanical property measurements by the Advanced Fuels Campaign (AFC). This report includes qualitative results of the cavity, dislocation loop and cluster/precipitate microstructures as well as radiation-induced segregation. Quantitative results are being shared with partner institutions and will be included in more detail in a future report in FY2026.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Solvent Effects on Catalytic Reactions and Related Phenomena at Liquid-Solid Interfaces

Catalytic reactions involve the direct interaction of reactants, intermediates and products with the catalyst surface. We not only need to control the atomic structure and electronic properties of the active site, but also explore the multiple molecular interactions that occur beyond the active site; they play an essential role in altering the binding and reactivity of surface species. In liquid-phase catalysis, solvents provide additional degrees of freedom in the design of the catalytic process for desirable activity and selectivity. Here, the multi-faceted effects of solvents have a profound impact on the catalyst performance by restricting the mass transfer to the site, tuning the chemical potential of the surface species, competing for active sites, stabilizing the initial and transition states, and causing mechanistic changes by participating in the kinetically relevant elementary steps. This review addresses the different aspects of solvent effects, using a few prototype solid-liquid interfaces to illustrate these fundamental features. Recent experimental and computational studies that provide new insight at the molecular level are examined. Solvent structures in the proximity of the catalyst surface are discussed along with their influence in molecular binding and reaction at the solid-liquid interfaces. Furthermore, opportunities to alter such a solid-liquid interaction by tuning the wettability of the catalyst surfaces are explored.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Unraveling Grain Boundary Instability in Dense Proton-Conducting Oxides

The long-term stability of protonic ceramic electrolysis cell (PCEC) materials under high-steam operating conditions remains a critical barrier to device commercialization. Here, we investigate the fundamental degradation mechanisms of dense BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3-δ (BCZYYb) electrolytes operated at 550 °C, 50% H 2 O in air. Over 1,000 h, the total electrolyte conductivity decreases by 11.1%, driven primarily by a >130% increase in grain-boundary resistivity. Post-mortem analyses reveal that damage is localized to near-surface grain boundaries extending ∼50 μm into the dense electrolyte pellet. This surface localization indicates that degradation is likely to be severe in thin, device-level electrolytes. Degradation is primarily attributed to chemo-mechanical grain-boundary weakening arising from hydration-induced chemical expansion, culminating in the formation of intergranular cracks oriented parallel to the pellet surface. These internal cracks subsequently react with steam and/or CO 2 , leading to the formation of nanoscale insulating phases, including Ba(OH) 2 , nanocrystalline BaCO 3 , and amorphous Ce/Zr/Y/Yb-containing oxides or hydroxycarbonates. After an initial degradation period of approximately 200 h, the overall conductivity stabilizes. Incorporating NiO sintering aids reduces grain-boundary density by an order of magnitude under identical sintering conditions. Although addition of NiO increases the initial resistivity by >160% at 550 °C, it substantially suppresses grain-boundary instability and mitigates chemical degradation. These findings underscore the urgent need for chemical and/or physical stabilization of BCZYYb electrolytes and offer design guidelines to enable durable, high-performance PCECs.

08 HYDROGEN↗
DOE OSTI2022

Summary of Savannah River Site FY21 Salt Waste Qualification Data and Characterization Data from Salt Waste Batches Processed Since 2008

The Savannah River National Laboratory analyzed samples from Savannah River Site Waste Tank 41H and 21H to support qualification of Salt Waste Processing Facility Waste Batches 3 and 4 (the FY21 Salt Batch Qualification samples) for processing. Neither of the samples displayed any unusual or unexpected characteristics such as large amounts of solids, floating solids, or unusual color. Characterization of these samples confirmed similar chemical composition and characteristics to previous salt waste batches. These results were initially provided to Savannah River Remediation (Liquid Waste Operations Sub-contractor) as External Sample Results LIMS Reports. The analytical results (both expedited and routine) are now summarized and discussed in this technical report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Summary of Savannah River Site FY22 Salt Waste Qualification Data

The Savannah River National Laboratory analyzed samples from Savannah River Site Waste Tank 41H, 21H and 42H to support qualification of Salt Waste Processing Facility (SWPF) Waste Batches 5, 6, and 7 for processing (the FY22 Salt Batch Qualification samples). None of the samples displayed any unusual or unexpected characteristics such as large amounts of solids, floating solids, or unusual color. Characterization of these samples confirmed similar chemical composition and characteristics to previous salt waste batches. These results were initially provided to Savannah River Mission Completion (Liquid Waste Operations subcontractor) as External Sample Results LIMS Reports and the Tank 41H (Salt Batch 5 sample) was also summarized in two separate technical reports. The analytical results (both rapid and long term) are now summarized and discussed in this technical report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1973

Hydrogen adsorption and diffusion, and subcritical-crack growth in high strength steels and nickel base alloys

Embrittlement, or the enhancement of crack growth by gaseous hydrogen in high strength alloys, is of primary interest in selecting alloys for various components in the space shuttle. Embrittlement is known to occur at hydrogen gas pressures ranging from fractions to several hundred atmospheres, and is most severe in the case of martensitic high strength steels. Kinetic information on subcritical crack growth in gaseous hydrogen is sparse at this time. Corroborative information on hydrogen adsorption and diffusion is inadequate to permit a clear determination of the rate controlling process and possible mechanism in hydrogen enhanced crack growth, and for estimating behavior over a range of temperatures and pressures. Therefore, coordinated studies of the kinetics of crack growth, and adsorption and diffusion of hydrogen, using identical materials, have been initiated. Comparable conditions of temperature and pressure will be used in the chemical and mechanical experiments. Inconel 718 alloy and 18Ni(200) maraging steel have been selected for these studies. Results from these studies are expected to provide not only a better understanding of the gaseous hydrogen embrittlement phenomenon itself, but also fundamental information on hydrogen adsorption and diffusion, and crack growth information that can be used directly for design.

Wei, R. P.↗
NASA NTRS1973

Theoretical model of vertical distributions of CO and CH4 in the mesosphere and upper stratosphere.

Carbon compounds (CO and CO2) and methane are added to an existing computer program (Shimazaki and Laird, 1970, 1972), which calculates diurnal variations in density distributions of pure oxygen compounds, hydrogen compounds, and nitrogen compounds. For the present model calculations, the original program was modified to calculate the steady-state profiles for all constituents, and then these solutions were used as the initial condition for solving the time-dependent equations. The effects of chemical reactions and vertical eddy diffusion transport on the CH4 and CO concentrations above 40 km, where one-dimensional treatment may not be too inappropriate, are discussed. Some of the 15 constituents included in the computer program interact with CO and CH4 chemically, and CH4 is more sensitive to the eddy diffusion coefficient than any other constituent. Therefore it is easier to determine the most appropriate profile of the eddy diffusion coefficient by comparing the model profile with the observations for CH4.

Shimazaki, T.↗