Search NASA⌕ Search

SEARCH · Search NASA

Results for “encounter geometry”

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.

280 records · Page 16

NASA Perspective on Requirements for Development of Advanced Methods Predicting Unsteady Aerodynamics and Aeroelasticity

Over the past three years, the National Aeronautics and Space Administration (NASA) has initiated design, development, and testing of a new human-rated space exploration system under the Constellation Program. Initial designs within the Constellation Program are scheduled to replace the present Space Shuttle, which is slated for retirement within the next three years. The development of vehicles for the Constellation system has encountered several unsteady aerodynamics challenges that have bearing on more traditional unsteady aerodynamic and aeroelastic analysis. This paper focuses on the synergy between the present NASA challenges and the ongoing challenges that have historically been the subject of research and method development. There are specific similarities in the flows required to be analyzed for the space exploration problems and those required for some of the more nonlinear unsteady aerodynamic and aeroelastic problems encountered on aircraft. The aggressive schedule, significant technical challenge, and high-priority status of the exploration system development is forcing engineers to implement existing tools and techniques in a design and application environment that is significantly stretching the capability of their methods. While these methods afford the users with the ability to rapidly turn around designs and analyses, their aggressive implementation comes at a price. The relative immaturity of the techniques for specific flow problems and the inexperience with their broad application to them, particularly on manned spacecraft flight system, has resulted in the implementation of an extensive wind tunnel and flight test program to reduce uncertainty and improve the experience base in the application of these methods. This provides a unique opportunity for unsteady aerodynamics and aeroelastic method developers to test and evaluate new analysis techniques on problems with high potential for acquisition of test and even flight data against which they can be evaluated. However, researchers may be required to alter the geometries typically used in their analyses, the types of flows analyzed, and even the techniques by which computational tools are verified and validated. This paper discusses these issues and provides some perspective on the potential for new and innovative approaches to the development of methods to attack problems in nonlinear unsteady aerodynamics.

Schuster, David M.↗

The Surface Compositions of Pluto and Charon

The surface of Pluto as it is understood on the eve of the encounter of the New Horizons spacecraft (mid-2015) consists of a spatially heterogeneous mix of solid N2, CH4, CO, C2H6, and an additional component that imparts color, and may not be an ice. The known molecular ices are detected by near-infrared spectroscopy. The N2 ice occurs in the hexagonal crystalline b-phase, stable at T greater than 35.6 K. Spectroscopic evidence for wavelength shifts in the CH4 bands attests to the complex mixing of CH4 and N2 in the solid state, in accordance with the phase diagram for N2 + CH4. Spectra obtained at several aspects of Plutos surface as the planet rotates over its 6.4-day period show variability in the distribution of CH4 and N2 ices, with stronger CH4 absorption bands associated with regions of higher albedo, in correlation with the visible rotational light curve. CO and N2 ice absorptions are also strongly modulated by the rotation period; the bands are strongest on the anti-Charon hemisphere of Pluto. Longer term changes in the strengths of Plutos absorption bands occur as the viewing geometry changes on seasonal time-scales, although a complete cycle has not been observed. The non-ice component of Plutos surface may be a relatively refractory material produced by the UV and cosmic-ray irradiation of the surface ices and gases in the atmosphere, although UV does not generally penetrate the atmospheric CH4 to interact with the surface. Laboratory simulations indicate that a rich chemistry ensues by the irradiation of mixtures of the ices known to occur on Pluto, but specific compounds have not yet been identified in spectra of the planet. Charons surface is characterized by spectral bands of crystalline H2O ice, and a band attributed to one or more hydrates of NH3. Amorphous H2O ice may also be present; the balance between the amorphization and crystallization processes on Charon remains to be clarified. The albedo of Charon and its generally spatially uniform neutral color indicate that a component, not yet identified, is mixed in some way with the H2O and NH3nH2O ices. Among the many known small bodies in the transneptunian region, several share characteristics with Pluto and Charon, including the presence of CH4, N2, C2H6, H2O ices, as well as components that yield a wide variety of surface albedo and color. The New Horizons investigation of the Pluto-Charon system will generate new insight into the physical properties of the broader transneptunian population, and eventually to the corresponding bodies expected in the numerous planetary systems currently being discovered elsewhere in the Galaxy.

surface composition↗

Unveiling the Journey of A Highly Inclined CME: Insights from the March 13, 2012, Event With 110° Longitudinal Separation

Context. A fast (∼2000 km s -1 ) and wide (>110°) coronal mass ejection (CME) erupted from the Sun on March 13, 2012. Its interplanetary counterpart was detected in situ two days later by STEREO-A and near-Earth spacecraft, such as ACE, Wind, and Cluster. We suggest that at 1 au the CME extended at least 110° in longitude, with Earth crossing its east flank and STEREO-A crossing its- west flank. Despite their separation, measurements from both positions showed very similar in situ CME signatures. The solar source region where the CME erupted was surrounded by three coronal holes (CHs). Their locations with respect to the CME launch site were east (negative polarity), southwest (positive polarity) and west (positive polarity). The solar magnetic field polarity of the area covered by each CH matches that observed at 1 au in situ. Suprathermal electrons at each location showed mixed signatures with only some intervals presenting clear counter streaming flows as the CME transits both locations. The strahl population coming from the shortest magnetic connection of the structure to the Sun showed more intensity. Aims. The aim of this work is to understand the propagation and evolution of the CME and its interaction with the surrounding CHs, to explain the similarities and differences between the observations at each spacecraft, and report what one of the most longitudinal expanded CME structures measured in situ would be. Methods. Known properties of the large-scale structures from a variety of catalogues and previous studies were used to have a better overview of this particular event. In addition, multipoint observations were used to reconstruct the 3D geometry of the CME and determine the context of the solar and heliospheric conditions before the CME eruption and during its propagation. The graduated cylindrical shell model (GCS) was used to reproduce the orientation, size and speed of the structure with a simple geometry. Also, the Drag-Based Model (DBM) was utilised to understand the conditions of the interplanetary medium better in terms of the drag undergone by the structure while propagating in different directions. Finally, a comparative analysis of the different regions of the structure through the different observatories was carried out in order to directly compare the in situ plasma and magnetic field properties at each location. Results. The study presents important findings regarding the in situ measured CME on March 15, 2012, detected at a longitudinal separation of 110° in the ecliptic plane despite its initial inclination being around 45° when erupted (March 13). This suggests that the CME may have deformed and/or rotated, allowing it to be observed near its legs with spacecraft at a separation angle greater than 100°. The CME structure interacted with high-speed streams generated by the surrounding CHs. The piled-up plasma in the sheath region exhibited an unexpected correlation in magnetic field strength despite the large separation in longitude. In situ observations reveal that at both locations there was a flank encounter – where the spacecraft crossed the first part of the CME – then encountered ambient solar wind, and finally passed near the legs of the structure. Conclusions. A scenario covering all evidence is proposed for both locations with a general view of the whole structure and solar wind conditions. Also, the study shows the necessity of having multipoint observations of large-scale structures in the heliosphere.

methods: data analysis↗

Progress in Grid Generation: From Chimera to DRAGON Grids

Hybrid grids, composed of structured and unstructured grids, combines the best features of both. The chimera method is a major stepstone toward a hybrid grid from which the present approach is evolved. The chimera grid composes a set of overlapped structured grids which are independently generated and body-fitted, yielding a high quality grid readily accessible for efficient solution schemes. The chimera method has been shown to be efficient to generate a grid about complex geometries and has been demonstrated to deliver accurate aerodynamic prediction of complex flows. While its geometrical flexibility is attractive, interpolation of data in the overlapped regions - which in today's practice in 3D is done in a nonconservative fashion, is not. In the present paper we propose a hybrid grid scheme that maximizes the advantages of the chimera scheme and adapts the strengths of the unstructured grid while at the same time keeps its weaknesses minimal. Like the chimera method, we first divide up the physical domain by a set of structured body-fitted grids which are separately generated and overlaid throughout a complex configuration. To eliminate any pure data manipulation which does not necessarily follow governing equations, we use non-structured grids only to directly replace the region of the arbitrarily overlapped grids. This new adaptation to the chimera thinking is coined the DRAGON grid. The nonstructured grid region sandwiched between the structured grids is limited in size, resulting in only a small increase in memory and computational effort. The DRAGON method has three important advantages: (1) preserving strengths of the chimera grid; (2) eliminating difficulties sometimes encountered in the chimera scheme, such as the orphan points and bad quality of interpolation stencils; and (3) making grid communication in a fully conservative and consistent manner insofar as the governing equations are concerned. To demonstrate its use, the governing equations are discretized using the newly proposed flux scheme, AUSM+, which will be briefly described herein. Numerical tests on representative 2D inviscid flows are given for demonstration. Finally, extension to 3D is underway, only paced by the availability of the 3D unstructured grid generator.

Liou, Meng-Sing↗

Adaptive Mesh Refinement for Microelectronic Device Design

Finite element and finite volume methods are used in a variety of design simulations when it is necessary to compute fields throughout regions that contain varying materials or geometry. Convergence of the simulation can be assessed by uniformly increasing the mesh density until an observable quantity stabilizes. Depending on the electrical size of the problem, uniform refinement of the mesh may be computationally infeasible due to memory limitations. Similarly, depending on the geometric complexity of the object being modeled, uniform refinement can be inefficient since regions that do not need refinement add to the computational expense. In either case, convergence to the correct (measured) solution is not guaranteed. Adaptive mesh refinement methods attempt to selectively refine the region of the mesh that is estimated to contain proportionally higher solution errors. The refinement may be obtained by decreasing the element size (h-refinement), by increasing the order of the element (p-refinement) or by a combination of the two (h-p refinement). A successful adaptive strategy refines the mesh to produce an accurate solution measured against the correct fields without undue computational expense. This is accomplished by the use of a) reliable a posteriori error estimates, b) hierarchal elements, and c) automatic adaptive mesh generation. Adaptive methods are also useful when problems with multi-scale field variations are encountered. These occur in active electronic devices that have thin doped layers and also when mixed physics is used in the calculation. The mesh needs to be fine at and near the thin layer to capture rapid field or charge variations, but can coarsen away from these layers where field variations smoothen and charge densities are uniform. This poster will present an adaptive mesh refinement package that runs on parallel computers and is applied to specific microelectronic device simulations. Passive sensors that operate in the infrared portion of the spectrum as well as active device simulations that model charge transport and Maxwell's equations will be presented.

Cwik, Tom↗

Source and Propagation of a Streamer Blowout Coronal Mass Ejection Observed by the Parker Solar Probe

In the first orbit of the Parker Solar Probe (PSP), in situ thermal plasma and magnetic field measurements were collected as close as 35RSun from the Sun, an environment that had not been previously explored. During the first orbit of PSP, the spacecraft flew through a streamer blowout coronal mass ejection (SBO-CME) on 2018 November 11 at 23:50 UT as it exited the science encounter. The SBO-CME on November 11 was directed away from the Earth and was not visible by L1 or Earth-based telescopes due to this geometric configuration. However, PSP and the STEREO-A spacecraft were able to make observations of this slow (v ≈ 380 kms−1) SBO-CME. Using the PSP data, STEREO-A images, and Wang–Sheeley–Arge model, the source region of the CME is found to be a helmet streamer formed between the northern polar coronal hole and a mid-latitude coronal hole. Using the YGUAZU-A model, the propagation of the CME is traced from the source at the Sun to PSP. This model predicts the travel time of the flux rope to the PSP spacecraft as 30 hr, which is within 0.33 hr of the actual measured arrival time. The in situ Solar Wind Electrons Alphas and Protons data were examined to determine that no shock was associated with this SBO-CME. Modeling of the SBO-CME shows that no shock was present at PSP; however, at other positions along the SBO-CME front, a shock could have formed. The geometry of the event requires in situ and remote sensing observations to characterize the SBO-CME and further understand its role in space weather.

Kelly Elizabeth Korreck↗

Particle signatures of magnetic topology at the magnetopause: AMPTE/CCE observations

Electron distributions at energies above 50 eV have been found to be a sensitive indicator of magnetic topology for magnetopause crossings of the AMPTE/CCE spacecraft. Progressing from the magnetosheath to the magnetosphere two abrupt transitions occur. First, the magnetosheath electron population directed either parallel or antiparallel to the magnetic field is replaced by a streaming, heated magnetosheath electron population. The other half of the distribution is unchanged. The region with unidirectional, heated magnetosheath electrons is identified as the magnetosheath boundary layer (MSBL). Second, the unheated magnetosheath electron population is replaced by a heated population nearly identical to the population encountered in the MSBL, resulting in a symmetric counterstreaming distribution. The region populated by the bidirectional heated magnetosheath electrons is identified as the low-latitude boundary layer (LLBL). The MSBL and LLBL identified by the electron transitions are the same as the regions identified using ion composition measurements. The magnetosheath-MSBL transition reflects a change in magnetic topology from a solar wind field line to one that threads the magnetopause, and the existence of a magnetosheath-MSBL transition implies that the magnetopause is open. When the current layer is easily identified, the MSBL-LLBL transition coincides with the magnetopause current layer, indicating that the magnetosheath electrons are heated in the current layer. Both magnetosheath-MSBL and MSBL-LLBL transitions are observed for low as well as high magnetic shears. Moreover, the transitions are particularly clear for low shear implying that magnetic topology boundaries are sharp even when abrupt changes in the field and other plasma parameters are absent. Furthermore, for low magnetic shear, solar wind ions with low parallel drift speeds make up the majority of the LLBL population indicating that the magnetosheath plasma has convected directly across the magnetosheath plasma has converted directly across the magnetopause. These observations are consistent with quasi-steady, high-latitude reconnection and indicate that the signatures of this reconnection geometry are commonly present in the subpolar region.

Fuselier, S. A.↗

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes

Cell-centered and node-centered approaches have been compared for unstructured finite-volume discretization of inviscid fluxes. The grids range from regular grids to irregular grids, including mixed-element grids and grids with random perturbations of nodes. Accuracy, complexity, and convergence rates of defect-correction iterations are studied for eight nominally second-order accurate schemes: two node-centered schemes with weighted and unweighted least-squares (LSQ) methods for gradient reconstruction and six cell-centered schemes two node-averaging with and without clipping and four schemes that employ different stencils for LSQ gradient reconstruction. The cell-centered nearest-neighbor (CC-NN) scheme has the lowest complexity; a version of the scheme that involves smart augmentation of the LSQ stencil (CC-SA) has only marginal complexity increase. All other schemes have larger complexity; complexity of node-centered (NC) schemes are somewhat lower than complexity of cell-centered node-averaging (CC-NA) and full-augmentation (CC-FA) schemes. On highly anisotropic grids typical of those encountered in grid adaptation, discretization errors of five of the six cell-centered schemes converge with second order on all tested grids; the CC-NA scheme with clipping degrades solution accuracy to first order. The NC schemes converge with second order on regular and/or triangular grids and with first order on perturbed quadrilaterals and mixed-element grids. All schemes may produce large relative errors in gradient reconstruction on grids with perturbed nodes. Defect-correction iterations for schemes employing weighted least-square gradient reconstruction diverge on perturbed stretched grids. Overall, the CC-NN and CC-SA schemes offer the best options of the lowest complexity and secondorder discretization errors. On anisotropic grids over a curved body typical of turbulent flow simulations, the discretization errors converge with second order and are small for the CC-NN, CC-SA, and CC-FA schemes on all grids and for NC schemes on triangular grids; the discretization errors of the CC-NA scheme without clipping do not converge on irregular grids. Accurate gradient reconstruction can be achieved by introducing a local approximate mapping; without approximate mapping, only the NC scheme with weighted LSQ method provides accurate gradients. Defect correction iterations for the CC-NA scheme without clipping diverge; for the NC scheme with weighted LSQ method, the iterations either diverge or converge very slowly. The best option in curved geometries is the CC-SA scheme that offers low complexity, second-order discretization errors, and fast convergence.

Diskin, Boris↗

Modern shelf ice, equatorial Aeolis Quadrangle, Mars

As part of a detailed study of the geological and geomorphological evolution of Aeolis Quadrangle, I have encountered evidence suggesting that near surface ice exists at low latitudes and was formed by partial or complete freezing of an inland sea. The area of interest is centered at approximately -2 deg, 196 deg. As seen in a suite of Viking Orbiter frames obtained at a range of approximately 600 km, the plains surface at this location is very lightly cratered or uncratered, and it is thus of late Amazonian age. Extant topographic data indicate that the Amazonian plains at this location occupy a trough whose surface lies at least 1000 m below the Mars datum. A reasonable hypothesis is that quite recent surface water releases, perhaps associated with final evolution of large 'outflow chasms' to the south, but possibly from other source areas, filled this trough, that ice floes formed almost immediately, and that either grounded ice or an ice-covered sea still persists. A reasonable hypothesis is that quite recent surface water releases, perhaps associated with final evolution of large 'outflow chasms' to the south, but possibly from other source areas, filled this trough, that ice floes formed almost immediately, and that either grounded ice or an ice-covered sea still persists. In either case, the thin (a few meters at most) high albedo, low thermal inertia cover of aeolian materials was instrumental in allowing ice preservation, and at least the lower portions of this dust cover may be cemented by water ice. Detailed mapping using Viking stereopairs and quantitative comparisons to terrestrial shelf ice geometries are underway.

Brakenridge, G. R.↗

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↗