Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model species”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Inferring the abundances of ClO and HO2 from Spacelab 3 atmospheric trace molecule spectroscopy observations

The vertical distributions of the important highly reactive stratospheric species, ClO and HO2, have been inferred from Spacelab 3 (May 1985) Atmospheric Trace Molecule Spectroscopy (ATMOS) measurements of more detectable radical and reservoir species. A simple steady state algebraic expression for ClO, utilizing the observed ClONO2/NO2 abundance ratio, approximates the ClO results of a time-dependent photochemical model at sunset (30 deg N). Balloon measurements of ClO and comparisons of time-dependent photochemical model calculations of ClONO2 and HCl with the corresponding ATMOS profiles suggest that the actual ClO values are less than the time-dependent model profiles for ClO at sunset and sunrise (47 deg S). Errors in the current model simulation of the partitioning among the principal free chlorine species (HCl, ClONO2, and ClO) are indicated. The accuracy of the time-dependent model calculations of HO2 is suggested by the agreement between model results and ATMOS observations for HO2NO2 and balloon measurements of HO2. This confirms for the first time the procedure suggested previously by a number of authors of deriving HO(x) abundances from observed fields of O3 and H2O.

Allen, Mark↗
DOE OSTI2021

A novel physics-based and data-supported microstructure model for part-scale simulation of laser powder bed fusion of Ti-6Al-4V

The elasto-plastic material behavior, material strength and failure modes of metals fabricated by additive manufacturing technologies are significantly determined by the underlying process-specific microstructure evolution. In this work a novel physics-based and data-supported phenomenological microstructure model for Ti-6Al-4V is proposed that is suitable for the part-scale simulation of laser powder bed fusion processes. The model predicts spatially homogenized phase fractions of the most relevant microstructural species, namely the stable β-phase, the stable α s -phase as well as the metastable Martensite α m -phase, in a physically consistent manner. In particular, the modeled microstructure evolution, in form of diffusion-based and non-diffusional transformations, is a pure consequence of energy and mobility competitions among the different species, without the need for heuristic transformation criteria as often applied in existing models. The mathematically consistent formulation of the evolution equations in rate form renders the model suitable for the practically relevant scenario of temperature- or time-dependent diffusion coefficients, arbitrary temperature profiles, and multiple coexisting phases. Due to its physically motivated foundation, the proposed model requires only a minimal number of free parameters, which are determined in an inverse identification process considering a broad experimental data basis in form of time-temperature transformation diagrams. Subsequently, the predictive ability of the model is demonstrated by means of continuous cooling transformation diagrams, showing that experimentally observed characteristics such as critical cooling rates emerge naturally from the proposed microstructure model, instead of being enforced as heuristic transformation criteria. Eventually, the proposed model is exploited to predict the microstructure evolution for a realistic selective laser melting application scenario and for the cooling/quenching process of a Ti-6Al-4V cube of practically relevant size. Numerical results confirm experimental observations that Martensite is the dominating microstructure species in regimes of high cooling rates, e.g., due to highly localized heat sources or in near-surface domains, while a proper manipulation of the temperature field, e.g., by preheating the base-plate in selective laser melting, can suppress the formation of this metastable phase.

Inverse parameter identification↗
DOE OSTI2021

Test of model of equivalence of tree height growth and transpiration rates in percolation-based phenomenology for root-soil interaction

We recently introduced a phenomenological model to represent root-soil interactions that accounts for variability in tree growth on local scales and within species (Hunt et al., 2020). Here, using publicly available datasets, we examine the developed model to provide additional empirical evidence, which also gives novel insight into two particular tree species, Eucalyptus regnans and Eucalyptus sieberri. In particular, we discovered two novel aspects of the model: a power-law time dependence with a different power-law exponent than usually heretofore applied, and apparently a direct equivalence of results for transpiration and growth rates. Also, although our original expectation was that growth and transpiration were merely proportional, comparison of the entire BAAD data set (Falster et al., 2015) for tree height with our model (Hunt et al., 2020; Fig. 1) already suggested approximate equality of these disparate measurements (volume/area, height). In this present manuscript we find further evidence for equivalence of the two measurements and provide additional theoretical basis for a model assumption regarding constancy of sap flow rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part I: m-cresol purple properties from 25-75 °C and Er hydrolysis

The speciation and mobility of rare earth elements (REE) strongly depends on pH which controls the formation of charged aqueous hydroxyl species. The latter potentially play an important role in controlling heavy REE adsorption on clay minerals in near-neutral to alkaline waters such as in regolith-hosted REE mineral deposits. However, accurate REE hydrolysis constants are needed for developing geochemical models that can predict the role of these charged species in natural systems. Here, we develop a robust experimental UV-Vis spectrophotometric method using m-cresol purple to determine in situ pH from 25 to 75 °C. This method is used to derive the average ligand number and hydrolysis constants of erbium (Er) at 25 °C in aqueous solutions with low ionic strength (≤ 0.001 mol/L) at pH from ~7 to 9.5 and in the presence of Er concentrations from 0 to 0.057 mM. The average ligand number ranges between 1 and 3 indicating that Er(OH) 2+ , Er(OH) 2 + and Er(OH) 3 0 control speciation in the experiments. The logarithm of the Er hydrolysis constants (log*β n ° , n= 1 to 3) derived at infinite dilution for the reaction Er 3+ + nH 2 O = Er(OH) n 3-n + nH + are: *β 1 ° = –7.22 ± 0.10, *β 2 ° = –14.52 ± 0.08, *β 3 ° = –23.24 ± 0.04. Implementation of these experimental data into a geochemical model indicates that the Er(OH) 2 + and Er(OH) 3 0 species are both stable in a much wider pH range than previously predicted. Consequently, the positively charged REE hydroxyl complexes can potentially control the fractionation of light vs. heavy REE via adsorption as observed in the formation of certain regolith-hosted REE deposits.

58 GEOSCIENCES↗
NASA NTRS2011

Large Eddy Simulation Study for Fluid Disintegration and Mixing

A new modeling approach is based on the concept of large eddy simulation (LES) within which the large scales are computed and the small scales are modeled. The new approach is expected to retain the fidelity of the physics while also being computationally efficient. Typically, only models for the small-scale fluxes of momentum, species, and enthalpy are used to reintroduce in the simulation the physics lost because the computation only resolves the large scales. These models are called subgrid (SGS) models because they operate at a scale smaller than the LES grid. In a previous study of thermodynamically supercritical fluid disintegration and mixing, additional small-scale terms, one in the momentum and one in the energy conservation equations, were identified as requiring modeling. These additional terms were due to the tight coupling between dynamics and real-gas thermodynamics. It was inferred that if these terms would not be modeled, the high density-gradient magnitude regions, experimentally identified as a characteristic feature of these flows, would not be accurately predicted without the additional term in the momentum equation; these high density-gradient magnitude regions were experimentally shown to redistribute turbulence in the flow. And it was also inferred that without the additional term in the energy equation, the heat flux magnitude could not be accurately predicted; the heat flux to the wall of combustion devices is a crucial quantity that determined necessary wall material properties. The present work involves situations where only the term in the momentum equation is important. Without this additional term in the momentum equation, neither the SGS-flux constant-coefficient Smagorinsky model nor the SGS-flux constant-coefficient Gradient model could reproduce in LES the pressure field or the high density-gradient magnitude regions; the SGS-flux constant- coefficient Scale-Similarity model was the most successful in this endeavor although not totally satisfactory. With a model for the additional term in the momentum equation, the predictions of the constant-coefficient Smagorinsky and constant-coefficient Scale-Similarity models were improved to a certain extent; however, most of the improvement was obtained for the Gradient model. The previously derived model and a newly developed model for the additional term in the momentum equation were both tested, with the new model proving even more successful than the previous model at reproducing the high density-gradient magnitude regions. Several dynamic SGS-flux models, in which the SGS-flux model coefficient is computed as part of the simulation, were tested in conjunction with the new model for this additional term in the momentum equation. The most successful dynamic model was a "mixed" model combining the Smagorinsky and Gradient models. This work is directly applicable to simulations of gas turbine engines (aeronautics) and rocket engines (astronautics).

Bellan, Josette↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

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

A model of Triton's role in Neptune's magnetosphere

Escape of neutral hydrogen (H) and nitrogen (N) from Triton's maintains a large neutral cloud, called the Triton torus, in Neptune's magnetosphere. We have developed the first detailed Monte Carlo simulation model of the Triton torus that includes the collisionality, the complex geometry, the injection of two neutral species from Triton (H and N), and the combined effects of photoionization, electron impact ionization, and charge exchange. Ionization in Neptune's plasma sheet was modeled using Voyager plasma observations. Collisions cause both the H and N neutral clouds to become more radially extended, both toward Neptune and out beyond the magnetopause, as well as more extended in latitude, when compared with collisionless models. Moreover, collisions of H with the much more massive N greatly enhance the collisional ejection of H from the system and into Neptune's atmosphere. This effect decreases the probability of H ionization within the magnetosphere relative to that for N, and furthermore causes model results for two-species injection from Triton to differ significantly from those for H injection alone. For a hydrogen escape rate from Triton of 5 x 10(exp 25)/s, as given by photo-chemical models of Triton's upper atmosphere, a nitrogen escape rate of 5 x 10(exp 24)/s gives proton and N(+) sources of 5.6 x 10(exp 24)/s and 3.3 x 10(exp 24)/s, respectively, whose ratio is close to the observed ratio of protons to heavies. A nitrogen escape rate of 2 x 10(exp 25)/s, yields an N(+) source more than twice that of protons, inconsistent with the Voyager data.

Decker, R. B.↗
NASA NTRS1997

A Three-Wave Model of the Stratosphere with Coupled Dynamics, Radiation and Photochemistry

A zonal mean chemistry transport model (2-D CTM) coupled with a semi-spectral dynamical model is used to simulate the distributions of trace gases in the present day atmosphere. The zonal-mean and eddy equations for the velocity and the geopotential height are solved in the semi-spectral dynamical model. The residual mean circulation is derived from these dynamical variables and used to advect the chemical species in the 2- D CTM. Based on a linearized wave transport equation, the eddy diffusion coefficients for chemical tracers are expressed in terms of the amplitude, frequency and growth rate of dynamical waves; local chemical loss rates; and a time constant parameterizing small scale mixing. The contributions to eddy flux are from the time varying wave amplitude (transient eddy), chemical reactions (chemical eddy) and small scale mixing. In spite of the high truncation in the dynamical module (only three longest waves are resolved), the model has simulated many observed characteristics of stratospheric dynamics and distribution of chemical species including ozone. Compared with the values commonly used in 2-D CTMs, the eddy diffusion coefficients for chemical species calculated in this model are smaller, especially in the subtropics. It is also found that the chemical eddy diffusion has only a small effects in determining the distribution of most slow species, including ozone in the stratosphere.

Shia, Run-Lie↗
DOE OSTI2023

Electrochemical Performance of Mixed Redox-Active Organic Molecules in Redox Flow Batteries

Designing electrolytes based on mixture of different organic redox active molecules brings the opportunity of enhancing the volumetric energy density of flow batteries and removes the requirement of high solubility for individual organic species in the mixture. Here, in the present work, we conduct computational and experimental analysis to investigate the electrochemical performance of mixed redox-active organic molecules. A zero-dimensional transient model is employed to investigate the changes in the half-cell potential and the concentrations and partial currents of individual redox reactions in a mixture of organic molecules over time. The model demonstrates the effects of individual properties of species such as kinetic rate constants, mass transfer coefficients, concentration ratios and standard redox potentials and reports the effect of energy-losing homogenous chemical redox reaction on the voltage efficiency and concentration ratios of the mixed species. Pairs of anthraquinone negolyte species were selected for an experimental case study. A mixture of 2,6-N-TSAQ and 2,6-DHAQ showed 40% increase in the volumetric energy density compared to the performance of 2,6-DHAQ alone. Based on the results of the experimental and computational analysis, we propose guidelines for the design of suitable mixed redox-active organic species.

25 ENERGY STORAGE↗
DOE OSTI2024

Molecular profiles of blood from numerous species that differ in sensitivity to acute inflammation

Abstract Vertebrates differ over 100,000-fold in responses to pro-inflammatory agonists such as bacterial lipopolysaccharide (LPS), complicating use of animal models to study human sepsis or inflammatory disorders. We compared transcriptomes of resting and LPS-exposed blood from six LPS-sensitive species (rabbit, pig, sheep, cow, chimpanzee, human) and four LPS-resilient species (mice, rats, baboon, rhesus), as well as plasma proteomes and lipidomes. Unexpectedly, at baseline, sensitive species already had enhanced expression of LPS-responsive genes relative to resilient species. After LPS stimulation, maximally different genes in resilient species included genes that detoxify LPS, diminish bacterial growth, discriminate sepsis from SIRS, and play roles in autophagy and apoptosis. The findings reveal the molecular landscape of species differences in inflammation. This may inform better selection of species for pre-clinical models and could lead to new therapeutic strategies that mimic mechanisms in inflammation-resilient species to limit inflammation without causing immunosuppression.

Biochemistry & Molecular Biology↗
DOE OSTI2023

Pyrolysis of Cyclohexane and 1-Hexene at High Temperatures and Pressures—A Photoionization Mass Spectrometry Study

Cycloalkanes are important components of a wide range of fuels. However, there are few experimental data at simultaneously high temperatures and pressures similar to those found in practical systems. Such data are necessary for developing and testing chemical kinetic models. In this study, data relevant to cycloalkane pyrolysis were obtained from high repetition rate shock tube experiments coupled with synchrotron-based photoionization mass spectrometry diagnostics. The pyrolysis of cyclohexane was studied over 1270–1550 K and ~9 bar, while the more reactive primary decomposition product, 1-hexene, was studied at 1160–1470 K and ~5 bar. Insights into the decomposition of the parent molecules, the formation of primary products and the production of aromatic species were gained. Simulations were performed with models for cyclohexane and 1-hexene that were based on literature models. The results indicate that over several hundred microseconds reaction time at high pressures and temperatures the pyrolysis of cyclohexane is largely dominated by reactions initiated by cyclohexyl radicals. Furthermore, good agreement between the simulations and the experiments were observed for cyclohexane and 1-hexene with a modified version of the cyclohexane model. Conversely, the 1-hexene model did not reproduce the experimental observations.

1-alkenes↗
DOE OSTI2023

cTULIP: application of a human-based RNA-seq primary tumor classification tool for cross-species primary tumor classification in canine

The domestic dog, Canis familiaris, is quickly gaining traction as an advantageous model for use in the study of cancer, one of the leading causes of death worldwide. Naturally occurring canine cancers share clinical, histological, and molecular characteristics with the corresponding human diseases. In this study, we take a deep-learning approach to test how similar the gene expression profile of canine glioma and bladder cancer (BLCA) tumors are to the corresponding human tumors. We likewise develop a tool for identifying misclassified or outlier samples in large canine oncological datasets, analogous to that which was developed for human datasets. We test a number of machine learning algorithms and found that a convolutional neural network outperformed logistic regression and random forest approaches. We use a recently developed RNA-seq-based convolutional neural network, TULIP, to test the robustness of a human-data-trained primary tumor classification tool on cross-species primary tumor prediction. Our study ultimately highlights the molecular similarities between canine and human BLCA and glioma tumors, showing that protein-coding one-to-one homologs shared between humans and canines, are sufficient to distinguish between BLCA and gliomas. The results of this study indicate that using protein-coding one-to-one homologs as the features in the input layer of TULIP performs good primary tumor prediction in both humans and canines. Furthermore, our analysis shows that our selected features also contain the majority of features with known clinical relevance in BLCA and gliomas. Our success in using a human-data-trained model for cross-species primary tumor prediction also sheds light on the conservation of oncological pathways in humans and canines, further underscoring the importance of the canine model system in the study of human disease.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1986

Tracer transport by the diabatic circulation deduced from satellite observations

Nimbus-7 sensor data were used to track the diabatic circulation in the stratosphere to study the advective transport of CH4 and N2O as tracer species. Advective transport by the mean circulation was found to be a function of the temperature field and associated deviations from radiative equilibrium. A photochemical model was applied to account for the disappearance of the tracer species from the stratosphere. Comparisons between the SAMS data and modeling on the basis of the chemical loss rates of the tracers and the LIMS circulation data showed that the model predictions underestimated the resident abundances, although the global distributions and circulations exhibited a good match.

Solomon, S.↗
DOE OSTI2021

Testing effects of invasive fire ants and disturbance on ant communities of the longleaf pine ecosystem

Invasive species and habitat disturbance are among the most important drivers of biodiversity loss and ecological change. Their individual effects, however, are difficult to disentangle because invasion and disturbance are often intimately linked. Here we test alternative hypotheses to determine if the invasive red imported fire ant, Solenopsis invicta, is a ‘driver’ or simply a ‘passenger’ of ecological change in a longleaf pine ecosystem. We randomly assigned treatments of (1) unmanipulated, (2) soil disturbance, (3) fire ant removal and (4) soil disturbance + fire ant removal to experimental blocks and measured how ant communities changed over 2 years in thirty-six 15-m 2 plots. Additionally, fire ant abundance in removal plots averaged 42% lower in pitfall traps and 95% lower on baits compared to unmanipulated, control plots. Species richness of co-occurring ants also decreased 42% in removal plots, with significant changes in community composition. Soil disturbance alone did not affect ant communities. Fire ant diet breadth—measured using carbon and nitrogen stable isotopes—increased up to 4.7-fold in soil disturbance + removal plots (i.e. 0.84‰ 2 to 3.94‰ 2 ). While non-target impacts of the fire ant removal treatment complicate interpretation, our results suggest fire ants follow an alternative ‘interacting drivers’ model in which partial recovery of some species occurs when populations of an invasive species are reduced. Further recovery of native ants may be limited by persisting, landscape-level effects of fire ants suppressing co-occurring ants below historical levels.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Zavitsas’ hydration model for electrolytes only stable in the presence of another: NaAl(OH) 4 in aqueous NaOH solution

Many aqueous electrolytes are only stable in the presence of another electrolyte, such as electrolytes that are only soluble in strong acid or base. An example is sodium aluminate [NaAl(OH) 4 ], which is only stable in aqueous NaOH. This complicates developing thermodynamic parameters because it is difficult to separate the contributions of individual electrolytes in multicomponent solutions to measured bulk thermodynamic properties. The present study develops a method to determine the liquid phase parameters from solubility data for Zavitsas’ Hydration model, a model that incorporates hydration parameters into the activities of dissolved species. Furthermore, this study uses gibbsite [Al(OH) 3 ] solubility data in aqueous NaOH solution to develop Zavitsas' model parameters for NaAl(OH) 4 simultaneously with the equilibrium constants. A good fit of the solubility data was found, showing that Zavitsas’ model is effective for this system.

Bayer process↗
NASA NTRS2001

Retrieval of Temperature and Species Distributions from Multispectral Image Data of Surface Flame Spread in Microgravity

Weight, size, and power constraints severely limit the ability of researchers to fully characterize temperature and species distributions in microgravity combustion experiments. A powerful diagnostic technique, infrared imaging spectrometry, has the potential to address the need for temperature and species distribution measurements in microgravity experiments. An infrared spectrum imaged along a line-of-sight contains information on the temperature and species distribution in the imaged path. With multiple lines-of-sight and approximate knowledge of the geometry of the combustion flowfield, a three-dimensional distribution of temperature and species can be obtained from one hyperspectral image of a flame. While infrared imaging spectrometers exist for collecting hyperspectral imagery, the remaining challenge is retrieving the temperature and species information from this data. An initial version of an infrared analysis software package, called CAMEO (Combustion Analysis Model et Optimizer), has been developed for retrieving temperature and species distributions from hyperspectral imaging data of combustion flowfields. CAMEO has been applied to the analysis of multispectral imaging data of flame spread over a PMMA surface in microgravity that was acquired in the DARTFire program. In the next section of this paper, a description of CAMEO and its operation is presented, followed by the results of the analysis of microgravity flame spread data.

Annen, K. D.↗
DOE OSTI2022

(Invited) Continuum Mathematical Modeling of Water Electrolysis: A Tutorial

Widespread use of hydrogen energy is contingent on the development of reliable and economical sources of hydrogen. Electrolysis from renewably-derived low-carbon electricity is a potentially viable method of hydrogen generation. Prime among the electrolysis technologies are those utilizing ion-conducting polymers (ionomers) including proton-exchange-membrane water electrolyzer (PEMWE). However, these technologies need to exhibit increased efficiency, performance, and durability to become commercially viable. Like most electrochemical devices, PEMWEs involve multiple components (e.g., catalyst, ionomer, transport layers, membrane, plates) and multiple phases, with phenomena occurring across different time and length scales. Furthermore, it is difficult to experimentally probe many of the species and phenomena during operation. Thus, mathematical modeling at the continuum level has been an invaluable aid in exploring, understanding, and optimizing PEMWE cell and components. Furthermore, this is especially true in the highly coupled and complex physics and chemistries that occur with the membrane-electrode assembly (MEA). The physics in a typical volume-averaged non-isothermal model include multiphase transport in porous media, concentrated-solution theory, Ohm's Law, and Butler-Volmer kinetics, and ion, gas, and water transport in the ionomer.

Dizon, Arthur↗
NASA NTRS2000

The Dynamical Evolution of A Tubular Leonid Persistent Train

The dynamical evolution of the persistent train of a bright Leonid meteor was examined for evidence of the source of the luminosity and the physical conditions in the meteor path. The train consisted of two parallel somewhat diffuse luminous tracks, interpreted as the walls of a tube. A general lack of wind shear along the trail allowed these structures to remain intact for nearly 200 s, from which it was possible to determine that the tubular structure expanded at a near constant 10.5 m/s, independent of altitude between 86 and 97 km. An initial fast decrease of train intensity below 90 km was followed by an increase in intensity and then a gradual decrease at longer times, whereas at high attitudes the integrated intensity was nearly constant with time. These results are compared to a model that describes the dynamical evolution of the train by diffusion, following an initial rapid expansion of the hot gaseous trail behind the meteoroid. The train luminosity is produced by O ((sup 1)S) emission at 557 nm, driven by elevated atomic O levels produced by the meteor impact, as well as chemiluminescent reactions of the ablated metals Na and Fe with O3. Ozone is rapidly removed within the train, both by thermal decomposition and catalytic destruction by the metallic species. Hence, the brightest emission occurs at the edge of the train between outwardly diffusing metallic species and inwardly diffusing O3. Although the model is able to account plausibly for a number of characteristic features of the train evolution, significant discrepancies remain that cannot easily be resolved.

Jenniskens, Peter↗