Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearity preservation”

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.

425 records · Page 24

Analog Techniques for Measuring the Frequency Response of Linear Physical Systems Excited by Frequency-Sweep Inputs

Data-reduction methods using general-purpose analog computer equipment and compatible testing techniques for determining the frequency response of linear physical systems are examined. The techniques considered may be classed as steady state or transient depending on the method of excitation. The relative merits of periodic, slow sweep, and transient (rapid sweep) forcing functions are discussed and applications are given that relate to dynamic-response tests of aeroelastic systems. Two frequency-sweep-input methods are considered in detail. one case the sweep rate is sufficiently slow that the response is approximately the same as that for steady-state conditions. With this input the frequency response can be evaluated and displayed in real time while the test is in progress. Errors due to treating sweep data as steady state can be eliminated, when desired, by reanalyzing tape-recorded time histories of the input and output as transient rather than as periodic data. In the second method the frequency-response function is deter- mined from the system's transient response to a very rapid sweep input. The purpose of frequency sweep in this case is to provide sufficient harmonic content in the input to overcome noise while keeping the test time as short as possible. tests and limited flight-test data presented herein, it appears that a transient-type rapid-sweep forcing function offers a considerable saving in test time while preserving the accuracy possible with steady-state sinusoidal inputs.

Wilmer H. Reed III↗

Ejecta patterns of Meteor Crater, Arizona derived from the linear un-mixing of TIMS data and laboratory thermal emission spectra

Accurate interpretation of thermal infrared data depends upon the understanding and removal of complicating effects. These effects may include physical mixing of various mineralogies and particle sizes, atmospheric absorption and emission, surficial coatings, geometry effects, and differential surface temperatures. The focus is the examination of the linear spectral mixing of individual mineral or endmember spectra. Linear addition of spectra, for particles larger than the wavelength, allows for a straight-forward method of deconvolving the observed spectra, predicting a volume percent of each endmember. The 'forward analysis' of linear mixing (comparing the spectra of physical mixtures to numerical mixtures) has received much attention. The reverse approach of un-mixing thermal emission spectra was examined with remotely sensed data, but no laboratory verification exists. Understanding of the effects of spectral mixing on high resolution laboratory spectra allows for the extrapolation to lower resolution, and often more complicated, remotely gathered data. Thermal Infrared Multispectral Scanner (TIMS) data for Meteor Crater, Arizona were acquired in Sep. 1987. The spectral un-mixing of these data gives a unique test of the laboratory results. Meteor Crater (1.2 km in diameter and 180 m deep) is located in north-central Arizona, west of Canyon Diablo. The arid environment, paucity of vegetation, and low relief make the region ideal for remote data acquisition. Within the horizontal sedimentary sequence that forms the upper Colorado Plateau, the oldest unit sampled by the impact crater was the Permian Coconino Sandstone. A thin bed of the Toroweap Formation, also of Permian age, conformably overlays the Coconino. Above the Toroweap lies the Permian Kiabab Limestone which, in turn, is covered by a thin veneer of the Moenkopi Formation. The Moenkopi is Triassic in age and has two distinct sub-units in the vicinity of the crater. The lower Wupatki member is a fine-grained sandstone, while the upper Moqui member is a fissile siltstone. Ejecta from these units are preserved as inverted stratigraphy up to 2 crater radii from the rim. The mineralogical contrast between the units, relative lack of post-emplacement erosion and ejecta mixing provide a unique site to apply the un-mixing model. Selection of the aforementioned units as endmembers reveals distinct patterns in the ejecta of the crater.

Ramsey, Michael S.↗

The ARADS Project as an Astrobiological Survey Mission Prototype

An astrobiological survey prior to human arrival will be important for preservation of signs of any existing or past evidence of Martian life, to reduce the risk of unintentional forward contamination of Martian environments in which terrestrial microbes might be able to gain a foothold (e.g., [1]), and to potentially protect the health of any human crews and enable their safe return to Earth. Evaluating for the presence of life and biosignatures may become a critical-path Mars exploration precursor in the not-so-far future, circa 2030-35. The near-subsurface, with its potential combination of volatiles, vapor exchange with the surface, and shelter from surface radiation, will be a primary initial target for in-situ astrobiological investigation. Sample acquisition, drilling and sample handling are key elements of robotic missions that target the subsurface [2], as well as near-surface ground ice deposits. In the mid 2000s, the Mars Astrobiology Field Laboratory (AFL) was proposed as a follow-on mission concept to the Mars Science Laboratory (MSL)/Curiosity rover mission [3]. The AFL would have included precision drilling to 2m, automated sample handling and remote processing of cores. However, the AFL concept needed further development and maturation. The Atacama Rover Astrobiology Drilling Studies (ARADS) project concept [4] was in part AFL-inspired. The 2014 ARADS concept included a medium sized rover, with habitability instruments, including the Phoenix Wet Chemistry Laboratory (WCL) [5], and an enhanced version of the ExoMars Mars Organic Molecule Analyzer (MOMA) Laser Desorption-Linear Ion Trap Mass Spectrometer (LITMS) [6]. ARADS likewise included two biosignature detection instruments -- the Signs of Life Detector (SOLID) immunoassay instrument [7,8] and the Chemical Laptop-derived [9] Planetary In-Situ Capillary Electrophoresis System (PISCES), while mitigating Planetary Protection-relevant contamination issues [10]. ARADS project objectives were to demonstrate the potential science return for such a mission in the context of the search for evidence of life, and technically to demonstrate the feasibility of roving with drilling missions to Mars. ARADS technology goals targeted analog-site demonstrations of automated drilling and sampling, integrated with its mobile habitability and biosignature detection instruments. And ARADS science operations goals developed remote science, contamination mitigation, and communications methods to enable simulating AFL-like science operations with time delays, directed by a remote science team[11]. Five subsequent Mars analog site field demonstrations (in Chile’s hyperarid Atacama Desert, 2016-19) provided a step-by-step development pathway for the ARADS platform’s science capabilities and systems integration, culminating in an integrated remote roving mission simulation in the final ARADS year (2019).

biosignatures↗

Total metals, carbon, nitrogen & anion concentration data; Slate River & East River floodplains, Crested Butte, CO; May 2022-October 2022

This data package includes processed and undiluted measurements for metal, total carbon, total nitrogen, and anion concentrations from pore water (groundwater) and surface water samples from the Slate River and East River floodplains of Crested Butte, CO, focus field sites for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and October of 2022. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS (metals):Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA. The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solutions (SPEX Certiprep, Metuchen, NJ). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace elements in water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. Lastly, a suitable internal standard (usually Rh, In, Ga or Ge) is added using on-line addition into the sample line and mixing tee.Analysis by Shimadzu TOC-L (TOC/TN):The TOC-L system is a combustion technique where liquid samples are injected and combusted into CO2 for carbon detection by non-dispersive infrared (NDIR) and NO for detection by chemiluminescence. A calibration curve using five standard solutions between 0.1 and 7 ppm for carbon and 0.05 and 3.5 ppm for nitrogen is made for each type of measurement with a linearity >0.99. All samples, standards, and QC’s are prepared in 24mL scintillation vials that have been baked for 4hrs at 475 Cº and made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.Analysis by Ion Chromatography (Anions):The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with sodium carbonate eluent. A calibration curve using five standard solutions between 5 and 250 umol/L is made with a linearity >0.99. Standards and QC’s are prepared in 15mL polypropylene conical tubes, pipetted along with the samples into 1.5mL polypropylene vials. Dilutions are made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗

Machine learning for nonintrusive model order reduction of the parametric inviscid transonic flow past an airfoil

Fluid flow in the transonic regime finds relevance in aerospace engineering, particularly in the design of commercial air transportation vehicles. Computational fluid dynamics models of transonic flow for aerospace applications are computationally expensive to solve because of the high degrees of freedom as well as the coupled nature of the conservation laws. While these issues pose a bottleneck for the use of such models in aerospace design, computational costs can be significantly minimized by constructing special, structure-preserving surrogate models called reduced-order models. In this work, we propose a machine learning method to construct reduced-order models via deep neural networks and we demonstrate its ability to preserve accuracy with a significantly lower computational cost. In addition, our machine learning methodology is physics-informed and constrained through the utilization of an interpretable encoding by way of proper orthogonal decomposition. Application to the inviscid transonic flow past the RAE2822 airfoil under varying freestream Mach numbers and angles of attack, as well as airfoil shape parameters with a deforming mesh, shows that the proposed approach adapts to high-dimensional parameter variation well. Notably, the proposed framework precludes the knowledge of numerical operators utilized in the data generation phase, thereby demonstrating its potential utility in the fast exploration of design space for diverse engineering applications. Comparison against a projection-based nonintrusive model order reduction method demonstrates that the proposed approach produces comparable accuracy and yet is orders of magnitude computationally cheap to evaluate, despite being agnostic to the physics of the problem.

97 MATHEMATICS AND COMPUTING↗

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ambient Mass Spectrometry Rapid Analysis Methods for Trace Inorganics

The ability to rapidly detect and characterize the chemistry of trace quantities of inorganic species has multiple potential applications, including nuclear safeguards, nuclear forensics, and environmental monitoring. Some mass spectrometric methods currently exist, but they are costly and time-intensive, often destroying chemical information in the process of ionization. Ambient mass spectrometry is a relatively new analytical method in which samples are introduced to the mass spectrometer from atmospheric pressure, rather than from high vacuum. Ambient MS can operate with gentler ionization methods such as electrospray ionization (ESI) and paper-spray ionization (PSI), reducing fragmentation and preserving the chemical information of the analyte. This may be useful for determining the intended purpose of material in a sample. Strontium, part of nuclear fallout, is analyzed here to test these new analytical methods for viability. Electrospray ionization mass spectrometry functions by spraying the sample as a solution through a needle charged to a few kV. A nebulizing gas flows out along with the sample to create the spray before droplets enter a drying chamber with heated drying gas to remove the solvent molecules. From there bare analyte ions enter the mass spectrometer. This method has the advantage of relatively fast sample preparation and stable ionization. More stable ionization produces higher-quality data. Paper-spray ionization mass spectrometry operates by applying sample to a paper triangle. The paper triangle is then electrified in a metal clamp to a few kV before spray solvent is applied. The spray solvent transports analyte in the sample to the tip, where it is ionized and enters the mass spectrometer. One advantage of this method is reduced sample preparation - an analyst can even swipe a surface and cut a triangle directly from the swipe. The data analysis serves to rapidly process tens of samples all at once. The mass spectrometer outputs a file containing a few hundred spectra for each sample. The pipeline processes each one to identify peaks of interest. A full palette of data analysis metrics are generated, including plots of the raw data and plots of calculated values such as isotopic ratios. The pipeline filters the raw data to remove noise and background peaks, as well as peaks representing detector saturation. A heat map is generated to show the overall spectrum in mass-time- intensity space. Post-processing, we obtained the desired calibration curves, demonstrating that both ESI-MS and PSI-MS are viable methods for quantifying trace amounts of inorganics while retaining their chemical speciation information. ESI-MS demonstrated a linear response region of 100 ppb to 20 ppm, with a limit-of-detection of around 100 ppb. At top left is an example spectrum with the {sup 88}SrNO{sub 3} peak identified to show the retention of chemical speciation. PSI-MS showed a linear response region of around 1 ppb to 1 ppm, with a limit-of-detection of near 1 ppb. PSI-MS produced more erratic results than ESI-MS, likely due to the variability in ionization resulting from manual positioning of the paper triangle near the MS inlet. Future efforts will focus on improving the repeatability of PSI-MS as well as analyzing other species of interest such as cerium and uranium. Improvements will also be made to create more reliable data extraction and intelligent data filtering.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of a Robust Reference Electrode in Aggressive Chemical and Radiation Environments in the Hanford Waste Tanks

The Hanford site stores more than 200 million liters of radioactive and chemically hazardous wastes from the production of weapons materials. The wastes are stored in 177 underground carbon-steel storage tanks, separated between 149 single shell tanks (SSTs) and 28 double shell tanks (DSTs). The DSTs provide critical retrieval and interim storage before the waste is vitrified in the Waste Treatment and Isolation Plant (WTP). The tanks will need to remain in-service far beyond the initial 40-year design life, and effective corrosion control practices must remain in force to extend the tanks’ lifespans. This effort includes direct measurements of corrosion rate (e.g., ultrasonic measurements and corrosion coupons) and electrochemical processes (e.g., linear polarization measurements and open circuit potential measurements). The Hanford site began monitoring the corrosion potential in select DSTs in 2008. Of the 45 reference electrodes that have been installed, 29 have failed and 6 others provided unreliable results. DOE-EM is supporting a 3-year program to develop a chemical and radiation resistant reference electrode for application in the Hanford tanks. The first year of the program focused on understanding the failure mechanism for the reference electrodes and identification of candidate construction materials that would mitigate degradation of the electrodes in the waste environment. During the second year of the program, the objectives were to: 1) test candidate materials under simulated waste conditions, 2) design components that will extend the service life of the electrode, 3) fabricate materials for prototype reference electrodes, and 4) assemble prototype reference electrodes for accelerated testing. The reference electrode is constructed of four principal parts: 1) junction, 2) casing, 3) inner chamber backfill materials, and 4) the sensing wire. Principally, improvements of the junction, casing, and inner chamber backfill materials are being pursued. The junction material at the interface between the waste and the inner chamber of the reference electrode was identified as a critical component in the failure of the reference electrodes. Nine candidate replacement junction materials were tested under simulated waste conditions to evaluate permeation rate. These materials included a variety of polymeric and ceramic materials, some of which were 3-D printed. Thus far, porous polyvinylidene fluoride materials have performed satisfactorily and are being considered for prototype development. The commercial electrode casing materials in general have performed well. Additionally, 3-D printing of chemically and mechanically stable materials is being investigated as a means for further improvement in fabrication consistency. SRNL has also investigated altering the reference electrode design to extend the service life. The new design of the interior of the reference electrode casing creates a longer, more tortuous path between the junction material and the electrode sensing wire. A finite element model was used to optimize the design without adversely impacting the circuit resistance of the electrode during the measurements, thus preserving the measurement accuracy while enhancing the service life. The inner chamber back fill materials are also critical to the performance of the reference electrode. Materials that are resistant to intruding tank waste and provide a conductive path to the sensing wire were investigated. Gel and powder materials that are interspersed with a conductive chloride bearing material were tested for their influence on diffusion and electrode resistance. All the investigated materials and components will be assembled, with collaboration from commercial vendors, to fabricate the initial prototypes. Accelerated testing of the prototypes will be initiated in Year 2 of the program and will be completed in Year 3. A recommendation on the materials of construction and the design of the new robust reference electrode will be presented to the Hanford tank farm facility.

Sykes, Kiana [Savannah River National Laboratory (↗

Zero-Gravity Locomotion Simulators: New Ground-Based Analogs for Microgravity Exercise Simulation

Maintaining health and fitness in crewmembers during space missions is essential for preserving performance for mission-critical tasks. NASA's Exercise Countermeasures Project (ECP) provides space exploration exercise hardware and monitoring requirements that lead to devices that are reliable, meet medical, vehicle, and habitat constraints, and use minimal vehicle and crew resources. ECP will also develop and validate efficient exercise prescriptions that minimize daily time needed for completion of exercise yet maximize performance for mission activities. In meeting these mission goals, NASA Glenn Research Center (Cleveland, OH, USA), in collaboration with the Cleveland Clinic (Cleveland, Ohio, USA), has developed a suite of zero-gravity locomotion simulators and associated technologies to address the need for ground-based test analog capability for simulating in-flight (microgravity) and surface (partial-gravity) exercise to advance the health and safety of astronaut crews and the next generation of space explorers. Various research areas can be explored. These include improving crew comfort during exercise, and understanding joint kinematics and muscle activation pattern differences relative to external loading mechanisms. In addition, exercise protocol and hardware optimization can be investigated, along with characterizing system dynamic response and the physiological demand associated with advanced exercise device concepts and performance of critical mission tasks for Exploration class missions. Three zero-gravity locomotion simulators are currently in use and the research focus for each will be presented. All of the devices are based on a supine subject suspension system, which simulates a reduced gravity environment by completely or partially offloading the weight of the exercising test subject s body. A platform for mounting treadmill is positioned perpendicularly to the test subject. The Cleveland Clinic Zero-g Locomotion Simulator (ZLS) utilizes a pneumatic subject load device to apply a near constant gravity-replacement load to the test subject during exercise, and is currently used in conjunction with the General Clinical Research Center for evaluating exercise protocols using a bedrest analog. The enhanced ZLS (eZLS) at NASA Glenn Research Center features an offloaded treadmill that floats on a thin film of air and interfaces to a force reaction frame via variably-compliant isolators, or vibration isolation system. The isolators can be configured to simulate compliant interfaces to the vehicle, which affects mechanical loading to crewmembers during exercise, and has been used to validate system dynamic models for new countermeasures equipment designs, such as the second International Space Station treadmill slated for use in 2010. In the eZLS, the test subject and exercise device can be pitched at the appropriate angle for partial gravity simulations, such as lunar gravity (1/6th earth gravity). On both the eZLS and the NASA-Johnson Space Center standalone ZLS installed at the University of Texas Medical Branch in Galveston, Texas, USA, the subject's body weight relative to the treadmill is controlled via a linear motor subject load device (LM-SLD). The LM-SLD employs a force-feedback closed-loop control system to provide a relatively constant force to the test subject during locomotion, and is set and verified for subject safety prior to each session. Locomotion data were collected during parabolic flight and on the eZLS. The purpose was to determine the similarities and differences between locomotion in actual and simulated microgravity. Subjects attained greater amounts of hip flexion during walking and running during parabolic flight. During running, subjects had greater hip range of motion. Trunk motion was significantly less on the eZLS than during parabolic flight. Peak impact forces, loading rate, and impulse were greater on the eZLS than during parabolic while walking with a low external load (EL) and rning with a high EL. Activation timing differences existed between locations in all muscles except for the rectus femoris. The tibialis anterior and gluteus maximus were active for longer durations on the eZLS than in parabolic flight during walking. Ground reaction forces were greater with the LM-SLD than with bungees during eZLS locomotion. While the eZLS serves as a ground-based analog, researchers should be aware that subtle, but measurable, differences in kinematics and leg musculature activities exist between the environments. Aside from space applications, zero-gravity locomotion simulators may help medical researchers in the future with development of rehabilitative or therapeutic protocols for injured or ill patients. Zero-gravity locomotion simulators may be used as a ground-based test bed to support future missions for space exploration, and eventually may be used to simulate planetary locomotion in partial gravity environments, including the Moon and Mars. Figure: Zero-gravity Locomotion Simulator at the Cleveland Clinic, Cleveland, Ohio, USA

Perusek, Gail P.↗

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 capability of the Position, Navigation, and Timing (PNT) system to support their scientific objectives. In addition, NASA's Commercial Lunar Payload Services (CLPS) program initiates the proliferation of public and private exploration partnerships using small satellites from commercial and private organizations, expanding traditionally confined low Earth orbit to be used for missions beyond geosynchronous orbit (Zucherman et al., 2022). Therefore, the Lunar PNT system is also required to provide navigation services compatible with the smaller platforms being sent by the public and private sectors, like CubeSats. However, traditional approaches to deep space missions’ navigation based on ground radio facilities have difficulties in providing sufficient support for the increasing number of users and communication at a distance from the Earth (Kaplev et al., 2022). In particular, the existing Lunar navigation technologies such as weak signal global positioning system (GPS) and deep space network (DSN) are not able to ensure operations of the upcoming small-scale Lunar missions due to their limitations in localization performance as well as capacity aspects. Another way to provide Lunar PNT service is to create a dedicated Lunar global navigation satellite system (GNSS) constellation, like GNSS systems on Earth. Space agencies like NASA, ESA, and JAXA are now developing the lunar communications relay and navigation systems (LCRNS) and Lunar navigation satellite systems (LNSS). In their systems, satellites will be deployed in moon orbits to provide the communication, positioning, navigation, and timing (CPNT) service at the lunar south pole region where the Artemis base camp will be expected (Murata et al., 2022). Meanwhile, common challenges considered in lunar PNT research arise from poor geometry of the terrestrial GNSS satellites when seen from the lunar user, highly perturbed lunar orbits, and limitations in power, size, and cost of the equipment on lunar satellites (Iiyama et al., 2023). It is also not clear if there will be enough Lunar users to support the cost and resources this would require as the Low-cost surface missions may not be able to support the large power, mass, and weight requirements that these navigation solutions entail (Niemoeller et al., 2022). As an alternative, existing Lunar science and exploration assets could 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 (Hagenau et al., 2021). Introducing the non-dedicated and ad-hoc Lunar navigation constellation gives a way to provide PNT services on-demand. The non-dedicated swarm assets of Lunar constellations are designed to localize themselves with minimal interaction with Earth by adding cooperative autonomous localization to lunar missions, freeing up valuable bandwidth and ground segment resources. An autonomous localization of Lunar constellations is based on the concept of the decentralized PNT system with a distributed extended Kalman filter (DEKF) approach to state estimation for minimal onboard operating costs. In the distributed data processing algorithm, computation is broken down and assigned to each satellite, resulting in a considerably decreased computational amount while maintaining the accuracy of the orbit ephemeris and clock offsets as the result of centralized data processing (Wen et al., 2019). The DEKF requires spacecraft to perform two-way ranging operations with each other to communicate simultaneously, leveraging neighbor two-way intersatellite link (ISL) measurements such as pseudoranges to, and relative velocities between, visible satellites as sensor values (Frank et al., 2021). The Lunar autonomous PNT simulation (LAPS) demonstrated the feasibility of orbital asset localization among ad-hoc Lunar small-sat constellations based on the DEKF in Hagenau et al. (2021) and evaluated the matching algorithm proposed by Frank et al. (2021) in scheduling position estimation updates. In previous papers, all assets and measurements are assumed to be always available without consideration of the impact of intermittent and permanent communication failure. This study presents localization performance with increasing levels of network degradation for swarm assets and users to demonstrate the robustness of the decentralized Lunar PNT service in more realistic scenarios. Main issues arising from communication failure include spacecraft permanent or transient loss, antenna failures, message delays, etc. We tested four possible reasons for network degradation for 7 days in 21 satellites frozen with an altitude of 5500 km, evenly spaced around 3 circular, 40 inclination orbital planes where each spacecraft has two directional antennas. As anchor nodes with an independent estimate of their position are required in the DEKF approach, two ground nodes in each pole and one node in the gateway were implemented in the simulation. First, the most probable failure scenario involves the loss of a single spacecraft due to solar interference and technical malfunctions of the assets. Losing the availability of a single spacecraft means losing the two-way ISL measurement of the asset in the DEKF update. In order to provide the best possible quality of PNT service with limited time and resources, the distributed Lunar constellations must schedule the communication activities. The scheduler leverages mixed-integer linear programming (MILP) for the coordination and scheduling of the desired “as-needed” localization service (Niemoeller et al., 2022). We assume the scheduler has completely excluded the spacecraft information before the DEKF update in the failure scenario. When a random spacecraft has been turned off at a specific time, the robustness of the autonomous Lunar PNT system is evaluated. The simulation results give an 11.5% degradation in median position accuracy compared to the idealized performance excluding the asset loss. Second, a large number of assets may vanish due to major hardware problems or meteor strikes around the moon. A multiple spacecraft loss can degrade the localization performance very fast by losing the communication ability to do cross-plane measurements and in-plane measurements in a 3-plane constellation. When the matching-based scheduler is aware of ISL availability, we investigate a large number of in-plane and cross-plane asset vanishments both in close proximity and equally spaced throughout the orbital plane. According to the simulations, the loss of in-plane measurements gives 40.2% degradation while cross-plane measurements degrade 50.5% of asset localization performance among available assets. Therefore, it is concluded that cross-plane measurements are more important in improving the position estimation accuracy. Third, spacecraft failure information can be lost due to the internal message delay, resulting in the DEKF update scheduler to solve the matching problem with unavailable assets. The DEKF update cycle is comprised of network setup, communication, and computations where a global broadcast network and a 2-way ISL network setup take 6 minutes in total (Frank et al., 2021). Once the broadcast network successfully transmits and receives information, a random spacecraft may lose its availability right before solving the matching problem. This means the matching solution is no longer optimal, resulting in degradation in the localization performance. A numerical assessment shows the matching-based scheduler with knowing failure holds 11.5% of position accuracy degradation, whereas the scheduler without knowing failure gives 34% degraded localization performance without asset loss. Fourth, a transient loss of a single or multiple spacecraft may occur due to their antenna outages. After losing the two-way ISL availability for a few DEKF update cycles, the availability of spacecraft can easily be recovered as their states have been independently updated using measurements from anchor nodes. It is likely that the longer failure will result in worse localization performance. We have tested the transient failure of a random single asset for 30 min in the simulation, which is losing 3 update cycles in the DEKF system. From the simulation results, the position accuracy has been degraded to 4.84% which is better than the degraded localization performance of 11.5% from the permanent loss scenario among available assets. In conclusion, the autonomous Lunar PNT system based on the DEKF approach shows the ability to maintain resilience and robustness in the possible communication failure scenarios, ensuring that localization accuracy is preserved across various network degradation and outages. Future studies on investigating user localization performance near the South Pole and the broadcast network system will be continued in the following months.

Yeji Kim↗

Analyses of Cometary Silicate Crystals: DDA Spectral Modeling of Forsterite

Comets are the Solar System's deep freezers of gases, ices, and particulates that were present in the outer protoplanetary disk. Where comet nuclei accreted was so cold that CO ice (approximately 50K) and other supervolatile ices like ethane (C2H2) were preserved. However, comets also accreted high temperature minerals: silicate crystals that either condensed (greater than or equal to 1400 K) or that were annealed from amorphous (glassy) silicates (greater than 850-1000 K). By their rarity in the interstellar medium, cometary crystalline silicates are thought to be grains that formed in the inner disk and were then radially transported out to the cold and ice-rich regimes near Neptune. The questions that comets can potentially address are: How fast, how far, and over what duration were crystals that formed in the inner disk transported out to the comet-forming region(s)? In comets, the mass fractions of silicates that are crystalline, f_cryst, translate to benchmarks for protoplanetary disk radial transport models. The infamous comet Hale-Bopp has crystalline fractions of over 55%. The values for cometary crystalline mass fractions, however, are derived assuming that the mineralogy assessed for the submicron to micron-sized portion of the size distribution represents the compositional makeup of all larger grains in the coma. Models for fitting cometary SEDs make this assumption because models can only fit the observed features with submicron to micron-sized discrete crystals. On the other hand, larger (0.1-100 micrometer radii) porous grains composed of amorphous silicates and amorphous carbon can be easily computed with mixed medium theory wherein vacuum mixed into a spherical particle mimics a porous aggregate. If crystalline silicates are mixed in, the models completely fail to match the observations. Moreover, models for a size distribution of discrete crystalline forsterite grains commonly employs the CDE computational method for ellipsoidal platelets (c:a:b=8.14x8.14xl in shape with geometrical factors of x:y:z=1:1:10, Fabian et al. 2001; Harker et al. 2007). Alternatively, models for forsterite employ statistical methods like the Distribution of Hollow Spheres (Min et al. 2008; Oliveira et al. 2011) or Gaussian Random Spheres (GRS) or RGF (Gielen et al. 200S). Pancakes, hollow spheres, or GRS shapes similar to wheat sheaf crystal habit (e.g., Volten et al. 2001; Veihelmann et al. 2006), however, do not have the sharp edges, flat faces, and vertices seen in images of cometary crystals in interplanetary dust particles (IDPs) or in Stardust samples. Cometary forsterite crystals often have equant or tabular crystal habit (J. Bradley). To simulate cometary crystals, we have computed absorption efficiencies of forsterite using the Discrete Dipole Approximation (DDA) DDSCAT code on NAS supercomputers. We compute thermal models that employ a size distribution of discrete irregularly shaped forsterite crystals (nonspherical shapes with faces and vertices) to explore how crystal shape affects the shape and wavelength positions of the forsterite spectral features and to explore whether cometary crystal shapes support either condensation or annealing scenarios (Lindsay et al. 2012a, b). We find forsterite crystal shapes that best-fit comet Hale-Bopp are tetrahedron, bricks or brick platelets, essentially equant or tabular (Lindsay et al. 2012a,b), commensurate with high temperature condensation experiments (Kobatake et al. 2008). We also have computed porous aggregates with crystal monomers and find that the crystal resonances are amplified. i.e., the crystalline fraction is lower in the aggregate than is derived by fitting a linear mix of spectral features from discrete subcomponents, and the crystal resonances 'appear' to be from larger crystals (Wooden et al. 2012). These results may indicate that the crystalline mass fraction in comets with comae dominated by aggregates may be lower than deduced by popular methods that only emoy ensembles of discrete crystals.

Wooden, Diane↗