Search NASA⌕ Search

SEARCH · Search NASA

Results for “best source selection”

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.

188 records · Page 11

Airborne Laser Polar Nephelometer

A polar nephelometer has been developed at NCAR to measure the angular variation of the intensity of light scattered by air molecules and particles. The system has been designed for airborne measurements using outside air ducted through a 5-cm diameter airflow tube; the sample volume is that which is common to the intersection of a collimated source beam and the detector field of view within the airflow tube. The source is a linearly polarized helium-neon laser beam. The optical system defines a collimated field-of-view (0.5deg half-angle) through a series of diaphragms located behind a I72-mm focal length objective lens. A photomultiplier tube is located immediately behind an aperture in the focal plane of the objective lens. The laser beam is mechanically chopped (on-off) at a rate of 5 Hz; a two-channel pulse counter, synchronized to the laser output, measures the photomultiplier pulse rate with the light beam both on and off. The difference in these measured pulse rates is directly proportional to the intensity of the scattered light from the volume common to the intersection of the laser beam and the detector field-of-view. Measurements can be made at scattering angles from 15deg to 165deg with reference to the direction of propagation of the light beam. Intermediate angles are obtained by selecting the angular increments desired between these extreme angles (any multiple of 0.1deg can be selected for the angular increment; 5deg is used in normal operation). Pulses provided by digital circuits control a stepping motor which sequentially rotates the detector by pre-selected angular increments. The synchronous photon-counting system automatically begins measurement of the scattered-light intensity immediately after the rotation to a new angle has been completed. The instrument has been flown on the NASA Convair 990 airborne laboratory to obtain data on the complex index of refraction of atmospheric aerosols. A particle impaction device is operated simultaneously to collect particles from the same airflow tube used to make the scattered-light measurements. A size distribution function is obtained by analysis of the particles collected by the impaction device. Calculated values of the angular variation of the scattered-light intensity are obtained by applying Mie scattering theory to the observed size distribution function and assuming different values of the complex index of refraction of the particles. The calculated values are then compared with data on the actual variation of the scattered-light intensity obtained with the polar nephelometer. The most probable value of the complex refractive index is that which provides the best fit between the experimental light scattering data and the parameters calculated from the observed size distribution function.

Grams, Gerald W.↗

AI-3000K Infrared Line List for Hot CO 2

AI-3000K is a semi-empirical IR line list constructed for hot CO 2 spectra analysis and simulation up to 3000 – 4000 K. Compared to previously published Ames-1000K and UCL-4000, it represents a major upgrade, utilizing a new algorithm for optimization and including the latest improvements in potential energy surface (PES), dipole moment surface (DMS), and room temperature IR line list (Ames-2021 296K). To maximize the success of introducing experimental based accurate line positions, a new PES (X01d) was refined with respect to >800 selected CDSD2019 [Tashkun et al. JQSRT (2019) 228, 124] energy levels in the range of 0 – 24,000 cm -1 , with σ rms = 0.5 – 0.7 cm -1 . Most differences between the X01d PES based levels and CDSD2019 energies are within ± 2 cm -1 . A new DMS is fitted from extrapolated CCSD(T)/aug-cc-pV(T,Q,5)Z dipole calculations, with σ rms = 5.1 × 10 -6 au for 11,155 geometries up to 40,000 cm -1 , denoted Ames-2021–40K. Compared to the best available Ames-2021 DMS and room temperature IR line list [Huang et al J. Phys. Chem. A (2022) 126, 5940], the relative intensity differences are expected to be ∼ 1 ‰. The line position accuracy of ''X01d + Ames-2021–40K'' IR line list is significantly improved by adopting CDSD2019 energy levels up to 24,000 cm -1 (J ≤ 150). The Einstein A21 coefficients for E' < 15,000 cm -1 transitions are replaced by more accurate values from the Ames-2021 296K IR line list. In short, the AI-3000K is the X01d PES and Ames-2021–40K DMS based line list enhanced with the A 21 of Ames-2021 296K line list and CDSD energy levels. It provides continuous coverage from 0 to 20,000 cm -1 for the four most abundant isotopologues: 12 C 16 O 2 (626), 13 C 16 O 2 (636), 16 O 12 C 18 O (628), and 16 O 12 C 17 O (627). The impacts of isotopologue and E'/E'' cutoffs have been examined. Intensity convergence (not accuracy) of AI-3000K line list is quantitatively estimated in 1 cm -1 bins. It is better than 99% in the whole range of 0 – 20,000 (or 10,000) cm -1 at 1000 K (or 2000 K), or better than 90% in the whole range of 0 – 15,000 (or 9000) cm -1 at 2000 K (or 3000 K), respectively. Convergence beyond 3000 K will require a new PES and DMS for E' > 40,000 cm -1 . The AI-3000K and HITEMP based IR simulations are compared to high resolution shock tube experiments for CO 2 -Ar mixture up to 2000 K. With line position accuracy comparable to that of HITEMP, AI-3000K IR line list yields better agreements at more frequencies. Potential sources of discrepancies with experiment are discussed.

AI-3000K↗

Laser Sounder for Global Measurement of CO2 Concentrations in the Troposphere from Space

Measurements of tropospheric CO2 abundance with global-coverage, a few hundred km spatial and monthly temporal resolution are needed to quantify processes that regulate CO2 storage by the land and oceans. The Orbiting Carbon Observatory (OCO) is the first space mission focused on atmospheric CO2 for measuring total column CO, and O2 by detecting the spectral absorption in reflected sunlight. The OCO mission is an essential step, and will yield important new information about atmospheric CO2 distributions. However there are unavoidable limitations imposed by its measurement approach. These include best accuracy only during daytime at moderate to high sun angles, interference by cloud and aerosol scattering, and limited signal from CO2 variability in the lower tropospheric CO2 column. We have been developing a new laser-based technique for the remote measurement of the tropospheric CO2 concentrations from orbit. Our initial goal is to demonstrate a lidar technique and instrument technology that will permit measurements of the CO2 column abundance in the lower troposphere from aircraft. Our final goal is to develop a space instrument and mission approach for active measurements of the CO2 mixing ratio at the 1-2 ppmv level. Our technique is much less sensitive to cloud and atmospheric scattering conditions and would allow continuous measurements of CO2 mixing ratio in the lower troposphere from orbit over land and ocean surfaces during day and night. Our approach is to use the 1570nm CO2 band and a 3-channel laser absorption spectrometer (i.e. lidar used an altimeter mode), which continuously measures at nadir from a near polar circular orbit. The approach directs the narrow co-aligned laser beams from the instrument's lasers toward nadir, and measures the energy of the laser echoes reflected from land and water surfaces. It uses several tunable fiber laser transmitters which allowing measurement of the extinction from a single selected CO2 absorption line in the 1570 nm band. This band is free from interference from other gases and has temperature insensitive absorption lines. During the measurement the lasers are tuned on- and off- a selected CO2 line near 1572 nm and a selected O2 line near 768 nm in the Oxygen A band at kHz rates. The lasers use tunable diode seed lasers followed by fiber amplifiers, and have spectral widths much narrower than the gas absorption lines. The receiver uses a 1-m diameter telescope and photon counting detectors and measures the background light and energies of the laser echoes from the surface. The extinction and column densities for the CO2 and O2 gases are estimated from the ratio of the on and offline surface echo via the differential optical absorption technique. Our technique rapidly alternates between several on-line wavelengths set to the sides of the selected gas absorption lines. It exploits the atmospheric pressure broadening of the lines to weight the measurement sensitivity to the atmospheric column below 5 km. This maximizes sensitivity to CO2 in the boundary layer, where variations caused by surface sources and sinks are largest. Simultaneous measurements of O2 column will use an identical approach with an O2 line. Thee laser frequencies are tunable and have narrow (MHz) line widths. In combination with sensitive photon counting detectors these enables much higher spectral resolution and precision than is possible with passive spectrometer. 1aser backscatter profiles are also measured, which permits identifying measurements made to cloud tops and through aerosol layers. The measurement approach using lasers in common-nadir-zenith path allows retrieving CO2 column mixing ratios in the lower troposphere irrespective of sun angle. Pulsed laser signals, time gated receiver and a narrow receiver field-of-view are used to isolate the surface laser echo signals and to exclude photons scattered from clouds and aerosols. Nonetheless, the optical absorption change due to a change of a few ppO2 is small, <1 % which makes achieving the needed measurement sensitivities and stabilities quite challenging. Measurement SNRs and stabilities of >600:1 are needed to estimate CO2 mixing ratio at the 1-2 ppm level. We have calculated characteristics of the technique and have demonstrated aspects of the laser, detector and receiver approaches in th e laboratory We have also measured O2 in an absorption cell, and made C02 measurements over a 400 m long (one way) horizontal path using a sensor breadboard. We will describe these and more details of our approach in the paper.

Abshire, James B.↗

Technology Development Roadmap: A Technology Development Roadmap for a Future Gravitational Wave Mission

Humankind will detect the first gravitational wave (GW) signals from the Universe in the current decade using ground-based detectors. But the richest trove of astrophysical information lies at lower frequencies in the spectrum only accessible from space. Signals are expected from merging massive black holes throughout cosmic history, from compact stellar remnants orbiting central galactic engines from thousands of close contact binary systems in the Milky Way, and possibly from exotic sources, some not yet imagined. These signals carry essential information not available from electromagnetic observations, and which can be extracted with extraordinary accuracy. For 20 years, NASA, the European Space Agency (ESA), and an international research community have put considerable effort into developing concepts and technologies for a GW mission. Both the 2000 and 2010 decadal surveys endorsed the science and mission concept of the Laser Interferometer Space Antenna (LISA). A partnership of the two agencies defined and analyzed the concept for a decade. The agencies partnered on LISA Pathfinder (LPF), and ESA-led technology demonstration mission, now preparing for a 2015 launch. Extensive technology development has been carried out on the ground. Currently, the evolved Laser Interferometer Space Antenna (eLISA) concept, a LISA-like concept with only two measurement arms, is competing for ESA's L2 opportunity. NASA's Astrophysics Division seeks to be a junior partner if eLISA is selected. If eLISA is not selected, then a LISA-like mission will be a strong contender in the 2020 decadal survey. This Technology Development Roadmap (TDR) builds on the LISA concept development, the LPF technology development, and the U.S. and European ground-based technology development. The eLISA architecture and the architecture of the Mid-sized Space-based Gravitational-wave Observatory (SGO Mid)-a competitive design with three measurement arms from the recent design study for a NASA-led mission after 2020-both use the same technologies. Further, NASA participation in an ESA-led mission would likely augment the eLISA architecture with a third arm to become the SGO Mid architecture. For these reasons, this TDR for a future GW mission applies to both designs and both programmatic paths forward. It is adaptable to the different timelines and roles for an ESA-led or a NASA-led mission, and it is adaptable to available resources. Based on a mature understanding of the interaction between technology and risk, the authors of this TDR have chosen a set of objectives that are more expansive than is usual. The objectives for this roadmap are: (1) reduce technical and development risks and costs; (2) understand and, where possible, relieve system requirements and consequences; (3) increase technical insight into critical technologies; and (4) validate the design at the subsystem level. The emphasis on these objectives, particularly the latter two, is driven by outstanding programmatic decisions, namely whether a future GW mission is ESA-led or NASA-led, and availability of resources. The relative emphasis is best understood in the context of prioritization.

Laser Interferometry Space Antenna↗

Impact of Powder Variability on the Microstructure and Mechanical Behavior of Selective Laser Melted Alloy 718

Powder-bed additive manufacturing processes use fine powders to build parts layer-by-layer. Alloy 718 powder feedstocks for selective laser melting (SLM) additive manufacturing are produced commercially by both gas and rotary atomization and are available typically in the 10-45 or 15-45 microns size ranges. A comprehensive investigation was conducted to understand the impact of powder variability on the microstructure and mechanical behavior of SLM 718 heat treated to Aerospace Material Specification (AMS) 5664. This study included sixteen virgin powders and three once-recycled powders within the 10-45 and 15-45 microns size ranges that were obtained from seven direct source suppliers and one reseller. Although alike as highly regular spheroids, these powders showed distinct differences in composition (especially Al, C and N contents), particle size distributions, and powder features such as degree of agglomeration, fusion and surface roughness. Compositional differences expectedly had the strongest impact on microstructure. High N and C contents formed TiN-nitrides and/or (Nb,Ti,Mo)-C carbides on the grain boundaries, prevented recrystallization during heat treatment, and resulted in retained (001)-scalloped shaped grains that ranged from 19 to 41 microns in average size. In the absence of this particle pinning, the average grain size of the heat treated SLM 718 ranged from 51 microns to 90 microns. Room temperature tensile and high cycle fatigue (HCF) testing compared as-fabricated (AF) and low stress ground (LSG) surface conditions. Tensile testing revealed consistent behavior between the two surface conditions and amongst the powder lots. The finer grained SLM 718 builds displayed the lowest tensile properties. A SLM 718 build fabricated from a powder with eight times lower C content showed statistically better tensile properties presumably due to enhanced coarsening of (delta)-Ni3Nb precipitates. The specimens from once-recycled powders had slightly higher tensile strengths and slightly higher ductility compared to their virgin equivalents; once-recycling also did not substantially degrade the mean HCF life. The LSG fatigue lives agreed with conventionally manufactured 718 data, while AF lives exhibited a knock-down due to surface roughness. The fatigue lives of AF specimens were statistically equivalent across powder lots except for one and failures typically initiated at stress concentrators associated with SLM surface asperities. Fatigue testing of low stress ground specimens result in both transgranular and within facet crack initiations. More than half of the cracks initiated from these facets for the machined condition; however, these facets appeared to be within grains that were larger-than-average in size. A nitrogen-atomized powder with fine prior particles of TiN-nitrides and M(Ti,Nb,Mo)C carbides from atomization on powder surfaces resulted in the best fatigue performance with segregation of these particles to the SLM 718 grain boundaries leading to higher resistance to early-stage crack propagation. Typically the fine-grained builds with minor phases along the grain boundaries did not perform well in fatigue, whereas a larger-grain build with lower carbon content and coarser delta-Ni3Nb precipitates showed the next best HCF response. Further details of the build microstructure and its impact on tensile and fatigue behavior was considered.

laser deposition↗

Robust Exploration and Commercial Missions to the Moon Using NTR LANTR Propulsion and Lunar-Derived Propellants

The nuclear thermal rocket (NTR) has frequently been identified as a key space asset required for the human exploration of Mars. This proven technology can also provide the affordable access through cislunar space necessary for commercial development and sustained human presence on the Moon. In his post-Apollo Integrated Space Program Plan (1970-1990), Wernher von Braun, proposed a reusable nuclear thermal propulsion stage (NTPS) to deliver cargo and crew to the Moon to establish a lunar base before undertaking human missions to Mars. The NTR option was selected by von Braun because it was a demonstrated technology capable of generating both high thrust and high specific impulse (Isp 900 s) twice that of todays best chemical rockets. In NASAs Mars Design Reference Architecture (DRA) 5.0 study, the crewed Mars transfer vehicle used three 25 klbf Pewee engines the smallest and highest performing engine tested in the Rover program along with graphite composite fuel. Smaller, lunar transfer vehicles consisting of a NTPS using three approximately 16.5 klbf Small Nuclear Rocket Engines (SNREs), an in-line propellant tank, plus the payload can enable a variety of reusable lunar missions. These include cargo delivery and crewed lunar landing missions. Even weeklong tourism missions carrying passengers into lunar orbit for a day of sightseeing and picture taking are possible. The NTR can play an important role in the next phase of lunar exploration and development by providing an affordable in-space lunar transportation system (LTS) that can allow initial outposts to evolve into settlements supported by a variety of commercial activities such as in-situ propellant production used to supply strategically located propellant depots and transportation nodes. The utilization of iron-rich volcanic glass or lunar polar ice (LPI) deposits (each estimated at billions of metric tons) for propellant production can significantly reduce the launch mass requirements from Earth and can enable reusable, surface-based lunar landing vehicles (LLVs) using liquid oxygen/hydrogen (LOX/LH2) chemical rocket engines. Afterwards, LOX/LH2 propellant depots can be established in lunar equatorial and polar orbits to supply the LTS. At this point a modified version of the conventional NTR called the LOX-augmented NTR, or LANTR would be introduced into the LTS allowing bipropellant operation and leveraging the mission benefits of refueling with lunar-derived propellants for Earth return. The bipropellant LANTR engine utilizes the large divergent section of its nozzle as an afterburner into which oxygen is injected and supersonically combusted with nuclear preheated hydrogen emerging from the engines choked sonic throat essentially scramjet propulsion in reverse. By varying the oxygen-to-hydrogen mixture ratio, LANTR engines can operate over a range of thrust and Isp values while the reactor core power level remains relatively constant. Eventually, a LANTR-based LTS can enable a rapid commuter shuttle with one-way trip times to and from the Moon ranging from 36 to 24 hours. Even if only 1 of the extracted propellant from identified volcanic glass and polar ice deposits were available for use in lunar orbit, such a supply could support daily commuter flights to the Moon for many thousands of years! An evolutionary mission architecture is outlined and a variety of lunar missions and transfer vehicle designs are examined, along with the increasing demands on propellant production as mission complexity increases. A comparison of vehicle features and engine operating characteristics, for both NTR and LANTR engines, is also provided along with a brief discussion on the propellant production issues associated with using volcanic glass and LPI as source material.

Spacecraft design↗

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↗

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↗