Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source release”

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,081 records · Page 60

Sewage Treatment Plant #1 Area, SWMU 117 Per- and Polyfluoroalkyl Substances (PFAS) Site Assessment Progress Report

This Per- and Polyfluoroalkyl Substances (PFAS) Site Assessment (SA) Progress Report (SAPR) discusses the investigation activities and findings for the Sewage Treatment Plant #1 (STP1) Area located at Kennedy Space Center (KSC), Florida (Figure 1-1). This site has been designated Solid Waste Management Unit (SWMU) 117 under KSC’s Resource Conservation and Recovery Act (RCRA) Corrective Action Program, as the sewage treatment plant and associated areas were identified as a potential source of PFAS to the environment. This PFAS SAPR was prepared by Tetra Tech, Inc., for the National Aeronautics and Space Administration (NASA) under Indefinite Delivery Indefinite Quantity Contract 80KSC019D0011-80KSC019F0070. This is the first progress report to document on-going SA activities; supplemental progress reports will be provided as additional data is collected. PFAS SA activities were conducted between April 2020 and March 2022 to collect additional data to supplement the existing datasets to better understand the extent of PFAS impacts to the environment in the STP1 Area. The SA for the STP1 Area covers an approximately 130-acre investigation area with multiple structures and buildings. The focus of the SA is the STP1 Complex and associated structures, including the former Polishing Pond, former Sludge Disposal8 Area, and former Spray Field. The STP1 Complex is located in the KSC Industrial Area, at the southwest corner of 4th Street SE and C Avenue SE. The STP1 Complex is located approximately ¼-mile south and downgradient of the Fire Station #1 site (SWMU 116), which is also currently undergoing a PFAS SA because of potential releases of PFAS-containing aqueous film-forming foam (AFFF). During the SA, a total of seven soil, 131 groundwater direct push technology (DPT), 24 groundwater monitoring well, and 11 surface water samples were collected between April 2020 and March 2022. Monitoring well samples were analyzed for 18 PFAS compounds, with all other samples analyzed for 28 PFAS compounds. The SA sample results were used along with historical results to evaluate the extent of PFAS impacts to the environment in the STP1 Area. Data generated to date and prior results were screened against the United States Environmental Protection Agency (USEPA) May 2022 Tap Water Regional Screening Levels (RSL) for groundwater and residential RSLs for soil (hazard quotient of 0.1). Surface water results were screened against the State of Florida Human Health Surface Water Screening Levels (SWSLs). Results from the SA showed exceedances of the applicable screening criteria for groundwater and surface water. Considering the current and historical dataset, PFOS is the prevalent PFAS compound. Based on these results, additional groundwater DPT and surface water sampling should be considered for PFAS analysis, focused on evaluating surface water bodies in the southeast portion of the Industrial Area, which discharge into the Banana River. Additionally, installation of monitoring wells should be considered to evaluate the interaction between the groundwater and surface water at the site. Collection of additional samples for TOC analysis should also be considered from representative groundwater (saturated soils) and surface water locations to further evaluate potential correlations between PFAS and TOC to provide a more comprehensive dataset to assist in fate and transport analyses.

Sewage Treatment Plant↗

Noachian Martian Volcanics a Water Source

H2O was supplied to the Noachian atmosphere by eruptions, or in association with large impacts. Most water outgassed into an extremely cold atmosphere, and condensate deposits were inevitable. High heat flow could lead to subglacial melting only if ice thicknesses were greater than 500-1000m, which is extremely unlikely. Subareal melting and flow is contingent upon temperatures periodically exceeding 273 K, and retarding evaporative loss of the flow. In still air, evaporation into a dry atmosphere is in the free convection regime, and a stream with 2 cu m/s discharge, flowing 1 m/s could persist for hundreds of days and cover distances greater than any valley reach. The zero-wind-shear condition is considered implausible however. We investigate the possibility that evaporation rates were suppressed because the atmosphere was regionally charged with H2O as it moved over snow/ice fields. Our initial concern is precipitation from volcanic plumes. A Kilauea-style eruption on the martian surface would cover a 10km circular deposit with 10cm of H2O, if all H2O could be precipitated near the vent. The characteristics of the eruption at the vent, (vent size, temperature, H2O content, etc.) are independent of the environmental conditions. The subsequent behavior of the plume, including precipitation of ash and H2O condensate depends strongly on the environment. Hence, the proximal fate of volcanic H2O is amenable to treatment in a model. A simple bulk thermodynamic model of the rise of an H2O plume through a stably stratified CO2 atmosphere, with only adiabatic cooling, produces runaway plume rise. A more complex treatment includes the effects of latent heat release, wind shear along the plume, divergence of ash and H2O, and will yield more realistic estimates of H2O transport in eruptive plumes. Results of these simulations will be presented.

Zent, A. P.↗

Faint blue counts from formation of dwarf galaxies at z approximately equals 1

The nature of faint blue objects (FBO's) has been a source of much speculation since their detection in deep CCD images of the sky. Their high surface density argues against them being progenitors of present-day bright galaxies and since they are only weakly clustered on small scales, they cannot be entities that merged together to form present-day galaxies. Babul & Rees (1992) have suggested that the observed faint blue counts may be due to dwarf elliptical galaxies undergoing their initial starburst at z is approximately equal to 1. In generic hierarchical clustering scenarios, however, dwarf galaxy halos (M is approximately 10(exp 9) solar mass) are expected to form at an earlier epoch; for example, typical 10(exp 9) solar mass halos will virialize at z is approximately equal to 2.3 if the power-spectrum for the density fluctuations is that of the standard b = 2 cold dark matter (CDM) model. Under 'ordinary conditions' the gas would rapidly cool, collect in the cores, and undergo star-formation. Conditions at high redshifts are far from 'ordinary'. The intense UV background will prevent the gas in the dwarf halos from cooling, the halos being released from their suspended state only when the UV flux has diminished sufficiently.

Babul, Arif↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

CFD Multiphysics Tool

The recent bold initiatives to expand the human presence in space require innovative approaches to the design of propulsion systems whose underlying technology is not yet mature. The space propulsion community has identified a number of candidate concepts. A short list includes solar sails, high-energy-density chemical propellants, electric and electromagnetic accelerators, solar-thermal and nuclear-thermal expanders. For each of these, the underlying physics are relatively well understood. One could easily cite authoritative texts, addressing both the governing equations, and practical solution methods for, e.g. electromagnetic fields, heat transfer, radiation, thermophysics, structural dynamics, particulate kinematics, nuclear energy, power conversion, and fluid dynamics. One could also easily cite scholarly works in which complete equation sets for any one of these physical processes have been accurately solved relative to complex engineered systems. The Advanced Concepts and Analysis Office (ACAO), Space Transportation Directorate, NASA Marshall Space Flight Center, has recently released the first alpha version of a set of computer utilities for performing the applicable physical analyses relative to candidate deep-space propulsion systems such as those listed above. PARSEC, Preliminary Analysis of Revolutionary in-Space Engineering Concepts, enables rapid iterative calculations using several physics tools developed in-house. A complete cycle of the entire tool set takes about twenty minutes. PARSEC is a level-zero/level-one design tool. For PARSEC s proof-of-concept, and preliminary design decision-making, assumptions that significantly simplify the governing equation sets are necessary. To proceed to level-two, one wishes to retain modeling of the underlying physics as close as practical to known applicable first principles. This report describes results of collaboration between ACAO, and Embry-Riddle Aeronautical University (ERAU), to begin building a set of level-two design tools for PARSEC. The "CFD Multiphysics Tool" will be the propulsive element of the tool set. The name acknowledges that space propulsion performance assessment is primarily a fluid mechanics problem. At the core of the CFD Multiphysics Tool is an open-source CFD code, HYP, under development at ERAU. ERAU is renowned for its undergraduate degree program in Aerospace Engineering the largest in the nation. The strength of the program is its applications-oriented curriculum, which culminates in one of three two-course Engineering Design sequences: Aerospace Propulsion, Spacecraft, or Aircraft. This same philosophy applies to the HYP Project, albeit with fluid physics modeling commensurate with graduate research. HYP s purpose, like the Multiphysics Tool s, is to enable calculations of real (three-dimensional; geometrically complex; intended for hardware development) applications of high speed and propulsive fluid flows.

Perrell, Eric R.↗

Non Radiation Hardened Microprocessors in Spaced Based Remote Sensing Systems

The CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) mission is a comprehensive suite of active and passive sensors including a 20Hz 230mj Nd:YAG lidar, a visible wavelength Earth-looking camera and an imaging infrared radiometer. CALIPSO flies in formation with the Earth Observing System Post-Meridian (EOS PM) train, provides continuous, near-simultaneous measurements and is a planned 3 year mission. CALIPSO was launched into a 98 degree sun synchronous Earth orbit in April of 2006 to study clouds and aerosols and acquires over 5 gigabytes of data every 24 hours. The ground track of one CALIPSO orbit as well as high and low intensity South Atlantic Anomaly outlines is shown. CALIPSO passes through the SAA several times each day. Spaced based remote sensing systems that include multiple instruments and/or instruments such as lidar generate large volumes of data and require robust real-time hardware and software mechanisms and high throughput processors. Due to onboard storage restrictions and telemetry downlink limitations these systems must pre-process and reduce the data before sending it to the ground. This onboard processing and realtime requirement load may mean that newer more powerful processors are needed even though acceptable radiation-hardened versions have not yet been released. CALIPSO's single board computer payload controller processor is actually a set of four (4) voting non-radiation hardened COTS Power PC 603r's built on a single width VME card by General Dynamics Advanced Information Systems (GDAIS). Significant radiation concerns for CALIPSO and other Low Earth Orbit (LEO) satellites include the South Atlantic Anomaly (SAA), the north and south poles and strong solar events. Over much of South America and extending into the South Atlantic Ocean the Van Allen radiation belts dip to just 200-800km and spacecraft entering this area are subjected to high energy protons and experience higher than normal Single Event Upset (SEU) and Single Event Latch-up (SEL) rates. Although less significant, spacecraft flying in the area around the poles experience similar upsets. Finally, powerful solar proton events in the range of 10MeV/10pfu to 100MeV/1pfu as are forecasted and tracked by NOAA's Space Environment Center in Colorado can result in Single Event Upset (SEU), Single Event Latch-up (SEL) and permanent failures such as Single Event Gate Rupture (SEGR) in some technologies. (Galactic Cosmic Rays (GCRs) are another source, especially for gate rupture) CALIPSO mitigates common radiation concerns in its data handling through the use of redundant processors, radiation-hardened Application Specific Integrated Circuits (ASIC), hardware-based Error Detection and Correction (EDAC), processor and memory scrubbing, redundant boot code and mirrored files. After presenting a system overview this paper will expand on each of these strategies. Where applicable, related on-orbit data collected since the CALIPSO initial boot on May 4, 2006 will be noted.

Decoursey, Robert J.↗

Is There a Cosmological Constant?

The grant contributed to the publication of 18 refereed papers and 5 conference proceedings. The primary uses of the funding have been for page charges, travel for invited talks related to the grant research, and the support of a graduate student, Charles Keeton. The refereed papers address four of the primary goals of the proposal: (1) the statistics of radio lenses as a probe of the cosmological model (#1), (2) the role of spiral galaxies as lenses (#3), (3) the effects of dust on statistics of lenses (#7, #8), and (4) the role of groups and clusters as lenses (#2, #6, #10, #13, #15, #16). Four papers (#4, #5, #11, #12) address general issues of lens models, calibrations, and the relationship between lens galaxies and nearby galaxies. One considered cosmological effects in lensing X-ray sources (#9), and two addressed issues related to the overall power spectrum and theories of gravity (#17, #18). Our theoretical studies combined with the explosion in the number of lenses and the quality of the data obtained for them is greatly increasing our ability to characterize and understand the lens population. We can now firmly conclude both from our study of the statistics of radio lenses and our survey of extinctions in individual lenses that the statistics of optically selected quasars were significantly affected by extinction. However, the limits on the cosmological constant remain at lambda < 0.65 at a 2-sigma confidence level, which is in mild conflict with the results of the Type la supernova surveys. We continue to find that neither spiral galaxies nor groups and clusters contribute significantly to the production of gravitational lenses. The lack of group and cluster lenses is strong evidence for the role of baryonic cooling in increasing the efficiency of galaxies as lenses compared to groups and clusters of higher mass but lower central density. Unfortunately for the ultimate objective of the proposal, improved constraints on the cosmological constant, the next large survey for gravitational lenses did not release its results during the term of the proposal. The research supported the career development. of six graduate students (polar, Fletcher, Herold, Keeton, Deng and Rusin) and two post-docs (Labor and Munoz).

Kochanek, Christopher↗

Dynamical Model for the Zodiacal Cloud and Sporadic Meteors

The solar system is dusty, and would become dustier over time as asteroids collide and comets disintegrate, except that small debris particles in interplanetary space do not last long. They can be ejected from the solar system by Jupiter, thermally destroyed near the Sun, or physically disrupted by collisions. Also, some are swept by the Earth (and other planets), producing meteors. Here we develop a dynamical model for the solar system meteoroids and use it to explain meteor radar observations. We find that the Jupiter Family Comets (JFCs) are the main source of the prominent concentrations of meteors arriving to the Earth from the helion and antihelion directions. To match the radiant and orbit distributions, as measured by the Canadian Meteor Orbit Radar (CMOR) and Advanced Meteor Orbit Radar (AMOR), our model implies that comets, and JFCs in particular, must frequently disintegrate when reaching orbits with low perihelion distance. Also, the collisional lifetimes of millimeter particles may be longer (approx. > 10(exp 5) yr at 1 AU) than postulated in the standard collisional models (approx 10(exp 4) yr at 1 AU), perhaps because these chondrule-sized meteoroids are stronger than thought before. Using observations of the Infrared Astronomical Satellite (IRAS) to calibrate the model, we find that the total cross section and mass of small meteoroids in the inner solar system are (1.7-3.5) 10(exp 11) sq km and approx. 4 10(exp 19) g, respectively, in a good agreement with previous studies. The mass input required to keep the Zodiacal Cloud (ZC) in a steady state is estimated to be approx. 10(exp 4)-10(exp 5) kg/s. The input is up to approx 10 times larger than found previously, mainly because particles released closer to the Sun have shorter collisional lifetimes, and need to be supplied at a faster rate. The total mass accreted by the Earth in particles between diameters D = 5 micron and 1 cm is found to be approx 15,000 tons/yr (factor of 2 uncertainty), which is a large share of the accretion flux measured by the Long Term Duration Facility (LDEF). Majority of JFC particles plunge into the upper atmosphere at <15 km/s speeds, should survive the atmospheric entry, and can produce micrometeorite falls. This could explain the compositional similarity of samples collected in the Antarctic ice and stratosphere, and those brought from comet Wild 2 by the Stardust spacecraft. Meteor radars such as CMOR and AMOR see only a fraction of the accretion flux (approx 1- 10% and approx 10-50%, respectively), because small particles impacting at low speeds produce ionization levels that are below these radars detection capabilities.

Nesvorny, David↗

Pre-Demonstration Development of Controlled-Release Corrosion Inhibitors and Healing Agents as Alternatives to Hexavalent Chromium

This Limited Scope Study was directed to achieve the following objectives: (1) Scale-up of materials that can meet MIL-PRF-23377 (solvent-based primer). (2) Provide evidence of resistance to aircraft alkaline cleaners and de-icing fluids. (3) Provide formulation for initial ecological and toxicity screening. (4) Submit an interim report that will provide the basis for a future ESTCP (the Dept. of Defense's (DoD's) Environmental Security Technology Certification Program) demonstration effort.The encapsulation procedure for two corrosion inhibitors was scaled up from lab scale to 2.0 kilograms. This scale will accommodate high volume production of coatings with encapsulated corrosion inhibitors for large surface areas during a follow-on ESTCP demonstration effort. Test results provided evidence of resistance to alkaline cleaners and aircraft de-icing fluids and compliance with MIL-PRF-23377. All the MIL-PRF-23377 requirements were met with the exception of the adhesion requirement, in top coated panels, and the flexibility requirement. Work on solving these two problems will be done prior to validation and demonstration of the technology. An initial ecological and toxicity screening identified several factors that were not critical to the acceptance of the encapsulated corrosion inhibitor. This new technology will lead to environmentally friendly alternatives to hexavalent chromium that will enable DoD to protect its assets. Field demonstration, licensing, qualification, and commercialization will allow its real world utilization. Abstract for Appendix A. Toxicology Report: Corrosion controls a significant problem for military equipment. Many surface treatment materials currently in use on items of military or aerospace materiel contain chromium (VI) as a component, a source of significant human health and environmental hazard. The formulation under evaluation contains no chromium. Most of the proposed components are low to moderate toxicity and are not a significant concern. Two of the formulation components are associated with some cancer hazard - 2-mercaptobenzothiazole and formaldehyde, but both find widespread use in industry and additional investigation is not required for this project.

Calle, Luz Marina↗

Time Variable Gravity modeling for Precise Orbits Across the TOPEX/Poseidon, Jason-l and Jason-2 Missions

Modeling of the Time Variable Gravity (TVG) is believed to constitute one of the the largest remaining source of orbit error for altimeter satellite POD. The GSFC operational TVG model consists of forward modeling the atmospheric gravity using ECMWF 6-hour pressure data, a GRACE derived 20x20 annual field to account for changes in the hydrology and ocean water mass, and linear rates for C20, C30, C40, based on 17 years of SLR data analysis (IERS 2003) using the EIGEN-GL04S1 (a GRACE+Lageos-based geopotential solution). Although the GSFC Operational model can be applied from 1987, there may be long-term variations not captured by these linear models, and more importantly the linear models may not be consistent with more recent surface mass trends due to global climate change, We have evaluated the impact of TVG in two different wavs: (1) by using the more recent EIGEN-6S gravity model developed by the GFZ/GRGS tearm, which consists of annual, semi-annual and secular changes in the coefficients to 50x50 determined over 8(?) years of GRACE+Lageos+GOCE data (2003-200?): (2) Application of 4x4 solutions developed from a multi satellite SLR+DORIS solution based on GGM03S that span the period from 1993 to 2011. We have evaluated the recently released EIGEN6s static and time-varying gravity field for Jason-2 (J2). Jason-I (J1), and TOPEX/Posiedon (TP) Precise Orbit Determination (POD) spanning 1993-2011. Although EIGEN6s shows significant improvement for J2POD spanning 2008 - 2011, it also shows significant degradation for TP POD from 1992. The GSFC 4x4 time SLR+DORIS-based series spans 1993 to mid 2011, and shows promise for POD. We evaluate the performance of the different TVG models based on analysis of tracking data residuals use of independent data such as altimeter crossovers, and through analysis of differences with internally-generated and externally generated orbits.

Zelensky, Nikita P.↗

Factors controlling fluxes of volatile sulfur compounds in Sphagnum peatlands

Exchange of DMS and OCS between the surface of Sphagnum peatlands and the atmosphere were measured with dynamic (S-free sweep air) and static enclosures. DMS emission rates determined by both methods were comparable. The dynamic method provided positive OCS flux rates (emission) for measurements performed at sites containing Sphagnum. Conversely, data from the static method indicated that OCS was consumed from the atmosphere. Short and long-term impacts of increased S deposition on fluxes of volatile S compounds (VSC's) from Sphagnum peatlands were investigated in a poor fen (Mire 239) at the Experimental Lakes Area, Ontario, Canada. Additional experiments were conducted in a poor fen (Sallie's Fen in Barrington, NH, USA). At Mire 239, emissions of VSC's were monitored, before and after acidification, at control and experimental sections within two major physiographic areas of the mire (oligotrophic and minerotrophic). DMS was the predominant VSC released from Mire 239 and varied largely with time and space. Sulfur addition did not affect DMS emissions in a period of hours to a few days. DMS emissions in the experimental oligotrophic area of the mire was approximately 3-fold greater than in the control oligotrophic area, and approximately 10-fold greater than in the minerotrophic zones. These differences could be due to a combination of differences in types of vegetation, nutritional status, and S input. At Sallie's Fen, DMS fluxes were not significantly affected by sulfate amendments, while DMS and MSH concentrations increased greatly with time in the top 10 cm of the peat column. The major environmental factors controlling fluxes of DMS in a Sphagnum-dominated peatland were investigated in Sallie's Fen, NH. DMS emissions from the surface of the peatland varied greatly over 24 hours and seasonally. Temperature seemed to be the major environmental factor controlling these variabilities. Concentrations of dissolved VSC's varied with time and space throughout the fen. Dissolved DMS, MSH, and OCS in the surface of the water table were supersaturated with respect to their concentrations in the atmosphere. Sphagnum mosses did not appear to be a direct source of VSC's, however they increase transport of DMS from the peat surface to the atmosphere.

Demello, William Zamboni↗

A coronal mass ejection model for the 1992 July 15 flare on AU Microscopii observed by the extreme ultraviolet explorer

The dM1e flare star AU Microscopii (AU Mic) was observed by the EUVE Deep Survey Instrument on 1992 July 14-18. A large flare was detected in the Deep Survey Lexan/Boron (DS Lex/B)(65-190 A) band and the SW (70-190 A) and MW (140-380 A) spectrometers. The flare consisted of a sharp impulsive peak lasting approximately 2 hours followed by a decaying tail lasting about a day. We present a simple, single temperature, dynamic model for the flare decay which is consistent with the DS Lex/B light curve and reproduces the strongest, high-temperature spectral lines in the released EUVE spectra. In this model, we assume the long decay time is due to an ejected, magnetically confined, low beta plasmoid expanding self similarly in the ambient medium in a manner reminiscent of solar coronal mass ejections. We demonstrate that the long tail of the DS Lex/B light curve can be explained by rapid expansion, causing the plasma to become tenuous sufficiently quickly that it avoids catastrophic radiative cooling. From this model, we estimate the mass of the plasmoid to be approximately = 10(exp 20) g and the total energy of the event to be approximately = 10(exp 36) ergs. These values are approximately 10(exp 4) times as large as those seen during the largest solar coronal mass ejection (CME) events. We argue that the results of our model are consistent with other measurements of stellar flare parameters. We also estimate a mass-loss rate of a few times 10(exp -13) Solar mass/yr and discuss the role of mass loss from dMe stars in the mass balance of the interstellar medium. We estimate the rotational braking timescale from these events to be less than 500 million years and suggest that CME's may be an important source of angular momentum loss from late-type stars.

Cully, Scott L.↗

HEASARC Software Archive

(1) Chandra Archive: SAO has maintained the interfaces through which HEASARC gains access to the Chandra Data Archive. At HEASARC's request, we have implemented an anonymous ftp copy of a major part of the public archive and we keep that archive up-to- date. SAO has participated in the ADEC interoperability working group, establishing guidelines or interoperability standards and prototyping such interfaces. We have provided an NVO-based prototype interface, intending to serve the HEASARC-led NVO demo project. HEASARC's Astrobrowse interface was maintained and updated. In addition, we have participated in design discussions surrounding HEASARC's Caldb project. We have attended the HEASARC Users Group meeting and presented CDA status and developments. (2) Chandra CALDB: SA0 has maintained and expanded the Chandra CALDB by including four new data file types, defining the corresponding CALDB keyword/identification structures. We have provided CALDB upgrades for the public (CIAO) and for Standard Data Processing. Approximately 40 new files have been added to the CALDB in these version releases. There have been in the past year ten of these CALDB upgrades, each with unique index configurations. In addition, with the inputs from software, archive, and calibration scientists, as well as CIAO/SDP software developers, we have defined a generalized expansion of the existing CALDB interface and indexing structure. The purpose of this is to make the CALDB more generally applicable and useful in new and future missions that will be supported archivally by HEASARC. The generalized interface will identify additional configurational keywords and permit more extensive calibration parameter and boundary condition specifications for unique file selection. HEASARC scientists and developers from SAO and GSFC have become involved in this work, which is expected to produce a new interface for general use within the current year. (3) DS9: One of the decisions that came from last year's HEADCC meeting was to make the ds9 image display program the primary vehicle for displaying line graphics (as well as images). The first step required to make this possible was to enhance the line graphics capabilities of ds9. SAO therefore spent considerable effort upgrading ds9 to use Tcl 8.4 so that the BLT line graphics package could be built and imported into ds9 from source code, rather than from a pre-built (and generally outdated) shared library. This task, which is nearly complete, allows us to extend BLT as needed for the HEAD community. Following HEADCC discussion concerning archiving and the display of archived data, we extended ds9 to support full access to many astronomical Web-based archives sites, including HEASARC, MAST, CHANDRA, SKYVIEW, ADS, NED, SIMBAD, IRAS, NVRO, SAO TDC, and FIRST. Using ds9's new internal Web access capabilities, these archives can be accessed via their Web page. FITS images, plots, spectra, and journal abstracts can be referenced, down-loaded, and displayed directly and easily in ds9. For more information, see: http://hea-www.harvard.edu/saord/ds9. Also after the HEADCC discussion concerning region filtering, we extended the Funtools sample implementation of region filtering as described in: http://hea-www.harvard.edu/saord/funtools/regions.html. In particular, we added several new composite regions for event and image filtering, including elliptical and box annuli. We also extended the panda (Pie AND Annulus) region support to include box pandas and elliptical pandas. These new composite regions are especially useful in programs that need to count photons in each separate region using only a single pass through the data. Support for these new regions was added to ds9. In the same vein, we developed new region support for filtering images using simple FITS image masks, i.e. 8-bit or 16-bit FITS images where the value of a pixel is the region id number for that pixel. Other important enhancements to DS9 this year, include supporor multiple world coordinate systems, three dimensional event file binning, image smoothing, region groups and tags, the ability to save images in a number of image formats (such as JPEG, TIFF, PNG, FITS), improvements in support for integrating external analysis tools, and support for the virtual observatory. In particular, a full-featured web browser has been implemented within D S 9 . This provides support for full access to HEASARC archive sites such as SKYVIEW and W3BROWSE, in addition to other astronomical archives sites such as MAST, CHANDRA, ADS, NED, SIMBAD, IRAS, NVRO, SA0 TDC, and FIRST. From within DS9, the archives can be searched, and FITS images, plots, spectra, and journal abstracts can be referenced, downloaded and displayed The web browser provides the basis for the built-in help facility. All DS9 documentation, including the reference manual, FAQ, Know Features, and contact information is now available to the user without the need for external display applications. New versions of DS9 maybe downloaded and installed using this facility. Two important features used in the analysis of high energy astronomical data have been implemented in the past year. The first is support for binning photon event data in three dimensions. By binning the third dimension in time or energy, users are easily able to detect variable x-ray sources and identify other physical properties of their data. Second, a number of fast smoothing algorithms have been implemented in DS9, which allow users to smooth their data in real time. Algorithms for boxcar, tophat, and gaussian smoothing are supported.

White, Nicholas↗

Tropical Diabatic Heating and the Role of Convective Processes as Represented in Several Contemporary Climate Models

One of the most fundamental properties of the global heat balance is the net heat input into the tropical atmosphere that helps drive the planetary atmospheric circulation. Although broadly understood in terms of its gross structure and balance of source / sink terms, incorporation of the relevant processes in predictive models is still rather poor. The work reported here examines the tropical radiative and water cycle behavior as produced by four contemporary climate models. Among these are the NSIPP-2 (NASA Seasonal to Interannual Prediction Project) which uses the RAS convective parameterization; the FVCCM, a code using finite volume numerics and the CCM3.6 physics; FVCCM-MCRAS again having the finite volume numerics, but MCRAS convective parameterization and a different radiation treatment; and, finally, the NCEP GSM which uses the RAS. Using multi-decadal integrations with specified SSTs we examine the statistics of radiative / convective processes and associated energy transports, and then estimate model energy flux sensitivities to SST changes. In particular the behavior of the convective parameterizations is investigated. Additional model integrations are performed specifically to assess the importance representing convective inhibition in regulating convective cloud-top structure and moisture detrainment as well as controlling surface energy fluxes. To evaluate the results of these experiments, a number of satellite retrievals are used: TRMM retrievals of vertical reflectivity structure, rainfall rate, and inferred diabatic heating are analyzed to show both seasonal and interannual variations in vertical structure of latent heat release. Top-of-atmosphere radiative fluxes from ERBS and CERES are used to examine shortwave and longwave cloud forcing and to deduce required seasonal energy transports. Retrievals of cloud properties from ISCCP and water vapor variations from SSM/T-2 are also used to understand behavior of the humidity fields. These observations are supplemented with output form the DOE Reanalysis-2.

Robertson, Franklin R.↗

Ground Operations of the ISS GNC Babb-Mueller Atmospheric Density Model

The ISS GNC system was updated recently with a new software release that provides onboard state determination capability. Prior to this release, only the Russian segment maintained and propagated the onboard state, which was periodically updated through Russian ground tracking. The new software gives the US segment the capability for maintaining the onboard state, and includes new GPS and state vector propagation capabilities. Part of this software package is an atmospheric density model based on the Babb-Mueller algorithm. Babb-Mueller efficiently mimics a full analytical density model, such as the Jacchia model. While lacchia is very robust and is used in the Mission Control Center, it is too computationally intensive for use onboard. Thus, Babb-Mueller was chosen as an alternative. The onboard model depends on a set of calibration coefficients that produce a curve fit to the lacchia model. The ISS GNC system only maintains one set of coefficients onboard, so a new set must be uplinked by controllers when the atmospheric conditions change. The onboard density model provides a real-time density value, which is used to calculate the drag experienced by the ISS. This drag value is then incorporated into the onboard propagation of the state vector. The propagation of the state vector, and therefore operation of the BabbMueller algorithm, will be most critical when GPS updates and secondary state vector sources fail. When GPS is active, the onboard state vector will be updated every ten seconds, so the propagation error is irrelevant. When GPS is inactive, the state vector must be updated at least every 24 hours, based on current protocol. Therefore, the Babb-Mueller coefficients must be accurate enough to fulfill the state vector accuracy requirements for at least one day. A ground operations concept was needed in order to manage both the on board Babb-Mueller density model and the onboard state quality. The Babb-Mueller coefficients can be determined operationally in two ways. The first method is to calibrate the coefficients in real-time, where a set of custom coefficients is generated for the real-time atmospheric conditions. The second approach is to generate pre-canned sets of coefficients that encompass the expected atmospheric conditions over the lifetime of the vehicle. These predetermined sets are known as occurrences. Even though a particular occurrence will not match the true atmospheric conditions, the error will be constrained by limiting the breadth of each occurrence. Both methods were investigated and the advantages and disadvantages of each were considered. The choice between these implementations was a trade-off between the additional accuracy of the real-time calibration and the simpler development for the approach using occurrences. The operations concept for the frequency of updates was also explored, and depends on the deviation in solar flux that still achieves the necessary accuracy of the coefficients. This was determined based on historical solar flux trends. This analysis resulted in an accurate and reliable implementation of the Babb-Mueller coefficients and how flight controllers use them during realtime operations.

Brogan, Jonathan↗

Trick Simulation Environment 07

The Trick Simulation Environment is a generic simulation toolkit used for constructing and running simulations. This release includes a Monte Carlo analysis simulation framework and a data analysis package. It produces all auto documentation in XML. Also, the software is capable of inserting a malfunction at any point during the simulation. Trick 07 adds variable server output options and error messaging and is capable of using and manipulating wide characters for international support. Wide character strings are available as a fundamental type for variables processed by Trick. A Trick Monte Carlo simulation uses a statistically generated, or predetermined, set of inputs to iteratively drive the simulation. Also, there is a framework in place for optimization and solution finding where developers may iteratively modify the inputs per run based on some analysis of the outputs. The data analysis package is capable of reading data from external simulation packages such as MATLAB and Octave, as well as the common comma-separated values (CSV) format used by Excel, without the use of external converters. The file formats for MATLAB and Octave were obtained from their documentation sets, and Trick maintains generic file readers for each format. XML tags store the fields in the Trick header comments. For header files, XML tags for structures and enumerations, and the members within are stored in the auto documentation. For source code files, XML tags for each function and the calling arguments are stored in the auto documentation. When a simulation is built, a top level XML file, which includes all of the header and source code XML auto documentation files, is created in the simulation directory. Trick 07 provides an XML to TeX converter. The converter reads in header and source code XML documentation files and converts the data to TeX labels and tables suitable for inclusion in TeX documents. A malfunction insertion capability allows users to override the value of any simulation variable, or call a malfunction job, at any time during the simulation. Users may specify conditions, use the return value of a malfunction trigger job, or manually activate a malfunction. The malfunction action may consist of executing a block of input file statements in an action block, setting simulation variable values, call a malfunction job, or turn on/off simulation jobs.

Lin, Alexander S.↗

Magnetic Field-Line Lengths in Interplanetary Coronal Mass Ejections Inferred from Energetic Electron Events

About one quarter of the observed interplanetary coronal mass ejections (ICMEs) are characterized by enhanced magnetic fields that smoothly rotate in direction over timescales of about 10-50 hr. These ICMEs have the appearance of magnetic flux ropes and are known as "magnetic clouds" (MCs). The total lengths of MC field lines can be determined using solar energetic particles of known speeds when the solar release times and the I AU onset times of the particles are known. A recent examination of about 30 near-relativistic (NR) electron events in and near 8 MCs showed no obvious indication that the field-line lengths were longest near the MC boundaries and shortest at the MC axes or outside the MCs, contrary to the expectations for a flux rope. Here we use the impulsive beamed NR electron events observed with the Electron Proton and Alpha Monitor instrument on the Advanced Composition Explorer spacecraft and type III radio bursts observed on the Wind spacecraft to determine the field-line lengths inside ICMEs included in the catalog of Richardson & Cane. In particular, we extend this technique to ICMEs that are not MCs and compare the field-line lengths inside MCs and non-MC ICMEs with those in the ambient solar wind outside the ICMEs. No significant differences of field-line lengths are found among MCs, ICMEs, and the ambient solar wind. The estimated number of ICME field-line turns is generally smaller than those deduced for flux-rope model fits to MCs. We also find cases in which the electron injections occur in solar active regions CARs) distant from the source ARs of the ICMEs, supporting CME models that require extensive coronal magnetic reconnection with surrounding fields. The field-line lengths are found to be statistically longer for the NR electron events classified as ramps and interpreted as shock injections somewhat delayed from the type III bursts. The path lengths of the remaining spike and pulse electron events are compared with model calculations of solar wind field-line lengths resulting from turbulence and found to be in good agreement.

Kahler, S. W.↗

The Cretaceous-Tertiary Impact Crater and the Cosmic Projectile that Produced it

Evidence gathered to date from topographic data, geophysical data, well logs, and drill-core samples indicates that the buried Chicxulub basin, the source crater for the approximately 65 Ma Cretaceous-Tertiary (K/T) boundary deposits, is approximately 300 km in diameter. A prominent topographic ridge and a ring of gravity anomalies mark the position of the basin rim at approximately 150 km from the center. Wells in this region recovered thick sequences of impact-generated breccias at 200-300 m below present sea level. Inside the rim, which has been severely modified by erosion following impact, the subsurface basin continues to deepen until near the center it is approximately 1 km deep. The best planetary analog for this crater appears to be the 270 km-diameter Mead basin on Venus. Seismic reflection data indicate that the central zone of downward displacement and excavation (the transient crater is approximately 130 km in diameter, consistent with previous studies of gravity anomaly data). Our analysis of projectile characteristics utilizes this information, coupled with conventional scaling relationships, and geochemical constraints on the mass of extraterrestrial material deposited within the K/T boundary layer. Results indicate that the Chicxulub crater would most likely be formed by a long-period comet composed primarily of nonsilicate materials (ice, hydrocarbons, etc.) and subordinate amounts (less than or equal to 50 percent) primitive chondritic material. This collision would have released the energy equivalent to between 4 x 10(exp 8) and 4 x 10(exp 9) megatons of TNT. Studies of terrestrial impact rates suggest that such an event would have a mean production rate of approximately 1.25 x 10(exp -9) y(exp -1). This rate is considerably lower than that of the major mass extinctions over the last 250 million years (approximately 5 x 10(exp -7) y(exp -1). Consequently, while there is substantial circumstantial evidence establishing the cause-effect link between the Chicxulub basin forming event and the K/T biological extinctions, the results of our analysis do not support models of impact as a common or singular causative agent of mass extinctions on Earth.

Sharpton, Virgil L.↗