Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter continuation method”

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.

478 records · Page 27

Development of a Tunable Electromechanical Acoustic Liner for Engine Nacelles

This report describes the development of a tunable electromechanical Helmholtz resonator (EMHR) for engine nacelles using smart materials technology. This effort addresses both near-term and long-term goals for tunable electromechanical acoustic liner technology for the Quiet Aircraft Technology (QAT) Program. Analytical models, i.e. lumped element model (LEM) and transfer matrix (TM) representation of the EMHR, have been developed to predict the acoustic behavior of the EMHR. The models have been implemented in a MATLAB program and used to compare with measurement results. Moreover, the prediction performance of models is further improved with the aid of parameter extraction of the piezoelectric backplate. The EMHR has been experimentally investigated using standard two-microphone method (TMM). The measurement results validated both the LEM and TM models of the EMHR. Good agreement between predicted and measured impedance is obtained. Short- and open circuit loads define the limits of the tuning range using resistive and capacitive loads. There is approximately a 9% tuning limit under these conditions for the non-optimized resonator configuration studied. Inductive shunt loads result in a 3 degree-of-freedom DOF) system and an enhanced tuning range of over 20% that is not restricted by the short- and open-circuit limits. Damping coefficient ' measurements for piezoelectric backplates in a vacuum chamber are also performed and indicate that the damping is dominated by the structural damping losses, such as compliant boundaries, and other intrinsic loss mechanisms. Based on models of the EMHR, a Pareto optimization design of the EMHR has been performed for the EMHR with non-inductive loads. The EMHR with non-inductive loads is a 2DOF system with two resonant fiequencies. The tuning ranges of the two resonant frequencies of the EMHR with non-inductive loads cannot be optimized simultaneously; a trade-off (i.e., a Pareto solution) must be reached. The Pareto solution provides the information for a designer that shows how design trade-offs can be used to satisfy specific design requirements. The optimization design of the EMHR with inductive loads aims at optimal tuning of these three resonant fiequencies. The results indicate that it is possible to keep the acoustic reactance of the resonator close to a constant over a given frequency range. An effort to mimic the second layer of the NASA 2DOF liner using a piezoelectric composite diaphragm has been made. The optimal acoustic reactance of the second layer of the NASA 2DOF liner is achieved using a thin PVDF composite diaphragm, but matching the acoustic resistance requires further investigation. Acoustic energy harvesting is achieved by connecting the EMHR to an energy reclamation circuit that converts the ac voltage signal across the piezoceramic to a conditioned dc signal. Energy harvesting experiment yields 16 m W continuous power for an incident SPL of 153 dB. Such a level is sufficient to power a variety of low power electronic devices. Finally, technology transfer has been achieved by converting the original NASA ZKTL FORTRAN code to a MATLAB code while incorporating the models of the EMHR. Initial studies indicate that the EMHR is a promising technology that may enable lowpower, light weight, tunable engine nacelle liners. This technology, however, is very immature, and additional developments are required. Recommendations for future work include testing of sample EMHR liner designs in NASA Langley s normal incidence dual-waveguide and the grazing-incidence flow facility to evaluating both the impedance characteristics as well as the energy reclamation abilities. Additional design work is required for more complex tuning circuits with greater performance. Poor electromechanical coupling limited the electromechanical tuning capabilities of the proof of concept EMHR. Different materials than those studies and perhaps novel composite material systems may dramatically improvehe electromechanical coupling. Such improvements are essential to improved mimicking of existing double layer liners.

Liu, Fei↗

Active Control of Rotating Stall Demonstrated for a Multistage Compressor With Inlet Distortion

Aircraft compressors can suffer debilitating consequences as a result of rotating stall and surge events caused by inlet distortions. This is particularly true of aircraft during takeoff, when the compressor is operating at peak performance close to the surge line. Significant research has been conducted by the NASA Glenn Research Center in the area of compressor stability enhancement through active and passive control methods. Most recently, an experiment was conducted at the Wright Patterson Air Force Base Research Laboratory on a two-stage fan with inlet guide vanes and inlet distortion. In this joint Small Business Innovation Research effort between Scientific Systems and Glenn, control of rotating stall was demonstrated in a multistage transonic fan. This twostage fan with inlet guide vanes was tested under clean and distorted inlet conditions. The compressor was also configured with a circumferential distortion screen capable of 180 of distortion and with 14 high-velocity injectors upstream of the first rotor. Twelve of these injectors could oscillate up to frequencies of 450 Hz. The additional two injectors were located next to each other and were used in concert with each other as a single, on/off, high-authority actuator. In a first test of injection in this multistage environment, 12 of the valves were opened 50 percent of their full stroke to assess steady injection through the compressor. This baseline injection is shown in the compressor characteristic of the following figure, and stall margin improvements are tracked from this baseline condition. The compressor was then tested with clean inlet conditions using 12 injectors and active control. Pressure disturbances were tracked before rotating stall, and a constant gain control scheme reduced the stalling mass flow by 10.8 percent over the baseline. With the distortion screen present in the inlet, a pole-zero cancellation control scheme was used to achieve a 6.4-percent decrease in stalling mass flow. These improvements also are shown in the figure. In a final experiment, actively controlled, high-frequency injection from the 12 valves was used in conjunction with the high-authority actuators. In this test, the stalling mass flow of the compressor was reduced by 27 percent as indicated in the graph. These results were obtained by injecting less than 2 percent of the total compressor throughflow into the rotor tip region via 14 injection ports. These results mark the first successful demonstration of actively controlled air injection as a stall-control strategy for multistage compressors operating at speeds typical of an actual gas turbine engine. A goal of continuing research is to determine the combination of air-injection parameters and control strategies that are most effective in providing stall control for both clean and distorted inlet flow conditions for multistage environments. Other goals include the demonstration of stall control at many locations along the core compressor and development and application of active stall control strategies that will be integral flightworthy components of onboard engine hardware.

VanSchalkwyk, Christian↗

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗

Intracranial Effects of Intermittent Lower Body Negative Pressure with Head Down Tilt Bed Rest: Comparison to Upright Posture

INTRODUCTION Spaceflight-associated neuro-ocular syndrome (SANS) affects a majority of astronauts during long-duration spaceflight. SANS is hypothesized to result from an unrelenting headward fluid shift that occurs in the microgravity environment. Altered intracranial structure and physiology documented in spaceflight and head-down tilt (HDT) experiments(1, 2) are also theorized to be related to chronic headward fluid shift and thereby can provide an independent quantitative assessment related to this mechanism. As a potential countermeasure, lower body negative pressure (LBNP) applied in the supine posture has shown efficacy in reducing headward fluid shift(3). However, it remains unknown if intermittent LBNP simulating daily upright posture fluid redistribution can mitigate SANS and intracranial changes associated with long-duration spaceflight. The goal of this study was to quantify the effects of daily application of LBNP or daily exposure to the upright posture on intracranial structure and physiology during long-term HDT in order to evaluate the potential efficacy of LBNP as a countermeasure to SANS. METHODS 22 healthy volunteers (11 men; 11 women; mean age = 35 years (SD=9.2)(age range = 24 to 52 years); mean BMI = 24.0 kg/m2 (SD=2.8)) completed an MRI study performed at the German Aerospace Facility in Cologne, Germany. Strict six-degree head-down tilt (HDT) bedrest was used as a spaceflight analog to induce a continuous headward fluid shift for 30 days. The subjects were divided equally into two groups of interventions: 1) LBNP for 3-hour sessions twice daily and 2) seated position for 3-hour sessions twice daily. Interventions were divided into morning and afternoon sessions for both groups. The LBNP intervention was maintained at 25 mmHg  2 mmHg. Pulse-gated MRI phase-contrast flow imaging was used to quantify cerebral artery stroke volume (CASV) and peak-to-peak cerebral spinal fluid (CSF) velocity (CSFVp-p) within the cerebral aqueduct. A 3D T1-MPRAGE sequence was used to quantify volumetric changes of the brain and intracranial CSF spaces using MRI Cloud software. MRI acquisitions were obtained at baseline (BDC, supine posture), 15 days into HDT (HDT15), 29 days into HDT (HDT29), and 12 days after recovery (R12, supine posture). The data were analyzed by a mixed model, which included intervention, time (four nominal levels BDC, HDT15, HDT29, R12), and intervention-time interaction as the fixed effects and included subject as a random effect. RESULTS Compared to BDC there was no statistically significant difference in CASV and CSFVp-p during HDT except for CASV at HDT29 (seated group) where there was a 1.7 mL (12%) decrease (P<0.01). Compared to BDC there was a statistically significant increase in intracranial volume (ICV; ICV = white matter + gray matter + CSF) for both interventions at HDT15 (Δ13 mL, 0.9 % (LBNP group); Δ15 mL, 1.0%, (seated group)) and HDT29 (Δ19 mL, 1.2%, (LBNP group); Δ23 mL, 1.5%) (seated group)) (All Ps<.001). Compared to baseline, lateral ventricular volume increased at HDT15 (Δ0.9 mL, 5.5%) and HDT29 (Δ1.8 mL, 10%)) for LBNP only (P=.001 for each). During HDT, white matter volume remained stable compared to baseline for both interventions. There were no significant intervention effects in the overall response to HDT (All Ps >.2). CONCLUSION LBNP had a similar response to the seated posture during 29 days of HDT. Although there was an increase in ICV and lateral ventricular volume with HDT there was no change in intracranial physiological parameters with LBNP suggesting an overall diminished response to the long-term effects of HDT. The lack of any significant increase in white matter volume during HDT with LBNP suggests maintenance of cerebral interstitial fluid transport. Final conclusions are pending the data collection for the control group (no intervention) which will be completed in 2023.

Larry A Kramer↗

Variation in Global Chemical Composition of PM2.5: Emerging Results from SPARTAN

The Surface PARTiculate mAtter Network (SPARTAN) is a long-term project that includes characterization of chemical and physical attributes of aerosols from filter samples collected worldwide. This paper discusses the ongoing efforts of SPARTAN to define and quantify major ions and trace metals found in fine particulate matter (PM (sub 2.5). Our methods infer the spatial and temporal variability of PM (sub 2.5) in a cost-effective manner. Gravimetrically weighed filters represent multi-day averages of PM (sub 2.5), with a collocated nephelometer sampling air continuously. SPARTAN instruments are paired with AErosol RObotic NETwork (AERONET) sun photometers to better understand the relationship between ground-level PM (sub 2.5) and columnar aerosol optical depth (AOD). We have examined the chemical composition of PM (sub 2.5) at 12 globally dispersed, densely populated urban locations and a site at Mammoth Cave (US) National Park used as a background comparison. So far, each SPARTAN location has been active between the years 2013 and 2016 over periods of 2-26 months, with an average period of 12 months per site. These sites have collectively gathered over 10 years of quality aerosol data. The major PM (sub 2.5) constituents across all sites (relative contribution plus or minus Standard Deviation) are ammoniated sulfate (20 percent plus or minus 11 percent), crustal material (13.4 percent plus or minus 9.9 percent), equivalent black carbon (11.9 percent plus or minus 8.4 percent), ammonium nitrate (4.7 percent plus or minus 3.0 percent), sea salt (2.3 percent plus or minus 1.6 percent), trace element oxides (1.0 percent plus or minus 1.1 percent), water (7.2 percent plus or minus 3.3 percent) at 35 percent relative humidity, and residual matter (40 percent plus or minus 24 percent). Analysis of filter samples reveals that several PM (sub 2.5) chemical components varied by more than an order of magnitude between sites. Ammoniated sulfate ranges from 1.1 microns per cubic meter (Buenos Aires, Argentina) to 17 microns per cubic meter (Kanpur, India in the dry season). Ammonium nitrate ranged from 0.2 microns per cubic meter (Mammoth Cave, in summer) to 6.8 microns per cubic meter (Kanpur, dry season). Equivalent black carbon ranged from 0.7 microns per cubic meter (Mammoth Cave) to over 8 microns per cubic meter (Dhaka, Bangladesh and Kanpur, India). Comparison of SPARTAN vs. coincident measurements from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network at Mammoth Cave yielded a high degree of consistency for daily PM (sub 2.5) (r squared equals 0.76, slope equals 1.12), daily sulfate (r squared equals 0.86, slope equals 1.03), and mean fractions of all major PM (sub 2.5) components (within 6 percent). Major ions generally agree well with previous studies at the same urban locations (e.g. sulfate fractions agree within 4 percent for 8 out of 11 collocation comparisons). Enhanced anthropogenic dust fractions in large urban areas (e.g. Singapore, Kanpur, Hanoi, and Dhaka) are apparent from high Zn to Al ratios. The expected water contribution to aerosols is calculated via the hygroscopicity parameter kappa (sub v (volume)) for each filter. Mean aggregate values ranged from 0.15 (Ilorin) to 0.28 (Rehovot). The all-site parameter mean is 0.20 plus or minus 0.04. Chemical composition and water retention in each filter measurement allows inference of hourly PM (sub 2.5) at 35 percent relative humidity by merging with nephelometer measurements. These hourly PM (sub 2.5) estimates compare favourably with a beta attenuation monitor (MetOne) at the nearby US embassy in Beijing, with a coefficient of variation r squared equals 0.67 (number equals 3167), compared to r squared equals 0.62 when v (volume) was not considered. SPARTAN continues to provide an open-access database of PM (sub 2.5) compositional filter information and hourly mass collected from a global federation of instruments.

Snider, Graydon↗

Field Test: Results of Tandem Walk Performance Following Long-Duration Spaceflight

BACKGROUND: Coordinated locomotion has proven to be challenging for many astronauts following long duration spaceflight. As NASA's vision for spaceflight points toward interplanetary travel, we must prepare for unassisted landings, where crewmembers may need to perform mission critical tasks within minutes of landing. Thus, it is vital to develop a knowledge base from which operational guidelines can be written that define when astronauts can be expected to safely perform certain tasks. Data obtained during the Field Test experiment (FT) will add important insight to this knowledge base. Specifically, we aim to develop a recovery timeline of functional sensorimotor performance during the first 24 hours and several days after landing. METHODS: FT is an ongoing study of 30 long-duration ISS crewmembers. Thus far, 9 have completed the full FT (5 U.S. Orbital Segment [USOS] astronauts and 4 Russian cosmonauts) and 4 more consented and launching within the next year. This is in addition to the eighteen crewmembers that participated in the pilot FT (11 USOS and 7 Russian crewmembers). The FT is conducted three times preflight and three times during the first 24 hours after landing. All crewmembers were tested in Kazakhstan in either the medical tent at the Soyuz landing site (~one hour post-landing), or at the airport (~four hours post-landing). The USOS crewmembers were also tested at the refueling stop (~12 hours post-landing) and at the NASA Johnson Space Center (~24 hours post-landing) and a final session 7 days post-landing. Crewmembers are instrumented with 9 inertial measurement unit sensors that measure acceleration and angular displacement (APDM's Emerald Sensors) and foot pressure-sensing insoles that measure force, acceleration, and center of pressure (Moticon GmbH, Munich, Germany) along with heart rate and blood pressure recording instrumentation. The FT consists of 12 tasks, but here we will focus on the most challenging task, the Tandem Walk, which was also performed as part of pilot FT. To perform the Tandem Walk, subjects begin with their feet together, their arms crossed at their chest and eyes closed. When ready, they brought one foot forward and touched the heel of their foot to their toe, repeating with the other foot, and continuing for about 10 steps. Three trials were collected with the eyes closed and a fourth trial was collected with eyes open. There are four metrics which are used to determine the performance level of the Tandem Walk. The first is percent correct steps. For a step to be counted as correct, the foot could not touch the ground while bringing it forward (no side stepping), eyes must stay closed during the eyes closed trials, the heel and toe should be touching, or almost touching (no large gaps) and there shouldn't be more than a three second pause between steps. Three judges score each step and the median of the three scores is kept. The second metric is the average step speed, or the number of steps/time to complete them. Thirdly, the root mean squared (RMS) error in the resultant trunk acceleration is used to determine the amount of upper body instability observed during the task. Finally, the RMS error of the mediolateral center of pressure as measured by the Moticon insoles is used to determine the mediolateral instability at the foot level. These four parameters are combined into a new overall Tandem Walk Parameter. RESULTS: Preliminary results show that crewmembers perform the Tandem Walk significantly worse the first 24 hours after landing as compared to their baseline performance. We find that each of the four performance metrics is significantly worse immediately after landing. We will present the results of tandem walk performance during the FT thus far. We will also combine these with the 18 crewmembers that participated in the pilot FT, concentrating on the level of performance and recovery rate. CONCLUSION: The Tandem Walk data collected as part of the FT experiment will provide invaluable information on the performance capabilities of astronauts during the first 24 hours after returning from long-duration spaceflight that can be used in planning future Mars, or other deep-space missions with unassisted landings. FT will determine the average sensorimotor recovery timeline and inform return-to-duty guidelines for unassisted landings.

Rosenberg, M. J. F.↗

Environmental Stress Testing of the Single Sample Cylinder: A Proven Consensus Standard for Internal Gas Analysis (IGA) or Residual Gas Analysis (RGA)

In August 2008, Schuessler Consulting was contracted by NASA GSFC in support of the NASA Electronic Parts and Packaging (NEPP) program to perform two separate studies on moisture laden air in a stainless steel cylinder that had been designed to become a consensus standard for Test Method 1018. This Test Method was originally released for hybrids under Mil. Std. 883 but was quickly utilized on other microelectronic devices under the auspice of Mil. Std. 750. The cylinder had subsequently been fabricated for the 750 community. It was back-filled with moist air and subsequently analyzed over a period of time under a previous NASA contract. It had been shown that moisture in the 4000 - 5000 ppm range could be analyzed rather precisely with a mass spectrometer, commonly referred to as a Residual Gas Analyzer (RGA). The scope of this study was to ascertain if the composition and precision varied as a function of thermal shock at sub-zero temperatures and whether there was consensus when the standard was submitted to other RGA units. It was demonstrated and published that the consensus standard would yield precise RGA data for moisture within +/- 1% when optimized for a given RGA unit. It has been subsequently shown in this study at Oneida Research Services, that sub-zero storage did not affect that precision when a well-defined protocol for the analysis was followed. The consensus standard was taken to a second facility for analysis where it was found that moisture adsorption on the transfer lines caused precision to drop to +/- 12%. The Single Sample Cylinder (SSC) is a one liter stainless steel cylinder with associated sampling valves and has considerable weight and volume. But this considerable size allows for approximately 300 gas samples of the same composition to be delivered to any RGA unit. Lastly, a smaller cylinder, approximately 75 cc, of a second consensus standard was fabricated and tested with a different mix of fixed gases where moisture was kept in the 100 ppm range. This second standard has the potential of providing 30 gaseous samples and can be readily shipped to any analytical facility that desires to generate comparison RGA data. A series of comparison residual gas analyses was performed at the Honeywell Federal Manufacturing & Technologies facility in the National Nuclear Facility Administration s plant in Kansas City to complete this project. It was shown that improvements in the precision of a given RGA unit can be done by controlling the cycle time for each analysis and increasing analysis temperatures to minimize moisture adsorption. It was also found that a "one time event" in the subzero storage of the large SSC did not effect the units ability to continuously supply precise samples of the same chemistry, however the "event" caused a permanent +8% shift in the reported value of the moisture content. Lastly, a set of SSC RGA results was plotted on a common graph with DSCC "correlation study" RGA data. The result demonstrates the ability of the SSC to remove many of the individual variances that single, individual samples introduce. The consensus standards are now in storage at Oneida Research Services, one of the DSCC certified houses that does RGA to Military Standards, where they await future studies. The analytical data and the operational parameters of the instruments used are provided in the following discussion. Limitations and suggested means for improvement of both precision and accuracy are provided.

Schuessler, Philipp WH↗

Lateral and Vertical Heterogeneity of Thorium in the Procellarum KREEP Terrane: As Reflected in the Ejecta Deposits of Post-Imbrium Craters

The Procellarum KREEP Terrane displays the highest concentrations of Th on the Moon. However, locations of elevated Th in this region appear to be random. As observed in the 5 deg per pixel equal-area Th data, and made more evident in the preliminary 2 deg data, Th is enhanced around the craters Aristillus, Aristarchus, Kepler, Mairan, the Apennine Bench formation, and the Fra Mauro region, while noticeably and unexpectedly lower in other locations (e.g., Archimedes, Copernicus, Eratosthenes, and Plato). We have examined the composition of the materials present in these regions with the goal of understanding the patchy nature to the distribution of Th and ultimately to decipher the geologic processes that have concentrated the Th. At present time, the published resolution of the Lunar Prospector Th gamma-ray data is low (5 deg per pixel), but this will soon be superceded by significantly higher-resolution data (2 deg per pixel). Even at this improved resolution, however, it is difficult to resolve the units that are the major source of Th. In an attempt to circumvent this problem, we employ the higher-resolution Clementine multispectral data for those regions mentioned above. We use the UV-VIS-derived compositional information and the spectral properties of craters, and their ejecta as drill holes through the mare-basalt surface to investigate the thickness and composition of underlying material. With this information we attempt to piece together the stratigraphy and geologic history of the Imbrium-Procellanim region. We processed the five-band multispectral data from the Clementine Mission (415, 750,900,950, and 1000nm) using ISIS software and calibration parameters developed by the USGS, Flagstaff, Arizona. Final image mosaics are in equal-area sinusoidal projection, and have a resolution of 250 m/pixel. Using the method of we produced maps of FeO and Ti02 composition. Here we examine the Th, FeO, and Ti02 composition and spectral properties of the craters discussed above and their ejecta, with the goal of describing the materials they excavate. One interpretation for the origin of the high-Th material is that subsurface KREEPy materials have been excavated by impact craters. The material excavated may be either volcanic KREEP (e,g., Apennine Bench Formation), KREEPy impact-melt breccia formed by the Imbrium impact (e.g., Fra Mauro Formation), or other KREEP-rich crustal material. Determining which type of material is responsible for the elevated Th and its extent is important to understanding the premare and possibly the prebasin stratigraphy of the Imbrium-Procellarum Region. Merging the 5 deg. Th data with the shaded relief map, we observe that the highest Th concentrations are not related to pre-Imbrium upper crustal materials. The Apennines, Alpes, and Caucasus Mountains represent the pre-Imbrian highlands material and do not express concentrations of Th, FeO, and TiO2 as high as the most Th-fich materials exposed within the Procellarum KREEP Terrane. We observe that, in general, these massifs contain 10-14 wt% FeO and 4-7 ppm Th. Determining whether the Th signal is from KREEP basalts or KREEPy impact-melt breccias cannot be done with the Clementine data because the two rock types are compositionally and mineralogically too similar (e.g., the Th-rich, mafic impact-melt breccias in the Apollo sample collection are dominated by a KREEP-basalt like component. Mapping-the distribution and sizes of craters and whether they display elevated Th concentrations or not, should reveal the depth and thickness of the KREEP-rich materials, and whether they are ubiquitous (i.e., impact-melt breccia) or more randomly distributed; this might be taken as an indicator of localized KREEP-basalt flows. Within the southeastern region of the Imbrium basin, there are two Th hot spots. The first is associated with the crater Aristillus, and the latter with the Apennine Bench Formation. Adjacent to these two hot spots are craters with a lower Th signature: Archimedes and Autolycus. We observe in the ejecta of Aristillus, a region of significantly lower FeO (10-14 wt%) relative to the surrounding mare basalt. The crater Autolycus, 50 km to the south, did not excavate similar low-FeO material. We suggest that the lower-FeO material in the ejecta of Aristillus corresponds to Th-rich material; the FeO content observed in Aristillus ejecta is comparable to that of KREEP basalt or mafic impact melt breccia (10-12 wt% FeO). We determine that this low FeO, Th-rich material is volcanic KREEP, as opposed to Imbrium impact melt, on the basis that the low-Fe material is exposed more prominently in ejecta in the northern portion of Aristillus. Our assumption is that if the layer underlying Aristillus was continuous, a more widespread and uniform low-Fe signature wouldbe observed in the ejecta deposit. Archimedes, 110 km southwest of Aristillus, impacted the northern portion of the Apennine Bench prior to the eruption of KREEP basalt. Archimedes rim material is not as enriched in Th as the Apennine Bench, and there are differences between the two in FeO concentration and in their continuum slope. Archimedes exhibits a much steeperor "redder" continuum slope than the Apennine Bench. This steepslope suggests the presence of glassy material. The glassy material is concentrated around an unnamed crater on the southern rim of Archimedes (4.5W, 28.2N) and along the northern rim of Archimedes. We suggest two possibilities, or a combination of the two, to explain the low-Th signal from Archimedes: (1) The Apennine Bench prior to KREEP basalt eruption was lower in Th (4-7 ppm, e.g., similar to the Apennine massifs) and KREEP basalts are absent in the rim of Archimedes; or (2) the glassy (possibly pyroclastic) material layering the rim of the Archimedes, dilutes any high-Th material present with low-Th material. (Additional information is contained in the original)

Gillis, J. J.↗

Towards an Enhanced Droplet Activation Scheme for Multi-Moment Bulk Microphysics Schemes

Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.

Loftus, Adrian M.↗

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↗