SEARCH · Engineering Papers

Results for “matrix multiplication”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

628 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Using Redundancy To Reduce Errors in Magnetometer Readings

A method of reducing errors in noisy magnetic-field measurements involves exploitation of redundancy in the readings of multiple magnetometers in a cluster. By "redundancy"is meant that the readings are not entirely independent of each other because the relationships among the magnetic-field components that one seeks to measure are governed by the fundamental laws of electromagnetism as expressed by Maxwell's equations. Assuming that the magnetometers are located outside a magnetic material, that the magnetic field is steady or quasi-steady, and that there are no electric currents flowing in or near the magnetometers, the applicable Maxwell 's equations are delta x B = 0 and delta(raised dot) B = 0, where B is the magnetic-flux-density vector. By suitable algebraic manipulation, these equations can be shown to impose three independent constraints on the values of the components of B at the various magnetometer positions. In general, the problem of reducing the errors in noisy measurements is one of finding a set of corrected values that minimize an error function. In the present method, the error function is formulated as (1) the sum of squares of the differences between the corrected and noisy measurement values plus (2) a sum of three terms, each comprising the product of a Lagrange multiplier and one of the three constraints. The partial derivatives of the error function with respect to the corrected magnetic-field component values and the Lagrange multipliers are set equal to zero, leading to a set of equations that can be put into matrix.vector form. The matrix can be inverted to solve for a vector that comprises the corrected magnetic-field component values and the Lagrange multipliers.

Kulikov, Igor↗

NASA NTRS2011

Post-Fragmentation Whole Genome Amplification-Based Method

This innovation is derived from a proprietary amplification scheme that is based upon random fragmentation of the genome into a series of short, overlapping templates. The resulting shorter DNA strands (<400 bp) constitute a library of DNA fragments with defined 3 and 5 termini. Specific primers to these termini are then used to isothermally amplify this library into potentially unlimited quantities that can be used immediately for multiple downstream applications including gel eletrophoresis, quantitative polymerase chain reaction (QPCR), comparative genomic hybridization microarray, SNP analysis, and sequencing. The standard reaction can be performed with minimal hands-on time, and can produce amplified DNA in as little as three hours. Post-fragmentation whole genome amplification-based technology provides a robust and accurate method of amplifying femtogram levels of starting material into microgram yields with no detectable allele bias. The amplified DNA also facilitates the preservation of samples (spacecraft samples) by amplifying scarce amounts of template DNA into microgram concentrations in just a few hours. Based on further optimization of this technology, this could be a feasible technology to use in sample preservation for potential future sample return missions. The research and technology development described here can be pivotal in dealing with backward/forward biological contamination from planetary missions. Such efforts rely heavily on an increasing understanding of the burden and diversity of microorganisms present on spacecraft surfaces throughout assembly and testing. The development and implementation of these technologies could significantly improve the comprehensiveness and resolving power of spacecraft-associated microbial population censuses, and are important to the continued evolution and advancement of planetary protection capabilities. Current molecular procedures for assaying spacecraft-associated microbial burden and diversity have inherent sample loss issues at practically every step, particularly nucleic acid extraction. In engineering a molecular means of amplifying nucleic acids directly from single cells in their native state within the sample matrix, this innovation has circumvented entirely the need for DNA extraction regimes in the sample processing scheme.

Benardini, James↗

NASA NTRS2008

Update on Development of SiC Multi-Chip Power Modules

Progress has been made in a continuing effort to develop multi-chip power modules (SiC MCPMs). This effort at an earlier stage was reported in 'SiC Multi-Chip Power Modules as Power-System Building Blocks' (LEW-18008-1), NASA Tech Briefs, Vol. 31, No. 2 (February 2007), page 28. The following recapitulation of information from the cited prior article is prerequisite to a meaningful summary of the progress made since then: 1) SiC MCPMs are, more specifically, electronic power-supply modules containing multiple silicon carbide power integrated-circuit chips and silicon-on-insulator (SOI) control integrated-circuit chips. SiC MCPMs are being developed as building blocks of advanced expandable, reconfigurable, fault-tolerant power-supply systems. Exploiting the ability of SiC semiconductor devices to operate at temperatures, breakdown voltages, and current densities significantly greater than those of conventional Si devices, the designs of SiC MCPMs and of systems comprising multiple SiC MCPMs are expected to afford a greater degree of miniaturization through stacking of modules with reduced requirements for heat sinking; 2) The stacked SiC MCPMs in a given system can be electrically connected in series, parallel, or a series/parallel combination to increase the overall power-handling capability of the system. In addition to power connections, the modules have communication connections. The SOI controllers in the modules communicate with each other as nodes of a decentralized control network, in which no single controller exerts overall command of the system. Control functions effected via the network include synchronization of switching of power devices and rapid reconfiguration of power connections to enable the power system to continue to supply power to a load in the event of failure of one of the modules; and, 3) In addition to serving as building blocks of reliable power-supply systems, SiC MCPMs could be augmented with external control circuitry to make them perform additional power-handling functions as needed for specific applications. Because identical SiC MCPM building blocks could be utilized in such a variety of ways, the cost and difficulty of designing new, highly reliable power systems would be reduced considerably. This concludes the information from the cited prior article. The main activity since the previously reported stage of development was the design, fabrication, and testing a 120- VDC-to-28-VDC modular power-converter system composed of eight SiC MCPMs in a 4 (parallel)-by-2 (series) matrix configuration, with normally-off controllable power switches. The SiC MCPM power modules include closed-loop control subsystems and are capable of operating at high power density or high temperature. The system was tested under various configurations, load conditions, load-transient conditions, and failure-recovery conditions. Planned future work includes refinement of the demonstrated modular system concept and development of a new converter hardware topology that would enable sharing of currents without the need for communication among modules. Toward these ends, it is also planned to develop a new converter control algorithm that would provide for improved sharing of current and power under all conditions, and to implement advanced packaging concepts that would enable operation at higher power density.

Lostetter, Alexander↗

NASA NTRS2012

Analysis of Direct Samples of Early Solar System Aqueous Fluids

Over the past three decades we have become increasingly aware of the fundamental importance of water, and aqueous alteration, on primitive solar-system bodies. Some carbonaceous and ordinary chondrites have been altered by interactions with liquid water within the first 10 million years after formation of their parent asteroids. Millimeter to centimeter-sized aggregates of purple halite containing aqueous fluid inclusions were found in the matrix of two freshly-fallen brecciated H chondrite falls, Monahans (1998, hereafter simply "Monahans") (H5) and Zag (H3-6) (Zolensky et al., 1999; Whitby et al., 2000; Bogard et al., 2001) In order to understand origin and evolution of the aqueous fluids inside these inclusions we much measure the actual fluid composition, and also learn the O and H isotopic composition of the water. It has taken a decade for laboratory analytical techniques to catch up to these particular nanomole-sized aqueous samples. We have recently been successful in (1) measuring the isotopic composition of H and O in the water in a few fluid inclusions from the Zag and Monahans halite, (2) mineralogical characterization of the solid mineral phases associated with the aqueous fluids within the halite, and (3) the first minor element analyses of the fluid itself. A Cameca ims-1270 equipped with a cryo-sample-stage of Hokkaido University was specially prepared for the O and H isotopic measurements. The cryo-sample-stage (Techno. I. S. Corp.) was cooled down to c.a. -190 C using liquid nitrogen at which the aqueous fluid in inclusions was frozen. We excavated the salt crystal surfaces to expose the frozen fluids using a 15 keV Cs+ beam and measured negative secondary ions. The secondary ions from deep craters of approximately 10 m in depth emitted stably but the intensities changed gradually during measurement cycles because of shifting states of charge compensation, resulting in rather poor reproducibility of multiple measurements of standard fluid inclusions of +/- 90 0/00(2 sigma) for delta D, and +/- 29 0/00 (2 sigma) for delta O-18. On the other hand, the reproducibility of Delta O-17 is plus or minus 8 /00 (2 sigma ) because the observed variations of isotope ratios follow a mass dependent fractionation law. Variations of delta D of the aqueous fluids range over sog,a 330(90; 2 sigma ) to +1200(90) 0/00 for Monahans and delta 300(96) 0/00 to +90(98)0/00 for Zag. Delta O-17 of aqueous fluids range over delta 16(22) 0/00 to +18(10) 0/00 for Monahans and +3(10) 0/00 to +27(11) 0/00 for Zag. These variations are larger than the reproducibility of standard analyses and suggest that isotope equilibria were under way in the fluids before trapping into halite. The mean values of delta D and Delta O-17 are +290 0/00 and +9 0/00, respectively. The mean values and the variations of these fluids are different from the representative values of ordinary chondrites, verifying our working hypothesis that the fluid inclusion-bearing halites were not indigenous to the H chondrite parent-asteroid but rather represent exogenous material delivered onto the asteroid from a separate cryovolcanically-active body. This initial isotopic work has demonstrated the feasibility of the measurements, but also revealed sample processing and analytical shortcomings that are now being addressed. Examination of solid mineral inclusions within Monahans and Zag halite grains by confocal Raman spectroscopy at the Carnegie Geophysical Laboratory has revealed them to be metal, magnetite, forsteritic olivine (Fo.98), macromolecular carbon (MMC), pyroxenes, feldspar with Raman spectral affinity to anorthoclase and, probably, fine-grained lepidocrocite (FeO(OH)). In addition, one inclusion features aliphatic material with Raman spectral features consistent with a mixture of short-chain aliphatic compounds. We have initiated analyses of the bulk composition of the fluids within the inclusions in Zag and Monahans halites at Virginia Tech by LA ICPMS using angilent 7500ce quadrupole ICPMS and a Lambda Physik GeoLas 193 nm excimer laser ablation system. Preliminary results reveal that the inclusion aqueous fluids contain highly charged cations of Ca, Mg and Fe. The minerals and compounds discovered thus far within Monahans/Zag halites are indicative of an originating body at least partly composed of unequilibrated anhydrous materials (high Fo olivine, pyroxenes, feldspars, possibly the metal) which were subjected to aqueous alteration (the halite parent brine) and containing a light organic component (the short-chain aliphatic compounds). This material was ejected from the originating body with little or no disruption, as evidenced with the presence of fluid inclusions. An actively geysering body similar to modern Enceladus (Postberg et al., 2011) may be a reasonable analogue in this respect. Also, the originating body should have been within close proximity to the H chondrite parent in order to generate the number of halite grains seen in Monahans and Zag. Other candidates for Monahans/Zag halite parent bodie(s) may include a young Ceres with its possible liquid ocean, or Main Belt comets.

Zolensky, Michael E.↗

NASA NTRS2019

High Velocity Impact Performance of a Dual Layer Thermal Protection System for the Mars Sample Return Earth Entry Vehicle

The Mars Sample Return (MSR) Earth Entry Vehicle (EEV) is currently planned on being released from its Micro-Meteorite/Orbital Debri (MM/OD) shielding housing about two days before the Earth entry phase. This leaves the EEV exposed to incoming MM/OD impacts, potentially damaging the heat shield and compromising its Entry, Decent, and Landing (EDL) integrity. Currently, two materials are proposed to comprise the MSR-EEV heat shield, a dual layer material Heat-shield for Extreme Entry Environment Technology (HEEET) and Phenolic Impregnated Carbon Ablator (PICA). PICA has been well characterized for OD class impacts, ~7km/s high mass impacts, from previous testing done in the Orion program, but hasn't under-gone extensive MM impact testing. HEEET is a rela-tively new material with minimal prior testing in re-gards to High Velocity Impacts (HVI). In order to in-form selectability of a material, it is crucial to under-stand the material performance when faced with an HVI, directly affecting mission success probability.The current measure of a Thermal Protection Sys-tem's (TPS) performance against an HVI is evaluating a thermally sized material against its derived Balistic Limit Equation (BLE). A BLE is generated empirically from multiple shots of HVI testing, and is used as a first order method in evaluating TPS's performance against the expected MM/OD environment. This method has proved useful for previous uniform densi-ty TPS materials, but has never been validated against a dual layer recession material such as HEEET. Testing in the MSR program for FY19 has a re-quirement to assess the effects of a dual layer TPS by testing various thicknesses of HEEET's Recession Layer, seen in Figure 1. From this data, BLE's will be derived for the individual thickness ratio samples, as well as the material as a whole to evaluate if a heritage form of the BLE can capture the complex physics associated with a dual layer system. Crater morpholo-gy will also be assessed with post-test Non-Destructive Evaluation (NDE) methods such as CT scanning to visualize if a BLE can well predict the associated pene-tration depths, since a BLE assumes full disinigration of the impacting particle and is generally only used to size a spherical crater ? disregarding any shrapenel effects from a high density impactor.To inform this analysis, expected MM environ-ments from the Meteoroid Engineering Model (MEM) and analytical equations for mass flux of incoming MM were evaluated against the notional MSR-EEV trajectory [1]. Using those dispersions, a monte-carlo was run to determine the most probable particle pa-rameters, as well as the riskiest in terms of full bondline penetration. From these probabilities, a test matrix was designed to test against bounding cases for the various parameters of the BLE: projectile density, projectile mass, projectile velocity, and the impact angle.This presentation will discuss the performance of the dual layer TPS material HEEET against a wide range of impact kinetic energies and densitites, as well as the comparison of HEEET to PICA in terms of MM/OD performance and selectability criteria.

Libben, B. J.↗

NASA NTRS1993

Relationships between pigment composition variation and reflectance for plant species from a coastal savannah in California

Advances in imaging spectroscopy have indicated that remotely sensed reflectance measurements of the plant canopy may be used to identify and qualify some classes of canopy biochemicals; however, the manner in which differences in biochemical compositions translate into differences is not well understood. Most frequently, multiple linear regression routines have been used to correlate narrow band reflectance values with measured biochemical concentrations. Although some success has been achieved with such methods for given data sets, the bands selected by multiple regression are not consistent between data sets, nor is it always clear what physical or biological basis underlies the correlation. To examine the relationship between biochemical concentration and leaf reflectance signal we chose to focus on the visible spectrum where the primary biochemical absorbances are due to photosynthetic pigments. Pigments provide a range of absorbance features, occur over a range of concentrations in natural samples, and are ecophysiologically important. Concentrations of chlorophyll, for example, have been strongly correlated to foliar nitrogen levels within a species and to photosynthetic capacity across many species. In addition pigments effectively absorb most of the photosynthetically active radiation between 400-700 nm, a spectral region for which silicon detectors have good signal/noise characteristics. Our strategy has been to sample a variety of naturally occurring species to measure leaf reflectance and pigment compositions. We hope to extend our understanding of pigment reflectance effects to interpret small overlapping absorbances of other biochemicals in the infrared region. For this reason, selected samples were also tested to determine total nitrogen, crude protein, cellulose, and lignin levels. Leaf reflectance spectra measured with AVIRIS bandwidths and wavelengths were compared between species and within species and for differences between seasons, for changes in the the shape of the spectra. We attempt to statistically correlate these shape changes with differences in pigment compositions. In parallel with our comparisons of pigment composition and leaf reflectance, we have modified the PROSPECT leaf reflectance model to test the contributions of pigments or pigment group concentrations. PROSPECT considers a leaf as a multi-layer dielectric plane with an uneven surface. Jacquemoud adapted the basic analysis of Allen for surface effects, a leaf thickness factor, and the absorption of water and chlorophyll (actually all pigments) and the plant matrix. Our modifications to PROSPECT in the forward direction include breaking out the pigment concentration parameter into separate components for chlorophyll a and b and a number of xanthophylls and carotenes, and introducing a shift and convolution function to model the spread and shift from their in vitro measurements to their in vivo state. Further, we have considered how the matrix elements (i.e., all biochemicals and structural effects not modeled explicity) vary with species.

Ustin, Susan L.↗

NASA NTRS2024

Qualitative Headspace GCMS Analysis of Lunar Regolith and Volatile Simulant Mixtures

Introduction: Future Artemis missions aim to return the volatile-bearing samples collected near lunar polar craters. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving the integrity of returned samples and the science value those samples contain. The extent to which that preservation is possible, and the trade-offs preservation requires (e.g. monetary costs, sample volume limitations) all must be considered. Even less-than pristine volatile-bearing samples will be of tremendous value to the scientific community seeking to unravel the history of lunar surface volatiles and, more broadly, volatiles in the solar system. A sample collected on the lunar surface will experience at least five distinctive periods during which any changes here referred to as “alteration” will certainly occur at some scale: collection on the lunar surface; transportation back to Earth; long-term storage; curatorial processing; and allocation/distribution. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on returned volatile samples by working with high-fidelity volatile-containing regolith simulants, setting the foundation for the future of cold curation. This abstract is focused on gas-surface interactions between LCROSS volatiles and readily available lunar regolith simulants. Experiments involved analyzing differences in headspace gas composition for various combinations of volatile and regolith simulants using gas chromatography/mass spectrometry (GC/MS). Experimental Procedure: The volatile simulants were chosen based on the molecules detected during the LCROSS mission.1 Stock solutions of condensed lunar volatile analytes were: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. The regolith simulants used were primarily JSC-1A and NU-LHT-4M. Additional regolith simulant control studies were conducted with <150 μm sieved sand and KBr. A nested vial sample preparation approach separated the liquid stock solutions from the regolith, eliminating potential matrix effects between liquid and solid phases. Fifteen microliter aliquots of liquid volatile simulants were added to a 2 mL liquid GC vial and capped in atmosphere. An 18G needle punctured the 2 mL GC vial immediately before being transferred and sealed in a 20 mL GC vial containing 0-0.3 g of regolith simulant, see Figure 1. Separating the analytes ensures any changes observed in the total headspace gases is a result of gas-surface and/or gas-gas interactions. Equipment and Method: Initial GC/MS method development for the separation and identification of relevant headspace gases can be found in Amick, et al. 2023.2 The only hardware change is a different column: a TG-1701MS 30 m × 0.25 mm × 1.00 μm column. The vials were sampled at 10°C, room temperature (~25°C) and 50°C. Low temperature samples were kept in a chilled autosampler stage for at least 1 hour prior to sampling. High temperature samples were agitated at 50°C for 5 minutes immediately prior to injection onto the column. Headspace chromatograms were collected for each combination of temperature, regolith, and volatile simulant, including controls without one or both types of analytes, in triplicate. The chromatogram elution window for each analyte or significant atmospheric gas was identified using the peak mass spectrum cross-referenced with Figure 1.Picture of the nested vial set-up. a NIST MS library search. Each analyte peak was integrated after filtering the mass spectrum trace for the parent or most unique mass-to-charge ratio. For example, the mass-to-charge ratio used to identify, filter for, and integrate the carbon dioxide peak was centered around 44 m/z. Results: Figure 2 shows the integrated peak area for hydrogen sulfide in all combinations of regolith simulants, temperatures, and the addition of ammonia in methanol solution. Multiple repeat experiments with H2S and regolith simulants have confirmed hydrogen sulfide is removed from the headspace within 1 hour when exposed to JSC-1A, NU-LHT-4M and sand (not pictured). The consumption of H2S by lunar regolith simulants at different temperatures indicates surface chemistry will be an integral component in sample integrity and preservation. It is important to note that while the effect of surface chemistry on gaseous hydrogen sulfide is intriguing, further investigation into more chemically accurate regolith simulants is necessary and ongoing. The sulfur in hydrogen sulfide is in its most reduced state while off-the-shelf, terrestrially sourced JSC-1A and NU-LHT-4M are more oxidized than most lunar materials,3which likely leads to different oxidation-reduction reactions than would be expected in lunar regolith. Figure 3 shows the integrated peak area of carbon dioxide for each sample combination that contained ammonia in methanol solution. The addition of ammonia to the GC vials results in a consistent and reproducible decrease in carbon dioxide gas, even at 50°C. The effect became more pronounced when JSC-1A or NU-LHT-4M were present. This set of experiments demonstrated that the sample composition will affect the chemical and physical state of each component present. Future spectroscopic and microscopy experiments will be geared towards identifying the cause for the CO2(g) concentration decrease and the consumption of hydrogen sulfide. As cold and volatile curation scientists, this information provides necessary insight on how to appropriately handle and analyze volatile bearing returned samples, as well as predict the effect chemical composition has on the various sample phases we will analyze upon return to Earth. Unlike traditional curation of geologic materials, the molecules in a sample cannot be identified or processed using the naked eye or even an optical microscope. Volatile curation will require a combination of analytical techniques, including but not limited to highly sensitive gas and solid/condensed phase spectroscopy. This set of experiments has demonstrated the need for more detailed studies of volatile mixtures with mineralogically and geochemically analogous lunar regolith simulants to prepare for the curation of volatile-rich lunar samples from the south polar region of the Moon. References: [1] Colaprete, A., et al. (2010) Science, 330, (463-468). [2] Amick, C. L., et al.(2023) Houston, Texas, [3] Heiken, G. H., et al.(1991) (778-778)

Cecilia L Amick↗

NASA NTRS2011

Lightweight, Ultra-High-Temperature, CMC-Lined Carbon/Carbon Structures

Carbon/carbon (C/C) is an established engineering material used extensively in aerospace. The beneficial properties of C/C include high strength, low density, and toughness. Its shortcoming is its limited usability at temperatures higher than the oxidation temperature of carbon . approximately 400 C. Ceramic matrix composites (CMCs) are used instead, but carry a weight penalty. Combining a thin laminate of CMC to a bulk structure of C/C retains all of the benefits of C/C with the high temperature oxidizing environment usability of CMCs. Ultramet demonstrated the feasibility of combining the light weight of C/C composites with the oxidation resistance of zirconium carbide (ZrC) and zirconium- silicon carbide (Zr-Si-C) CMCs in a unique system composed of a C/C primary structure with an integral CMC liner with temperature capability up to 4,200 F (.2,315 C). The system effectively bridged the gap in weight and performance between coated C/C and bulk CMCs. Fabrication was demonstrated through an innovative variant of Ultramet fs rapid, pressureless melt infiltration processing technology. The fully developed material system has strength that is comparable with that of C/C, lower density than Cf/SiC, and ultra-high-temperature oxidation stability. Application of the reinforced ceramic casing to a predominantly C/C structure creates a highly innovative material with the potential to achieve the long-sought goal of long-term, cyclic high-temperature use of C/C in an oxidizing environment. The C/C substructure provided most of the mechanical integrity, and the CMC strengths achieved appeared to be sufficient to allow the CMC to perform its primary function of protecting the C/C. Nozzle extension components were fabricated and successfully hot-fire tested. Test results showed good thermochemical and thermomechanical stability of the CMC, as well as excellent interfacial bonding between the CMC liner and the underlying C/C structure. In particular, hafnium-containing CMCs on C/C were shown to perform well at temperatures exceeding 3,500 F (.1,925 C). The melt-infiltrated CMC-lined C/C composites offered a lower density than Cf/SiC. The melt-infiltrated composites offer greater use temperature than Cf/SiC because of the more refractory ceramic matrices and the C/C substructure provides greater high-temperature strength. The progress made in this work will allow multiple high-temperature components used in oxidizing environments to take advantage of the low density and high strength of C/C combined with the high-temperature oxidation resistance of melt-infiltrated CMCs.

Wright, Matthew J.↗

NASA NTRS2017

Constraints on Mechanisms of Chondrule Formation from Chondrule Precursors and Chronology of Transient Heating Events in the Protoplanetary Disk

The mineralogy, petrography, and oxygen-isotope compositions of porphyritic chondrules—dominant chondrule type in most chondrite groups—suggest formation by incomplete melting of isotopically diverse precursors during localized transient heating events in dust-rich regions of the protoplanetary disk characterized by 16O-poor compositions (D17Odust+gas~ –7‰ to +4‰) relative to the inferred Sun’s value (D17O ~ –28 ± 2‰). The chondrule precursors included Ca,Al-rich inclusions (CAIs), amoeboid olivine aggregates (AOAs), chondrules of earlier generations, fine-grained matrix-like material, and possibly fragments of pre-existing planetesimals. Like porphyritic chondrules, igneous CAIs formed by melting of isotopically diverse precursors during transient heating events, but in an isotopically distinct, solar-like reservoir of the protoplanetary disk (D17O dust + gas ~ –24‰), probably near the proto Sun. Based on a narrow range of the initial 26Al/27Alratios inferred from the internal Al-Mg isochrons in igneous CAIs, their melting started at the very beginning of Solar System formation (t0), defined by the CV CAIs with U-corrected Pb-Pb age of 4567.3 ± 0.16 Ma and the canonical 26Al/27Al ratio of (5.25 ± 0.02) ¥ 10–5, and lasted at least 0.3 Ma. The U-corrected Pb-Pb absolute and 26Al-26Mg relative ages of porphyritic chondrules from type 3 ordinary, CO, CV, and CR carbonaceous chondrites (assuming uniform distribution of 26Al in the disk at the canonical level) suggest chondrule formation started at t0 and lasted for about 4 Ma. These observations may preclude formation of the majority of porphyritic chondrules by splashing of differentiated planetesimals and by collisions between planetesimals; instead, they are consistent with melting of dust balls by bow shocks or magnetized turbulence in the disk. Some porphyritic chondrules in equilibrated (petrologic type 4–6) ordinary chondrites contain relict fragments of coarse-grained chromite, ilmenite, phosphates, and albitic plagioclase. The similar mineral assemblage is commonly observed in type 4–6 ordinary chondrites, but is absent in type 3 chondrites, suggesting these chondrules formed by incomplete melting of thermally metamorphosed ordinary chondrite material, possibly by impacts. The CB metal-rich carbonaceous chondrites contain exclusively magnesian non-porphyritic chondrules crystallized from complete melts. These chondrules formed in a gas-melt plume generated by a hypervelocity (≥20 km/s) collision between planetesimals ~4.8 Ma after t0 in a transition or a debris disk. One of the colliding bodies was probably differentiated. The CH metal-rich carbonaceous chondrites contain chondrules formed by different mechanisms. The magnesian non-porphyritic chondrules formed in the CB impact plume ~4.8 Ma after t0. The chemically diverse (magnesian, ferroan, and Al-rich) porphyritic chondrules formed by incomplete melting of isotopically diverse precursors in the protoplanetary disk, most likely prior the CB impact plume event. We conclude that there are multiple mechanisms of chondrule formation that operated over the entire life-time of the disk.

Chondrules↗

NASA NTRSDate not supplied

NASA Delay Tolerant Networks: Operational, Evolving, an Ready for Expansion

The future of humanity’s presence beyond Earth depends on the successful commercialization of space. For commercialization to succeed, companies need cost-efficient architectures to support their business models and minimize risks for human capital, design, development, and operations. An ongoing challenge to any space enterprise is the reality that terrestrial network technologies are insufficient to provide reliable communications between assets in space. Whether you need to ensure your valuable data is safely transmitted to the ground or reliably delivered between platforms in orbit, ensuring data integrity over intermittent communication links is a necessity. Current solutions to space communications rely heavily on manual recording, storing, and retrieval of data from spacecraft. The current standard in space communication protocols, Consultative Committee for Space Data Systems (CCSDS) Space Packet standard, is reliant on inflexible network architectures based around mission-critical infrastructure to ensure data delivery. However, by automating the recording, storing, retrieval, and verification of data with Delay Tolerant Networks (DTN), the operator is freed from the dependence on manual data management and expensive mission critical infrastructure. NASA has been developing delay tolerant systems since the late 1990’s. Multiple DTN implementations have been established during that time, each suited to different use cases. Most notably, the DTN deployment for the International Space Station (ISS) includes demonstration of two DTN technologies: Interplanetary Overlay Network (ION) and Delay Tolerant Network Marshall Enterprise (DTNME). Beyond ISS, there are even more NASA DTN deployments being considered. Now that DTN implementations are maturing, it is appropriate to reflect upon these decades of work, review the integration and performance of the existing ISS deployment, and explore the future possibilities for DTN deployment industry-wide. The ISS DTN deployment is a complex architecture consisting of different DTN implementations for the onboard and ground network environments. The ION DTN implementation is being used in the on-board network. The Huntsville Operations Support Center (HOSC) DTN implementation, DTNME, is used by the ground network supporting ISS and will soon be a second onboard gateway too. The two implementations work cooperatively to provide high fidelity data services to flight operations users and payload developers across the globe. Though the two implementations yield a quality service, limitations are evident. Data rate, data storage, and device management are constrained by the services themselves and the complex nature of the deployment. Evolution of operations concepts will improve system capabilities and stability, but significant improvement will require additional development to the implementations themselves and to the overall deployment architecture. Taking advantage of the ongoing development and operation of the ISS DTN service will be central to the success of the future evolutions of NASA DTN deployments while demonstrating the benefits of DTN’s low-cost reliable data communication protocols for the growing commercial space industry. A broad effort on DTN integration and support is necessary to promote expansion beyond existing applications. NASA is developing several useful DTN implementations across a number of different systems: ION, DTNME, High-Rate DTN (HDTN), Bundle Protocol Library (BPLib), and others. To prevent fragmentation, DTN implementation teams need to communicate, collaborate, and integrate with one another to build a solid operational foundation for new DTN deployments. The establishment of a group that can assist new DTN users with understanding the purpose of each DTN implementation, provide best practices, and serve as a general knowledge base is paramount. Potential use of DTN on Gateway and other future NASA missions further drives the need for streamlined communication between DTN implementation teams. A well-integrated and highly engaged NASA DTN working group should help provide system architects the best DTN solutions for future commercial space efforts. This paper will first review the history of DTN implementations, explore the shortcoming of current space networking solutions given available limits in technology, and therefore establish the need for Delay Tolerant Networking in space communications. Secondly, the authors will explore NASA’s array of DTN implementations and highlight their usefulness to space applications. Thirdly, this paper will establish general DTN implementation distinguishing factors. Fourthly, the authors will discuss attempts to create a generic DTN comparison matrix, and the authors will review potential future topics in DTN innovation and collaboration, highlighting several key future efforts. Finally, this paper will describe how the institution of a NASA DTN Working Group will benefit DTN adoption across the governmental and commercial space sector. The goal of this paper is to encourage enthusiasm for DTN, share strategies for improving DTN on both current and future applications, promote the collaboration of DTN implementation groups within the international space operations community, and open the conversations about DTN, priorities, complexities, and innovation to the wider spaceflight industry.

DTN↗

NASA NTRSDate not supplied

NASA Delay Tolerant Networks: Operational, Evolving, an Ready for Expansion

DTN↗

NASA NTRS2024

Hybrid Thermally Efficient Core (HyTEC) HyTEC Phase 1 – Combustor Final Report

The objective of the HyTEC – Combustor Technology project is to develop technology for a compact, low emissions, rich-burn combustor that maintains a high-level of durability. To accomplish this, the combustor will incorporate CMC liners and a CMC Dome to enable improved air utilization for mixing and NOx reduction. In particular, three technology areas were matured: metallic swirler attachment to the CMC dome, Nickel diffusion in CMC, and EBC advancements with improved durability. The maturation of these technology areas is critical to successfully designing a combustor for an engine intended to demonstrate increased thermal efficiency with integrated high-power density-core engine technologies. As part of the project, multiple swirler attachment architectures were matured through the design process, manufactured, and tested in simulated engine conditions. These architectures utilized two different approaches to attachment: a pressure loaded approach and a clamped approach. Each approach required testing at simulated engine conditions to mature the technology. These tests were developed to ensure that the probable failure mechanisms for each architecture were assessed relative to its long-term durability. After successfully testing both approaches, the pressure loaded approach was chosen based on its simplicity and its similarity to the attachment approach used on lean burn combustors. Finally, this successful testing resulted in the swirler attachment technology achieving TRL 4 status. Another aspect of concern with the metallic swirler attachment is the potential for nickel diffusion into the CMC to cause durability shortfalls. A series of tests were completed to understand both the amount of nickel that is expected to diffuse into the CMC and the effect this diffusion will have on mechanical properties of the CMC. After successfully completing this testing, the reduction in local mechanical properties was found to not have an impact on the overall durability of the CMC Dome. To further improve the overall durability, three different EBC coating architectures were also studied. One of these coatings being the current state of the art for CMC liners with the other two being proposed improvements. Coupons were produced using each coating and subjected to a series of laboratory testing. One coating was found to have an overall superior durability in laboratory testing and survived rig testing without indication of degradation. This successful testing resulted in the EBC technology achieving TRL 5 status. Overall, a design approach was successfully down-selected that achieves both the TPM requirements and the need for TRL 4+ status. This design incorporates a pressure loaded swirler attachment with an improved EBC coating (Architecture 1). This combination allows the continued use of a CMC dome and results in an improved durability for the proposed compact core combustor.

CMC↗

NASA NTRS2016

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗

NASA NTRSDate not supplied

Thermodynamic constraints on the textural evolution of eucrite, EET 90020

Basaltic eucrites, which are possibly parts of the Vestan crust, formed as either lava flows or intrusions, and can offer insights into early crust formation. However most eucrites have experienced thermal metamorphism, which exacerbates the challenges of understanding these samples. A variety of petrogenetic models such as, partial melting of a primitive source, fractional crystallization of magmas emplaced in a primitive crust, and/or partial melting of a eucritic source coupled with melt mixing and assimilation have been proposed in order to explain observed petrologic and chemical characteristics of eucrites. EET 90020 is an unbrecciated eucrite that has experienced significant thermal metamorphism, with temperatures of metamorphic equilibration calculated to be ~840-1040°C. Its petrologic history remains contentious, in part due to differences in trace element analyses [10,17, this work], which has led to multiple metamorphic interpretations [8,10-11]. Here, we combined petrologic observations, chemical analyses, and thermodynamic modeling, to interpret micro-domain textures identified in EET 90020 and better constrain its petrologic history. Ultimately, the development of such textures is a direct result of the geologic processes operating on a young Vesta or similar asteroid. Sample description: We identified three textural domains in EET 90020 (Figure 1); (I) A coarse grain domain with granoblastic plagioclase and pyroxene, (II) a medium grain domain with curved grain boundaries between plagioclase and pyroxene, and (III) a fine grain domain dominated by tridymite that fills in interstitial space between spherical plagioclase. Major element chemistry in phases is consistent across domains. Pyroxene have distinctive high Ca lamellae (Wo40.92En28.1Fs31.0) forming from the low-Ca pigeonite host (Wo3.5En33.7Fs62.8). There are small amounts of fayalitic olivine along pyroxene/oxide grain boundaries, where pyroxenes have reacted with ilmenite/chromite clasts. Plagioclase is anorthitic (Avg ~An88), with little variation (An86-92.5). Trace element bulk rock analyses were collected via ICP-MS and individual mineral analyses with LA-ICP-MS (Figure 2). The bulk sample has an observable enrichment in light rare earth elements, with a notable depletion in Eu. Plagioclase is enriched in LREEs, depleted in HREEs, and has a positive Eu anomaly while pyroxene is enriched in HREEs with a negative Eu anomaly. Variations in trace element analyses is likely influenced by the presence of phosphates because are highly enriched in REEs and have a dramatic effect on the bulk domain compositions. Thermodynamic modeling: The software package, Perple_X, was used to calculate mineral phase equilibria over a range of conditions using a Gibbs free energy minimization approach [12]. Models used all major/minor elements except P2O5. Thermodynamic properties from [13] were used to constrain endmember phase stabilities. Activity models were used to describe mixing in phases with solid solution. Isochemical P-T phase diagrams were constructed for the bulk thin section, and the coarse, medium and fine grain domain compositions. Domain compositions were determined using microprobe analyses of phases identified in each domain and observed modal mineral abundancies. Thermodynamically stable pyroxene compositions were extracted from model results and compared to pyroxene compositions collected via EMPA, allowing us to calculate temperatures for metamorphic equilibrium. Results: Results that use the coarse domain composition (Fig. 3a) demonstrate overlap at T ~ 1020°C indicating high T equilibration. There was no overlap for the medium and fine grain domains, or the bulk composition (e.g., Fig. 3b) indicating that these are not equilibrium assemblages. These results were replicated using data from [8,20]. Discussion: The disequilibrium results from the bulk rock, fine grain and medium grain models imply that EET 90020 was modified during and/or after peak metamorphism (open system behavior). However, this is in contrast to the coarse domain where it appears as if equilibrium was maintained. We suggest that this conflict is explained by highly localized equilibrium occurring at millimeter length-scales, which has been observed in terrestrial metamorphism where aqueous fluid infiltrate [14] or melt loss occurs [15]. It is unlikely that aqueous alteration was significant in modifying the bulk composition of EET 90020, because it lacks hydrated minerals. However, the presence of melt is supported by the presence of an abundance of curved grain boundaries and spherical grains in the medium and fine grain domains This contrasts the coarse domain, where grain boundaries are angular. Given the evidence for disequilibrium in the medium and fine grain domains, it would be inappropriate to draw additional conclusions about textural evolution using thermodynamic models of those domains. However, modeling results of the coarse grain domain, where equilibrium was maintained, can provide useful insights into the geologic evolution of EET 90020. We approximate maximum metamorphic temperatures, using the Ca component in pyroxene, at ~1020°C (Fig. 4), which is consistent with previous two-pyroxene thermometry [11]. At this temperature, up to 10 vol % melt can be produced in the coarse grain domain, which is at the boundary of minimum melt needed in order to segregate and form a melt network [16]. Tridymite and plagioclase would be the first minerals to melt out during heating starting ~990°C. This is consistent with the observed melt textures formed by tridymite and plagioclase in the fine grain domain. We suggest that a partial melting model best explains the textural development of EET 90020. Upon heating, melt generated was either trapped, resulting in no net change in the localized rock composition, or migrated, resulting in disequilibrium between pyroxene and the surrounding matrix. A positive Eu anomaly in plagioclase indicates that all of the plagioclase was melted previously. This is consistent with previous studies that have suggested that EET 90020 represents a residual eucrite that experienced partial melting and subsequent melt loss [8]. We speculate that partial melting could have occurred when EET 90020 was heated in the asteroid’s lower crust or adjacent to an intrusion [7]. Future work will study other, texturally heterogeneous samples to see if EET 90020 presents a unique circumstance or whether crustal partial melting was ubiquitous on the eucrite asteroid.

Eucrites↗

NASA NTRSDate not supplied

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

NASA NTRSDate not supplied

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗