Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pattern formation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Weather Research and Forecasting Model Sensitivity Comparisons for Warm Season Convective Initiation

Mesoscale weather conditions can significantly affect the space launch and landing operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). During the summer months, land-sea interactions that occur across KSC and CCAFS lead to the formation of a sea breeze, which can then spawn deep convection. These convective processes often last 60 minutes or less and pose a significant challenge to the forecasters at the National Weather Service (NWS) Spaceflight Meteorology Group (SMG). The main challenge is that a "GO" forecast for thunderstorms and precipitation at the Shuttle Landing Facility is required at the 90 minute deorbit decision for End Of Mission (EOM) and at the 30 minute Return To Launch Site (RTLS) decision. Convective initiation, timing, and mode also present a forecast challenge for the NWS in Melbourne, FL (MLB). The NWS MLB issues such tactical forecast information as Terminal Aerodrome Forecasts (TAF5), Spot Forecasts for fire weather and hazardous materials incident support, and severe/hazardous weather Watches, Warnings, and Advisories. Lastly, these forecasting challenges can also affect the 45th Weather Squadron (45 WS), which provides comprehensive weather forecasts for shuttle launch, as well as ground operations, at KSC and CCAFS. The need for accurate mesoscale model forecasts to aid in their decision making is crucial. This study specifically addresses the skill of different model configurations in forecasting warm season convective initiation. Numerous factors influence the development of convection over the Florida peninsula. These factors include sea breezes, river and lake breezes, the prevailing low-level flow, and convergent flow due to convex coastlines that enhance the sea breeze. The interaction of these processes produces the warm season convective patterns seen over the Florida peninsula. However, warm season convection remains one of the most poorly forecast meteorological parameters. To determine which configuration options are best to address this specific forecast concern, the Weather Research and Forecasting (WRF) model, which has two dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM) was employed. In addition to the two dynamical cores, there are also two options for a "hot-start" initialization of the WRF model - the Local Analysis and Prediction System (LAPS; McGinley 1995) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS; Brewster 1996). Both LAPS and ADAS are 3- dimensional weather analysis systems that integrate multiple meteorological data sources into one consistent analysis over the user's domain of interest. This allows mesoscale models to benefit from the addition of highresolution data sources. Having a series of initialization options and WRF cores, as well as many options within each core, provides SMG and MLB with considerable flexibility as well as challenges. It is the goal of this study to assess the different configurations available and to determine which configuration will best predict warm season convective initiation.

Watson, Leela R.↗

Multispectral observations of the Jovian aurora

The upper atmospheres of the Earth and the outer planets form a screen on which precipitating charged particles, like the electron beam in a television, trace fleeting, but revealing patterns of visible, ultraviolet, infrared, and x ray emissions that offer valuable clues to processes occurring within the planetary magnetospheres. At Earth, years of in situ measurements, as well as ground based observations, have yielded a picture (still fuzzy) where the interaction of the solar wind with the magnetosphere of the Earth provides a complex path for the storage and release of energy during magnetic substorms; the ultimate manifestation of terrestrial auroral processes. More recent global imaging of substorm events from high above the Earth (greater than 3.5 R(sub e)) by Dynamics Explorer have made a unique contribution towards understanding the global and temporal evolution of such auroral events by providing a morphological perspective and by providing the crucial observational link that allows the separation of spatial and temporal variations inherent in the interpretation of in situ data. A similar role was played by the Hubble Space Telescope (HST) during the recent encounter of Ulysses with Jupiter in helping to define a new paradigm in Jovian auroral physics. The old paradigm portrayed Jupiter's magnetosphere as totally dominated by internal processes (i.e. Io related tori, heavy ions, etc.) where energetic heavy ion precipitation in the inner magnetosphere was solely responsible for the observed auroral phenomena. Ulysses and HST portray a more Earth-like paradigm where electron acceleration in the outer magnetosphere near the boundary with the solar wind plays a distinct role in the formation of auroral hot spots, yet energetic heavy ions also enter into the picture (similar to the role of the energetic ions from the terrestrial ring current during magnetic substorms). These heavy ions as a result of excitation during their transit through the atmosphere produce the x ray emissions observed in Roentgensatellit (ROSAT) x ray energy spectra.

Source record↗

Strategies and methodologies to develop techniques for computer-assisted analysis of gas phase formation during altitude decompression

It would be of operational significance if one possessed a device that would indicate the presence of gas phase formation in the body during hypobaric decompression. Automated analysis of Doppler gas bubble signals has been attempted for 2 decades but with generally unfavorable results, except with surgically implanted transducers. Recently, efforts have intensified with the introduction of low-cost computer programs. Current NASA work is directed towards the development of a computer-assisted method specifically targeted to EVA, and we are most interested in Spencer Grade 4. We note that Spencer Doppler Grades 1 to 3 have increased in the FFT sonogram and spectrogram in the amplitude domain, and the frequency domain is sometimes increased over that created by the normal blood flow envelope. The amplitude perturbations are of very short duration, in both systole and diastole and at random temporal positions. Grade 4 is characteristic in the amplitude domain but with modest increases in the FFT sonogram and spectral frequency power from 2K to 4K over all of the cardiac cycle. Heart valve motion appears to characteristic display signals: (1) the demodulated Doppler signal amplitude is considerably above the Doppler-shifted blow flow signal (even Grade 4); and (2) demodulated Doppler frequency shifts are considerably greater (often several kHz) than the upper edge of the blood flow envelope. Knowledge of these facts will aid in the construction of a real-time, computer-assisted discriminator to eliminate cardiac motion artifacts. There could also exist perturbations in the following: (1) modifications of the pattern of blood flow in accordance with Poiseuille's Law, (2) flow changes with a change in the Reynolds number, (3) an increase in the pulsatility index, and/or (4) diminished diastolic flow or 'runoff.' Doppler ultrasound devices have been constructed with a three-transducer array and a pulsed frequency generator.

Powell, Michael R.↗

The Large Impact Process Inferred from the Geology of Lunar Multiring Basins

The study of the geology of multiring impact basins on the Moon over the past ten years has given us a rudimentary understanding of how these large structures have formed and evolved on the Moon and other bodies. Two-ring basins on the Moon begin to form at diameters of about 300 km; the transition diameter at which more than two rings appear is uncertain, but it appears to be between 400 and 500 km in diameter. Inner rings tend to be made up of clusters or aligned segments of massifs and are arranged into a crudely concentric pattern; scarp-like elements may or may not be present. Outer rings are much more scarp-like and massifs are rare to absent. Basins display textured deposits, interpreted as ejecta, extending roughly an apparent basin radius exterior to the main topographic rim. Ejecta may have various morphologies, ranging from wormy and hummocky deposits to knobby surfaces; the causes of these variations are not known, but may be related to the energy regime in which the ejecta are deposited. Outside the limits of the textured ejecta are found both fields of satellitic craters (secondaries) and light plains deposits. Impact melt sheets are observed on the floors of relatively unflooded basins. Samples of impact melts from lunar basins have basaltic major-element chemistry, characterized by K, rare-earth elements (REE), P, and other trace elements of varying concentration (KREEP); ages are between 3.8 and 3.9 Ga. These lithologies cannot be produced through the fusion of known pristine (plutonic) rock types, suggesting the occurrence of unknown lithologies within the Moon. These melts were probably generated at middle to lower crustal levels. Ejecta compositions, preservation of pre-basin topography, and deposit morphologies all indicate that the excavation cavity of multiring basins is between about 0.4 and 0.6 times the diameter of the apparent crater diameter. Basin depths of excavation can be inferred from the composition of basin ejecta. A variety of mechanisms has been proposed to account for the formation of basin rings but none of them are entirely plausible. Mechanisms can be divided into two broad groups: (1) forcible uplift due to fluidization of the target; (2) concentric, brittle, fracturing and failure of the target, on regional (megaterraces) to global scales (lithospheric fracturing). Most basin rings are spaced at a constant factor on all planets. Evidence supports divergent ringforming models, so it may be that the ring-locating mechanism differs from the ring-forming mechanism. Thus, large-scale crustal foundering (megaterracing) could occur along concentric zones of weakness created by some type of resonant wave mechanism (fluidization and uplift); such immediate crustal adjustment could then be followed by long-term adjustment of the fractured lithosphere.

Spudis, Paul D.↗

Dissolution Rates of Allophane with Variable Fe Contents: Implications for Aqueous Alteration and the Preservation of X-ray Amorphous Materials on Mars

Recent measurements from Mars document X-ray amorphous/nano-crystalline materials in multiple locations across the planet. However, despite their prevalence, little is known about these materials or what their presence implies for the history of Mars. The X-ray amorphous component of the martian soil in Gale crater has an X-ray diffraction pattern that can be partially fit with allophane (approximately Al2O3•(SiO2)1.3-2•(H2O)2.5-3), and the low-temperature water release data are consistent with allophane. The chemical data from Gale crater suggest that other silicate materials similar to allophane, such as Fe-substituted allophane (approximately (Fe2O3)0.01-0.5(Al2O3)0.5-0.99•(SiO2)2•3H2O), may also be present. In order to investigate the properties of these potential poorly crystalline components of the martian soil, Fe-free allophane (Fe:Al = 0), Fe-poor allophane (Fe:Al = 1:99), and Fe-rich allophane (Fe:Al = 1:1) were synthesized and then characterized using electron microscopy and Mars-relevant techniques, including infrared spectroscopy, X-ray diffraction, and evolved gas analysis. Dissolution experiments were performed at acidic (initial pH values pH0 = 3.01, pH0 = 5.04), near-neutral (pH0 = 6.99), and alkaline (pH0 = 10.4) conditions in order to determine dissolution kinetics and alteration phases for these poorly crystalline materials. Dissolution rates (rdiss), based on the rate of Si release into solution, show that these poorly crystalline materials dissolve approximately an order of magnitude faster than crystalline phases with similar compositions at all pH conditions. For Fe-free allophane, logrdiss = -10.65 – 0.15 × pH; for Fe-poor allophane, logrdiss = -10.35 – 0.22 × pH; and for Fe-rich allophane, logrdiss = -11.46 – 0.042 × pH at 25°C, where rdiss has the units of mol m-2 s-1. The formation of incipient phyllosilicate-like phases was detected in Fe-free and Fe-rich allophane reacted in aqueous solutions with pH0 = 10.4 (steady-state pH ≈ 8). Mars-analog instrument analyses demonstrate that Fe-free allophane, Fe-poor allophane, and Fe-rich allophane are appropriate analogs for silicate phases in the martian amorphous soil component. Therefore, similar materials on Mars must have had limited interaction with liquid water since their formation. Combined with chemical changes expected from weathering, such as phyllosilicate formation, the rapid alteration of these poorly crystalline materials may be a useful tool for evaluating the extent of aqueous alteration in returned samples of martian soils.

S J Ralston↗

Micro and Macro Segregation in Alloys Solidifying with Equiaxed Morphology

To understand macro segregation formation in Al-Cu alloys, experiments were run under terrestrial gravity (1g) and under low gravity during parabolic flights (10(exp -2) g). Alloys of two different compositions (2% and 5% Cu) were solidified at two different cooling rates. Systematic microscopic and SEM observations produced microstructural and segregation maps for all samples. These maps may be used as benchmark experiments for validation of microstructure evolution and segregation models. As expected, the macro segregation maps are very complex. When segregation was measured along the central axis of the sample, the highest macro segregation for samples solidified at 1g was obtained for the lowest cooling rate. This behavior is attributed to the longer time available for natural convection and shrinkage flow to affect solute redistribution. In samples solidified under low-g, the highest macro-segregation was obtained at the highest cooling rate. In general, low-gravity solidification resulted in less segregation. To explain the experimental findings, an analytical (Flemings-Nereo) and a numerical model were used. For the numerical model, the continuum formulation was employed to describe the macroscopic transports of mass, energy, and momentum, associated with the microscopic transport phenomena, for a two-phase system. The model proposed considers that liquid flow is driven by thermal and solutal buoyancy, and by solidification shrinkage. The Flemings-Nereo model explains well macro segregation in the initial stages of low-gravity segregation. The numerical model can describe the complex macro segregation pattern and the differences between low- and high-gravity solidification.

Stefanescu, Doru M.↗

Improved Rainfall Data in the Philippines through Concurrent Use of GPM IMERG and Ground-Based Measurements

The availability of accurate and reliable rainfall data that are applicable to various phenomenological, climatological, and modeling studies is important, especially in the Philippines, which is considered to be highly vulnerable to natural hazards and a changing climate. The presented strategy involved constructing a dataset consisting of synoptic data, automatic rain gauge (ARG) measurements, and satellite data that are co-registered, consistent, and formatted in the same manner. Although sparse in number, the synoptic stations provide the most accurate rainfall information and were used as the baseline for creating the dataset. The ARGs that are within a distance of 1 km to the synoptic stations were used to determine the correction factors needed to make the synoptic and ARG data consistent. Subsequently, the corrected ARGs were used to make the satellite IMERG data consistent with both ARG and synoptic data. In case of the latter, only IMERG pixels with at least 10 ARGs within the relatively large footprint of the satellite sensor were used in estimating the required correction parameters derived from a combination of a power transform and linear regression correction techniques. The final results show good agreement of synoptic and corrected ARG data with correlation coefficients of 0.94 and 0.97 for the 10 day and monthly data, respectively, and improvement in the linear regression slope from 0.67 to 0.90 for 10 day data, and 0.70 to 0.94 for monthly data. In addition, the corrected ARG data agree well with the corrected IMERG data, with correlation coefficients of 0.88 and 0.93 for the 10 day and monthly data, respectively, and an improvement in slope from 0.66 to 0.87 for 10 day data, and 0.74 to 0.99 for monthly data. The merit of using a combined dataset is illustrated through comparative analyses of the IMERG data and spatially interpolated synoptic and ARG data. The results show general agreements in spatial patterns of rainfall across the datasets, especially in areas where in situ measurements are recorded. The observed discrepancy when ground data is limited emphasizes the need for satellite IMERG data to obtain the true spatial patterns of rainfall distribution.

rainfall↗

NASA Tech Briefs, September 2006

Topics covered include: Improving Thermomechanical Properties of SiC/SiC Composites; Aerogel/Particle Composites for Thermoelectric Devices; Patches for Repairing Ceramics and Ceramic- Matrix Composites; Lower-Conductivity Ceramic Materials for Thermal-Barrier Coatings; An Alternative for Emergency Preemption of Traffic Lights; Vehicle Transponder for Preemption of Traffic Lights; Automated Announcements of Approaching Emergency Vehicles; Intersection Monitor for Traffic-Light-Preemption System; Full-Duplex Digital Communication on a Single Laser Beam; Stabilizing Microwave Frequency of a Photonic Oscillator; Microwave Oscillators Based on Nonlinear WGM Resonators; Pointing Reference Scheme for Free-Space Optical Communications Systems; High-Level Performance Modeling of SAR Systems; Spectral Analysis Tool 6.2 for Windows; Multi-Platform Avionics Simulator; Silicon-Based Optical Modulator with Ferroelectric Layer; Multiplexing Transducers Based on Tunnel-Diode Oscillators; Scheduling with Automated Resolution of Conflicts; Symbolic Constraint Maintenance Grid; Discerning Trends in Performance Across Multiple Events; Magnetic Field Solver; Computing for Aiming a Spaceborne Bistatic- Radar Transmitter; 4-Vinyl-1,3-Dioxolane-2-One as an Additive for Li-Ion Cells; Probabilistic Prediction of Lifetimes of Ceramic Parts; STRANAL-PMC Version 2.0; Micromechanics and Piezo Enhancements of HyperSizer; Single-Phase Rare-Earth Oxide/Aluminum Oxide Glasses; Tilt/Tip/Piston Manipulator with Base-Mounted Actuators; Measurement of Model Noise in a Hard-Wall Wind Tunnel; Loci-STREAM Version 0.9; The Synergistic Engineering Environment; Reconfigurable Software for Controlling Formation Flying; More About the Tetrahedral Unstructured Software System; Computing Flows Using Chimera and Unstructured Grids; Avoiding Obstructions in Aiming a High-Gain Antenna; Analyzing Aeroelastic Stability of a Tilt-Rotor Aircraft; Tracking Positions and Attitudes of Mars Rovers; Stochastic Evolutionary Algorithms for Planning Robot Paths; Compressible Flow Toolbox; Rapid Aeroelastic Analysis of Blade Flutter in Turbomachines; General Flow-Solver Code for Turbomachinery Applications; Code for Multiblock CFD and Heat-Transfer Computations; Rotating-Pump Design Code; Covering a Crucible with Metal Containing Channels; Repairing Fractured Bones by Use of Bioabsorbable Composites; Kalman Filter for Calibrating a Telescope Focal Plane; Electronic Absolute Cartesian Autocollimator; Fiber-Optic Gratings for Lidar Measurements of Water Vapor; Simulating Responses of Gravitational-Wave Instrumentation; SOFTC: A Software Correlator for VLBI; Progress in Computational Simulation of Earthquakes; Database of Properties of Meteors; Computing Spacecraft Solar-Cell Damage by Charged Particles; Thermal Model of a Current-Carrying Wire in a Vacuum; Program for Analyzing Flows in a Complex Network; Program Predicts Performance of Optical Parametric Oscillators; Processing TES Level-1B Data; Automated Camera Calibration; Tracking the Martian CO2 Polar Ice Caps in Infrared Images; Processing TES Level-2 Data; SmaggIce Version 1.8; Solving the Swath Segment Selection Problem; The Spatial Standard Observer; Less-Complex Method of Classifying MPSK; Improvement in Recursive Hierarchical Segmentation of Data; Using Heaps in Recursive Hierarchical Segmentation of Data; Tool for Statistical Analysis and Display of Landing Sites; Automated Assignment of Proposals to Reviewers; Array-Pattern-Match Compiler for Opportunistic Data Analysis; Pre-Processor for Compression of Multispectral Image Data; Compressing Image Data While Limiting the Effects of Data Losses; Flight Operations Analysis Tool; Improvement in Visual Target Tracking for a Mobile Robot; Software for Simulating Air Traffic; Automated Vectorization of Decision-Based Algorithms; Grayscale Optical Correlator Workbench; "One-Stop Shopping" for Ocean Remote-Sensing and Model Data; State Analysis Database Tool; Generating CAHV and CAHVOmages with Shadows in ROAMS; Improving UDP/IP Transmission Without Increasing Congestion; FORTRAN Versions of Reformulated HFGMC Codes; Program for Editing Spacecraft Command Sequences; Flight-Tested Prototype of BEAM Software; Mission Scenario Development Workbench; Marsviewer; Tool for Analysis and Reduction of Scientific Data; ASPEN Version 3.0; Secure Display of Space-Exploration Images; Digital Front End for Wide-Band VLBI Science Receiver; Multifunctional Tanks for Spacecraft; Lightweight, Segmented, Mostly Silicon Telescope Mirror; Assistant for Analyzing Tropical-Rain-Mapping Radar Data; and Anion-Intercalating Cathodes for High-Energy- Density Cells.

Source record↗

A Layered Debris Disk around M Star TWA 7 in Scattered Light

We have obtained Hubble Space Telescope (HST) coronagraphic observations of the circumstellar disk around M star TWA 7 using the Space Telescope Imaging Spectrograph (STIS) instrument in visible light. Together with archival observations, including HST/NICMOS using the F160W filter and Very Large Telescope/SPHERE at the H-band in polarized light, we investigate the system in scattered light. By studying this nearly face-on system using geometric disk models and Henyey–Greenstein phase functions, we report a new discovery of a tertiary ring and a clump. We identify a layered architecture: three rings, a spiral, and an ≈150 au 2 elliptical clump. The most extended ring peaks at 28 au, and the other components are on its outskirts. Our point-source detection-limit calculations demonstrate the necessity of disk modeling in imaging fainter planets. Morphologically, we witness a clockwise spiral motion, and the motion pattern is consistent with both solid body motion and local Keplerian motion; we also observe underdensity regions for the secondary ring that might result from mean-motion resonance or moving shadows: both call for re-observations to determine their nature. Comparing multi-instrument observations, we obtain blue STIS-NICMOS color, a STIS-SPHERE radial distribution peak difference for the tertiary ring, and a high SPHERE-NICMOS polarization fraction; these aspects indicate that TWA 7 could retain small dust particles. By viewing the debris disk around M star TWA 7 at a nearly face-on vantage point, our study allows for the understanding of such disks in scattered light in both system architecture and dust property.

Debris disks↗

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↗

In the Pursuit of Regolithic Howardites

The HED (Howardite, Eucrite and Diogenite) meteorite clan likely originate from the asteroid 4-Vesta [1]. Howardites (polymict breccias of eucritic and diogenitic material) are believed to originate from the vestan surface, and many contain regolith-like features (impact and/or melt clasts, fragmental breccia clasts, carbonaceous chondrite fragments), which may relate to regolith-formation processes. Noble gas analysis can help determine true regolithic nature, as Solar Wind (SW) is im-planted into grains at the upper-surfaces of solar system bodies lacking an atmosphere or magnetic field. Howardites from Ves-ta s true regolith would thus show evidence for SW components. Thus far, we have identified 5 regolithic howardites: LEW 85313; MET 00423; PRA 04401; SCO 06040 and EET 87513; based on our noble gas analyses, with the latter 3 showing some evidence for a planetary(+SW) component, likely related to CM-like material present in the assemblage [2-5]. However, we did not find a good correlation between SW content and other petro-logic regolithic features. Nor did we find an obvious correlation between SW and high siderophile element contents (Ni greater than 300 micrograms/gram), or an Al2O3 range of 8-9 wt% and eucrite/diogenite (E/D) ratio of 2:1 as suggested by [6] to be further regolith indicators. Here, we report our latest noble gas data for two howardites GRO 95535, GRO 95602 and a polymict eucrite EET 87518 in continuing research aimed at better understanding the vestan regolith. Results: Noble gas analysis was performed on an MAP 215-50 noble gas mass spectrometer using furnace step-heating. Our results, shown in Table 1, are compared with SW- (LEW 85313) and planetary-dominated (PRA 04401 ~60% CM) howardites [2]. EET 87518 is dominated by cosmogenic components. By comparison, both howardites show strong evidence for SW, with total Ne-20/Ne-22 ~8.7-8.8 (SW: Ne-20/Ne-22 13.78 [7]), and identical release patterns to our other CM-poor SW-rich samples. This suggests that these samples are from the vestan regolith. As they have lower Ni contents than suggested by [6], this further illustrates that these parameters may show some bias [2].

Cartwright, J. A.↗

Decision Aid for Conjunction Risk Mitigation by Differential Drag

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

conjunction assessment↗

Spatially Resolved, In Situ Carbon Isotope Analysis of Archean Organic Matter

Spatiotemporal variability in the carbon isotope composition of sedimentary organic matter (OM) preserves information about the evolution of the biosphere and of the exogenic carbon cycle as a whole. Primary compositions, and imprints of the post-depositional processes that obscure them, exist at the scale of individual sedimentary grains (mm to micron). Secondary ion mass spectrometry (SIMS) (1) enables analysis at these scales and in petrographic context, (2) permits morphological and compositional characterization of the analyte and associated minerals prior to isotopic analysis, and (3) reveals patterns of variability homogenized by bulk techniques. Here we present new methods for in situ organic carbon isotope analysis with sub-permil precision and spatial resolution to 1 micron using SIMS, as well as new data acquired from a suite of Archean rocks. Three analytical protocols were developed for the CAMECA ims1280 at WiscSIMS to analyze domains of varying size and carbon concentration. Average reproducibility (at 2SD) using a 6 micron spot size with two Faraday cup detectors was 0.4 %, and 0.8 % for analyses using 1 micron and 3 micron spot sizes with a Faraday cup (for C-12) and an electron multiplier (for C-13). Eight coals, two ambers, a shungite, and a graphite were evaluated for micron-scale isotopic heterogeneity, and LCNN anthracite (delta C-13 = -23.56 +/- 0.1 %, 2SD) was chosen as the working standard. Correlation between instrumental bias and H/C was observed and calibrated for each analytical session using organic materials with H/C between 0.1 and 1.5 (atomic), allowing a correction based upon a C-13H/C-13 measurement included in every analysis. Matrix effects of variable C/SiO2 were evaluated by measuring mm to sub-micron graphite domains in quartzite from Bogala mine, Sri Lanka. Apparent instrumental bias and C-12 count rate are correlated in this case, but this may be related to a crystal orientation effect in graphite. Analyses of amorphous Archean OM suggest that instrumental bias is consistent for 12C count rates as low as 10% relative to anthracite. Samples from the ABDP-9 (n=3; Mount McRae Shale, approximately 2.5 Ga), RHDH2a (n=2; Carrawine Dolomite and Jeerinah Fm, approximately 2.6 Ga), WRL1 (n=3; Wittenoom Fm, Marra Mamba Iron Formation, and Jeerinah Fm, approximately 2.6 Ga), and SV1 (n=1; Tumbiana Fm, approximately 2.7 Ga) drill cores, each previously analyzed for bulk organic carbon isotope composition, yielded 100 new, in situ data from Neoarchean sedimentary OM. In these samples, delta C-13 varies between -53.1 and -28.3 % and offsets between in situ and bulk compositions range from -8.3 to 18.8%. In some cases, isotopic composition and mode of occurrence (e.g. morphology and mineral associations) are statistically correlated, enabling the identification of distinct reservoirs of OM. Our results support previous evidence for gradients of oxidation with depth in Neoarchean environments driven by photosynthesis and methane metabolism. The relevance of these findings to questions of bio- and syngenicity as well as the alteration history of previously reported Archean OM will be discussed.

Williford, Kenneth H.↗

Formation and evolution of radial fracture systems on Venus

A survey of approximately 90 percent of the surface of Venus using Magellan data has been carried out to locate all radial fracture systems and to assess their association with other features such as volcanic edifices and coronae. Squyres et al. and Stofan et al. have discussed the association of radial fracture features in relation to coronae features, our approach was to assess the associations of all of the fracture systems. These fracture systems have two broad types of form - some fracture systems are associated with updomed topography, radiate from a point and have relatively uniform fracture lengths while others have a wider range of fracture lengths and radiate from the outer edge of a central caldera. Squyres et al. and Stofan et al. have interpreted both types of feature as reflecting tectonic fracturing resulting from uplift of the surface as a mantle plume impinges upon the crust. While it is true that a number of features are related to uplift and that such uplift will induce stresses consistent with radial fracturing, we explore the possibility that these fractures are not exclusively of tectonic origin. Purely tectonic fracturing will tend to generate a few main fractures/faults along which most of the stresses due to uplift will be accommodated leading to the triple-junction form common for terrestrial updoming. Though this type of feature is observed on Venus (e.g., feature located at 34S86), the majority of radial fracture systems display much more intensive fracturing than this through a full 360 degrees; this is difficult to explain by purely tectonic processes. The association of many of the fractures with radial lava flows leads us to interpret these fractures as reflecting dike emplacement: the form of the fractures being consistent with primarily vertical propagation from the head of a mantle plume. In the case of the second type of fracture system (those radiating from a central caldera), an even stronger case can be made that the fractures are not of tectonic origin. These features are not as commonly associated with updoming of the surface and where they are, the fractures extend out well beyond the edge of the topographic rise - an observation which is not consistent with the fractures being of tectonic uplift origin. Furthermore the fractures have a distribution of lengths (many short, fewer long) which is characteristic of dike swarms, and show direct associations with calderas and lava flows consistent with a volcanic origin. In addition, the longest fractures have a radial pattern only close to the center of the system but bend with distance to align themselves with the regional stress field - this behavior is very difficult to explain on purely tectonic grounds but is a pattern commonly seen for terrestrial dikes. For these reasons, we argue that many, if not the majority, of radial fracture systems found on Venus are the surface reflection of dike swarms, those associated with positive topography reflecting vertical emplacement and those radiating from calderas reflecting lateral propagation.

Parfitt, E. A.↗

Tool for Forecasting Cool-Season Peak Winds Across Kennedy Space Center and Cape Canaveral Air Force Station

The expected peak wind speed for the day is an important element in the daily morning forecast for ground and space launch operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45th Weather Squadron (45 WS) must issue forecast advisories for KSC/CCAFS when they expect peak gusts for >= 25, >= 35, and >= 50 kt thresholds at any level from the surface to 300 ft. In Phase I of this task, the 45 WS tasked the Applied Meteorology Unit (AMU) to develop a cool-season (October - April) tool to help forecast the non-convective peak wind from the surface to 300 ft at KSC/CCAFS. During the warm season, these wind speeds are rarely exceeded except during convective winds or under the influence of tropical cyclones, for which other techniques are already in use. The tool used single and multiple linear regression equations to predict the peak wind from the morning sounding. The forecaster manually entered several observed sounding parameters into a Microsoft Excel graphical user interface (GUI), and then the tool displayed the forecast peak wind speed, average wind speed at the time of the peak wind, the timing of the peak wind and the probability the peak wind will meet or exceed 35, 50 and 60 kt. The 45 WS customers later dropped the requirement for >= 60 kt wind warnings. During Phase II of this task, the AMU expanded the period of record (POR) by six years to increase the number of observations used to create the forecast equations. A large number of possible predictors were evaluated from archived soundings, including inversion depth and strength, low-level wind shear, mixing height, temperature lapse rate and winds from the surface to 3000 ft. Each day in the POR was stratified in a number of ways, such as by low-level wind direction, synoptic weather pattern, precipitation and Bulk Richardson number. The most accurate Phase II equations were then selected for an independent verification. The Phase I and II forecast methods were compared using an independent verification data set. The two methods were compared to climatology, wind warnings and advisories issued by the 45 WS, and North American Mesoscale (NAM) model (MesoNAM) forecast winds. The performance of the Phase I and II methods were similar with respect to mean absolute error. Since the Phase I data were not stratified by precipitation, this method's peak wind forecasts had a large negative bias on days with precipitation and a small positive bias on days with no precipitation. Overall, the climatology methods performed the worst while the MesoNAM performed the best. Since the MesoNAM winds were the most accurate in the comparison, the final version of the tool was based on the MesoNAM winds. The probability the peak wind will meet or exceed the warning thresholds were based on the one standard deviation error bars from the linear regression. For example, the linear regression might forecast the most likely peak speed to be 35 kt and the error bars used to calculate that the probability of >= 25 kt = 76%, the probability of >= 35 kt = 50%, and the probability of >= 50 kt = 19%. The authors have not seen this application of linear regression error bars in any other meteorological applications. Although probability forecast tools should usually be developed with logistic regression, this technique could be easily generalized to any linear regression forecast tool to estimate the probability of exceeding any desired threshold . This could be useful for previously developed linear regression forecast tools or new forecast applications where statistical analysis software to perform logistic regression is not available. The tool was delivered in two formats - a Microsoft Excel GUI and a Tool Command Language/Tool Kit (Tcl/Tk) GUI in the Meteorological Interactive Data Display System (MIDDS). The Microsoft Excel GUI reads a MesoNAM text file containing hourly forecasts from 0 to 84 hours, from one model run (00 or 12 UTC). The GUI then displays e peak wind speed, average wind speed, and the probability the peak wind will meet or exceed the 25-, 35- and 50-kt thresholds. The user can display the Day-1 through Day-3 peak wind forecasts, and separate forecasts are made for precipitation and non-precipitation days. The MIDDS GUI uses data from the NAM and Global Forecast System (GFS), instead of the MesoNAM. It can display Day-1 and Day-2 forecasts using NAM data, and Day-1 through Day-5 forecasts using GFS data. The timing of the peak wind is not displayed, since the independent verification showed that none of the forecast methods performed significantly better than climatology. The forecaster should use the climatological timing of the peak wind (2248 UTC) as a first guess and then adjust it based on the movement of weather features.

Barrett, Joe H., III↗

Predictions of a population of cataclysmic variables in globular clusters

We have studied the number of cataclysmic variables (CVs) that should be active in globular clusters during the present epoch as a result of binary formation via two-body tidal capture. We predict the orbital period and luminosity distributions of CVs in globular clusters. The results arebased on Monte Carlo simulations combined with evolution calculations appropriate to each system formed during the lifetime of two specific globular clusters, omega Cen and 47 Tuc. From our study of these two clusters, which represent the range of core densities and states of mass segregation that are likely to be interesting, we extrapolate our results to the Galactic globlular cluster system. Although there is at present little direct observational evidence of CVs in globular clusters, we find that there should be a large number of active systems. We predict that there should be more than approximately 100 CVs in both 47 Tuc and omega Cen and several thousand in the Galactic globular cluster system. These numbers are based on two-body processes alone and represent a lower bound on the number of systems that may have been formed as a result of stellar interaction within globular clusters. The relation between these calculations and the paucity of optically detected CVs in globular clusters is discussed. Should future observations fail to find convincing evidence of a substantial population of cluster CVs, then the two-body tidal capture scenario is likely to be seriously constrained. Of the CVs we espect in 47 Tuc and omega Cen, approximately 45 and 20, respectively, should have accretion luminosities above 10(exp 33) ergs/s. If one utilizes a relation for converting accretion luminosity to hard X-ray luminosity that is based on observations of Galactic plane CVs, even these sources will not exhibit X-ray luminosities above 10(exp 33) ergs/s. While we cannot account directly for the most luminous subset of the low-luminosity globular cluster X-ray sources without assuming an evolutionary pattern that is different from that of the majority of CVs in the disk, we are able to account for all of the observed lower luminosity subset of these sources, many of which have been recently discovered through ROSAT observations. In order for our predicted integrated cluster X-ray luminosities to be consistent with observational upper limits, the relation between accretion and X-ray luminosities should be something like that inferred from the Galactic plane population of CVs. Our calculations predict a large number of systems with L(sub acc) is less than 10(exp 32) ergs/s. Although our calculations imply that globular clusters should have an enhancement of CVs relative to the number thought to be present in the Galactic disk, this enhancement is at most roughly an order of magnitude, not comparable to the factor of approximately 100 for low-mass X-ray binaries (LMXBs).

Di Stefano, R.↗

Tessera Terrain on Venus: A Survey of the Global Distribution, Characteristics, and Relation to Surrounding Units from Magellan Data

The tessera terrain on Venus, comprised of areas of high radar backscatter, complex deformation patterns relative to other units, and topography standing higher than surrounding plains, covers approximately 35.33 x 10(exp 6)sq km, about 8% of the surface of Venus, and is nonrandomly distributed, being preferentially located at equatorial and higher northern latitudes with a distinct paucity below about 30degS. Individual tessera occurrences range in area from the lower limits of our measurements (about 200 sq km) up to the largest tessera, Ovda, with an area of about 8.6 x 10(exp 6)sq km, or about 2% of the surface area of Venus. The size-frequency distribution of tessera patches is strongly unimodal and skewed toward smaller sizes, reflecting the great abundance of small tessera fragments. Modes of occurrence include (1) large clusters (e.g., Aphrodite Terra and Ishtar Terra); (2) arc-like segments which may extend for thousands of kilometers and are either concave inward toward the major tessera cluster development or away from it; (3) areas where tesserae are rare or absent which occur both as low-lying plains (e.g., Guinevere Planitia), and as elevated regions (e.g., Ata Regio). Tessera terrain has a bimodal elevation-frequency distribution, with the main peak at about 0- 1 km and an additional peak at about 3 km above mean planetary radius. In terms of number of occurrences, however, tesserae do not display a correlation with elevation at the global scale, since small tessera patches commonly occupy low-lying regions. Although tessera exhibit a range of gravity signatures, many occurrences are interpreted to represent relatively shallow (crustal) levels of compensation. Tessera boundaries include Type I (sinuous/embayed, dominated by adjacent lava plains embaying tessera massifs; 73% of total tesserae boundaries) and Type II (linear/tectonic). Only a small percentage of the length of all boundary types show no lava embayment and could be interpreted as tectonically active for long periods subsequent to initial tessera formation. Occurrence of broad slopes of post-tessem embayed plains away from tessera boundaries suggests that regional tilting occurred subsequent to final temra deformation in some places.

Ivanov, Mikhail↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗