Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust estimation”

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.

At least 451 records · Page 25

Viral Enumeration Using Cost-Effective Wet-Mount Epifluorescence Microscopy for Aquatic Ecosystems and Modern Microbialites

Enumeration is a fundamental measure of community ecology in which viruses represent the most numerous biological identities. Epifluorescence microscopy (EFM) has been the gold standard method for environmental viral enumeration for over 25 years. Currently, standard EFM methods using the Anodisc filters are no longer cost-effective (>$15 per slide) and have yet to be applied to modern microbialites. Microbialites are microbially driven benthic organosedimentary deposits that have been present for most of Earth’s history. We present a cost-effective method for environ mental viral enumeration from aquatic samples, microbial mats, and exopolymeric substances (EPSs) within modern microbialites using EFM. Our integrated approach, which includes filtration, differential centrifugation, chloroform treatment, glutaralde hyde fixation, benzonase nuclease treatment, probe sonication (EPS and mat only), SYBR Gold staining, wet mounting, and imaging, provides a robust method for modern microbialites and aquatic samples. Viral abundances of modern microbialites and aquatic samples collected from Fayetteville Green Lake (FGL) and Great Salt Lake (GSL) did not differ across ecosystems by sample type. EPS and microbial mat samples had an order of magnitude higher viral-like particle abundance when compared to water regardless of the ecosystem (107 vs 106). Viral enumeration allows for estimates of total viral numbers and weights. The entire weight of all the viruses in FGL and GSL are ~598 g and ~2.2 kg, respectively. Further development of EFM methods and software is needed for viral enumeration. Our method provides a robust and cost-effective (~$0.75 per sample) viral enumeration within modern microbialites and aquatic ecosystems.

Madeline Bellanger↗

Nonlinear control of magnetic bearings

In this paper we present a variety of nonlinear controllers for the magnetic bearing that ensure both stability and robustness. We utilize techniques of discontinuous control to design novel control laws for the magnetic bearing. We present in particular sliding mode controllers, time optimal controllers, winding algorithm based controllers, nested switching controllers, fractional controllers, and synchronous switching controllers for the magnetic bearing. We show existence of solutions to systems governed by discontinuous control laws, and prove stability and robustness of the chosen control laws in a rigorous setting. We design sliding mode observers for the magnetic bearing and prove the convergence of the state estimates to their true values. We present simulation results of the performance of the magnetic bearing subject to the aforementioned control laws, and conclude with comments on design.

Pradeep, A. K.↗

An Efficient and Robust Singular Value Method for Star Pattern Recognition and Attitude Determination

A new star pattern recognition method is developed using singular value decomposition of a measured unit column vector matrix in a measurement frame and the corresponding cataloged vector matrix in a reference frame. It is shown that singular values and right singular vectors are invariant with respect to coordinate transformation and robust under uncertainty. One advantage of singular value comparison is that a pairing process for individual measured and cataloged stars is not necessary, and the attitude estimation and pattern recognition process are not separated. An associated method for mission catalog design is introduced and simulation results are presented.

Juang, Jer-Nan↗

Algorithms for High-Speed Noninvasive Eye-Tracking System

Two image-data-processing algorithms are essential to the successful operation of a system of electronic hardware and software that noninvasively tracks the direction of a person s gaze in real time. The system was described in High-Speed Noninvasive Eye-Tracking System (NPO-30700) NASA Tech Briefs, Vol. 31, No. 8 (August 2007), page 51. To recapitulate from the cited article: Like prior commercial noninvasive eyetracking systems, this system is based on (1) illumination of an eye by a low-power infrared light-emitting diode (LED); (2) acquisition of video images of the pupil, iris, and cornea in the reflected infrared light; (3) digitization of the images; and (4) processing the digital image data to determine the direction of gaze from the centroids of the pupil and cornea in the images. Most of the prior commercial noninvasive eyetracking systems rely on standard video cameras, which operate at frame rates of about 30 Hz. Such systems are limited to slow, full-frame operation. The video camera in the present system includes a charge-coupled-device (CCD) image detector plus electronic circuitry capable of implementing an advanced control scheme that effects readout from a small region of interest (ROI), or subwindow, of the full image. Inasmuch as the image features of interest (the cornea and pupil) typically occupy a small part of the camera frame, this ROI capability can be exploited to determine the direction of gaze at a high frame rate by reading out from the ROI that contains the cornea and pupil (but not from the rest of the image) repeatedly. One of the present algorithms exploits the ROI capability. The algorithm takes horizontal row slices and takes advantage of the symmetry of the pupil and cornea circles and of the gray-scale contrasts of the pupil and cornea with respect to other parts of the eye. The algorithm determines which horizontal image slices contain the pupil and cornea, and, on each valid slice, the end coordinates of the pupil and cornea. Information from multiple slices is then combined to robustly locate the centroids of the pupil and cornea images. The other of the two present algorithms is a modified version of an older algorithm for estimating the direction of gaze from the centroids of the pupil and cornea. The modification lies in the use of the coordinates of the centroids, rather than differences between the coordinates of the centroids, in a gaze-mapping equation. The equation locates a gaze point, defined as the intersection of the gaze axis with a surface of interest, which is typically a computer display screen (see figure). The expected advantage of the modification is to make the gaze computation less dependent on some simplifying assumptions that are sometimes not accurate

Talukder, Ashit↗

Deriving Mercury Geodetic Parameters with Altimetric Crossovers from the Mercury Laser Altimeter (MLA)

Based on the previous applications of laser altimetry to planetary geodesy at GSFC, we use the recently developed PyXover software package to analyze altimetric crossovers from the Mercury Laser Altimeter (MLA). Using PyXover, we place new constraints on Mercury's geodetic parameters via least squares minimization of crossover discrepancies. We simultaneously solve for orbital corrections for each MLA ground track, for the geodetic parameters of the International Astronomical Union-recommended orientation model for Mercury (pole right-ascension and declination coordinates, prime meridian rotation rate, and librations), and for the Mercury's Love number h2. We calibrate the formal errors of our solution based on closed-loop simulations and on the level of robustness against a priori values, data selection, and parametrization. Our solution of the Mercury's rotational parameters is consistent with published values. In particular, our new estimate for the orientation of the pole places Mercury in a Cassini state, with an obliquity ϵ = 2.031 ± 0.03 arcmin compatible with previous “surface” related measurements. Moreover, we provide a first data-based estimate of the Love number h_(2) = 1.55 ± 0.65. The latter is consistent with expectations from models of Mercury's interior, although its precision does not enable their refinement.

S. Bertone↗

Using Space Weather Variability in Evaluating the Radiation Environment Design Specifications for NASA's Constellation Program

NASA's Constellation program, initiated to fulfill the Vision for Space Exploration, will create a new generation of vehicles for servicing low Earth orbit, the Moon, and beyond. Space radiation specifications for space system hardware are necessarily conservative to assure system robustness for a wide range of space environments. Spectral models of solar particle events and trapped radiation belt environments are used to develop the design requirements for estimating total ionizing radiation dose, displacement damage, and single event effects for Constellation hardware. We first describe the rationale using the spectra chosen to establish the total dose and single event design environmental specifications for Constellation systems. We then compare variability of the space environment to the spectral design models to evaluate their applicability as conservative design environments and potential vulnerabilities to extreme space weather events

Coffey, Victoria N.↗

Ozone Production from the 2004 North American Boreal Fires

We examine the ozone production from boreal forest fires based on a case study of wildfires in Alaska and Canada in summer 2004. The model simulations were performed with the chemistry transport model, MOZART-4, and were evaluated by comparison with a comprehensive set of aircraft measurements. In the analysis we use measurements and model simulations of carbon monoxide (CO) and ozone (O3) at the PICO-NARE station located in the Azores within the pathway of North American outflow. The modeled mixing ratios were used to test the robustness of the enhancement ratio deltaO3/deltaCO (defined as the excess O3 mixing ratio normalized by the increase in CO) and the feasibility for using this ratio in estimating the O3 production from the wildfires. Modeled and observed enhancement ratios are about 0.25 ppbv/ppbv which is in the range of values found in the literature, and results in a global net O3 production of 12.9 2 Tg O3 during summer 2004. This matches the net O3 production calculated in the model for a region extending from Alaska to the East Atlantic (9-11 Tg O3) indicating that observations at PICO-NARE representing photochemically well-aged plumes provide a good measure of the O3 production of North American boreal fires. However, net chemical loss of fire related O3 dominates in regions far downwind from the fires (e.g. Europe and Asia) resulting in a global net O3 production of 6 Tg O3 during the same time period. On average, the fires increased the O3 burden (surface-300 mbar) over Alaska and Canada during summer 2004 by about 7-9%, and over Europe by about 2-3%.

Pfister, G. G.↗

Localization of Ad-Hoc Lunar Constellations in Communication Failure Modes for Distributed Spacecraft Autonomy

As Lunar missions increase in complexity, inspired by NASA’s Artemis Program, they will require reliable and sufficient Position, Navigation, and Timing (PNT) capability to support the upcoming Lunar users. The navigation service should also be compatible with the smaller platforms, like CubeSats, being sent by the public and private sectors. A non-dedicated, ad-hoc Lunar navigation constellation can provide PNT services on-demand using the non-dedicated swarm assets. Swarm members cooperatively and autonomously localize themselves with minimal interaction from Earth, freeing up valuable bandwidth and ground segment resources. The autonomous localization of Lunar constellations utilizes neighbor two-way intersatellite link (ISL) measurements in a distributed extended Kalman filter (DEKF) system to minimize operating costs. Because the decentralized Lunar PNT system relies on relay communication amongst the agents, network failures or loss of assets among ad-hoc Lunar constellations may impact localization performance. This study presents an evaluation of localization performance under increasing levels of network degradation. A simulation of an ad-hoc Lunar PNT swarm is augmented to include system faults and the impacts of intermittent and permanent failures on localization performance are evaluated. We investigate three potential causes of network degradation: single spacecraft loss, multiple spacecraft loss, and antenna failure. The numerical assessments from the simulation show that the LPNT system under study, based on an autonomous decentralized concept of operation, is highly robust and resilient to communication failures. Minor faults, such as single spacecraft loss, solar interference, technical malfunctions, message delays, and antenna outages, have minimal impact on state estimation, with only a 4.47% and 3.75% degradation in median position error for assets and a representative ground user, respectively, compared to an ideal communication scenario. However, major faults, such as hardware failures or meteor strikes leading to the loss of multiple spacecrafts, are more concerning. The permanent loss of three spacecraft results in a more severe performance degradation, with median position error increasing by 23.3% for assets and 11.7% for a representative ground user, despite the Lunar PNT system remaining functional.

Yeji Kim↗

Interaction Between Eddies and Mean Flow in Jupiter's Atmosphere: Analysis of Cassini Imaging Data

Beebe et al. [Beebe, R.F., et al., 1980. Geophys. Res. Lett. 17, 1-4] and Ingersoll et al. [Ingersoll, A.P., et al., 1981. J. Geophys. Res. 86, 8733-8743] used images from Voyagers 1 and 2 to analyze the interaction between zonal winds and eddies in Jupiter's atmosphere. They reported a high positive correlation between Jupiter's eddy momentum flux, pu'v', and the variation of zonal velocity with latitude, du/dy. This correlation implied a surprisingly high rate of conversion of energy from eddies to zonal flow: approx. 1.5-3.0 W/sq m, a value more than 10% of Jupiter s thermal flux emission. However, Sromovsky et al. [Sromovsky, L.A., et al., 1982. J. Atmos. Sci. 39,1413-1432] argued that possible biases in the analysis could have caused an artificially high correlation. In addition, significant differences in the derived eddy flux between datasets put into question the robustness of any one result. We return to this long-standing puzzle using images of Jupiter from the Cassini flyby of December 2000. Our method is similar to previous analyses, but utilizes an automatic feature tracker instead of the human eye. The number of velocity vectors used in this analysis is over 200,000, compared to the 14,000 vectors used by Ingersoll et al. We also find a positive correlation between u'v' and du/dy and derive a global average power per unit mass, u'v' du/dy, ranging from (7.1-12.3) x 10(exp -5)W/kg. Utilizing Ingersoll et al.'s estimate of the mass per unit area involved in the transport, this would imply a rate of energy conversion of approx.0.7-1.2 W/sq m. We discuss the implications of this result and employ several tests to demonstrate its robustness.

Salyk, Colette↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗

Dissipative compensators for flexible spacecraft control

The problem of controller design for flexible spacecraft is addressed. Model-based compensators, which rely on the knowledge of the system parameters to tune the state estimator, are considered. The instability mechanisms resulting from high sensitivity to parameter uncertainties are investigated. Dissipative controllers, which use collocated actuators and sensors, are also considered, and the robustness properties of constant-gain dissipative controllers in the presence of unmodeled elastic-mode dynamics, sensor/actuator nonlinearities, and actuator dynamics are summarized. In order to improve the performance without sacrificing robustness, a class of dissipative dynamic compensators is proposed and is shown to retain robust stability in the presence of second-order actuator dynamics if acceleration feedback is employed. A class of dissipative dynamic controllers is proposed which consists of a low-authority, constant-gain controller and a high-authority dynamic compensator. A procedure for designing an optimal dissipative dynamic compensator is given which minimizes a quadratic performance criterion. Such compensators offer the promise of better performance while still retaining robust stability.

Joshi, S. M.↗

A nonparametric software-reliability growth model

The authors (1985) previously introduced a nonparametric model for software-reliability growth which is based on complete monotonicity of the failure rate. The authors extend the completely monotone software model by developing a method for providing long-range predictions of reliability growth, based on the model. They derive upper and lower bounds on extrapolation of the failure rate and the mean function. These are then used to obtain estimates for the future software failure rate and the mean future number of failures. Preliminary evaluation indicates that the method is competitive with parametric approaches, while being more robust.

Sofer, Ariela↗

Plasma-Based Detector of Outer-Space Dust Particles

A report presents a concept for an instrument to be flown in outer space, where it would detect dust particles - especially those associated with comets. The instrument would include a flat plate that would intercept the dust particles. The anticipated spacecraft/dust-particle relative speeds are so high that the impingement of a dust particle on the plate would generate a plasma cloud. Simple electric dipole sensors located equidistantly along the circumference of the plate would detect the dust particle indirectly by detecting the plasma cloud. The location of the dust hit could be estimated from the timing of the detection pulses of the different dipoles. The mass and composition of the dust particle could be estimated from the shapes and durations of the pulses from the dipoles. In comparison with other instruments for detecting hypervelocity dust particles, the proposed instrument offers advantages of robustness, large collection area, and simplicity.

Tsurutani, Bruce↗

Piloted Simulator Evaluation of Maneuvering Envelope Information for Flight Crew Awareness

The implementation and evaluation of an efficient method for estimating safe aircraft maneuvering envelopes are discussed. A Bayesian approach is used to produce a deterministic algorithm for estimating aerodynamic system parameters from existing noisy sensor measurements, which are then used to estimate the trim envelope through efficient high- fidelity model-based computations of attainable equilibrium sets. The safe maneuverability limitations are extended beyond the trim envelope through a robust reachability analysis derived from an optimal control formulation. The trim and maneuvering envelope limits are then conveyed to pilots through three axes on the primary flight display. To evaluate the new display features, commercial airline crews flew multiple challenging approach and landing scenarios in the full motion Advanced Concepts Flight Simulator at NASA Ames Research Center, as part of a larger research initiative to investigate the impact on the energy state awareness of the crew. Results show that the additional display features have the potential to significantly improve situational awareness of the flight crew.

maneuvers↗

Advanced Background Subtraction Applied to Aeroacoustic Wind Tunnel Testing

An advanced form of background subtraction is presented and applied to aeroacoustic wind tunnel data. A variant of this method has seen use in other fields such as climatology and medical imaging. The technique, based on an eigenvalue decomposition of the background noise cross-spectral matrix, is robust against situations where isolated background auto-spectral levels are measured to be higher than levels of combined source and background signals. It also provides an alternate estimate of the cross-spectrum, which previously might have poor definition for low signal-to-noise ratio measurements. Simulated results indicate similar performance to conventional background subtraction when the subtracted spectra are weaker than the true contaminating background levels. Superior performance is observed when the subtracted spectra are stronger than the true contaminating background levels. Experimental results show limited success in recovering signal behavior for data where conventional background subtraction fails. They also demonstrate the new subtraction technique's ability to maintain a proper coherence relationship in the modified cross-spectral matrix. Beam-forming and de-convolution results indicate the method can successfully separate sources. Results also show a reduced need for the use of diagonal removal in phased array processing, at least for the limited data sets considered.

Bahr, Christopher J.↗

150 Shades of Green: Using the Full Spectrum of Remote Sensing Reflectance to Elucidate Color Shifts in the Ocean

This article proposes a simple and intuitive classification system by which to define full spectral remote sensing reflectance (Rrs(λ)) data with a quantitative output that enables a more manageable handling of spectral information for aquatic science applications. The weighted harmonic mean of the Rrs(λ) wavelengths outputs an Apparent Visible Wavelength (in units of nanometers), representing a one-dimensional geophysical metric of color that is inherently correlated to spectral shape. This dimensionality reduction of spectral information combined with the output along a continuum of wavelength values offers a robust and user-friendly means to describe and analyze spectral Rrs(λ) in terms of spatial and temporal trends and variability. The uncertainty in the algorithm's estimation of spectral shape is demonstrated on a global scale, in addition to the utility of the algorithm to discern spectral-spatial-temporal trends in the ocean, on a per-pixel basis for the entire 22 year continuous ocean color (SeaWiFS and MODIS-Aqua) time-series. This technique can be applied to datasets of varying multi- and hyper-spectral resolutions, providing continuity between heritage and future satellite sensors, and further enabling an effective means of elucidating similarities or differences in complex spectral signatures within the constraints of two dimensions. This straightforward means of conceptualizing multi-dimensional variability can help maximize the potential of the spectral information embedded in remote sensing data.

ocean color↗

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↗

Science Yield Estimates and Sensitivities With the AstroPIC Integrated Photonic Coronagraph for the Habitable Worlds Observatory

The Habitable Worlds Observatory (HWO) flagship has a goal to survey ~100 of the nearest stellar systems and their habitable zones in order to detector and spectroscopically characterize ~25 potentially “Earth-like planets” (or “Exoearths”). The HWO telescope will feature a >6-m diameter aperture and a coronagraph that will suppress starlight by ten orders of magnitude to enable the detection and characterization of Earth analogues in the habitable zone. Three Exploratory Analytic Cases (EACs) are currently considering different aperture configurations for HWO including on-axis and off-axis apertures and both different hexagonal and keystone segmentation configurations. Photonic Integrated Circuits (PICs) offer a path to miniaturize traditional, bulk optics coronagraph functionality reducing risk margins and increasing mission science exoearth yields for the HWO mission. We have developed a simulation pipeline that enables evaluation of the AstroPIC photonic integrated coronagraph through to estimates mission yields. We present sensitivities to aberrations for different types of segmentation and low-order aberrations for HWO EACs demonstrating that AstroPIC is suitable for both on-axis and off-axis configurations. We perform a parametric sweep of the exoearth science yield sensitivities in terms of key coronagraph parameters including coronagraphic end-to-end throughput, inner working angle, robustness to aberrations, and bandwidth. We consider both a hybrid and a pure photonic architecture; a pure integrated photonic architecture provides the greatest potential in science yield improvements, although a hybrid PIC limited in its numbers of modes can provide complementary science to a traditional, bulk-optic coronagraph and increasing the estimated science yield when operated jointly.

Photonics↗