Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Raman excitation profiles”

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.

87 records · Page 4

Lunar Electrical Conductivity, Permeability and Temperature From Apollo Magnetometer Experiments

Magnetometers have been deployed at four Apollo sites on the Moon to measure remanent and induced lunar magnetic fields. Measurements from this network of instruments have been used to calculate the electrical conductivity, temperature, magnetic permeability, and iron abundance of the lunar interior. The measured lunar remanent fields range from 3 gammas (γ) minimum at the Apollo 15 site to 327 γ maximum at the Apollo 16 site. Simultaneous magnetic field and solar plasma pressure measurements show that the remanent fields at the Apollo 12 and 16 sites interact with, and are compressed by, the solar wind. Remanent fields at Apollo 12 and Apollo 16 are increased 16 γ and 32 γ, respectively, by a solar plasma bulk pressure increase of 1.5 x 10 -7 dynes/cm 2 . Global lunar fields due to eddy currents, induced in the lunar interior by magnetic transients, have been analyzed to calculate an electrical conductivity profile for the Moon. From nightside magnetometer data in the solar wind it has been found that deeper than 170 km into the Moon the conductivity rises from 3 x 10 -4 mhos/m to 10 -2 mhos/m at 1000-km depth. Recent analysis of data obtained in the geomagnetic tail, in regions free of complicating plasma effects, yields results consistent with nightside values. Conductivity profiles have been used to calculate the lunar temperature for an assumed lunar material of olivine. In the outer layer (~ 170 km thick) the temperature rises to 110° C, after which it gradually increases with depth to 1500° C at a depth of ~ 1000 km. Simultaneous measurements by magnetometers on the lunar surface and in orbit around the moon are used to construct a whole-moon hysteresis curve, from which the global lunar magnetic permeability is determined to be μ = 1.012 ± 0.006. The corresponding global induced dipole moment is 2 x 10 18 gauss-cm 3 for typical inducing fields of 10 -4 gauss in the lunar environment. Lunar free iron abundance corresponding to the global permeability is determined to be 2.5 ± 2.0 wt. %. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed compositional models of the lunar interior. For a free iron/orthopyroxene lunar composition the total iron content is calculated to be 12.8 ± 1.0 wt. %; for a free iron/olivine composition, total iron content is 5.5 ± 1.2 wt. %. Other lunar models with an iron core and with a shallow iron-rich layer are also discussed in light of the measured global lunar permeability. Velocities and thicknesses of the Earth's magnetopause and bow shock have been estimated from simultaneous magnetometer measurements. Average speeds are determined to be about 50 km/s for the magnetopause and 70 km/s for the bow shock, although there are large variations in the measurement for any particular boundary crossing. Corresponding measured boundary thicknesses average about 2300 km for the magnetopause and 1400 km for the bow shock.

P Dyal

Behavior of Ullage Bubbles During Blowdown in Low-g Experiment (BUBBLE): Overview of a Cryogenic Tank Depressurization Test

Behavior of Ullage Bubbles during Blowdown in Low-g Experiment (BUBBLE) is a cryogenic tank depressurization drop tower experiment currently being designed at NASA Glenn Research Center. There is a need for liquid level rise data while venting a cryogenic propellant tank below the liquid saturation pressure in a reduced gravity environment. When a cryogenic tank is vented to a vapor pressure below the liquid saturation pressure, bubbles become entrained in the liquid causing the bulk liquid-vapor interface to rise. In reduced gravity, bubble rise velocity is diminished leading to a larger liquid level rise compared to a 1g environment. The purpose of the experiment is to gain further understanding of tank depressurization fluid physics and obtain high-fidelity data for model validation in 1g and reduced gravity environments. On-orbit cryogenic propellant tanks operating at high fill levels must efficiently manage venting operations to avoid the risk of liquid entrainment in the vent line, which could lead to asymmetric control thruster loads, freezing and clogging, and loss of liquid propellant. Validated models could be used to design settling and venting profiles to reduce risk and increase efficiency for cryogenic storage and transfer operations.

tank venting

Properties of a Ni19.5Pd30Ti50.5 High-Temperature Shape Memory Alloy in Tension and Compression

Potential applications involving high-temperature shape memory alloys have been growing in recent years. Even in those cases where promising new alloys have been identified, the knowledge base for such materials contains gaps crucial to their maturation and implementation in actuator and other applications. We begin to address this issue by characterizing the mechanical behavior of a Ni 19.5 Pd 30 Ti 50.5 high-temperature shape memory alloy in both uniaxial tension and compression at various temperatures. Differences in the isothermal uniaxial deformation behavior were most notable at test temperatures below the martensite finish temperature. The elastic modulus of the material was very dependent on strain level; therefore, dynamic Young’s Modulus was determined as a function of temperature by an impulse excitation technique. More importantly, the performance of a thermally activated actuator material is dependent on the work output of the alloy. Consequently, the strain-temperature response of the Ni 19.5 Pd 30 Ti 50.5 alloy under various loads was determined in both tension and compression and the specific work output calculated and compared in both loading conditions. It was found that the transformation strain and thus, the specific work output were similar regardless of the loading condition. Also, in both tension and compression, the strain-temperature loops determined under constant load conditions did not close due to the fact that the transformation strain during cooling was always larger than the transformation strain during heating. This was apparently the result of permanent plastic deformation of the martensite phase with each cycle. Consequently, before this alloy can be used under cyclic actuation conditions, modification of the microstructure or composition would be required to increase the resistance of the alloy to plastic deformation by slip.

High-Temperature Shape Memory Alloy

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle

Advanced nickel-hydrogen spacecraft battery development

Eagle-Picher currently has several advanced nickel-hydrogen (NiH2) cell component and battery designs under development including common pressure vessel (CPV), single pressure vessel (SPV), and dependent pressure vessel (DPV) designs. A CPV NiH2 battery, utilizing low-cost 64 mm (2.5 in.) cell diameter technology, has been designed and built for multiple smallsat programs, including the TUBSAT B spacecraft which is currently scheduled (24 Nov. 93) for launch aboard a Russian Proton rocket. An advanced 90 mm (3.5 in.) NiH2 cell design is currently being manufactured for the Space Station Freedom program. Prototype 254 mm (10 in.) diameter SPV batteries are currently under construction and initial boilerplate testing has shown excellent results. NiH2 cycle life testing is being continued at Eagle-Picher and IPV cells have currently completed more than 89,000 accelerated LEO cycles at 15% DOD, 49,000 real-time LEO cycles at 30 percent DOD, 37,800 cycles under a real-time LEO profile, 30 eclipse seasons in accelerated GEO, and 6 eclipse seasons in real-time GEO testing at 75 percent DOD maximum. Nickel-metal hydride battery development is continuing for both aerospace and electric vehicle applications. Eagle-Picher has also developed an extensive range of battery evaluation, test, and analysis (BETA) measurement and control equipment and software, based on Hewlett-Packard computerized data acquisition/control hardware.

Coates, Dwaine K.

Li-ion Battery Aging with Hybrid Physics-Informed Neural Networks and Fleet-wide Data

In this work, we propose a hybrid model for Li-ion battery discharge and aging prediction that leverages fleet-wide data to predict future capacity drops.The model is built upon an hybrid approach merging physics-based and empirical equations, as well as neural network models in a recurrent neural network cell. The hybrid physics-informed neural network can predict voltage discharge cycles given the loading profile, and estimate the used capacity of the battery under random-loading conditions by tracking aging parameters connected to the residual capacity of the battery. By merging information on the battery aging parameters with existing fleet-wide aging data, the model can predict the future residual capacity of the battery that is being monitored, and therefore enable predictions of voltage discharge curves far ahead in the battery life cycle. We validated the approach using the NASA Prognostics Data Repository Battery data-set, which contains experimental data on Li-ion batteries discharged at random loading conditions in a controlled environment. The approach also allows the identification of discrepancies between the battery aging trend and the trend observed at the fleet level, so that batteries behaving differently from the rest of the fleet can be subject to closer monitoring and further testing to refine predictions.

PINN

Radiation-Conditioned Ti–hBN Coatings for Space Mechanisms: Bridging Tribology and Irradiation Effects from Laboratory to Low-Earth Orbit

The longevity of tribological components in extraterrestrial environments is challenged by abrasive regolith, extreme temperatures, and ionizing radiation. This study evaluates the performance of vacuum plasma sprayed (VPS) Ti-2vol.% hexagonal boron nitride (hBN) coatings applied to joint mechanisms such as rod and slot, ball and socket, and hinge assemblies, fabricated from Al6061 and Ti6Al4V for planetary structures. Coated and uncoated samples were tested under vacuum with lunar regolith simulant (JSC-1A), followed by environmental exposures including electron radiation and thermal cycling. Results show that coated samples consistently outperformed uncoated counterparts across all configurations. Wear depth in coated components was reduced by over 80%, as confirmed by 3D profilometry and surface imaging. SEM analysis revealed severe abrasion and localized plastic deformation in uncoated surfaces, while coated samples retained structural integrity due to synergy of hard protective phases formed during deposition and retained hBN providing solid lubrication in the sliding surface. Coated and exposed samples exhibited smoother actuation profiles and reduced friction, indicating a beneficial "conditioning" effect from environmental exposure. The coatings-maintained performance even on complex geometries with minimal delamination. These findings demonstrate the potential of VPS Ti-hBN coatings to extend the operational life of lunar mechanisms, bridging laboratory-scale tribology with flight-relevant mechanical applications.

tribology

Development of global/chemistry model for jet-fuel thermal stability based on observations from static and flowing experiments

Two global-chemistry models for oxidative deposition of jet fuels are evaluated by integrating them into a Computational Fluid Dynamics with Chemistry (CFDC) code. A previously developed two-step global-chemistry model was found to be insufficient to describe the thermal-oxidation and -deposition rates associated with a Jet-A fuel. A new global-chemistry model has been developed systematically based on observations from flowing and static experiments. The global-autoxidation reaction is modified such that the reaction rate becomes zeroth-order with respect to the dissolved oxygen concentration. The generation of deposit-forming precursor is coupled with the autoxidation reaction by introducing a radical species ROO. A formulation for the sticking probability has also been developed. Deposition profiles are well represented by this new model under a variety of temperature and flow conditions. The model correctly predicts the changes in magnitude and spatial location of the deposition peak due to changes in flow. The CFDC model, which is designed for flowing systems, has been extended to static experiments. The model incorporates a non-depleting species F(sub s) representing all non-oxygen compounds responsible for deposition. Static experiments were found to provide a useful and inexpensive method for estimating the concentration of F(sub s) in the fuel.

V. R. Katta

Risk Reduction Project in Pad Abort-1 Launch Vehicle Loads & Dynamics

The Pad Abort (PA-1) system includes the Launch Abort System (LAS) and the Command Module (CM). The PA-1 abort flight test will launch from the White Sands Missile Range. Prior to the ignition, the vehicle will be resting, without restraints, on a short support structure. The Launch Abort Motor will ignite and burn for less than 5 seconds before a 20 second coast. Then the LAS will separate from the CM and the CM will descend via parachute to a soil landing. The static firing of the Abort Motor (ST-1) for PA-1 Launch Vehicle resulted in unexpected higher levels that superseded the environment predictions and all the design and test loads of all subassemblies and components. A rapid project was put in place to develop the flight environments, loads and associated uncertainties for the Verification Load Cycle (VLC) of PA-1. An acoustic test and shaker test of the CM were planned and executed in support of the developments of damping and transfer functions. Instead of using the traditional structure-borne envelops, the actual ST-1 pressures and forces were used to develop the internal accelerations. The PA-1 motor was cooled to 60 Deg F using thermal conditioning to reduce the thrust profile. The plumes of the abort motor were studied to determine if the airborne environments would be reduced. The combination of all items, noted above, allowed a reduction in loads for VLC.

Sasan C Armand

Design and Optimization of an Open Rotor Stage

Open rotor propulsors enable greatly increased effective bypass ratios over current aircraft engines, resulting in significantly improved efficiencies. One of the most promising of these geometries is the “open fan”, which consists of a single rotor followed by a single stator. This configuration will be referred to as an “open rotor stage” for design work at NASA. In this work we consider the multidisciplinary design and optimization of an open rotor stage using a variety of 2D and 3D techniques. The 2D aerodynamic techniques used for design include a 2D Trefftz plane method, axisymmetric throughflow, and quasi-3D blade-to-blade CFD. Aerodynamic performance is prioritized, while structural requirements are met and some acoustic considerations are included. Between the rotor, stator, and hub profiles, a total of 101 design variables are used. Differences in the process of designing an open rotor stage vs. designing other high-speed propellers are also discussed. The current design, the NASA Optimized Propulsor for Efficiency and Noise (N-OPEN-1), achieves 84.7% net efficiency as predicted by 3D CFD, with further improvements expected to be possible. This design makes effective use of a slower rotor and a less tip-loaded work distribution than typical single-rotation propellers, allowing the stator to produce 14.8% of the thrust. Laminar flow on the blades was responsible for a 1.7% efficiency benefit and warrants further study.

Open rotors

Design and Optimization of an Open Rotor Stage

Open rotor propulsors enable greatly increased effective bypass ratios over current aircraft engines, resulting in significantly improved efficiencies. One of the most promising of these geometries is the “open fan”, which consists of a single rotor followed by a single stator. This configuration will be referred to as an “open rotor stage” for design work at NASA. In this work we consider the multidisciplinary design and optimization of an open rotor stage using a variety of 2D and 3D techniques. The 2D aerodynamic techniques used for design include a 2D Trefftz plane method, axisymmetric throughflow, and quasi-3D blade-to-blade CFD. Aerodynamic performance is prioritized, while structural requirements are met and some acoustic considerations are included. Between the rotor, stator, and hub profiles, a total of 101 design variables are used. Differences in the process of designing an open rotor stage vs. designing other high-speed propellers are also discussed. The current design, the NASA Optimized Propulsor for Efficiency and Noise (N-OPEN-1), achieves 84.7% net efficiency as predicted by 3D CFD, with further improvements expected to be possible. This design makes effective use of a slower rotor and a less tip-loaded work distribution than typical single-rotation propellers, allowing the stator to produce 14.8% of the thrust. Laminar flow on the blades was responsible for a 1.7% efficiency benefit and warrants further study.

turbomachinery

Comparison Testing and Analysis of Braided Cords Constructed from Varying Kevlar Fiber Types

Material advancements move quickly, resulting in material obsolescence issues that require parachute engineers to understand how the material properties change with a material upgrade and if there are expected impacts to overall performance of the parachute. The objective of this test was to measure the elongation of different Kevlar types at both the cord-level and at the yarn-level to determine if a change in Kevlar type impacted parachute performance. An Instron machine with a pin-to-pin connection was used to measure load and displacement time histories for various Kevlar cord and yarn samples in a cyclic load profile. Twenty-four cord samples and twenty-four yarn samples were tested; each sample was cycled five times to the relevant maximum load (either to breaking strength or to the 30% peak load expected in flight). In addition to the Instron measurements, a photogrammetric technique was also implemented to measure elongation of the samples using button targets. Modulus test data shows approximately a 10% difference in elongation between Kevlar 29AP and Kevlar 129, and another 10% difference in elongation between Kevlar 29 and Kevlar 29AP. Although this is a measurable change in modulus, the effect on parachute performance is likely small, particularly for non-primary structure elements of relatively small length. This test campaign provides an example of the type of testing and analysis that parachute engineers could do to understand how material upgrades affect performance of their systems.

Modulus Testing

Comparison Testing and Analysis of Braided Cords Constructed from Varying Kevlar Fiber Types

Material advancements move quickly, resulting in material obsolescence issues that require parachute engineers to understand how the material properties change with a material upgrade and if there are expected impacts to overall performance of the parachute. The objective of this test was to measure the elongation of different Kevlar types at both the cord-level and at the yarn-level to determine if a change in Kevlar type impacted parachute performance. An Instron machine with a pin-to-pin connection was used to measure load and displacement time histories for various Kevlar cord and yarn samples in a cyclic load profile. Twenty-four cord samples and twenty-four yarn samples were tested; each sample was cycled five times to the relevant maximum load (either to breaking strength or to the 30% peak load expected in flight). In addition to the Instron measurements, a photogrammetric technique was also implemented to measure elongation of the samples using button targets. Modulus test data shows approximately a 10% difference in elongation between Kevlar 29AP and Kevlar 129, and another 10% difference in elongation between Kevlar 29 and Kevlar 29AP. Although this is a measurable change in modulus, the effect on parachute performance is likely small, particularly for non-primary structure elements of relatively small length. This test campaign provides an example of the type of testing and analysis that parachute engineers could do to understand how material upgrades affect performance of their systems.

Parachutes

Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities

Satellite radar retrievals of clouds and precipitation rely on assumptions about the microphysical properties of hydrometeors such as particle size distributions (PSDs) and massdimensional relationships, which have been shown to vary in distinct environments. A deeper understanding of linkages between mass-dimensional relationships, ice crystal shape, and radar reflectivities across cloud and precipitation regimes is necessary for developing new retrievals for upcoming satellite radar missions, including the INvestigation of Convective UpdraftS (INCUS), as well as future missions that will feature space-borne radars observing clouds, convection, and precipitation. This study investigates ice microphysical characteristics through aircraft in-situ and remote sensing observations from two convective events during the Mid-latitude Continental Convective Clouds Experiment (MC3E) in the Southern Great Plains, and tests implications for ice particle habit in forward modeled radar observations. During MC3E, the Ka/Ku- band High Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) was onboard the NASA ER-2 aircraft, while in-situ measurements were taken onboard the University of North Dakota Citation. In-situ PSDs and particle shape information from optical array probe data are used to forward model radar reflectivities from particle scattering databases and are compared against HIWRAP reflectivities. Various particle combinations can be used to find agreement with HIWRAP observations, and the forward modeled reflectivities were sensitive to particle type. Varying the riming on the aggregate and graupel particles resulted in larger than 10 dB differences in reflectivities. The in-situ and radar observations imply the occurrence of aggregation, size sorting and lofting particles.

Julia A. Shates

Physical Parameters of Space Mission Asteroid Targets

Ground-based characterization of asteroids that are planned or potential targets of space missions provides important data on their physical parameters and properties. Knowledge of the properties of mission targets is important especially during mission preparation and planning, as it serves to select best or suitable targets for specific purpose of a given mission, to prepare mission plans, to constrain possible mission scenarios, and to design mission experiments. Such characterization efforts may be particularly critical for flyby missions that take only limited data during the high-speed flybys of their target asteroids, but they also provide very crucial data for targets of rendezvous missions. Moreover, long-term observations taken from Earth also allows the modelling of, or constraining, processes acting on the asteroids over extended time scales. Over the past years and decades we have obtained rich data on physical parameters of 82 asteroids that are planned or potential targets of space missions. Our primary observing technique is time-resolved (lightcurve) photometry, but we use also data obtained with other techniques, such as spectroscopy, thermal or radar observations. 15 of the 82 characterized asteroids are planned or possible targets of several space missions that are in flight or preparation, such as ESA’s Hera, RAMSES and PRIAMOS, NASA’s OSIRIS-APEX, JAXA’s Hayabusa2#, DESTINY+ and Next Generation Sample Return (NGSR), the Emirates Mission to Asteroids (EMA), and Karman+’s High Frontier, but we have also characterized 67 asteroids that are potential targets of space missions for their low delta-V’s and were announced as being “of interest to NASA” in the Small-Bodies-Observations- NASA mailing list or listed on the NHATS page of “Accessible NEAs”. The sample of asteroid targets span 3 orders of magnitude in size, with absolute magnitudes H from 12.57 to 26.8, corresponding to diameters from about 10 m to about 10 km. The sample contains a variety of taxonomic types and physically or dynamically interesting objects. Among them, we have identified 5 binary asteroids and 20 tumblers (i.e., asteroids in excited, non-principal axis rotation states). Rotation periods of the characterized asteroids range from 1.45 min to 280 h, reflecting diversity of their properties and formation/evolution paths. We will present an overview of the data set and highlight several representative cases.

Petr Pravec

Performance of Two Battery Prognostic Applications used by Two Octocopters for Safe Low Altitude Autonomous Flight Operations

This paper addresses the problem of building trust in online predictions of the remaining available flying time for two different electric Unmanned Aerial Vehicles (eUAVs) powered by lithium-ion-polymer batteries. Flight tests for various automation research missions for the two vehicles were monitored using two on-board battery health management applications to make predictions of the remaining flying time (RFT) for each eUAV and to predict the state of the battery. Playback of the voltage, current and temperature profiles of the battery discharge were used to assess the accuracy of the estimation of the voltage and the charge states of the models as well as the estimate of the RFT. The reference ground truth values were the observed landing time and the measured battery pack resting pack voltage 20 minutes after the flight. The predicted RFT, state of charge (SoC), and state of energy (SoE) were compared with the observed results. Noise values of one standard deviation from the mean values of the internal charge states of the battery model during a reference run were used to vary the states during simulation. One application used an equivalent circuit model of the electrical dynamics of the battery pack, and the other application used a reduced-order electrochemistry model. The variation of the model state components was compared to the variation in the estimate of the RFT and the variation in the SoE to estimate a confidence factor. Variation in the estimates caused by factors affecting the off-line laboratory parameter identification experiments is considered. Variation in the estimates due to environmental factors are discussed.

Assurance

The Effect of Gravity on the Combustion of Bulk Metals

In recent years, metal combustion studies at the University of Colorado have focused on the effects of gravity (g) on the ignition and burning behavior of bulk metals. The impetus behind this effort is the understanding of the ignition conditions and flammability properties of structural metals found in oxygen (O,sub>2 ) systems for space applications. Since spacecraft are subjected to higher-than-1g loads during launch and reentry and to a zero-gravity environment while in orbit, the study of ignition and combustion of bulk metals at different gravitational accelerations is of great practical concern. From the scientific standpoint, studies conducted under low gravity conditions provide simplified boundary conditions, since buoyancy is removed, and make possible the identification of fundamental ignition and combustion mechanisms. This investigation is intended to provide experimental verification of the influence of natural convection on the burning behavior of metals. In addition, the study offers the first findings of the influence of gravity on ignition of bulk metals and on the combustion mechanism and structure of metal-oxygen, vapor-phase diffusion flames in a buoyancy-free environment. Titanium (Ti) and magnesium (Mg) metals were chosen because of their importance as elements of structural materials and their simple chemical composition-pure metals instead of multicomponent alloys to simplify chemical and spectroscopic analyses. In addition, these elements present the two different combustion modes observed in metals: heterogeneous or surface burning (for Ti) and homogeneous or gas-phase reaction (for Mg). Finally, Mg, Ti, and their oxides exhibit a wide range of thermophysical and chemical properties. Metal surface temperature profiles, critical and ignition temperatures, propagation rates, burning times, and spectroscopic measurements are obtained under normal and reduced gravity. Visual evidence of all phenomena is provided by high-speed photography.

Melvyn C Branch

Sulfur outgassing and in-gassing in lunar orange glass beads and implications for 33S “Anomaly” in the Moon

Our recent investigations have discovered inward diffusion (in-gassing) of moderately volatile elements (MVEs; e.g., Na, K and Cu) from volcanic gas into volcanic beads/droplets. In this work, we examine the distribution of sulfur in lunar orange glass beads. Our analyses reveal that sulfur exhibits a non-uniform distribution across the beads, forming “U” or “W” shaped profiles typical of in-gassing. A model developed to assess sulfur contributions from different sources (original magmatic sulfur versus atmospheric in-gassed sulfur) in the orange beads indicates that atmospheric sulfur in-gassed during eruption contributes approximately 9–24 % to the total sulfur content of an orange bead, averaging around 16 %. This in-gassed sulfur is derived from the eruption plume, where atmospheric sulfur could undergo photochemical reactions induced by UV light, leading to mass independent fractionation and a distinct sulfur isotope signature. Interestingly, a recent study discovered a small mass independent isotope fractionation of sulfur in lunar orange glass beads in drive tube 74002/1 and a lack of such mass independent isotope fractionation in black glass beads in the same lunar sample. This finding contrasts with sulfur in lunar basalts, which typically exhibit mass dependent fractionation. With our work, the observed mass independent fractionation signal in sulfur isotopes of orange beads can be attributed to the in-gassing of photolytic sulfur in the optically thin part of the eruption plume where UV light can penetrate. Using the sulfur isotope data of lunar orange beads, we estimate that the Δ33S value of atmospheric sulfur is approximately −0.18 ‰. Our study provides new insights into the complex dynamics of volatile elements in lunar volcanic processes, highlighting the role of in-gassing in shaping sulfur isotope signatures in volcanic glass beads.

Xue Su