Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Assessment of model assumptions and budget terms of the unsteady flamelet equations for a turbulent reacting jet-in-cross-flow

An a priori analysis of the flamelet model for diffusion flames is conducted to systematically assess model assumptions that are associated with the asymptotic expansion, the omission of higher-order expansion terms, the consideration of preferential diffusion effects, and the one-dimensional flamelet representation. For this, a recent direct numerical simulation database of a reacting hydrogen/air jet-in-cross-flow (JICF) by Grout et al. [15,16] is used. The full flamelet equation for temperature, exact to the order of the Eulerian transport equation and general to different definitions of the mixture fraction, is derived. Analysis of mixture fraction conditioned profiles of temperature and scalar dissipation rate along the jet trajectory identified different ignition and flame stabilization scenarios on the windward and leeward sides of the reacting JICF. A balance analysis of the temperature flamelet equation is conducted to quantify contributions in flame-aligned and flame-orthogonal directions. Consistent with the flamelet assumption, it is shown that terms arising from scalar diffusion, heat-release, and species-diffusion-induced enthalpy flux are the dominant contributions, while the flame-aligned preferential diffusion promotes entrainment of heat into the flamelet structure. Finally, for the current JICF-configuration, it is found that contributions along the flame-orthogonal direction are on average negligible.

42 ENGINEERING↗
DOE OSTI2025

Engineering quorum-sensing circuits in Synechococcus elongatus PCC 7942 towards self-inducible systems

Despite significant potential for cyanobacteria as sustainable bioproduction chases, there are limited examples of scaled cyanobacterial bioproduction. In part, this is because most cyanobacterial species are poorly adapted to bioreactor cultivation conditions and lack features that facilitate biomass growth and harvesting at scale. We explored quorum sensing (QS) pathways derived from heterotrophic microbes as a method for autoinduction of gene expression circuits coordinated to population density in cyanobacteria. Here, we integrated genetic modules designed to produce and detect the diffusible QS signal, acyl-homoserine lactones (AHLs), in the cyanobacterial model, Synechococcus elongatus PCC 7942 (S. elongatus). We demonstrate that S. elongatus heterologously produces sufficient AHL signals to activate gene expression in a dose-dependent and population density-responsive manner. A hybrid combination of AHL synthesis enzyme from Vibrio fischeri (Lux system) with the transcription factor receiver from Pseudomonas aeruginosa (Las system) provides an ideal activation ratio and mitigates toxicity observed with some AHL systems. As a proof of concept, we coupled the QS pathway to the expression of a cell division inhibitory gene, cdv3, facilitating late-phase cell elongation, cell sedimentation, and improved biomass recovery. Our findings provide a foundation for the development of auto-induction systems leverageable to improve cyanobacterial biotechnology applications.

Acyl homoserine lactones (AHL)↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Searching for synchrotron emission from the geminga TeV halo using the planck satellite

Pulsars convert a significant fraction of their total spin-down power into very high-energy electrons, leading to the formation of TeV halos. While these halos are well characterized at TeV energies, it remains unclear whether pulsars also accelerate electrons efficiently at lower energies and how these particles propagate through their surrounding environments. We aim to test whether synchrotron emission from $\sim 50$–$300 \, \textrm{GeV}$ electrons around the Geminga pulsar can be detected in the frequency range observed by the Planck satellite. This would help constrain low-energy particle acceleration and diffusion in the vicinity of pulsars. We model the expected synchrotron emission from Geminga’s TeV halo based on various diffusion and injection spectrum scenarios and compare these predictions to publicly available multi-frequency Planck data. We find no conclusive evidence of spatially extended synchrotron emission associated with Geminga in any of Planck’s frequency bands. Our calculations show that even under favorable diffusion and injection conditions, the predicted synchrotron flux lies well below Planck’s measured background levels.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Interfacial Effects on Transport Coefficient Measurements in Li-ion Battery Electrolytes

Development of Li + -containing electrolytes with improved transport properties requires reliable, reproducible, and ideally low volume techniques to rigorously understand ion-transport with varying composition. Precisely measuring the complete set of transport coefficients in liquid electrolytes under battery-relevant operating conditions is difficult and the reliability of these methods are sparsely described in electrolyte transport literature. In this work, we apply a potentiostatic polarization-based transport characterization approach typically used for polymer electrolytes to liquid electrolyte systems in an attempt to fully measure all transport coefficients (conductivity, total salt diffusion coefficient, thermodynamic factor and transference number) for the model system of LiPF 6 in an ethylene carbonate—ethyl methyl carbonate (EC:EMC) mixture. Using systematic timescale and statistical analyses, we find that transport coefficients measured using potentiostatic polarization of Li-Li symmetric cells exhibit strong correlation to Li electrode interfacial resistance, indicating that such methods are probing both bulk and interfacial phenomena. This reveals a major roadblock in characterizing electrolyte systems where the interfacial resistance is significantly larger than ohmic electrolyte resistance. As a result, we find that methods that rely on potentiostatic Li metal stripping/plating do not readily result in reliable liquid electrolyte transport coefficients, unlike similar methods for solid polymer electrolytes, where interfacial resistances are typically smaller than electrolyte resistances at the elevated temperatures typically of interest for such electrolytes.

25 ENERGY STORAGE↗
NASA NTRS1991

Silsesquioxane-derived ceramic fibres

Fibers formed from blends of silsesquioxane polymers were characterized to study the pyrolytic conversion of these precursors to ceramics. The morphology of fibers pyrolyzed to 1400 C revealed primarily amorphous glasses whose conversion to beta-SiC is a function of both blend composition and pyrolysis conditions. Formation of beta-SiC crystallites within the glassy phase is favored by higher than stoichiometric C/Si ratios, while carbothermal reduction of Si-O bonds to form SiC with loss of SiO and CO occurs at higher methyl/phenylpropyl silsesquioxane (lower C/Si) ratios. As the carbothermal reduction is assumed to be diffusion controlled, the fibers can serve as model systems to gain understanding of the silsesquioxane pyrolysis behavior, and therefore are useful in the development of polysilsesquioxane-derived ceramic matrices and coatings as well.

Hurwitz, F. I.↗
NASA NTRS2003

Transient Interfacial Phenomena in Miscible Polymer Systems (TIPMPS)

Almost one hundred years ago Korteweg published a theory of how stresses could be induced in miscible fluids by concentration gradients, causing phenomena that would appear to be the same as with immiscible fluids. Miscible fluids could manifest a transient or effective interfacial tension (EIT). To this day, there has been no definitive experiment to confirm Korteweg's model but numerous fascinating and suggestive experiments have been reported. The goal of TIPMPS is to answer the question: Can concentration and temperature gradients in miscible materials induce stresses that cause convection? Many polymer processes involving miscible monomer and polymer systems could be affected by fluid flow and so this work could help understand miscible polymer processing, not only in microgravity, but also on earth. Demonstrating the existence of this phenomenon in miscible fluids will open up a new area of study for materials science. The science objectives of TIPMPS are: (1) Determine if convection can be induced by variation of the width of a miscible interface; (2) Determine if convection can be induced by variation of temperature along a miscible interface; (3) Determine if convection can be induced by variation of conversion along a miscible interface An interface between two miscible fluids can best be created via a spatially-selective photopolymerization of dodecyl acrylate with a photoinitiator, which allows the creation of precise and accurate concentration gradients between polymer and monomer. Optical techniques will be used to measure the refractive index variation caused by the resultant temperature and concentration fields. The viscosity of the polymer will be measured from the increase in the fluorescence of pyrene. Because the large concentration and temperature gradients cause buoyancy-driven convection that prevents the observation of the predicted flows, the experiment must be done in microgravity. In this report, we will consider our efforts to estimate the square gradient parameter, k, and our use of the estimates in modeling of the planned TIPMPS experiments. We developed a model consisting of the heat and diffusion equations with convective terms and of the Navier-Stokes equations with an additional volume force written in the form of the Korteweg stresses arising from nonlocal interaction in the fluid. The fluid's viscosity dependence on polymer conversion and temperature was taken from measurements of poly(dodecyl acrylate). Numerical modeling demonstrated that significant flows would arise for conditions corresponding to the planned experiments.

Pojman, John A.↗
DOE OSTI2018

Measurement and modeling of the multiwavelength optical properties of uncoated flame-generated soot

Optical properties of flame-generated black carbon (BC) containing soot particles were quantified at multiple wavelengths for particles produced using two different flames: a methane diffusion flame and an ethylene premixed flame. Measurements were made for (i) nascent soot particles,(ii) thermally denuded nascent particles, and (iii) particles that were coated and then thermally denuded, leading to the collapse of the initiallylacy, fractal-like morphology. The measured mass absorption coefficients(MACs) depended on soot maturity and generation but were similar between flames for similar conditions. For mature soot, here corresponding to particles with volume-equivalent diameters > ~ 160 nm, the MAC and absorption Ångström exponent (AAE) values were independent of particle collapse while the single-scatter albedo increased. The MAC values for these larger particles were also size-independent. The mean MAC value at 532 nm for larger particles was 9.1±1.1 m 2 g -1 , about 17 % higher than that recommended by Bond and Bergstrom (2006), and the AAE was close to unity. Effective, theory-specific complex refractive index(RI) values are derived from the observations with two widely used methods:Lorenz–Mie theory and the Rayleigh–Debye–Gans (RDG) approximation. Mietheory systematically underpredicts the observed absorption cross sections at all wavelengths for larger particles (with x >0.9 ) independent of the complex RI used, while RDG provides good agreement. (The dimensionless size parameter x = π d p / λ , where d p is particle diameter and λ is wavelength.) Importantly, this implies that the use of Mie theory within air quality and climate models, as is common, likely leads to under predictions in the absorption by BC, with the extent of underprediction depending on the assumed BC size distribution and complex RI used. We suggest that it is more appropriate to assume a constant,size-independent (but wavelength-specific) MAC to represent absorption by uncoated BC particles within models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Artificial ground reflector size and position effects on energy yield and economics of single‐axis‐tracked bifacial photovoltaics

Abstract Artificial ground reflectors improve bifacial energy yield by increasing both front and rear‐incident irradiance. Studies have demonstrated an increase in energy yield due to the addition of artificial reflectors; however, they have not addressed the effect of varying reflector dimensions and placement on system performance and the impact of these parameters on the reflectors' financial viability. We studied the effect of high albedo (70% reflective) artificial reflectors on single‐axis‐tracked bifacial photovoltaic systems through ray‐trace modeling and field measurements. In the field, we tested a range of reflector configurations by varying reflector size and placement and demonstrated that reflectors increased daily energy yield up to 6.2% relative to natural albedo for PERC modules. To confirm the accuracy of our model, we compared modeled and measured power and found a root mean square error (RMSE) of 5.4% on an hourly basis. We modeled a typical meteorological year in Golden, Colorado, to demonstrate the effects of artificial reflectors under a wide range of operating conditions. Seventy percent reflective material can increase total incident irradiance by 1.9%–8.6% and total energy yield by 0.9%–4.5% annually after clipping is considered with a DC–AC ratio of 1.2. Clipping has a significant effect on reflector impact and must be included when assessing reflector viability because it reduces reflector energy gain. We calculated a maximum viable cost for these improvements of up to $2.50–4.60/m 2 , including both material and installation, in Golden. We expanded our analysis to cover a latitude range of 32–48°N and demonstrated that higher‐latitude installations with lower energy yield and higher diffuse irradiance content can support higher reflector costs. In both modeling and field tests, and for all locations, the ideal placement of the reflectors was found to be directly underneath the module due to the optimized rear irradiance increase.

14 SOLAR ENERGY↗
DOE OSTI2023

Effects of surface vacancy-adatom pair formation on PFC tungsten surface morphological response

Abstract We report a simulation study of the effect of He-irradiation-induced surface vacancy-adatom pair formation on the surface morphological evolution of plasma-facing component (PFC) tungsten and examine a number of factors that impact such evolution. Our analysis is based on self-consistent dynamical simulations according to an atomistically-informed, continuum-scale surface evolution model that has been developed following a hierarchical multiscale modeling strategy and can access the spatiotemporal scales of relevance to fuzz formation. The model accounts for the flux of surface adatoms generated as a result of the surface vacancy-adatom pair formation effect upon He implantation, which contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. We have carried out atomic-scale computations of optimal diffusion pathways along and across island step edges on the W(110) surface and calculated Ehrlich–Schwoebel (ES) barriers in adatom diffusion along and across such step edges. This aspect of surface adatom diffusion contributes to anisotropic surface atomic fluxes, terrace and step diffusive currents, and has been incorporated into our PFC surface evolution model, which predicts the formation of preferentially aligned nanoridge stripe patterns on the PFC surface. We establish that these anisotropic diffusive currents accelerate nanotendril growth on the PFC surface and the onset of surface nanostructure pattern formation. We also explore systematically the dependence of the PFC surface morphological response on the surface temperature and He ion incident flux, characterize in detail the resulting surface topographies and growth kinetics, and compare the predicted surface morphologies with experimental observations. Our simulation predictions for the emerging surface nanostructure patterns under certain plasma exposure conditions are consistent with experimental findings in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

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

Entropy modeling on serrated flows in carburized steels

Samples comprised of carburized steel alloys, 4615, 4720, and 10B22, were tension tested at room temperature and at a strain rate of 2 × 10 -4 s -1 . The serrated-flow behavior was modeled and analyzed, using the refined composite multiscale entropy (RCMSE) method. High energy X-ray diffraction revealed the body-centered-cubic (BCC) structure of the three alloys in the as-received condition. Moreover, Vickers-hardness experiments were performed, which showed an increase in the hardness of the alloy that corresponded to increasing impurity concentrations. Here the diffusion coefficient of carbon in the three alloys were estimated, respectively, to be 2.5 × 10 -11 m 2 /s, 5.1 × 10 -11 m 2 /s, and 2.8 × 10 -11 m 2 /s for alloys, 10B22, 4615, and 4720, respectively, which agree with the values found in the literature for α-iron. Moreover, the sample entropy, on average, increased with respect to the heating time, and, consequently, the amount of carbon that had diffused into the alloy. Finally, the above results show that the complexity of the serration behavior in the alloys studied increases with respect to the amount of carbon that has diffused into the material.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Novel vertical silicon photodiodes based on salicided polysilicon trenched contacts

The classical concept of silicon photodiodes comprises of a planar design characterized by heavily doped emitters. Such geometry has low collection efficiency of the photons absorbed close to the surface. An alternative, promising, approach is to use a vertical design. Nevertheless, realization of such design is technologically challenged, hence hardly explored. Herein, a novel type of silicon photodiodes, based on salicided polysilicon trenched contacts, is presented. These contacts can be prepared up to 10 μm in depth, without showing any leakage current associated with the increase in the contact area. Consequently, the trenched photodiodes revealed better performance than no-trench photodiodes. A simple two dimensional model was developed, allowing to estimate the conditions under which a vertical design has the potential to have better performance than that of a planar design. At large, the deeper the trench is, the better is the vertical design relative to the planar (up to 10 μm for silicon). The vertical design is more advantageous for materials characterized by short diffusion lengths of the carriers. Salicided polysilicon trenched contacts open new opportunities for the design of solar cells and image sensors. For example, these contacts may passivate high contact area buried contacts, by virtue of the conformity of polysilicon interlayer, thus lowering the via resistance induced recombination enhancement effect.

14 SOLAR ENERGY↗
DOE OSTI2019

Formulation analysis and computation of an optimization-based local-to-nonlocal coupling method

We present an optimization-based coupling method for local and nonlocal continuum models. Our approach couches the coupling of the models into a control problem where the states are the solutions of the nonlocal and local equations, the objective is to minimize their mismatch on the overlap of the local and nonlocal problem domains, and the virtual controls are the nonlocal volume constraint and the local boundary condition. We present the method in the context of Local-to-Nonlocal diffusion coupling. Numerical examples illustrate the theoretical properties of the approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Mass transfer coefficients, kL, and air-CO2 ingassing rates in 3.4 m2 and 1-acre raceway ponds.

The US DOE 2016 Billion Ton Update (Langholtz 2016) projected that CO2 in flue-gas from power plants or similar sources would limit U.S. algal biofuel potential to under 5 billion gallons gasoline equivalent (gge) per year, due to lack of sufficient land, water and other requirements near such flue gas sources. 2nd generation carbon capture technologies are proposed to overcome this constraint (Davis 2017), and could expand this resource potential nearly 10-fold, assuming CO2 costs near the future, 2025 DOE NETL, goal of $40/metric ton CO2. An alternative pathway is the direct uptake of air CO2 into the algal ponds through an increase in the air-CO2 transfer rate resulting from the reaction of CO2 with hydroxide ions. However, such a ‘chemical enhancement’ in mass transfer depends on a high pH in the culture, and, by extension, the ability of microalgae to maintain high rates of carbon assimilation under such conditions. Abiotic experiments characterized the air-CO2 mass-transfer rate in 1-acre and 3.4 m2 raceway ponds as a function of pH, with rates approaching 10 g C/m2-day at pH 12 in the brackish water tested, equivalent to a biomass productivity potential of near 20 g AFDW/m2-day (@ 0.47 g C/g AFDW). The mass transfer rate was found to be nearly independent of the mass-transfer coefficient, indicating that at turbulence levels achievable in typical raceway ponds, the system is reaction rate, rather than diffusion limited. In biotic trials with an unknown green microalga, biomass productivity in ponds receiving CO2 from only direct surface air-CO2 exchange averaged 5 +/- 0.5 g AFDW/m2-day in August - September central California conditions, about half-that of experimental controls fertilized with supplemental CO2. pH in the air-only treatment reached a maximum of 10.5, supporting a model predicted air-CO2 ingassing rate between 2-3 g C/m2-day, consistent with the observed biomass carbon assimilation rate. Results suggest that unique alkaliphilic strains are required if reliant on direct in-pond air-CO2 transfer for algal inorganic carbon supply, and that reaching an economically viable biomass productivity will require strains that thrive at pH 11 and above.

09 BIOMASS FUELS↗
NASA NTRS1978

The use of airborne and ground-based measurements of launch vehicle HCl puffs in a finite difference parameter estimation procedure

Model 4 of the multilayer diffusion models, version 5 (MDMV5) is discussed. MDMV5 is designed to help factor considerations of the potential environment impact of a given launch operation into the actual launch decision process. The model supports this decision process by providing quantitative estimates of particular air quality indices at a set of spatial positions in the launch area and surrounding communities, for specified meteorological conditions. Modeling concepts, the finite difference estimation algorithm, and data application and results are examined.

Englar, T. S., Jr.↗
NASA NTRS1999

A Process-based, Climate-Sensitive Model to Derive Methane Emissions from Natural Wetlands: Application to 5 Wetland Sites, Sensitivity to Model Parameters and Climate

Methane emissions from natural wetlands constitutes the largest methane source at present and depends highly on the climate. In order to investigate the response of methane emissions from natural wetlands to climate variations, a 1-dimensional process-based climate-sensitive model to derive methane emissions from natural wetlands is developed. In the model the processes leading to methane emission are simulated within a 1-dimensional soil column and the three different transport mechanisms diffusion, plant-mediated transport and ebullition are modeled explicitly. The model forcing consists of daily values of soil temperature, water table and Net Primary Productivity, and at permafrost sites the thaw depth is included. The methane model is tested using observational data obtained at 5 wetland sites located in North America, Europe and Central America, representing a large variety of environmental conditions. It can be shown that in most cases seasonal variations in methane emissions can be explained by the combined effect of changes in soil temperature and the position of the water table. Our results also show that a process-based approach is needed, because there is no simple relationship between these controlling factors and methane emissions that applies to a variety of wetland sites. The sensitivity of the model to the choice of key model parameters is tested and further sensitivity tests are performed to demonstrate how methane emissions from wetlands respond to climate variations.

Walter, Bernadette P.↗
NASA NTRS2004

Current State of Modeling the Photochemistry of Titan's Mutually Dependent Atmosphere and Ionosphere

In the context of recent observations, microphysical models, and laboratory data, a photochemical model of Titan's atmosphere, including updated chemistry focusing on rate coefficients and cross sections measured under appropriate conditions, has been developed to increase understanding of these processes and improve upon previous Titan photochemical models. The model employs a two-stream discrete ordinates method to characterize the transfer of solar radiation, and the effects of electron-impact, cosmic-ray deposition, and aerosol opacities from fractal and Mie particles are analyzed. Sensitivity studies demonstrate that an eddy diffusion profile with a homopause level of 850 km and a methane stratospheric mole fraction of 2.2% provides the best fit of stratospheric and upper atmosphere observations and an improved fit over previous Titan photochemical models. Lack of fits for C3H8, HC3N, and possibly C2H3CN can be resolved with adjustments in aerosol opacity. The model presents a benzene profile consistent with its detection in Titan's stratosphere [Coustenis et al., 2003], which may play an important role in the formation of Titan hazes. An electron peak concentration of 4200 cm(exp -3) is calculated, which exceeds observations by 20%, considerably lower than previous ionosphere models. With adjustments in aerosol opacities and surface fluxes the model illustrates that reasonable fits to existing observations are possible with a single eddy diffusion profile, contrary to the conclusions of previous Titan models. These results will aid in the receipt and interpretation of data from Cassini-Huygens, which will arrive at Titan in 2004 and deploy a probe into Titan's atmosphere in January 2005.

ion chemistry↗