Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Saturation problem”

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

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

At least 91 records · Page 5

Control strategies for systems with limited actuators

This work investigates the effects of actuator saturation in multi-input, multi-output (MIMO) control systems. The adverse system behavior introduced by the saturation nonlinearity is viewed here as resulting from two mechanisms: controller windup - a problem caused by the discrepancy between the limited actuator commands and the corresponding control signals, and directionality - the problem of how to use nonlimited actuators when a limited condition exists. The tracking mode and Hanus methods are two common strategies for dealing with the windup problem. It is seen that while these methods alleviate windup, performance problems remain due to plant directionality. Though high gain conventional antiwindup as well as more general linear methods have the potential to address both windup and directionality, no systematic design method for these schemes has emerged; most approaches used in practice are application driven. An alternative method of addressing the directionality problem is presented which involves the introduction of a control direction preserving nonlinearity to the Hanus antiwindup system. A nonlinearity is subsequently proposed which reduces the conservation inherent in the former direction-preserving approach, improving performance. The concept of multivariable sensitivity is seen to play a key role in the success of the new method.

Marcopoli, Vincent R.↗

Missile impact craters (White Sands Missile Range, New Mexico) and applications to lunar research: Contributions to astrogeology

Craters in natural materials at White Sands Missile Range, N. Mex., were produced by the impact of high-velocity to hypervelocity missiles traveling along oblique trajectories with kinetic energies between 2.1 and 81 × 1014 ergs. The oblique impacts produce craters 2 to 10 m across with morphologies and ejecta that are bilaterally symmetrical with respect to the plane of the missile trajectory. Rims are high and the amount of ejecta large in down-trajectory and lateral directions, whereas rims are low to nonexistent and ejecta thin to absent up-trajectory. Symmetry development and modifications of the symmetry are a function of target material, local topography, and angle of impact. Seven mappable units can be recognized in and around the craters. Three of these are ejecta: thick ejecta near the crater, thin to discontinuous ejecta at greater distances, and scattered ejecta at the greatest distances to the limit of throwout. These ejecta units may be absent on the up-trajectory side; if present, they are rarely as thick or continuous as on other sides of the crater. Three units are target materials: undeformed target material exposed in local patches through thin to discontinuous ejecta and everywhere between the fragments of scattered ejecta, tilted and broken target material exposed in upper crater walls, and shattered and fractured target material exposed on the up-trajectory crater wall. The seventh unit is slope material composed of talus and fallback within the crater. Development, character, and exposure of these units varies chiefly with the target material. Ejecta from the craters is chiefly broken but relatively undeformed target material that may range in size from very fine grained debris to large blocks. Where the target is porous, significant amounts of the ejecta are composed of sheared and compressed fragments, some coated with dark layers of mixed projectile pieces, powder, and fused metal mixed with crushed target material. For layered targets, the original stratigraphic sequence is crudely preserved and in inverted order in thick ejecta. Secondary impact craters are produced by the impact of ejected fragments when the surrounding surface materials are sufficiently weak. A wide variety of secondary impact crater relations may result. Secondary craters nearest the primary crater have blocks in them that are larger than or the same size as the crater they produced. Farther from the primary crater, the fragments are generally smaller than the secondary crater and are ejected from it. Excavation of four craters revealed a mixed breccia beneath the crater floor composed of missile pieces, sheared and compressed target material, and crushed debris. Banded, disaggregated target material and nonmixed breccia surrounded the mixed breccia, and these breccias were surrounded by a zone of conjugate fractures. Beneath the ejecta on the lateral and down-trajectory crater flanks, the target materials were tilted upward and broken. Up-trajectory, open fractures and downward displacement occurred in two of the craters. No displacement was observed for the other two. Beneath the down-trajectory rims of craters with distinct layering, overturned synclines were observed. Missile breakup and behavior during cratering are a function of target and missile properties. Missile breakup depends on missile velocity and is most extensive at high velocities, where the missile is fragmented, powdered, and partly fused. Burial of missile or its fragmented, powdered, and fused remains is greatest for porous targets and least for dense cohesive targets. For very porous targets, camouflet structures containing the fragmented missile may form. Least squares fit to the data on craters in dry to moist targets indicate V(a) = 10^(-11.433)E(p)^(1.205) where V(a) is the volume of the apparent crater and E(p) is the kinetic energy of the missile. This equation is consistent with expectations of the equations relating apparent depth and radius to kinetic energy. Extrapolation of displaced masses and kinetic energies for laboratory impacts with sand and rock converge near 10^(15) to 10^(16) ergs, where the extrapolations are near the data on missile impact craters, corrected for impact angle. Displaced masses of craters produced by missile impacts and by chemical explosives with small scaled depths of burial are about the same when the kinetic energies of the missiles (corrected for angle of impact) are equal to the TNT equivalent energy of the explosive. The problem of equivalent scaled depth of burst for an impact crater is complicated and not entirely resolved, however. Both missile impact craters and chemical explosive craters in water-saturated targets are larger than their counterparts in dry to moist materials. Data collected during the study of missile impact craters have helped resolve a number of problems in lunar research: (1) the soillike nature of lunar surface materials was predicted, (2) sizes of craters produced by artificial impacts were correctly predicted, (3) certain features imaged by Surveyor were found to be analogous to features associated with missile impact craters, {4) missile impacts were used in support of the Apollo passive seismic experiment, (5) craters seen in Apollo orbital photographs were found to be similar to some missile impact craters, (6) missile impact craters supplied data on sample collection and crater phenomenology used in training astronauts, and (7) some returned lunar samples are similar to coated, sheared, and compressed fragments ejected from missile impact craters.

H. J. Moore↗

Reduction of plasma electron density in a gas ionized by an electron beam - Use of a gaseous dielectric

Propagation of an electron beam through a gas creates a secondary electron/ion plasma which can have subsequent deleterious effects on the propagation of the beam. In the case of pulsed electron beams with short micropulse durations, these effects can be greatly reduced through the use of a small doping fraction of an electron attachment gas. We present a model which allows the calculation of the reduction in unbound plasma electron density attainable with a gaseous dielectric dopant. Potential problems with a dopant, including increased ionization, increased scattering, altered refractive index, and dopant saturation and fragmentation, are discussed.

Reid, Max B.↗

Optimal Output Trajectory Redesign for Invertible Systems

Given a desired output trajectory, inversion-based techniques find input-state trajectories required to exactly track the output. These inversion-based techniques have been successfully applied to the endpoint tracking control of multijoint flexible manipulators and to aircraft control. The specified output trajectory uniquely determines the required input and state trajectories that are found through inversion. These input-state trajectories exactly track the desired output; however, they might not meet acceptable performance requirements. For example, during slewing maneuvers of flexible structures, the structural deformations, which depend on the required state trajectories, may be unacceptably large. Further, the required inputs might cause actuator saturation during an exact tracking maneuver, for example, in the flight control of conventional takeoff and landing aircraft. In such situations, a compromise is desired between the tracking requirement and other goals such as reduction of internal vibrations and prevention of actuator saturation; the desired output trajectory needs to redesigned. Here, we pose the trajectory redesign problem as an optimization of a general quadratic cost function and solve it in the context of linear systems. The solution is obtained as an off-line prefilter of the desired output trajectory. An advantage of our technique is that the prefilter is independent of the particular trajectory. The prefilter can therefore be precomputed, which is a major advantage over other optimization approaches. Previous works have addressed the issue of preshaping inputs to minimize residual and in-maneuver vibrations for flexible structures; Since the command preshaping is computed off-line. Further minimization of optimal quadratic cost functions has also been previously use to preshape command inputs for disturbance rejection. All of these approaches are applicable when the inputs to the system are known a priori. Typically, outputs (not inputs) are specified in tracking problems, and hence the input trajectories have to be computed. The inputs to the system are however, difficult to determine for non-minimum phase systems like flexible structures. One approach to solve this problem is to (1) choose a tracking controller (the desired output trajectory is now an input to the closed-loop system and (2) redesign this input to the closed-loop system. Thus we effectively perform output redesign. These redesigns are however, dependent on the choice of the tracking controllers. Thus the controller optimization and trajectory redesign problems become coupled; this coupled optimization is still an open problem. In contrast, we decouple the trajectory redesign problem from the choice of feedback-based tracking controller. It is noted that our approach remains valid when a particular tracking controller is chosen. In addition, the formulation of our problem not only allows for the minimization of residual vibration as in available techniques but also allows for the optimal reduction fo vibrations during the maneuver, e.g., the altitude control of flexible spacecraft. We begin by formulating the optimal output trajectory redesign problem and then solve it in the context of general linear systems. This theory is then applied to an example flexible structure, and simulation results are provided.

Devasia, S.↗

Design and Performance Improvements of the Prototype Open Core Flywheel Energy Storage System

A prototype magnetically suspended composite flywheel energy storage (FES) system is operating at the University of Maryland. This system, designed for spacecraft applications, incorporates recent advances in the technologies of composite materials, magnetic suspension, and permanent magnet brushless motor/generator. The current system is referred to as an Open Core Composite Flywheel (OCCF) energy storage system. This paper will present design improvements for enhanced and robust performance. Initially, when the OCCF prototype was spun above its first critical frequency of 4,500 RPM, the rotor movement would exceed the space available in the magnetic suspension gap and touchdown on the backup mechanical bearings would occur. On some occasions it was observed that, after touchdown, the rotor was unable to re-suspend as the speed decreased. Additionally, it was observed that the rotor would exhibit unstable oscillations when the control system was initially turned on. Our analysis suggested that the following problems existed: (1) The linear operating range of the magnetic bearings was limited due to electrical and magnetic saturation; (2) The inductance of the magnetic bearings was affecting the transient response of the system; (3) The flywheel was confined to a small movement because mechanical components could not be held to a tight tolerance; and (4) The location of the touchdown bearing magnifies the motion at the pole faces of the magnetic bearings when the linear range is crucial. In order to correct these problems an improved design of the flywheel energy storage system was undertaken. The magnetic bearings were re-designed to achieve a large linear operating range and to withstand load disturbances of at least 1 g. The external position transducers were replaced by a unique design which were resistant to magnetic field noise and allowed cancellation of the radial growth of the flywheel at high speeds. A central rod was utilized to ensure the concentricity of the magnetic bearings, the motor/generator, and the mechanical touchdown bearings. In addition, the mechanical touchdown bearings were placed at two ends of the magnetic bearing stack to restrict the motion at pole faces. A composite flywheel was made using a multi-ring interference assembled design for a high specific energy density. To achieve a higher speed and better efficiency, a permanent magnet DC brushless motor was specially designed and fabricated. A vacuum enclosure was constructed to eliminate windage losses for testing at high speeds. With the new improvements the OCCF system was tested to 20,000 RPM with a total stored energy of 15.9 WH and an angular momentum of 54.8 N-m-s (40.4 lb-ft-s). Motor current limitation, caused by power loss in the magnetic bearings, was identified as causing the limit in upper operating speed.

Pang, D.↗

Method of Distributions for Two‐Phase Flow in Heterogeneous Porous Media

Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.

Yang, Hyung Jun↗

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗

Recursive estimator for OSO-8 attitude

Modifications and enhancements that have been made to the Recursive Estimation Attitude Program (REAP) are discussed. Continuous attitudes for OSO-8 to + or - 0.05 degree accuracy are determined from Sun and star slit sensors mounted on the spinning portion of the spacecraft. The bulk of the attitude production is performed by a Weighted Least Squares (WLS) batch processor, but REAP is used for problem passes such as those involving gas jet maneuvers, sparse star fields, or star sensor saturation by high energy particles in the South Atlantic Anomaly.

Headrick, R. D.↗

Application of Modified Meshgraphnets for Subsurface Prediction during CO2 Sequestration

In the face of the increasingly dire consequences of anthropogenic climate change, capturing and storing carbon dioxide is paramount. However, several impediments exist to the safe and effective subsurface storage of CO2, such as cost of transport, identification of suitable sites for subsurface storage, and assessment of long-term risk from storage in subsurface aquifers. Accurate subsurface modeling is necessary to ensure that CO2 storage is both safe and effective. Still, such modeling has traditionally required either substantial time and computational power (numerical simulation) or a substantial amount of pre-existing data for training (machine learning models). Additionally, these models lack flexibility in dealing with both changes in discretization of the input data and generalizability beyond the data on which they are trained. In order to address these issues, this research applies graph neural networks (GNNs) to predict subsurface saturation and pressure during CO₂ injection in a model of the Illinois Basin-Decatur Project (IBDP). GNNs provide a flexible, intuitive method for representing and manipulating complex unstructured data, which is often found in many practical domain problems such as fluid flow and subsurface characterization. These unstructured grids are easily represented in GNNs by representing spatially-localized features such as permeability, porosity, saturation, and pressure as nodes in a graph and relationships between these properties as edges connecting these nodes. This research applies a specific GNN model called MeshGraphNets (MGN) to model the change in CO2 saturation and pressure over a 50-month time period (36 months of injection, 14 months post-injection). The MGN model leverages a message passing process that allows the network to learn both the spatial and temporal dynamics of this system simultaneously. Additionally, training on a limited dataset (64 realizations, 20 time points each) resulted in a high degree of accuracy in saturation prediction both within the same timeframe as the training (20 months, 0.039 average RMSE) and when projecting out to the end of injection (36 months, 0.053 average RMSE). Temporal predictions such as those generated by MGNs and other similar models are prone to accumulated error over time; in order to address this, a multi-step rollout (MSR) training process was applied to calculate training loss. This method mimics the forward prediction during inference by “rolling out” multiple time points in a single training step using the previous prediction as input to the MGN model. By calculating the loss several time steps forward from the current prediction, the model is forced to find a more stable state over time. Application of MSR to the MGN model resulted in an average 15% reduction in inference error over time during forward prediction. This study showcases the immense potential of GNNs as a game-changing methodology for predicting pressure and saturation evolution in CCS projects, ultimately paving the way for more sustainable and effective carbon storage solutions. Presentation prepared for the 2024 AiChE Annual Meeting, October 27 to November 1 2024, San Diego, CA.

Holcomb, Paul↗

Optimal Control Allocation for Distributed Electric Propulsion in A Series/Parallel Partial Hybrid Powertrain

The SUbsonic Single Aft eNgine (SUSAN) Electrofan is a NASA concept transport aircraft representative of technology anticipated for a 2040 entry-into-service date. The powertrain consists of a single thrust-producing geared turbofan engine with generators driving a series/parallel partial hybrid power/propulsion system. The architecture includes 16 underwing contrarotating fans, eight on each side. The distributed fans can be used by the flight control system to augment or replace the rudder function. This paper sets up the optimal control problem of setpoint determination for individual wingfans in the distributed propulsion system, accounting for electrical string efficiencies, saturations, and failures. The solution minimizes power consumption while maintaining thrust and torque on the airframe for maneuvering. Additionally, thrust that would have been lost due to temporary fan speed or power saturation is optimally redistributed to maintain overall desired thrust and torque on the aircraft. A simulation of a coordinated turn utilizing the distributed electric propulsion for yaw rate control in a multiple wingfan failure scenario demonstrates the robustness of the powertrain design to failures and helps define its limitations.

Distributed Electric Propulsion↗

Optimal Control Allocation for Distributed Electric Propulsion in A Series/Parallel Partial Hybrid Powertrain

The SUbsonic Single Aft eNgine (SUSAN) Electrofan is a NASA concept transport aircraft representative of technology anticipated for a 2040 entry-into-service date. The powertrain consists of a single thrust-producing geared turbofan engine with generators driving a series/parallel partial hybrid power/propulsion system. The architecture includes 16 underwing contrarotating fans, eight on each side. The distributed fans can be used by the flight control system to augment or replace the rudder function. This paper sets up the optimal control problem of setpoint determination for individual wingfans in the distributed propulsion system, accounting for electrical string efficiencies, saturations, and failures. The solution minimizes power consumption while maintaining thrust and torque on the airframe for maneuvering. Additionally, thrust that would have been lost due to temporary fan speed or power saturation is optimally redistributed to maintain overall desired thrust and torque on the aircraft. A simulation of a coordinated turn utilizing the distributed electric propulsion for yaw rate control in a multiple wingfan failure scenario demonstrates the robustness of the powertrain design to failures and helps define its limitations.

Distributed Electric Propulsion↗

Super Resolving Unrolled Neural Networks for Remote Sensing

In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.

47 OTHER INSTRUMENTATION↗

Cognitive consequences of clumsy automation on high workload, high consequence human performance

The growth of computational power has fueled attempts to automate more of the human role in complex problem solving domains, especially those where system faults have high consequences and where periods of high workload may saturate the performance capacity of human operators. Examples of these domains include flightdecks, space stations, air traffic control, nuclear power operation, ground satellite control rooms, and surgical operating rooms. Automation efforts may have unanticipated effects on human performance, particularly if they increase the workload at peak workload times or change the practitioners' strategies for coping with workload. Smooth and effective changes in automation requires detailed understanding of the congnitive tasks confronting the user: it has been called user centered automation. The introduction of a new computerized technology in a group of hospital operating rooms used for heart surgery was observed. The study revealed how automation, especially 'clumsy automation', effects practitioner work patterns and suggest that clumsy automation constrains users in specific and significant ways. Users tailor both the new system and their tasks in order to accommodate the needs of process and production. The study of this tailoring may prove a powerful tool for exposing previously hidden patterns of user data processing, integration, and decision making which may, in turn, be useful in the design of more effective human-machine systems.

Cook, Richard I.↗

A New Technique for Compensating Joint Limits in a Robot Manipulator

A new robust, optimal, adaptive technique for compensating rate and position limits in the joints of a six degree-of-freedom elbow manipulator is presented. In this new algorithm, the unmet demand as a result of actuator saturation is redistributed among the remaining unsaturated joints. The scheme is used to compensate for inadequate path planning, problems such as joint limiting, joint freezing, or even obstacle avoidance, where a desired position and orientation are not attainable due to an unrealizable joint command. Once a joint encounters a limit, supplemental commands are sent to other joints to best track, according to a selected criterion, the desired trajectory.

Litt, Jonathan↗

Numerical comparison of strong Langmuir turbulence models

Two models of Langmuir turbulence, the nonlinear Schroedinger equation and the Zakharov equations, are solved numerically for an initial value problem in which the electric field evolves from an almost flat initial condition via the modulational instability and finally saturates into a set of solitons. The two models agree well with each other only when the initial dimensionless electric field has an amplitude less than unity. An analytic soliton gas model consisting of equal-amplitude, randomly spaced, zero-speed solitons is remarkably good at reproducing the time-averaged Fourier spectra in both cases.

Shen, Mei-Mei↗

Integration of airborne Thematic Mapper Simulator (TMS) data and digitized aerial photography via an ISH transformation

A simple method for enhancing the spatial and spectral resolution of disparate data sets is presented. Two data sets, digitized aerial photography at a nominal spatial resolution 3,7 meters and TMS digital data at 24.6 meters, were coregistered through a bilinear interpolation to solve the problem of blocky pixel groups resulting from rectification expansion. The two data sets were then subjected to intensity-saturation-hue (ISH) transformations in order to 'blend' the high-spatial-resolution (3.7 m) digitized RC-10 photography with the high spectral (12-bands) and lower spatial (24.6 m) resolution TMS digital data. The resultant merged products make it possible to perform large-scale mapping, ease photointerpretation, and can be derived for any of the 12 available TMS spectral bands.

Ambrosia, Vincent G.↗

Removing Spectral Diagnostics of Galactic and Stellar X-Ray Emission from Charged Exchange Recombination

Our research uses the electron beam ion trap (EBIT) at the Lawrence Livermore National Laboratory to study X-ray emission from the charge exchange (CX) of highly charged ions with neutral gases. The resulting data help to fill a void in existing experimental and theoretical understanding of this atomic physics process, and are needed to explain all or part of the observed X-ray emission from the soft X-ray background, stellar winds, the Galactic Center and Galactic Ridge, supernova ejecta, and photoionized nebulae. Appreciation of the astrophysical relevance of our work continues to grow with the publication of roughly a dozen papers in the past four years describing Chandra and XMM observations of geocoronal and heliospheric CX emission, the temporal variation of such emission and correlation with X-ray emission enhancements observed by ROSAT, the theoretical spatial distribution of that emission, and CX emission around other stars. A similar number of papers were also published during that time describing CX emission from planets and comets. We expect that the launch of ASTRSE2, with its second-generation XRS microcalo- (with 6-eV resolution), will reveal even more clearly the contributions of CX to astrophysical emission. In our EBIT work we collected CX spectra from such ions as H-like and He-like Ne, Ar, and Fe. Our early measurements were made with a high-purity Ge detector, but during the second year we began operation of the first-generation XRS microcalorimeter (a twin of the XRS on ASTRO-E) and greatly improved the resolution of our measurements from roughly 150 eV (FWHM) with the Ge detectors to 10 eV with the XRS. We found that saturation of the XRS counting apparatus, which we described in our proposal as a potential concern, is not a problem for studying CX. During the course of our research, we expanded the number of injection gases permitted by the LLNL safety team, purchased and eventually operated an atomic H source, and clearly demonstrated the feasibility of our longer-range plan. For example, we successfully injected He into EBIT (not a small feat because of the difficulty of maintaining a good vacuum with He and avoiding electrical breakdown) to collect a H-like oxygen CX spectrum. The highest energy CX spectrum recorded with the XRS to date is that of the Ar K-shell emission. These measurements provided the first observation of the relative intensity ratios of resolved He-like singlet and triplet n=2->1 lines. We also carried out measurements of He-like Ne as a function of collision energy (i.e., ion temperature). Significant differences in the resulting x-ray spectra were noted. In all cases, the intensity of high-n H-like Lyman lines is significantly higher than current theoretical CX models predict.

Wargelin, Brad↗

Lunar Impact Basins: Stratigraphy, Sequence and Ages from Superposed Impact Crater Populations Measured from Lunar Orbiter Laser Altimeter (LOLA) Data

Impact basin formation is a fundamental process in the evolution of the Moon and records the history of impactors in the early solar system. In order to assess the stratigraphy, sequence, and ages of impact basins and the impactor population as a function of time, we have used topography from the Lunar Orbiter Laser Altimeter (LOLA) on the Lunar Reconnaissance Orbiter (LRO) to measure the superposed impact crater size-frequency distributions for 30 lunar basins (D ≥ 300 km). These data generally support the widely used Wilhelms sequence of lunar basins, although we find significantly higher densities of superposed craters on many lunar basins than derived by Wilhelms (50% higher densities). Our data also provide new insight into the timing of the transition between distinct crater populations characteristic of ancient and young lunar terrains. The transition from a lunar impact flux dominated by Population 1 to Population 2 occurred before the mid-Nectarian. This is before the end of the period of rapid cratering, and potentially before the end of the hypothesized Late Heavy Bombardment. LOLA-derived crater densities also suggest that many Pre-Nectarian basins, such as South Pole-Aitken, have been cratered to saturation equilibrium. Finally, both crater counts and stratigraphic observations based on LOLA data are applicable to specific basin stratigraphic problems of interest; for example, using these data, we suggest that Serenitatis is older than Nectaris, and Humboldtianum is younger than Crisium. Sample return missions to specific basins can anchor these measurements to a Pre-Imbrian absolute chronology.

Fassett, C. I.↗