Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched bound”

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 91 records · Page 5

Sensitivity of Simulated Radiative Emission from Nitrogen Flows to Chemical Parameters

Numerical models of hypersonic flows rely on chemistry data inferred from experiments and quantum calculations. Experiments from NASA's Electric Arc Shock Tube (EAST) measure the spectral emission of strong shocks in pure nitrogen flows. These spectral profiles are only approximately captured by numerical simulations. To address these discrepancies, one wishes to calibrate some of the chemical parameters to better match experiment. Because there are hundreds of parameters in the numerical model, a sensitivity analysis was carried out on the entire chemistry database to determine which parameters are most important to analyze in the future. In order to perform this work it was necessary to generate a consistent chemical database that interfaced with both the CFD and line-by-line radiation solver. Using Monte Carlo methods, the total Sobol index was calculated for important parameters. Twenty parameters were identified as important across the nitrogen system. However, it was noted that in some spectral regions, the parametric uncertainty was unable to bound EAST measurements. Further improvements could be made by increasing the parameter uncertainties or addressing other sources of modeling error.

Kaelan B Hansson↗

Sensitivity of Simulated Radiative Emission from Nitrogen Flows to Chemical Parameters

Numerical models of hypersonic flows rely on chemistry data inferred from experiments and quantum calculations. Experiments from NASA's Electric Arc Shock Tube (EAST) measure the spectral emission of strong shocks in pure nitrogen flows. These spectral profiles are only approximately captured by numerical simulations. To address these discrepancies, one wishes to calibrate some of the chemical parameters to better match experiment. Because there are hundreds of parameters in the numerical model, a sensitivity analysis was carried out on the entire chemistry database to determine which parameters are most important to analyze in the future. In order to perform this work it was necessary to generate a consistent chemical database that interfaced with both the CFD and line-by-line radiation solver. Using Monte Carlo methods, the total Sobol index was calculated for important parameters. Twenty parameters were identified as important across the nitrogen system. However, it was noted that in some spectral regions, the parametric uncertainty was unable to bound EAST measurements. Further improvements could be made by increasing the parameter uncertainties or addressing other sources of modeling error.

Kaelan Hansson↗

Orbital Trade Study for the PREFIRE Mission

The Polar Radiant Energy in the Far-Infrared Experiment (PREFIRE) is a low-cost cubesat-based Earth Science mission selected for rapid implementation through NASAs Earth Ventures-Instrument program. Two cubesats carrying thermal infrared spectrometers (TIRS) in polar orbits will test long-held assumptions regarding radiative balance at long wavelengths where the cold-dry Arctic and Antarctic surfaces expel significant amounts of the earth's radiant energy. Thermal radiometry with better than one micron sampling at wavelengths beyond 15 microns will provide critical observational bounds for inputs into glacier melt and climate change models through improvements in surface emissivities, longwave greenhouse effect parameterizations, and cloud presence/types. Additional information on potential far-infrared surface/atmosphere feedbacks during rapid melt processes may be captured with sub-diurnal revisits of targeted areas. A baseline mission with two cubesats in different orbits provides a means for optimization of sub-diurnal sampling for science purposes. In this presentation we provide details of orbit studies that inform expected sampling metrics. Candidate orbits are modeled in SGP and then TIRS sampling is projected into surface footprints with appropriate geo-location. Matching algorithms then compile statistics for time-delayed revisits between and within the two orbits. Binning of results by latitude allows for abstraction of precession issues as well as for correlations with surface types as defined by radiometric classification schemes. Multiple acceptable pairs of science-driven orbits were found that optimize sampling in areas known to experience rapid melt events. The final choice of PREFIRE orbital parameters includes the systems trade-space, particularly the power budget, as well as the availability of launch opportunities.

"L'Ecuyer, Tristan"↗

Equations of state and impact-induced shock-wave attenuation on the moon

Equation of state formulations are considered in a framework that permits comparison with one-dimensional impedance match solutions. The problem considered is the peak pressures attained along the impact symmetry axis when a sphere impacts with a half-space. The regimes of melting and vaporization - in particular, incipiently melted, completely melted, incipiently vaporized, and completely vaporized states - are examined, and the pressures at which critical isentropes intersect the Hugoniots of iron and gabbroic anorthosite are considered. A means of representing the spatial attenuation of shock pressure along the impact axis by two regimes is introduced, and results for the near-field and far-field regime are presented. It is thought that the treatment can be used to obtain quantitative bounds on the impact velocity of the meteorite.

Ahrens, T. J.↗

Spatial learning and memory is preserved in rats after early development in a microgravity environment

This study evaluated the cognitive mapping abilities of rats that spent part of their early development in a microgravity environment. Litters of male and female Sprague-Dawley rat pups were launched into space aboard the National Aeronautics and Space Administration space shuttle Columbia on postnatal day 8 or 14 and remained in space for 16 days. These animals were designated as FLT groups. Two age-matched control groups remained on Earth: those in standard vivarium housing (VIV) and those in housing identical to that aboard the shuttle (AGC). On return to Earth, animals were tested in three different tasks that measure spatial learning ability, the Morris water maze (MWM), and a modified version of the radial arm maze (RAM). Animals were also tested in an open field apparatus to measure general activity and exploratory activity. Performance and search strategies were evaluated in each of these tasks using an automated tracking system. Despite the dramatic differences in early experience, there were remarkably few differences between the FLT groups and their Earth-bound controls in these tasks. FLT animals learned the MWM and RAM as quickly as did controls. Evaluation of search patterns suggested subtle differences in patterns of exploration and in the strategies used to solve the tasks during the first few days of testing, but these differences normalized rapidly. Together, these data suggest that development in an environment without gravity has minimal long-term impact on spatial learning and memory abilities. Any differences due to development in microgravity are quickly reversed after return to earth normal gravity.

NASA Discipline Neuroscience↗

Dense gas in high-latitude molecular clouds

The nearby molecular clouds MBM 7, 12, 30, 32, 40, 41, and 55 were surveyed for tracers of dense gas, including the (1-0), (2-1), and (3-2) rotational lines of CS and the (1-0) lines of HCO(+) and HCN. MBM 7 and MBM 12 contain dense cores, while the other clouds contain little or no traces of dense gas. Comparison of the emission from dense gas tracers to that of (13)CO reveals that the former are more compact in angular size as well as line width. An extensive CS(2-1) survey of part of MBM 12 reveals that the emission is characterized by clumps on approximately 3 min scales as well as extended emission. Observations of the CS(1-0) and (3-2) lines using telescopes with matched beam sizes reveal that the volume density must be at least approximately 10(exp 4.5)/cc within the (3-2) emitting regions, which are approximately 0.03 pc in radius. Electron excitation of the CS rotational levels is ruled out (in the cores) by comparing the (3-2)/(1-0) line ratios with models including H2 and electron collisions. The volume density in the cores is substantially larger than in the portions of the cloud traced by CO emission. The density increases into the cores as r(exp -2), suggesting dynamical collapse. The masses of the cores are close to the virial mass, suggesting they are dynamically bound. The cores in MBM 7 and MBM 12 are thus likely to form stars; they are the nearest sites of star formation.

Reach, William T.↗

RCC Plug Repair Thermal Tools for Shuttle Mission Support

A thermal math model for the Space Shuttle Reinforced Carbon-Carbon (RCC) Plug Repair was developed to increase the confidence in the repair entry performance and provide a real-time mission support tool. The thermal response of the plug cover plate, local RCC, and metallic attach hardware can be assessed with this model for any location on the wing leading edge. The geometry and spatial location of the thermal mesh also matches the structural mesh which allows for the direct mapping of temperature loads and computation of the thermoelastic stresses. The thermal model was correlated to a full scale plug repair radiant test. To utilize the thermal model for flight analyses, accurate predictions of protuberance heating were required. Wind tunnel testing was performed at CUBRC to characterize the heat flux in both the radial and angular directions. Due to the complexity of the implementation of the protuberance heating, an intermediate program was developed to output the heating per nodal location for all OML surfaces in SINDA format. Three Design Reference Cases (DRC) were evaluated with the correlated plug thermal math model to bound the environments which the plug repair would potentially be used.

Rodriguez, Alvaro C.↗

JWST ISIM Distortion Analysis Challenge

Very tight distortion requirements are imposed on the JWST's ISM structure due to the sensitivity of the telescope's mirror segment and science instrument positioning. The ISIM structure is a three dimensional truss with asymmetric gusseting and metal fittings. One of the primary challenges for ISIM's analysis team is predicting the thermal distortion of the structure both from the bulk cooldown from ambient to cryo, and the smaller temperature changes within the cryogenic operating environment. As a first cut to estimate thermal distortions, a finite element model of bar elements was created. Elements representing joint areas and metal fittings use effective properties that match the behavior of the stack-up of the composite tube, gusset and adhesive under mechanical and thermal loads. These properties were derived by matching tip deflections of a solid model simplified T-joint. Because of the structure s asymmetric gusseting, this effective property model is a first attempt at predicting rotations that cannot be captured with a smeared CTE approach. In addition to the finite element analysis, several first order calculations have been performed to gauge the feasibility of the material design. Because of the stringent thermal distortion requirements at cryogenic temperatures, a composite tube material with near zero or negative CTE is required. A preliminary hand analysis of the contribution of the various components along the distortion path between FGS and the other instruments, neglecting second order effects were examined. A plot of bounding tube longitudinal and transverse CTEs for thermal stability requirements was generated to help determine the feasibility of meeting these requirements. This analysis is a work in progress en route to a large degree of freedom hi-fidelity FEA model for distortion analysis. Methods of model reduction, such as superelements, are currently being investigated.

Cifie, Emmanuel↗

On Target to Mars

This viewgraph presentation reviews the use of Descent Image Motion Estimation System (DIMES) for the descent of a spacecraft onto the surface of Mars. In the past this system was used to assist in the landing of the MER spacecraft. The overall algorithm is reviewed, and views of the hardware, and views from Spirit's descent are shown. On Spirit, had DIMES not been used, the impact velocity would have been at the limit of the airbag capability and Spirit may have bounced into Endurance Crater. By using DIMES, the velocity was reduced to well within the bounds of the airbag performance and Spirit arrived safely at Mars. Views from Oppurtunity's descent are also shown. The system to avoid and detect hazards is reviewed next. Landmark Based Spacecraft Pinpoint Landing is also reviewed. A cartoon version of a pinpoint landing and the various points is shown. Mars s surface has a large amount of craters, which are ideal landmarks . According to literatures on Martian cratering, 60 % of Martian surface is heavily cratered. The ideal (craters) landmarks for pinpoint landing will be between 1000 to 50 meters in diagonal The ideal altitude for position estimation should greater than 2 km above the ground. The algorithms used to detect and match craters are reviewed.

Mars↗

Frequency-Modulated, Continuous-Wave Laser Ranging Using Photon-Counting Detectors

Optical ranging is a problem of estimating the round-trip flight time of a phase- or amplitude-modulated optical beam that reflects off of a target. Frequency- modulated, continuous-wave (FMCW) ranging systems obtain this estimate by performing an interferometric measurement between a local frequency- modulated laser beam and a delayed copy returning from the target. The range estimate is formed by mixing the target-return field with the local reference field on a beamsplitter and detecting the resultant beat modulation. In conventional FMCW ranging, the source modulation is linear in instantaneous frequency, the reference-arm field has many more photons than the target-return field, and the time-of-flight estimate is generated by balanced difference- detection of the beamsplitter output, followed by a frequency-domain peak search. This work focused on determining the maximum-likelihood (ML) estimation algorithm when continuous-time photoncounting detectors are used. It is founded on a rigorous statistical characterization of the (random) photoelectron emission times as a function of the incident optical field, including the deleterious effects caused by dark current and dead time. These statistics enable derivation of the Cramér-Rao lower bound (CRB) on the accuracy of FMCW ranging, and derivation of the ML estimator, whose performance approaches this bound at high photon flux. The estimation algorithm was developed, and its optimality properties were shown in simulation. Experimental data show that it performs better than the conventional estimation algorithms used. The demonstrated improvement is a factor of 1.414 over frequency-domainbased estimation. If the target interrogating photons and the local reference field photons are costed equally, the optimal allocation of photons between these two arms is to have them equally distributed. This is different than the state of the art, in which the local field is stronger than the target return. The optimal processing of the photocurrent processes at the outputs of the two detectors is to perform log-matched filtering followed by a summation and peak detection. This implies that neither difference detection, nor Fourier-domain peak detection, which are the staples of the state-of-the-art systems, is optimal when a weak local oscillator is employed.

Erkmen, Baris I.↗

Kriging as a Means of Improving WAAS Availability

The Wide Area Augmentation System (WAAS), an augmentation of the Global Positioning System (GPS), provides safe and reliable use of GPS signals for airline navigation over much of North America. Currently the largest source of positioning error in the system is signal delay caused by the ionosphere. To allow the user to take account of such error, WAAS computes and broadcasts ionospheric vertical delays at a set of regularly-spaced grid points. In addition, WAAS computes and broadcasts a safety-critical integrity bound at each ionospheric grid point (IGP) called the Grid Ionospheric Vertical Error (GIVE). GIVEs are constructed to be sufficiently large to protect the user against positioning error due to the presence of ionospheric irregularity. In the initial operating capability (IOC) of WAAS, the vertical delay estimate at each IGP is determined from a planar fit of neighboring slant delay measurements, projected to vertical using an obliquity factor specified by the standard thin-shell model of the ionosphere. In WAAS Follow-On (WFO) Release 3, however, the vertical delay will be estimated by an established, geo-statistical technique known as kriging. Compared to the planar fit model, the kriging model is found, in general, to match better the observed random structure of the vertical delay. This paper presents the kriging methodology that will be used to estimate the vertical delay and its uncertainty at each IGP, and it assesses the subsequent improvement in WAAS availability enabled by kriging.

GPS↗

Data compression using adaptive transform coding. Appendix 1: Item 1

Adaptive low-rate source coders are described in this dissertation. These coders adapt by adjusting the complexity of the coder to match the local coding difficulty of the image. This is accomplished by using a threshold driven maximum distortion criterion to select the specific coder used. The different coders are built using variable blocksized transform techniques, and the threshold criterion selects small transform blocks to code the more difficult regions and larger blocks to code the less complex regions. A theoretical framework is constructed from which the study of these coders can be explored. An algorithm for selecting the optimal bit allocation for the quantization of transform coefficients is developed. The bit allocation algorithm is more fully developed, and can be used to achieve more accurate bit assignments than the algorithms currently used in the literature. Some upper and lower bounds for the bit-allocation distortion-rate function are developed. An obtainable distortion-rate function is developed for a particular scalar quantizer mixing method that can be used to code transform coefficients at any rate.

Rost, Martin Christopher↗

Atmospheric Energy Deposition Modeling and Inference for Varied Meteoroid Structures

Asteroids populations are highly diverse, ranging from coherent monoliths to loosely-bound rubble piles with a broad range of material and compositional properties. These different structures and properties could significantly affect how an asteroid breaks up and deposits energy in the atmosphere, and how much ground damage may occur from resulting blast waves. We have previously developed a fragment-cloud model (FCM) for assessing the atmospheric breakup and energy deposition of asteroids striking Earth. The approach represents ranges of breakup characteristics by combining progressive fragmentation with releases of variable fractions of debris and larger discrete fragments. In this work, we have extended the FCM to also represent asteroids with varied initial structures, such as rubble piles or fractured bodies. We have used the extended FCM to model the Chelyabinsk, Benesov, Kosice, and Tagish Lake meteors, and have obtained excellent matches to energy deposition profiles derived from their light curves. These matches provide validation for the FCM approach, help guide further model refinements, and enable inferences about pre-entry structure and breakup behavior. Results highlight differences in the amount of small debris vs. discrete fragments in matching the various flare characteristics of each meteor. The Chelyabinsk flares were best represented using relatively high debris fractions, while Kosice and Benesov cases were more notably driven by their discrete fragmentation characteristics, perhaps indicating more cohesive initial structures. Tagish Lake exhibited a combination of these characteristics, with lower-debris fragmentation at high altitudes followed by sudden disintegration into small debris in the lower flares. Results from all cases also suggest that lower ablation coefficients and debris spread rates may be more appropriate for the way in which debris clouds are represented in FCM, offering an avenue for future model refinement.

Deposition↗

A Resonator for Low-Threshold Frequency Conversion

A proposed toroidal or disklike dielectric optical resonator (dielectric optical cavity) would be made of an optically nonlinear material and would be optimized for use in parametric frequency conversion by imposition of a spatially periodic permanent electric polarization. The poling (see figure) would suppress dispersions caused by both the material and the geometry of the optical cavity, thereby effecting quasi-matching of the phases of high-resonance-quality (high-Q) whispering-gallery electromagnetic modes. The quasi-phase-matching of the modes would serve to maximize the interactions among them. Such a resonator might be a prototype of a family of compact, efficient nonlinear devices for operation over a broad range of optical wavelengths. A little background information is prerequisite to a meaningful description of this proposal: (1) Described in several prior NASA Tech Briefs articles, the whispering-gallery modes in a component of spheroidal, disklike, or toroidal shape are waveguide modes that propagate circumferentially and are concentrated in a narrow toroidal region centered on the equatorial plane and located near the outermost edge. (2) For the sake of completeness, it must be stated that even though optical resonators of the type considered here are solid dielectric objects and light is confined within them by total internal reflection at dielectric interfaces without need for mirrors, such components are sometimes traditionally called cavities because their effects upon the light propagating within them are similar to those of true cavities bounded by mirrors. (3) For a given set of electromagnetic modes interacting with each other in an optically nonlinear material (e.g., modes associated with the frequencies involved in a frequency-conversion scheme), the threshold power for oscillation depends on the mode volumes and the mode-overlap integral. (4) Whispering-gallery modes are attractive in nonlinear optics because they maximize the effects of nonlinearities by occupying small volumes and affording high Q values

Iltchenko, Vladimir↗

Venusian Applications of 3D Convection Modeling

This study models mantle convection on Venus using the 'cubed sphere' code OEDIPUS, which models one-sixth of the planet in spherical geometry. We are attempting to balance internal heating, bottom mantle viscosity, and temperature difference across Venus' mantle, in order to create a realistic model that matches with current planetary observations. We also have begun to run both lower and upper mantle simulations to determine whether layered (as opposed to whole-mantle) convection might produce more efficient heat transfer, as well as to model coronae formation in the upper mantle. Upper mantle simulations are completed using OEDIPUS' Cartesian counterpart, JOCASTA. This summer's central question has been how to define a mantle plume. Traditionally, we have defined a hot plume the region with temperature at or above 40% of the difference between the maximum and horizontally averaged temperature, and a cold plume as the region with 40% of the difference between the minimum and average temperature. For less viscous cases (1020 Pa?s), the plumes generated by that definition lacked vigor, displaying buoyancies 1/100th of those found in previous, higher viscosity simulations (1021 Pa?s). As the mantle plumes with large buoyancy flux are most likely to produce topographic uplift and volcanism, the low viscosity cases' plumes may not produce observable deformation. In an effort to eliminate the smallest plumes, we experimented with different lower bound parameters and temperature percentages.

Venus↗

Flow field measurements in the cell culture unit

The cell culture unit (CCU) is being designed to support cell growth for long-duration life science experiments on the International Space Station (ISS). The CCU is a perfused loop system that provides a fluid environment for controlled cell growth experiments within cell specimen chambers (CSCs), and is intended to accommodate diverse cell specimen types. Many of the functional requirements depend on the fluid flow field within the CSC (e.g., feeding and gas management). A design goal of the CCU is to match, within experimental limits, all environmental conditions, other than the effects of gravity on the cells, whether the hardware is in microgravity ( micro g), normal Earth gravity, or up to 2g on the ISS centrifuge. In order to achieve this goal, two steps are being taken. The first step is to characterize the environmental conditions of current 1g cell biology experiments being performed in laboratories using ground-based hardware. The second step is to ensure that the design of the CCU allows the fluid flow conditions found in 1g to be replicated from microgravity up to 2g. The techniques that are being used to take these steps include flow visualization, particle image velocimetry (PIV), and computational fluid dynamics (CFD). Flow visualization using the injection of dye has been used to gain a global perspective of the characteristics of the CSC flow field. To characterize laboratory cell culture conditions, PIV is being used to determine the flow field parameters of cell suspension cultures grown in Erlenmeyer flasks on orbital shakers. These measured parameters will be compared to PIV measurements in the CSCs to ensure that the flow field that cells encounter in CSCs is within the bounds determined for typical laboratory experiments. Using CFD, a detailed simulation is being developed to predict the flow field within the CSC for a wide variety of flow conditions, including microgravity environments. Results from all these measurements and analyses of the CSC flow environment are presented and discussed. The final configuration of the CSC employs magnetic stir bars with angled paddles to achieve the necessary flow requirements within the CSC.

Cell Culture↗

The Balmer-dominated northeast limb of the Cygnus loop supernova remnant

We present a comprehensive investigation of the Balmar-dominated northeast limb of the Cygnus Loop supernova remnant. Data presented include H alpha (O III), and X-ray images, UV and visible spectrophotometry, and high-resolution spectroscophy. The two relatively bright Balmer-dominated filaments visible on the POSS prints are seen to be part of a very smooth and regular complex of filaments. These filaments mark the current location of the blast wave and are seen to bound the sharply limb-brightened X-ray emission, including the previously reported X-ray, 'halo.' The (O III)/h beta ratio throughout the region is approximately 0.1, except for regions in which the shock is undergoing a transition from nonradiative to incomplete radiative to incomplete radiative conditions. At these locations (O III) emission from the cooling region is quite strong, while collisionally excited Balmer-line emission can be weak because of photoionization of the preshock medium by UV from the nascent cooling region. As a result (O III)/H beta is greater than 100 in some locations. The nonradiative/radiative transition is best studied along the length of the northwestern of the two brightest filaments, where the shock velocity and swept-up column go from approximately 180 km/s and 10(exp 17)/sq cm at one end to approximately 140 km/s and 8 x 10(exp 17)/cm at the other. There are also a number of locations of such incomplete radiative emission where the shock has recently encountered denser regions with characteristic sizes of approximately 10(exp 18) cm. There is a considerable amount of evidence that the shock has decelerated from approximately 400 km/s to less than 200 km/s in the last 1000 yr. We interpret this as the result of the blast wave hitting the wall of a cavity which surround the supernova precursor and succeed in matching a wide range of data with a reflected shock model in which the density ofthe cavity wall is approximately 1.2/cu cm and the density in the interior of the cavity is about 0.09/cu cm. This analysis supports previous suggestions that the Cygnus Loop was a cavity explosion. We opredict that the X-ray shell will continue to brighten for the next 1000 years or so, during which time the shock front will become radiative throughout the northeast region. The similarity of the time scales for post-shock cooling, shock deceleration, and the 1ressure evolution of the region violates assumptions of steady flow, and greatly complicated interpretation of shock specptra in the region. This may be a common and important effect in many supernova remnants. High-resolutoin spectroscopy of the H alpha line profile shows that the broad component is surprisingly narrow (approximately 130 km/s). Consistency between this line width and the shock velocity inferred from spectrophotometry requires either rapid postshock equilibration of T(sub i) and T(sub e), or that the shock front has declerated from 180 km/s to 150 km/s in the last approximately 200 yr. While we cannot unambiguously distinquish between these models, the required deceleration is consistent with extrapolation of the shock decleration rate inferred over the last 1000 years. The width of the narrow component of the h alpha line profile is approximately 3 km/s, which indicates that the preshock medium has been heated to approximately 25,000 K by a shock precursor. The most promising types of precursor are heating either by fast neutrals which overtake the shock or by cosmic rays. In the latter case the observed properties of the precursor constrain the cosmic-ray diffusion coefficient to be between about 10(exp 23) sq cm/s. The characteristics of the precursor seem to be about right to explain the presence of weak (S II) and (N II) emission previously reported by Fesen & Itoh (1985).

Hester, J. Jeff↗

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↗