Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pattern matching”

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 397 records · Page 22

Analysis of a disk-on-rod surface wave element inside a corrugated horn using the mode-matching technique

A disk-on-rod inside a corrugated horn is one of the horn configurations for dual-frequency or wide-band operation. A mode-matching analysis method is described. A disk-on-rod inside a corrugated horn is represented as a series of coaxial waveguide sections and circular waveguide sections connected to each other. Three kinds of junctions need to be considered: coaxial-to-coaxial, coaxial-to-circular, and circular-to-circular. A computer program was developed to calculate the scattering matrix and the radiation pattern of a disk-on-rod inside a corrugated horn. The software as verified by experiment, and good agreement between calculation and measurement was obtained. The disk-on-rod inside a corrugated horn design gives an option to the Deep Space Network dual-frequency operation system, which currently is a two-horn/one-dichroic plate system.

Chen, J. C.↗

Reconstruction of the 1994 Pittsburgh Airplane Accident Using a Computer Simulation

On September 8, 1994, a Boeing 737-300 passenger airplane was on a downwind approach to the Pittsburgh International Airport at an altitude of 5000 feet above ground level (6000 feet MSL). While in a shallow left turn onto a downwind approach heading, the airplane crossed into the vortex trail of a Boeing 727 flying in the same approach pattern about 4 miles ahead. The B-737 airplane rolled and turned sharply to the left, exited the vortex wake and plunged into the ground. Weather was not a factor in the accident. The airplane was equipped with a 11+ channel digital Flight Data Recorder (FDR) and a multiple channel Cockpit Voice Recorder (CVR). Both recorders were recovered from the crash site and provided excellent data for the development of an accident scenario. Radar tracking of the two airplanes as well as the indicated air speed (IAS) perturbations clearly visible on the B-737 FDR recordings indicate that the upset was apparently initiated by the airplane's crossing into the wake of the B-727 flying ahead in the same traffic pattern. A 6 degree-of-freedom simulation program for the B-737 airplane using MATLAB and SIMULINK was constructed. The simulation was initialized at the stabilized flight conditions of the airplane about 13 seconds prior to its entry into the vortex trail of the B-727 airplane. By assuming a certain combination of control inputs, it was possible to produce a simulated motion that closely matched that recorded on the FDR.

Parks, Edwin K.↗

Phase Field Modeling of Microstructure Development in Microgravity

This newly funded project seeks to extend our NASA-sponsored project on modeling of dendritic microstructures to facilitate collaboration between our research group and those of other NASA investigators. In our ongoing program, we have applied advanced computational techniques to study microstructural evolution in dendritic solidification, for both pure isolated dendrites and directionally solidified alloys. This work has enabled us to compute dendritic microstructures using both realistic material parameters and experimentally relevant processing conditions, thus allowing for the first time direct comparison of phase field computations with laboratory observations. This work has been well received by the materials science and physics communities, and has led to several opportunities for collaboration with scientists working on experimental investigations of pattern selection and segregation in solidification. While we have been able to pursue these collaborations to a limited extent, with some important findings, this project focuses specifically on those collaborations. We have two target collaborations: with Prof. Glicksman's group working on the Isothermal Dendritic Growth Experiment (IDGE), and with Prof. Poirier's group studying directional solidification in Pb-Sb alloys. These two space experiments match well with our two thrusts in modeling, one for pure materials, as in the IDGE, and the other directional solidification. Such collaboration will benefit all of the research groups involved, and will provide for rapid dissemination of the results of our work where it will have significant impact.

Dantzig, Jonathan A.↗

Sub-Audible Speech Recognition Based upon Electromyographic Signals

Method and system for processing and identifying a sub-audible signal formed by a source of sub-audible sounds. Sequences of samples of sub-audible sound patterns ("SASPs") for known words/phrases in a selected database are received for overlapping time intervals, and Signal Processing Transforms ("SPTs") are formed for each sample, as part of a matrix of entry values. The matrix is decomposed into contiguous, non-overlapping two-dimensional cells of entries, and neural net analysis is applied to estimate reference sets of weight coefficients that provide sums with optimal matches to reference sets of values. The reference sets of weight coefficients are used to determine a correspondence between a new (unknown) word/phrase and a word/phrase in the database.

Jorgensen, Charles C.↗

Synthetic estimation filters for a control system

Hybrid optical correlators, which use spatial light modulators and solid state imagers to interface computers and optical Fourier correlators, will allow the rapid changing of filters and digital pre- or post-correlation processing in vision systems. However, optical correlation (particulary with phase-only filtration) is notoriously sensitive to mismatches of rotation and magnification between the reference image and the input image. Such sharpness is not always desirable if the correlation pattern is to be used in a control system. Two methods of conforming the response of an inherently sharp-response optical correlation to an in plane rotation between reference and object images are discussed. The objective is to condition the correlation for better suitability to a vision-controlled system. The first method is controlled convolutional blurring of a reference image followed by construction of a matched filter (perhaps implemented as a phase-only filter). An alternative is discussed in more detail: it is a synthesis from sharp filters, so weighted as to produce a specified off-center behavior.

Juday, Richard D.↗

A Short-Range Distance Sensor with Exceptional Linearity

A sensor has been demonstrated that can measure distance over a total range of about 300 microns to an accuracy of about 0.1 nm (resolution of about 0.01 nm). This represents an exceptionally large dynamic range of operation - over 1,000,000. The sensor is optical in nature, and requires the attachment of a mirror to the object whose distance is being measured. This work resulted from actively developing a white light interferometric system to be used to measure the depths of defects in the Space Shuttle Orbiter windows. The concept was then applied to measuring distance. The concept later expanded to include spectrometer calibration. In summary, broadband (i.e., white) light is launched into a Michelson interferometer, one mirror of which is fixed and one of which is attached to the object whose distance is to be measured. The light emerging from the interferometer has traveled one of two distances: either the distance to the fixed mirror and back, or the distance to the moving mirror and back. These two light beams mix and produce an interference pattern where some wavelengths interfere constructively and some destructively. Sending this light into a spectrometer allows this interference pattern to be analyzed, yielding the net distance difference between the two paths. The unique feature of this distance sensor is its ability to measure accurately distance over a dynamic range of more than one million, the ratio of its range (about 300 microns) to its accuracy (about 0.1 nanometer). Such a large linear operating range is rare and arises here because both amplitude and phase-matching algorithms contribute to the performance. The sensor is limited by the need to attach a mirror of some kind to the object being tracked, and by the fairly small total range, but the exceptional dynamic range should make it of interest.

Simmons, Steven↗

Decision Aid for Conjunction Risk Mitigation by Differential Drag

In the previous five years, the rate of conjunctions that the NASA Conjunction Assessment Risk Analysis (CARA) team processed and analyzed has more than tripled. (NASA CARA, 2024) New missions in the early development phases are now required to plan for dealing with conjunctions under the present space environment, and also projecting forward into a future likely with even further increased utilization of the space environment. Some missions are investigating the possibility of using differential drag to remediate conjunctions without expending limited fuel or for missions without propulsive capabilities. The NASA CARA team studied the historical record of conjunctions to evaluate the circumstances under which differential drag may be successfully applied and have developed a series of tables to use as a decision aid for missions considering differential drag. Currently, if a CARA-protected mission with maneuvering capabilities is predicted to have a conjunction with probability of collision (Pc) greater than 7E-5 (the default value of the ‘yellow threshold’, which may have some other value agreed by CARA and the mission during the Orbital Collison Avoidance Planning (OCAP) process), CARA will use its Maneuver Trade Space (MTS) tool to evaluate and recommend options for the timing and magnitude of a risk mitigation maneuver (RMM), based on the mission’s capabilities. If the conjunction’s Pc is greater than 1E-4, the ‘red threshold’, then an RMM must be executed per NASA Procedural Requirements (NPR) 8079.1 (NASA, 2023), although missions may execute an RMM even if the Pc is lower. The magnitude of the maneuver is typically a few cm/s, and CARA estimates how many will be required for the mission’s nominal lifetime – typically a few per year – during the OCAP process, to inform the mission’s delta-V requirement. For traditional satellites, this is usually smaller than other requirements for orbit insertion, maintenance, and disposal, but for CubeSats or other small satellite missions, a propulsion system may not provide much more than a few cm/s of delta-V or may not fit at all within the available budget of money, time, size, weight, and/or power (SWaP). Conversely, CubeSats often have deployable solar panels, which offer the capacity to have much higher areas facing some directions than others. Such a mission can instead use ‘differential drag’ to remediate a conjunction -- in other words, change its drag area (usually increasing) to deviate from the predicted collision course. This is how Planet’s Dove spacecraft maintain their formations and remediate conjunction risks without having on-board propulsion (Foster, et al., 2017) (Griffith, et al., 2021). CARA has been developing improvements to MTS to support differential-drag for NASA's missions -- where it is effective. CARA records all conjunctions of their protected payloads, with historical records starting in 2005 (with significant conjunction events starting to occur on or after 2013). From this record, approximately 7,300 had a Pc greater than 1E-4 at 3 days prior to the time of closest approach, the time analyzed for differential drag efficacy. Of those, approximately 4,300 had fully-defined covariance matrices stored for both the primary and secondary objects; this set of conjunctions is the basis for the analysis of this work. CARA’s MTS tool was used to propagate the primary satellite forward from that decision point with varying degrees of increase to ballistic coefficient (BC). Because these conjunctions came from multiple missions, the nondimensional ‘delta-BC’ factor was used to quantify and normalize the increase in ballistic coefficient, defined as follows: Delta-BC = BC_new / BC_old - 1 Positive delta-BC factors represent an increase in drag compared to the nominal attitude, while negative delta-BC factors (to a minimum of -1) represent a decrease in drag. At the conclusion of the differential-drag ‘maneuver’, the Pc was recalculated to evaluate whether or not the conjunction was mitigated (Pc < 3E-6). These results were then binned and sorted along several axes, including altitude, amount of delta-BC, and (pre-maneuver) rate of energy dissipation (EDR), to identify underlying patterns. The altitude plot is shown in Figure 1. To validate this analysis, we consulted the record of a NASA mission which uses differential drag to maintain its orbit and remediate conjunction risk. CARA’s empirical record of the mission’s orbit history suggests it achieves a delta-BC of 2.2. Of the twenty-one RMM plans that were submitted by this mission, eighteen were matched with conjunctions in the historical record; of those, twelve were successfully remediated (final measured Pc < 3E-6), and six were not. This is consistent with the expected efficacy for missions orbiting at that altitude. We are presently simulating this mission’s RMMs with MTS; this work is ongoing, but so far, the MTS results are qualitatively in agreement with the empirical results -- correctly predicting that a maneuver would or would not remediate a conjunction, if not exactly matching the final post-remediation Pc value. We found that differential drag was most successful for satellites with perigees below 560 km, and which could adopt an average delta-BC of 2 or greater (that is, increasing their ballistic coefficient by a factor of 3). However, this is a difficult threshold for a mission to clear; very few spacecraft are capable of adopting a high-drag configuration for 72 hours continuously. Planet’s Dove spacecraft use differential drag to remediate conjunctions (Griffith, et al., 2021), and they have a maximum delta-BC factor of 9, but in practice (with mission and charging constraints) they achieve a time-averaged delta-BC that is closer to 2 (Foster, et al., 2017). A mission’s differential drag utility strongly depends on the operational constraints that has the capacity to limit the time-averaged delta-BC. A mission with a high maximum delta-BC of 5 or more can have an effective delta-BC of less than 1 due to operational constraints such as instrument and solar panel pointing, especially if this constraint results in holding an intermediate drag value for most of its orbit. CARA has developed tables that can be used as decision aids to advise missions-in-development about the best way to utilize their differential drag capabilities. For missions below 560 km with the operational flexibility to devote multiple days to holding a high-drag configuration (or a sufficiently high drag ratio to compensate for limitations on that time), they are -- more likely than not -- able to successfully remediate high-risk conjunctions. Conversely, missions that do not meet these exacting criteria -- most missions -- can instead be advised to use on-board propulsion systems to perform RMMs, or to turn their minimum-area face towards the approach vector, thereby reducing Pc at the moment of conjunction due to the decreased Hard-Body Radius (HBR), that is a strongly correlated variable in the Pc calculations. (NASA, 2023)

conjunction assessment↗

MOVPE of GaSb/InGaAsSb Multilayers and Fabrication of Dual Band Photodetectors

Metalorganic vapor phase epitaxy (MOVPE) of GaSb/InGaAsSb multilayer thin films and fabrication of bias-selectable dual band photodetectors are reported. For the dual band photodetectors the short wavelength detector, or the upper p- GaSb/n-GaSb junction photodiode, is placed optically ahead of the long wavelength one, or the lower photodiode. The latter is based on latticed-matched In0.13Ga0.87As0.11Sb0.89 with bandgap near 0.6 eV. Specifically, high quality multilayer thin films are grown sequentially from top to bottom as p+-GaSb/p-GaSb/n-GaSb/n-InGaAsSb/p-InGaAsSb/p-GaSb on undoped p-type GaSb substrate, and as n-GaSb/p-GaSb/p-InGaAsSb/n-InGaAsSb/n-GaSb on Te-doped n-type GaSb substrate respectively. The multilayer thin films are characterized by optical microscope, atomic force microscope (AFM), electron microprobe analyses etc. The photodiode mesa steps are patterned by photolithography with wet chemical etching and the front metallization is carried out by e-beam evaporation with Pd/Ge/Au/Ti/Au to give ohmic contact on both n- and p-type Sb based layer surfaces. Dark I-V measurements show typical diode behavior for both the upper and lower photodiodes. The photoresponsivity measurements indicate that both the upper and lower photodiodes can sense the infrared illumination corresponding to their cutoff wavelengths respectively, comparable with the simulation results. More work is underway to bring the long wavelength band to the medium infrared wavelength region near 4 micrometers.

Xiao, Ye-Gao↗

Variability of particulate organic carbon concentration in the north polar Atlantic based on ocean color observations with Sea-viewing Wide Field-of-view Sensor (SeaWiFS)

We use satellite data from Sea-viewing Wide Field-of-view Sensor (SeaWiFS) to investigate distributions of particulate organic carbon (POC) concentration in surface waters of the north polar Atlantic Ocean during the spring summer season (April through August) over a 6-year period from 1998 through 2003. By use of field data collected at sea, we developed regional relationships for the purpose of estimating POC from remote-sensing observations of ocean color. Analysis of several approaches used in the POC algorithm development and match-up analysis of coincident in situ derived and satellite-derived estimates of POC resulted in selection of an algorithm that is based on the blue-to-green ratio of remote-sensing reflectance R(sub rs) (or normalized water-leaving radiance L(sub wn)). The application of the selected algorithm to a 6-year record of SeaWiFS monthly composite data of L(sub wn) revealed patterns of seasonal and interannual variability of POC in the study region. For example, the results show a clear increase of POC throughout the season. The lowest values, generally less than 200 mg per cubic meters, and at some locations often less than 50 mg per cubic meters, were observed in April. In May and June, POC can exceed 300 or even 400 mg per cubic meters in some parts of the study region. Patterns of interannual variability are intricate, as they depend on the geographic location within the study region and particular time of year (month) considered. By comparing the results averaged over the entire study region and the entire season (April through August) for each year separately, we found that the lowest POC occurred in 2001 and the highest POC occurred in 2002 and 1999.

Stramska, Malgorzata↗

Neural-net Processed Electronic Holography for Rotating Machines

This report presents the results of an R&D effort to apply neural-net processed electronic holography to NDE of rotors. Electronic holography was used to generate characteristic patterns or mode shapes of vibrating rotors and rotor components. Artificial neural networks were trained to identify damage-induced changes in the characteristic patterns. The development and optimization of a neural-net training method were the most significant contributions of this work, and the training method and its optimization are discussed in detail. A second positive result was the assembly and testing of a fiber-optic holocamera. A major disappointment was the inadequacy of the high-speed-holography hardware selected for this effort, but the use of scaled holograms to match the low effective resolution of an image intensifier was one interesting attempt to compensate. This report also discusses in some detail the physics and environmental requirements for rotor electronic holography. The major conclusions were that neural-net and electronic-holography inspections of stationary components in the laboratory and the field are quite practical and worthy of continuing development, but that electronic holography of moving rotors is still an expensive high-risk endeavor.

Decker, Arthur J.↗

Optimization of an oligonucleotide microchip for microbial identification studies: a non-equilibrium dissociation approach

The utility of a high-density oligonucleotide microarray (microchip) for identifying strains of five closely related bacilli (Bacillus anthracis, Bacillus cereus, Bacillus mycoides, Bacillus medusa and Bacillus subtilis) was demonstrated using an approach that compares the non-equilibrium dissociation rates ('melting curves') of all probe-target duplexes simultaneously. For this study, a hierarchical set of 30 oligonucleotide probes targeting the 16S ribosomal RNA of these bacilli at multiple levels of specificity (approximate taxonomic ranks of domain, kingdom, order, genus and species) was designed and immobilized in a high-density matrix of gel pads on a glass slide. Reproducible melting curves for probes with different levels of specificity were obtained using an optimized salt concentration. Clear discrimination between perfect match (PM) and mismatch (MM) duplexes was achieved. By normalizing the signals to an internal standard (a universal probe), a more than twofold discrimination (> 2.4x) was achieved between PM and 1-MM duplexes at the dissociation temperature at which 50% of the probe-target duplexes remained intact. This provided excellent differentiation among representatives of different Bacillus species, both individually and in mixtures of two or three. The overall pattern of hybridization derived from this hierarchical probe set also provided a clear 'chip fingerprint' for each of these closely related Bacillus species.

Non-NASA Center↗

Residual Study: Testing Jupiter Atmosphere Models Against Juno MWR Observations

The Juno spacecraft provides unique close-up views of Jupiter underneath the synchrotron radiation belts while circling Jupiter in its 53-day orbits. The microwave radiometer (MWR) onboard measures Jupiter thermal radiation at wavelengths between 1.37 and 50 cm, penetrating the atmosphere to a pressure of a few hundred bars and greater. The mission provides the first measurements of Jupiter's deep atmosphere, down to ~250 bars in pressure, constraining the vertical distributions of its kinetic temperature and constituents. As a result, vertical structure models of Jupiter's atmosphere may now be tested by comparison with MWR data. Taking into account the MWR beam patterns and observation geometries, we test several published Jupiter atmospheric models against MWR data. Our residual analysis confirms Li et al.'s (2017, https://doi.org/10.1002/2017GL073159) result that ammonia depletion persists down to 50–60 bars where ground-based Very Large Array was not able to observe. We also present an extension of the study that iteratively improves the input model and generates Jupiter brightness temperature maps which best match the MWR data. A feature of Juno's north-to-south scanning approach is that latitudinal structure is more easily obtained than longitudinal, and the creation of optimum two-dimensional maps is addressed in this approach.

Zhimeng Zhang↗

An Approach for Practical Grid-Resolved Roughness Aerothermodynamic Simulations for Woven Heatshield Surfaces

An approach for simulating the aerothermal environment over a patterned roughness surface within an otherwise axisymmetric flow is developed. The patterned roughness is simplified into a sinusoidal shape that captures the basic features of a charred 3MDCP ablator. Two orientations of this pattern, which are equally flight relevant due to the forming process, enable periodic boundary conditions along the edges of a grid spanning a single roughness element. This simplified three-dimensional grid represents the minimum problem size for a patterned roughness simulation, which enables computationally efficient grid-resolved roughness simulations. This approach is validated by simulating the Langley Mach 6 measurements made on a similar sinusoidal surface, which results in convective heating within 5% of the experimental data. Applying this approach to Mars Sample Return (MSR) Earth Entry System (EES) flight cases, considering simulations with both 11 species air and 30 species air with ablation products, results in heating augmentation that is slightly lower than the widely-used Dahm correlation approach, with the difference being dependent on pattern orientation. This result provides the only insight into the behavior of the roughness augmentation at the high temperature reacting flow conditions present for EES flight cases, which are not captured by available ground test measurements. For shear augmentation, the simulations revealed that the surface-parallel pressure component is dominant. For arc-jet tests targeted to flight values for shear, the higher arc-jet pressures required to match flight values for shear result in up to a 400% increase in the rough wall shear value, which is not captured by the Dahm correlation. This indicates that, to avoid significant over-testing, grid-resolved roughness simulations are required to determine arc-jet conditions that produce flight shear levels in rough-wall scenarios.

Christopher O Johnston↗

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↗

Response of Global Lightning Activity Observed by the TRMM/LIS During Warm and Cold ENSO Phases

This paper investigates the response of global lightning activity to the transition from the warm (January February March-JFM 1998) to the cold (JFM 1999) ENSO phase. The nine-year global lightning climatology for these months from the Tropical Rainfall Measuring Mission (TRMM) Lightning Imaging Sensor (LIS) provides the observational baseline. Flash rate density is computed on a 5.0x5.0 degree lat/lon grid within the LIS coverage area (between approx.37.5 N and S) for each three month period. The flash rate density anomalies from this climatology are examined for these months in 1998 and 1999. The observed lightning anomalies spatially match the documented general circulation features that accompany the warm and cold ENSO events. During the warm ENSO phase the dominant positive lightning anomalies are located mostly over the Western Hemisphere and more specifically over Gulf of Mexico, Caribbean and Northern Mid-Atlantic. We further investigate specifically the Northern Mid-Atlantic related anomaly features since these show strong relation to the North Atlantic Oscillation (NAO). Furthermore these observed anomaly patterns show strong spatial agreement with anomalous upper level (200 mb) cold core cyclonic circulations. Positive sea surface temperature anomalies during the warm ENSO phase also affect the lightning activity, but this is mostly observed near coastal environments. Over the open tropical oceans, there is climatologically less lightning and the anomalies are less pronounced. Warm ENSO related anomalies over the Eastern Hemisphere are most prominent over the South China coast. The transition to the cold ENSO phase illustrates the detected lightning anomalies to be more pronounced over East and West Pacific. A comparison of total global lightning between warm and cold ENSO phase reveals no significant difference, although prominent regional anomalies are located over mostly oceanic environments. All three tropical "chimneys" (Maritime Continent, Central Africa, and Amazon Basin) do not show any particular response to this transition.

Chronis, Themis G.↗

Spectroscopic characterization of collagen cross-links in bone

Collagen is the most abundant protein of the organic matrix in mineralizing tissues. One of its most critical properties is its cross-linking pattern. The intermolecular cross-linking provides the fibrillar matrices with mechanical properties such as tensile strength and viscoelasticity. In this study, Fourier transform infrared (FTIR) spectroscopy and FTIR imaging (FTIRI) analyses were performed in a series of biochemically characterized samples including purified collagen cross-linked peptides, demineralized bovine bone collagen from animals of different ages, collagen from vitamin B6-deficient chick homogenized bone and their age- and sex-matched controls, and histologically stained thin sections from normal human iliac crest biopsy specimens. One region of the FTIR spectrum of particular interest (the amide I spectral region) was resolved into its underlying components. Of these components, the relative percent area ratio of two subbands at approximately 1660 cm(-1) and approximately 1690 cm(-1) was related to collagen cross-links that are abundant in mineralized tissues (i.e., pyridinoline [Pyr] and dehydrodihydroxylysinonorleucine [deH-DHLNL]). This study shows that it is feasible to monitor Pyr and DHLNL collagen cross-links spatial distribution in mineralized tissues. The spectroscopic parameter established in this study may be used in FTIRI analyses, thus enabling the calculation of relative Pyr/DHLNL amounts in thin (approximately 5 microm) calcified tissue sections with a spatial resolution of approximately 7 microm.

Non-NASA Center↗

A Fourier approach to cloud motion estimation

A Fourier phase-difference technique for cloud motion estimation from pairs of pictures is described, and results obtained using this technique are compared with the results of a Fourier-domain cross-correlation scheme. The phase-difference technique makes use of the phase of the cross-spectral density and allows motion estimates to be made for individual spatial frequencies, which are related to cloud pattern dimensions. When objects being tracked do not change their shape, size, and orientation to more than a limited degree, both techniques are effective. The phase difference technique is relatively sensitive to the presence of mixtures of motions, changes in cloud shape, and edge effects; in these circumstances, the cross-correlation scheme is preferable. It is suggested that the Fourier transform phase difference estimation methods can be applied in problems such as landmark matching.

Arking, A.↗

A technique for interpretation of multispectral remote sensor data

The author has identified the following significant results. The U.S. Army Engineer Waterways Experiment Station is engaged in a study to detect from ERTS-1 satellite data alterations to the absorption and scattering properties caused by movement of suspended particles and solutes in selected areas of the Chesapeake Bay and to correlate the data to determine the feasibility of delineating flow patterns, flushing action of the estuary, and sediment and pollutant dispersion. As a part of this study, ADP techniques have been developed that permit automatic interpretation of data from any multispectral remote sensor with computer systems which have limited memory capacity and computing speed. The multispectral remote sensor is considered as a reflectance spectrophotometer. The data which define the spectral reflectance characteristics of a scene are scanned pixel by pixel. Each pixel whose spectral reflectance matches a reference spectrum is identified, and the results are shown in a map that identifies the locations where spectrum matches were detected and spectrum that was matched. The interpretation technique is described and an example of interpreted data from ERTS-1 is presented.

Williamson, A. N.↗