Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density estimation”

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 937 records · Page 52

Deep VLA imaging of twelve extended 3CR quasars

A representative sample of 12 extended quasars from the 3CR catalog has been imaged at 4.9 GHz using the Very Large Array (VLA). These full synthesis observations typically achieve an rms noise of 20 micro-Jy per beam. Jets are detected on at least one side of every source. The jets are well collimated compared with those in less powerfull sources, but spreading is detected in most of them. The opening angles of several jets are not constant, but show recollimation after an initial regime of rapid spreading. Many of the jets contain quasiperiodic strings of knots, of which the knot closest to the central feature is usually the brightest (until the jet nears its hot spot). The degrees of linear polarization at the jet knots range from less than 5% to approx. 50%, but show no common trend with distance along the jets. In knots that are elongated in directions close to that of the jet, the E vectors tend to be orthogonal to the jet axis. The prominence of the inner, straighter jet segments relative to the extended lobes correlates significantly with the prominence of the milliarcsecond-scale central features, but the prominence of the more bent jet segments does not. Candidates for counterjet emission are detected in seven sources, but there is no unambiguous, continuous counterjet in any of them. Estimates of the flux density ratios between the straighter jet segments and the counterjets based on these tentative detections range from 1.2:1 to greater than 175:1. There is no evidence in this sample that counterjet detectability correlates with such putative inclination indicators as central feature prominence or projected linear size. There is also no evidence that the prominence of the counterjets anticorrelates with that of the jets as predicted by simple relativistic-beaming models for the jet/counterjet asymmetry. There is, however, strong evidence that large bends in the main jet favor counterjet detection, and there are no counterjet candidates opposite long, uninterrupted straight segments of the main jets. The detectability of the counterjets in these quasars may therefore be strongly influenced by interactions between the underlying beams and inhomogeneities in the surrounding material. We offer a new empirical definition of the term 'hot spot' that is intended to improve the distinction between such features and 'jet knots.' Both the compactness of hot spots and their position in the lobe are affected by whether they are fed by a detectable jet. When the hot spots differ significantly in compactness, the more compact one is always on the jetted side. Jetted hot spots are also more likely to be recessed deeply from the outer edge of their lobes than are their counterjetted counterparts. The jetted hot spot is less prominent relative to other extended emission if the jet bends through a large angle, particularly if a large bend occurs abruptly. The counterjetted hot spot is also less well defined if the jet is more bent. The lobes of several sources show considerable inhomogeneity, including filamentation.

Bridle, Alan H.↗

What we have learned about Mars from SNC meteorites

The SNC meteorites are thought to be igneous martian rocks, based on their young crystallization ages and a close match between the composition of gases implanted in them during shock and the atmosphere of Mars. A related meteorite, ALH84001, may be older and thus may represent ancient martian crust. These petrologically diverse basalts and ultramafic rocks are mostly cumulates, but their parent magmas share geochemical and radiogenic isotopic characteristics that suggest they may have formed by remelting the same mantle source region at different times. Information and inferences about martian geology drawn from these samples include the following: Planetary differentiation occured early at approximately 4.5 GA, probably concurrently with accretion. The martian mantle contains different abundances of moderately volatile and siderophile elements and is more Fe-rich than that of the Earth, which has implications for its mineralogy, density, and origin. The estimated core composition has a S abundance near the threshold value for inner core solidification. The former presence of a core dynamo may be suggested by remanent magnetization in Shergottite-Nakhlite-Chassignite (SNC) meteorites, although these rocks may have been magnetized during shock. The mineralogy of martian surface units, inferred from reflectance spectra, matches that of basaltic shergottites, but SNC lithologies thought to have crystallized in the subsurface are not presently recognized. The rheological properties of martian magmas are more accurately derived form these metorites than from observations of martian flow morphology, although the sampled range of magma compositions islimited. Estimates of planetary water abundance and the amount of outgassed water based on these meteorites are contridictory but overlap estimates based on geological observations and atmospheric measurements. Stable isotope measurements indicate that the martian hydrosphere experienced only limited exchange with the lithosphere, but it is in isotopic equilibrium with the atmosphere and has been since 1.3 Ga. The isotopically heavy atmosphere/hydrosphere composition deduced from these rocks reflects a loss process more severe than current atmospheric evolution models, and the occurence of carbonates in SNC meteorites suggest that they, rather than scapolite or hydrous carbonates, are the major crustal sink for CO2. Weathering products in SNC meteorites support the idea of limited alteration of the lithosphere by small volumes of saline, CO2-bearing water. Atmospheric composition and evolution are further constrained by noble gases in these meteorites, although Xe and Kr isotopes suggest different origins for the atmosphere. Planetary ejection of these rocks has promoted an advance in the understanding of impact physics, which has been accomplished by a model involving spallation during large cratering events. Ejection of all the SNC meteorites (except ALH84001) in one or two events may provide a plausible solution to most constraints imposed by chronology, geochemistry, and cosmic ray exposure, although problems remain with this scenario; ALH84001 may represent older martian crust sampled during a separate impact.

Mcsween, Harry Y., Jr.↗

Properties of a B0 I stellar wind and interstellar grains derived from Ginga observations of the binary X-ray pulsar 4U 1538-52

From measurements of the X-ray eclipse phenomena of the binary X-ray pulsar 4U 1538-52, we derive properties of the stellar wind of its B0 I companion, QV Nor, and a constraint on models of interstellar grains. Estimates of the wind density as a function of the distance from the center of QV Nor are obtained from an analysis of the variation of X-ray attenuation during an eclipse egress. A Monte Carlo computation of the absorption and scattering of X-rays in the X-ray-ionized wind accounts for approximately two-thirds of the spectrum of X-rays with energies above 4.5 keV observed during the eclipse. Our upper limit on R(sub XV) is 0.06/mag, which implies that the X-ray scattering efficiency of interstellar dust is less than expected for solid grains with a size distribution of the form n(sub g)(a) approximately a(exp -3.5) in the range from 0.005 to 0.25 microns and composed of silicate (R(sub XV) = 0.22/mag) or a silicate-graphite mixture (R(sub XV) = 0.11/mag) as derived from the calculations of Martin & Rouleau (1991). This lends support to the idea (Mathis & Whiffen 1989) that interstellar grains are 'fluffy' aggregates with an average bulk density less than that of their constitutent particles. Such aggregates would have a smaller ratio of X-ray scattering efficiency to optical extinction efficiency compared with solid grains of the same material.

Clark, George W.↗

IR excitation of contaminant water by oxygen for the space shuttle at low Earth orbit altitude

As the water outgas of a space shuttle passes through the rarefied atmosphere at orbital altitude, collisions occur between the gases with sufficient energy to excite infrared-active water molecules to various vibrational and rotational states. An infrared contaminant model (IR model) has been developed to study the shuttle-induced excitation and emission of water molecules outgassed from the space shuttle. The focus of the first application of the model is translation-to-vibration (T-V) energy transfer since estimates suggest that this process should dominate the production of vibrationally excited H2O under typical low Earth orbit conditions. Using the velocity and position distribution functions of interacting neutral gases obtained from a neutral gases interaction model, the spatial distributions of excitation and IR radiation from contaminant water are computed, and typical results are presented. Infrared spectral data (450 - 2500/cm), measured by the Cryogenic Infrared Radiance Instrumentation for Shuttle (CIRRIS-1A) sensor on STS-39 (April 28 to May 6, 1991) at an altitude near 265 km, are used to test model predictions. The dependence of the radiant emission structure and brightness on outgassing rates and altitudes is discussed. The time history of the contaminant water outgassing rate is inferred for STS-39, and it is compared with the mass-spectrometer-based results for STS-4 (June 26 to July 4, 1982). Also, estimates of H2O column density at mission elapsed time (MET) 50 hours are compared for missions STS-2, STS-3, STS-4, and STS-39.

Zhou, D. K.↗

A high spectral resolution VLA search for H I absorption towards A496, A1795, and A2584

In this paper, we present the results of a Very Large Array (VLA) search for H I absorption with high spectral resolution (1.6 km/s) towards A496, A1795, A2584, and A2597. These observations are well matched to the properties of cold, optically thick H I clouds, where the line width is given by the width of an individual cloud rather than the dispersion in an ensemble of clouds. We do not detect any H I absorption with narrow linewidths in these clusters. Our limits mainly apply to clouds which are larger than a few tenths parsec-i.e., if the clouds are much smaller than the background radio source and have a low covering factor in velocity space, they could still escape detection. The estimated limits on column density (for clouds in this regime of parameter space) are 2-3 orders of magnitude less than the 10(exp 21)/sq cm required to explain the x-ray absorption seen in some cooling flow clusters. The combination of our high spectral resolution H I absorption searches with the existing lower spectral resolution H I absorption searches and the searches for H I emission makes it unlikely that atomic hydrogen is the dominant component of the cold x-ray absorbing gas in the inter-cloud medium (ICM).

O'Dea, Christopher P.↗

Thermal instability in the inner coma of a comet

The spacecraft and ground based observations of comet Halley inner coma showed a localized ion density depletion region whose origin is not well understood. Although it has been linked to a thermal instability associated with negative ions, the photodetachment lifetime of negative ions (approximately 1 sec) is too short compared to the electron attachment time scale (approximately 100 sec) for this process to have a significant effect. A mechanism for the ion density depletion based on the thermal instability of the cometary plasma due to the excitation of rotational and vibrational levels of water molecules is proposed. The electron energy losses due to these processes peak near 4000 K (0.36 eV) and at temperatures higher than this value a localized cooling leads to further cooling (thermal instability) due to the increased radiation loss. The resulting increase in recombination leads to an ion density depletion and the estimates for this depletion at comet Halley agree with the observations.

Milikh, G. M.↗

Combustion Zone Characterization of GO2/GH2 Rocket Using Laser-Induced Fluorescence of OH

A series of studies were undertaken to measure various parameters such as velocity, species concentration, and temperature downstream of a shear coaxial injector in an optically accessible uni-element rocket chamber. Techniques applied include: laser Doppler velocimetry (LDV) for velocity; laser light scattering (LLS) for flow visualization and estimating mixture fraction and density; laser induced fluorescence (LIF) of hydroxyl radials (OH) to determine the characteristics and extent of the reaction zone; and Raman spectroscopy to measure major species concentration and temperature.

Moser, M. D.↗

Implications of Depleted flux Tubes in the Jovian Magnetosphere

A rare but persistent phenomenon in the jovian magnetosphere is the occurrence of apparently depleted flux tubes, whose magnetic pressures are significantly above ambient levels. These flux tubes occur about 0.25% of the observing time in the region of the Io torus in the Galileo high resolution data. The importance of these tubes is that they can return to the inner magnetosphere the magnetic flux that has been convected radially outward with the iogenic plasma to the tail. The paucity of these tubes is consistent with the expected flux return rates if the tubes are moving inward at an average rate of about 5-10 km/s in the torus. Depleted flux tubes have yet to be observed inside of the lo orbit where the plasma beta is lower than in the hot torus. Estimates of the plasma density outside the tube from plasma wave measurements enable the average perpendicular temperature to be obtained from the magnetic field change. Extrapolating this temperature back to lo, we obtain an average ion temperature of approximately 60 eV. These values are generally consistent with earlier Voyager observations but on the low side of their range of uncertainty, and agree quite well with contemporaneous Galileo measurements where these are available.

Russell, C. T.↗

The Intercluster Medium in Clusters

With its large collecting area and adequate angular resolution XMM will enable large samples of well determined masses and baryon fractions for groups and clusters of galaxies. This will allow an estimate of the mass density of collapsed objects at the present epoch. Rich clusters can be detected and spatially resolved at r greater than 2 if they exist. There are numerous observing proposals to determine the cluster evolution with redshift which will have strong cosmological implications. The Fe abundance of rich clusters out to z-1 can be determined and thus strongly constrain the star formation history of the universe. Astro-E high spectral resolution will allow detailed determination of the elemental composition of galaxies groups and clusters and determine the origin of these elements. Astro-E data will also allow measurement of the velocity field in clusters for comparison with detailed numerical simulations.

Mushotzky, Richard↗

A Multi-Wavelength Study of the Hot Component of the Interstellar Medium

This research focuses on the kinematics and evolution of the hot phase of the interstellar medium in the Galaxy. The plan was to measure the UV spectra of all hot stars observed with IUE, in order to identify and measure the main component and any high velocity components to the interstellar lines. Collection of data from higher resolution instruments on HST has been proposed for some of the interesting lines of sight. IUE spectra of 240 stars up to 8 kpc in 2 quadrants of the galactic plane have been examined to (1) estimate the total column density per kpc as a function of direction and distance, and (2) to obtain a lower limit to the number of high velocity components to the interstellar lines, thus giving an approximation of the number of conductive interfaces encountered per line of sight. By determining an approximation to the number of components per unit distance we aim to derive statistics on interfaces between hot and cold gas in the Galaxy. We find that 20% of the stars in this sample show at least one high velocity component in the C IV interstellar line. Two successful FUSE programs address this research and collected data for several of the lines of sight identified as locations of hot, expanding gas with the IUE data. One FUSE program is complete for the Vela SNR region. Data from another FUSE program to investigate the Cygnus superbubble region are being analyzed.

Oliversen, Ronald J.↗

Anisotropic Piezocomposite Actuator Incorporating Machined PMN-PT Single Crystal Fibers

The design, fabrication, and testing of a flexible, planar, anisotropic piezoelectric composite actuator utilizing machined PMN-32%PT single crystal fibers is presented. The device consists of a layer of rectangular single crystal piezoelectric fibers in an epoxy matrix, packaged between interdigitated electrode polyimide films. Quasistatic free-strain measurements of the single crystal device are compared with measurements from geometrically identical specimens incorporating polycrystalline PZT-5A and PZT-5H piezoceramic fibers. Free-strain actuation of the single crystal actuator at low bipolar electric fields (+/- 250 V/mm) is approximately 400% greater than that of the baseline PZT-5A piezoceramic device, and 200% greater than that of the PZT-5H device. Free-strain actuation under high unipolar electric fields (0-4kV/mm) is approximately 200% of the PZT-5A baseline device, and 150% of the PZT-5H alternate piezoceramic device. Performance increases at low field are qualitatively consistent with predicted increases based on scaling the low-field d(sub 33) piezoelectric constants of the respective piezoelectric materials. High-field increases are much less than scaled d(sub 33) estimates, but appear consistent with high-field freestrain measurements reported for similar bulk single-crystal and piezoceramic compositions. Measurements of single crystal actuator capacitance and coupling coefficient are also provided. These properties were poorly predicted using scaled bulk material dielectric and coupling coefficient data. Rules-of-mixtures calculations of the effective elastic properties of the single crystal device and estimated actuation work energy densities are also presented. Results indicate longitudinal stiffnesses significantly lower (50% less) than either piezoceramic device. This suggests that single-crystal piezocomposite actuators will be best suited to low induced-stress, high strain and deflection applications.

Wilkie, W. Keats↗

Gravitational evidence for an undifferentiated Callisto

Before the arrival of the Galileo spacecraft at Jupiter, models for the interior structure of the four galilean satellites--Io, Europa, Ganymede and Callisto-ranged from uniform mixtures of rock and ice (that is, undifferentiated objects) or rocky cores surrounded by a mantle of water ice. Now it appears that Io has a large metallic core and that Ganymede is strongly differentiated, most probably into a three-layer structure consisting of a metallic core, a silicate mantle and a deep outer layer of ice. Direct information on the interior structure of Callisto determined from previous spacecraft fly-bys was essentially limited to an estimate of the mean density being intermediate between pure ice and pure rock. Here we report measurements of Callisto's gravitational field which reveal that, in contrast to Io and Ganymede, this galilean satellite is most probably a homogeneous object consisting of a solar mixture of 40% compressed ice and 60% rock (including iron and iron sulphide). Callisto's undifferentiated state is consistent with the apparent lack of an intrinsic magnetic field, and indicates that the outermost galilean satellite has not experienced a heating phase sufficiently high to separate its rock and metal components from the lighter ices.

unmanned↗

Thermospheric Studies with Mars Global Surveyor

The Mars Global Surveyor spacecraft has been located in a near-circular, polar, and low-altitude mapping orbit about Mars for six years, since February 1999. The spacecraft is tracked routinely by the antennae of the Deep Space Network (DSN), using the X Band radio system of the spacecraft. These tracking data have been used for routine spacecraft navigation, and for radio science studies, such as the estimation of the static and time-varying gravity field of Mars. In this paper we describe the methodology for reduction of these data in order to estimate the Mars atmospheric density (normalized to an altitude 380 km) over half a solar cycle, where we discern the correlation of the density with the incident solar flux, and the 27-day solar rotation. The results show that the density at the MGS altitude varies from a mean of 0.7 x 10(exp -17) grams/cu cm near aphelion to a mean of 3.0 x 10(exp -17)grams/cu cm near perihelion.

Lemoine, F. G.↗

Coronas-F Orbit Monitoring and Re-Entry Prediction

Russian scientific satellite CORONAS-F was launched on July, 31, 2001. The object was inserted in near-circular orbit with the inclination 82.5deg and a mean altitude approx. 520 km. Due to the upper atmosphere drag CORONAS-F was permanently descended and as a result on December, 6, 2005 it has finished the earth-orbital flight, having lifetime in space approx. 4.5 years. The satellite structural features and its flight attitude control led to the significant variations of its ballistic coefficient during the flight. It was a cause of some specific difficulties in the fulfillment of the ballistic and navigation support of this space vehicle flight. Besides the main mission objective CORONAS-F also has been selected by the Inter-Agency Space Debris Coordination Committee (IADC) as a target object for the next regular international re-entry test campaign on a program of surveillance and re-entry prediction for the hazard space objects within their de-orbiting phases. Spacecraft (S/C) CORONAS-F kept its working state right up to the end of the flight - down to the atmosphere entry. This fact enabled to realization of the additional research experiments, concerning with an estimation of the atmospheric density within the low earth orbits (LEO) of the artificial satellites, and made possible to continue track the S/C during final phase of its flight by means of Russian regular command & tracking system, used for it control. Thus there appeared a unique possibility of using for tracking S/C at its de-orbiting phase not only passive radar facilities, belonging to the space surveillance systems and traditionally used for support of the IADC re-entry test campaigns, but also more precise active trajectory radio-tracking facilities from the ground control complex (GCC) applied for this object. Under the corresponding decision of the Russian side such capability of additional high-precise tracking control of the CORONAS-F flight in this period of time has been implemented. The organizing of the CORONAS-F ballistic and navigational support (BNS) and solving its main tasks (such as S/C orbit determination (OD) and its motion prediction and connected with them) both for regular mission stage and for additional flight program were realized by the group of specialists from the Mission Control Center (MCC). MCC was also assigned as a principal organization from the Russian side for participation in the 7th IADC re-entry test campaign on CORONAS-F. The CORONAS-F flight features and space environments circumstances during its flight as well as a methodology and technology of spacecraft ballistic and navigational support are given below. The BNS results for different phases of S/C flight, including the results of its re-entry predictions, obtained during the realization of the 7th IADC test campaign are submitted. The accuracy of space vehicle re-entry prediction and its dependence on various factors are analyzed in more details.

Ivanov, N. M.↗

Anisotropic Laminar Piezocomposite Actuator Incorporating Machined PMN-PT Single Crystal Fibers

The design, fabrication, and testing of a flexible, laminar, anisotropic piezoelectric composite actuator utilizing machined PMN-32%PT single crystal fibers is presented. The device consists of a layer of rectangular single crystal piezoelectric fibers in an epoxy matrix, packaged between interdigitated electrode polyimide films. Quasistatic free-strain measurements of the single crystal device are compared with measurements from geometrically identical specimens incorporating polycrystalline PZT-5A and PZT-5H piezoceramic fibers. Free-strain actuation of the single crystal actuator at low bipolar electric fields (+/- 250 V/mm) is approximately 400% greater than that of the baseline PZT-5A piezoceramic device, and 200% greater than that of the PZT-5H device. Free-strain actuation under high unipolar electric fields (0-4kV/mm) is approximately 200% of the PZT-5A baseline device, and 150% of the PZT-5H alternate piezoceramic device. Performance increases at low field are qualitatively consistent with predicted increases based on scaling the low-field d33 piezoelectric constants of the respective piezoelectric materials. High-field increases are much less than scaled d33 estimates, but appear consistent with high-field freestrain measurements reported for similar bulk single-crystal and piezoceramic compositions. Measurements of single crystal actuator capacitance and coupling coefficient are also provided. These properties were poorly predicted using scaled bulk material dielectric and coupling coefficient data. Rules-of-mixtures calculations of the effective elastic properties of the single crystal device and estimated actuation work energy densities are also presented. Results indicate longitudinal stiffnesses significantly lower (50% less) than either piezoceramic device. This suggests that single-crystal piezocomposite actuators will be best suited to low induced-stress, high strain and deflection applications.

Wilkie, W. Keats↗

Multiple-Beam Detection of Fast Transient Radio Sources

A method has been designed for using multiple independent stations to discriminate fast transient radio sources from local anomalies, such as antenna noise or radio frequency interference (RFI). This can improve the sensitivity of incoherent detection for geographically separated stations such as the very long baseline array (VLBA), the future square kilometer array (SKA), or any other coincident observations by multiple separated receivers. The transients are short, broadband pulses of radio energy, often just a few milliseconds long, emitted by a variety of exotic astronomical phenomena. They generally represent rare, high-energy events making them of great scientific value. For RFI-robust adaptive detection of transients, using multiple stations, a family of algorithms has been developed. The technique exploits the fact that the separated stations constitute statistically independent samples of the target. This can be used to adaptively ignore RFI events for superior sensitivity. If the antenna signals are independent and identically distributed (IID), then RFI events are simply outlier data points that can be removed through robust estimation such as a trimmed or Winsorized estimator. The alternative "trimmed" estimator is considered, which excises the strongest n signals from the list of short-beamed intensities. Because local RFI is independent at each antenna, this interference is unlikely to occur at many antennas on the same step. Trimming the strongest signals provides robustness to RFI that can theoretically outperform even the detection performance of the same number of antennas at a single site. This algorithm requires sorting the signals at each time step and dispersion measure, an operation that is computationally tractable for existing array sizes. An alternative uses the various stations to form an ensemble estimate of the conditional density function (CDF) evaluated at each time step. Both methods outperform standard detection strategies on a test sequence of VLBA data, and both are efficient enough for deployment in real-time, online transient detection applications.

Thompson, David R.↗

A Fast and Sensitive New Satellite SO2 Retrieval Algorithm based on Principal Component Analysis: Application to the Ozone Monitoring Instrument

We describe a new algorithm to retrieve SO2 from satellite-measured hyperspectral radiances. We employ the principal component analysis technique in regions with no significant SO2 to capture radiance variability caused by both physical processes (e.g., Rayleigh and Raman scattering and ozone absorption) and measurement artifacts. We use the resulting principal components and SO2 Jacobians calculated with a radiative transfer model to directly estimate SO2 vertical column density in one step. Application to the Ozone Monitoring Instrument (OMI) radiance spectra in 310.5-340 nm demonstrates that this approach can greatly reduce biases in the operational OMI product and decrease the noise by a factor of 2, providing greater sensitivity to anthropogenic emissions. The new algorithm is fast, eliminates the need for instrument-specific radiance correction schemes, and can be easily adapted to other sensors. These attributes make it a promising technique for producing longterm, consistent SO2 records for air quality and climate research.

Principal component analysis↗

Next Generation Aura-OMI SO2 Retrieval Algorithm: Introduction and Implementation Status

We introduce our next generation algorithm to retrieve SO2 using radiance measurements from the Aura Ozone Monitoring Instrument (OMI). We employ a principal component analysis technique to analyze OMI radiance spectral in 310.5-340 nm acquired over regions with no significant SO2. The resulting principal components (PCs) capture radiance variability caused by both physical processes (e.g., Rayleigh and Raman scattering, and ozone absorption) and measurement artifacts, enabling us to account for these various interferences in SO2 retrievals. By fitting these PCs along with SO2 Jacobians calculated with a radiative transfer model to OMI-measured radiance spectra, we directly estimate SO2 vertical column density in one step. As compared with the previous generation operational OMSO2 PBL (Planetary Boundary Layer) SO2 product, our new algorithm greatly reduces unphysical biases and decreases the noise by a factor of two, providing greater sensitivity to anthropogenic emissions. The new algorithm is fast, eliminates the need for instrument-specific radiance correction schemes, and can be easily adapted to other sensors. These attributes make it a promising technique for producing long-term, consistent SO2 records for air quality and climate research. We have operationally implemented this new algorithm on OMI SIPS for producing the new generation standard OMI SO2 products.

Sulfur dioxide↗