Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “feasible power flow solution”

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.

A zonal method for modeling powered-lift aircraft flow fields

A zonal method for modeling powered-lift aircraft flow fields is based on the coupling of a three-dimensional Navier-Stokes code to a potential flow code. By minimizing the extent of the viscous Navier-Stokes zones the zonal method can be a cost effective flow analysis tool. The successful coupling of the zonal solutions provides the viscous/inviscid interations that are necessary to achieve convergent and unique overall solutions. The feasibility of coupling the two vastly different codes is demonstrated. The interzone boundaries were overlapped to facilitate the passing of boundary condition information between the codes. Routines were developed to extract the normal velocity boundary conditions for the potential flow zone from the viscous zone solution. Similarly, the velocity vector direction along with the total conditions were obtained from the potential flow solution to provide boundary conditions for the Navier-Stokes solution. Studies were conducted to determine the influence of the overlap of the interzone boundaries and the convergence of the zonal solutions on the convergence of the overall solution. The zonal method was applied to a jet impingement problem to model the suckdown effect that results from the entrainment of the inviscid zone flow by the viscous zone jet. The resultant potential flow solution created a lower pressure on the base of the vehicle which produces the suckdown load. The feasibility of the zonal method was demonstrated. By enhancing the Navier-Stokes code for powered-lift flow fields and optimizing the convergence of the coupled analysis a practical flow analysis tool will result.

Roberts, D. W.↗

Lunar Power Transmission for Fission Surface Power

This work reviews key considerations and challenges in the implementation of lunar power transmission for the 40 kW Fission Surface Power (FSP) System. The overall FSP electric power flow is presented, and potential bulk power transmission strategies are discussed. Metallic conductors operating at elevated voltage are identified as the most efficient and feasible solution, and hardware implementation challenges for AC and DC voltage step up/down are discussed. This work is intended to drive progress on identifying the most feasible power transmission strategy for landed fission power and does not represent consensus or preference by NASA regarding a particular course of action.

Lunar↗

Lunar Power Transmission for Fission Surface Power

This work reviews key considerations and challenges in the implementation of lunar power transmission for the 40 kW Fission Surface Power (FSP) System. The overall FSP electric power flow is presented, and potential bulk power transmission strategies are discussed. Metallic conductors operate at elevated voltage are identified as the most efficient and feasible solution, and hardware implementation challenges for AC and DC voltage boost are discussed. This work is intended to drive progress on identifying the most feasible power transmission strategy for landed fission power and does not represent consensus or preference by NASA regarding a particular course of action

Lunar↗

Navier-Stokes Solutions of Rotor and Rotor-Body Flows

This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.

Srinivasan, G. R.↗

Squeeze film dampers with oil hole feed

To improve the damping capability of squeeze film dampers, oil hole feed rather than circumferential groove feed is a practical proposition. However, circular orbit response can no longer be assumed, significantly complicating the design analysis. This paper details a feasible transient solution procedure for such dampers, with particular emphasis on the additional difficulties due to the introduction of oil holes. It is shown how a cosine power series solution may be utilized to evaluate the oil hole pressure contributions, enabling appropriate tabular data to be compiled. The solution procedure is shown to be applicable even in the presence of flow restrictors, albeit at the expense of introducing an iteration at each time step. Though not of primary interest, the procedure is also applicable to dynamically loaded journal bearings with oil hole feed.

Chen, P. Y. P.↗

Multiobjective Optimization of Rocket Engine Pumps Using Evolutionary Algorithm

A design optimization method for turbopumps of cryogenic rocket engines has been developed. Multiobjective Evolutionary Algorithm (MOEA) is used for multiobjective pump design optimizations. Performances of design candidates are evaluated by using the meanline pump flow modeling method based on the Euler turbine equation coupled with empirical correlations for rotor efficiency. To demonstrate the feasibility of the present approach, a single stage centrifugal pump design and multistage pump design optimizations are presented. In both cases, the present method obtains very reasonable Pareto-optimal solutions that include some designs outperforming the original design in total head while reducing input power by one percent. Detailed observation of the design results also reveals some important design criteria for turbopumps in cryogenic rocket engines. These results demonstrate the feasibility of the EA-based design optimization method in this field.

Oyama, Akira↗

Description of concept and first feasibility test results of a life support subsystem of the Botany Facility based on water reclamation

The Botany Facility allows the growth of higher plants and fungi over a period of 6 months maximum. It is a payload planned for the second flight of the Eureca platform around 1990. Major tasks of the Life Support Subsystem (LSS) of the Botany Facility include the control of the pressure and composition of the atmosphere within the plant/fungi growth chambers, control of the temperature and humidity of the air and the regulation of the soil water content within specified limits. Previous studies have shown that various LSS concepts are feasible ranging from heavy, simple and cheap to light, complex and expensive solutions. A summary of those concepts is given. A new approach to accomplish control of the temperature and humidity of the air within the growth chambers based on water reclamation is discussed. This reclamation is achieved by condensation with a heat pump and capillary transport of the condensate back into the soil of the individual growth chamber. Some analytical estimates are given in order to obtain guidelines for circulation flow rates and to determine the specific power consumption.

Loeser, H. R.↗

NaK Plugging Meter Design for the Feasibility Test Loops

The design and predicted performance of a plugging meter for use in the measurement of NaK impurity levels are presented. The plugging meter is incorporated into a Feasibility Test Loop (FTL), which is a small pumped-NaK loop designed to enable the rapid, small-scale evaluation of techniques such as in situ purification methods and to permit the measurement of bulk material transport effects (not mechanisms) under flow conditions that are representative of a fission surface power reactor. The FTL operates at temperatures similar to those found in a reactor, with a maximum hot side temperature of 900 K and a corresponding cold side temperature of 860 K. In the plugging meter a low flow rate bypass loop is cooled until various impurities (primarily oxides) precipitate out of solution. The temperatures at which these impurities precipitate are indicative of the level of impurities in the NaK. The precipitates incrementally plug a small orifice in the bypass loop, which is detected by monitoring changes in the liquid metal flow rate.

Pearson, J. Boise↗

The shape of the future long-haul transport airplane

Long-haul transport plane configurations are discussed that are considered technically feasible within a practical level of near-term technology. Supersonic and transonic aircraft are considered, as well as different types of subsonic transports. It is indicated that superior second-generation supersonic transports are quite feasible technically; in particular, it is shown that problems of fuel consumption, efficient overland flight, and questions on noise and nitric oxides have straightforward engineering solutions. Transonic aircraft are practical but would have to have a different geometry than subsonic aircraft to provide the versatility necessary to compete at one end of their capability with supersonic transports and at the other end, for shorter ranges, with subsonic transports. The shape of subsonic transports is not expected to change noticeably, although recent advances of airfoil technology imply that new subsonic transports would have wings with higher aspect ratios and somewhat less sweep. But large specialized freighter aircraft and very long range laminar flow control or nuclear-powered aircraft with radically new shapes are possible.

Clay, C. W.↗

Unsteady flow through compressor stages

The application of the nonsteady Lax-Wendroff technique to problems with asymptotically periodic solution which offers a potentially powerful method for the investigation of the interaction of rotating and stationary blade rows in turbomachinery is reported. A technique for specifying boundary conditions with phase lag was developed to accomplish this. A complete nonlinear analysis is carried out numerically to determine the entire flow field without recourse to the assumption of small disturbances of linear equations which underlie the previous acoustic theories. The result, obtained for the case of equal number of rotor and stator blades shows that transonic flow can be handled without difficulty. In addition, the program is not limited with regard to blade thickness, camber or loading. Extension of this method to incorporate viscous wakes and to analysis of fully three dimensional configuration is feasible, and would greatly expand its utility in practical applications.

Alzner, E.↗

Aqueous, QUick-Charging Battery Integration For Electric Flight Research

AQUIFER establishes technical feasibility of an early-stage technology, a high-energy density, aqueous-based, flow battery, resulting in a near-term increase of 1.7 times range over an all-electric battery, while retiring fire and explosion hazards associated with lithium-based chemistries. The Nano-electrofuel (NEF) flow battery will be integrated with a rim-driven motor (RDM) as a multi-functional design to eliminate conductive EMI and weight from long cable runs, and provide liquid cooling from the aqueous fuel. When successful, the technology provides an improved safety energy storage solution for emission-free electric propulsion in commercial aviation.

power systems↗

A Simple Method for Determining Heat Transfer, Skin Friction, and Boundary-Layer Thickness for Hypersonic Laminar Boundary-Layer Flows in a Pressure Gradient

A procedure based on the method of similar solutions is presented by which the skin friction, heat transfer, and boundary-layer thickness in a laminar hypersonic flow with pressure gradient may be rapidly evaluated if the pressure distribution is known. This solution, which at present is. restricted to power-law variations of pressure with surface distance, is presented for a wide range of exponents in the power law corresponding to both favorable and adverse pressure gradients. This theory has been compared to results from heat-transfer experiments on blunt-nose flat plates and a hemisphere cylinder at free-stream Mach numbers of 4 and 6.8. The flat-plate experiments included tests made at a Mach number of 6.8 over a range of angle of attack of +/- 10 deg. Reasonable agreement of the experimental and theoretical heat-transfer coefficients has been obtained as well as good correlation of the experimental results over the entire range of angle of attack studied. A similar comparison of theory with experiment was not feasible for boundary-layer-thickness data; however, the hypersonic similarity theory was found to account satisfactorily for the variation in boundary-layer thickness due to local pressure distribution for several sets of measurements.

Bertram, Mitchel H.↗

Towards Real-time Global Assessment of Post-fire Debris Flow Hazard

As the risk of wildfires increases worldwide, burned steeplands are vulnerable to the secondary hazard of widespread sediment mobilization through debris flows. Following an initial burn, sediment and soil previously restrained by vegetation are no longer consolidated, allowing for easy mobilization into channels and along steep hillslopes through runoff. Sufficiently powerful rainfall incorporates entrained material into turbulent flows and serves as the primary trigger for debris flow initiation. There is thus an ongoing need to establish the relationship between rainfall and debris flow initiation based on a variety of spatiotemporal preconditions. Previous work establishes regional and local thresholds to constrain the effect of rainfall in recently burned areas, but no empirical or numerical solution has worldwide application. Building from regionally-based efforts in the U.S., this work considers how remote sensing data can be applied to better approximate the post-fire debris flow hazards worldwide using freely available global datasets and software. Our work assesses the utility of remote sensing resources for analyzing burn characteristics, topography, rainfall intensity/duration, and, thus, debris flow initiation. Early results show that global observations are sufficient to delineate background rainfall rates from storms likely to cause debris flows across a variety of burn severity and topographic conditions. However, the dearth of publicly-available post-fire debris flow inventories globally limit the ability to test how the model framework performs within different climatologic and morphologic areas. This work will present preliminary analysis over the Western United States and demonstrate the feasibility of a global, near-real time model to provide situational awareness of potential hazards within recently burned areas worldwide. Future work will also consider how global or regional precipitation forecasts may increase the lead time for improved early warning of these hazards.

Elijah Orland↗

Forced Mixer Nozzle Optimization

Computational fluid dynamic (CFD) and computational acoustic analyses (CAA) were performed for a TFE731-40 compound nozzle, a TFE731-60 mixer nozzle and an Energy Efficient Engine (E(sup 3)) mixer nozzle for comparison with available data. The CFD analyses were performed with a three dimensional, Navier-Stokes solution of the flowfield on an unstructured grid using the RAMPANT program. The CAA analyses were performed with the NASA Glenn MGB program using a structured grid. A successful aerodynamic solution for the TFE731-40 compound nozzle operating statically was obtained, simulating an engine operating on a test stand. Analysis of the CFD results of the TFE731-40 with the MGB program produced predicted sound power levels that agree quite well with the measured data front full-scale static engine tests. Comparison of the predicted sound pressure with the data show good agreement near the jet axis, but the noise levels are overpredicted at angles closer to the inlet. The predicted sound power level for the TFE731-60 did not agree as well with measured static engine data as the TFE731-40. Although a reduction in the predicted noise level due to the mixed flow was observed, the reduction was not as significant as the measured data. The analysis of the V2 mixer from the E(sup 3) study showed that peak temperatures predicted in the mixer exit flowfield were within 5 percent of the values measured by the exit probes. The noise predictions of the V2 mixer nozzle tended to be 3-5 dB higher in peak noise level than the measurements. In addition, the maximum frequency of the noise was also overpredicted. An analysis of the 3 candidate mixer nozzle configurations demonstrated the feasibility of using centerbody lobes and porosity to improve mixing efficiency. A final configuration was designed with a predicted thermal mixing efficiency that was 5 percent higher than the 3 candidate mixers. The results of the MGB noise calculations show that the final design will exceed the design goal of a 3 dB reduction in noise as compared to the baseline TFE731-40.

Sheoran, Yogi↗

Automated Euler and Navier-Stokes Database Generation for a Glide-Back Booster

The past two decades have seen a sustained increase in the use of high fidelity Computational Fluid Dynamics (CFD) in basic research, aircraft design, and the analysis of post-design issues. As the fidelity of a CFD method increases, the number of cases that can be readily and affordably computed greatly diminishes. However, computer speeds now exceed 2 GHz, hundreds of processors are currently available and more affordable, and advances in parallel CFD algorithms scale more readily with large numbers of processors. All of these factors make it feasible to compute thousands of high fidelity cases. However, there still remains the overwhelming task of monitoring the solution process. This paper presents an approach to automate the CFD solution process. A new software tool, AeroDB, is used to compute thousands of Euler and Navier-Stokes solutions for a 2nd generation glide-back booster in one week. The solution process exploits a common job-submission grid environment, the NASA Information Power Grid (IPG), using 13 computers located at 4 different geographical sites. Process automation and web-based access to a MySql database greatly reduces the user workload, removing much of the tedium and tendency for user input errors. The AeroDB framework is shown. The user submits/deletes jobs, monitors AeroDB's progress, and retrieves data and plots via a web portal. Once a job is in the database, a job launcher uses an IPG resource broker to decide which computers are best suited to run the job. Job/code requirements, the number of CPUs free on a remote system, and queue lengths are some of the parameters the broker takes into account. The Globus software provides secure services for user authentication, remote shell execution, and secure file transfers over an open network. AeroDB automatically decides when a job is completed. Currently, the Cart3D unstructured flow solver is used for the Euler equations, and the Overflow structured overset flow solver is used for the Navier-Stokes equations. Other codes can be readily included into the AeroDB framework.

Chaderjian, Neal M.↗

Numerical investigation of internal high-speed viscous flows using a parabolic technique

A feasibility study has been conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves are present. A set of test problems with weak shock and expansion waves have been analyzed with this modified ADD method and stable and accurate solutions were demonstrated provided the streamwise step size was maintained at levels larger than the boundary layer displacement thickness. Calculations made with further reductions in step size encountered departure solutions consistent with strong interaction theory. Calculations were also performed for a flow field with a flame front in which a specific heat release was imposed to simulate a SCRAMJET combustor. In this case the flame front generated relatively thick shear layers which aggravated the departure solution problem. Qualitatively correct results were obtained for these cases using a marching technique with the convective terms in the normal momentum equation suppressed. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

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↗