Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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,045 records · Page 58

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗

Modeling and Simulation of Tank Pressure Control using Zero-Boiloff Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗

Basalt generation at the Apollo 12 site. Part 1: New data, classification, and re-evaluation

New data are reported from five previously unanalyzed Apollo 12 mare basalts that are incorporated into an evaluation of previous petrogenetic models and classification schemes for these basalts. This paper proposes a classification for Apollo 12 mare basalts on the basis of whole-rock Mg# (molar 100*(Mg/(Mg+Fe))) and Rb/Sr ratio (analyzed by isotope dilution), whereby the ilmenite, olivine, and pigeonite basalt groups are readily distinguished from each other. Scrutiny of the Apollo 12 feldspathic 'suite' demonstrates that two of the three basalts previously assigned to this group (12031, 12038, 12072) can be reclassified: 12031 is a plagioclase-rich pigeonite basalt; and 12072 is an olivine basalt. Only basalt 12038 stands out as a unique sample to the Apollo 12 site, but whether this represents a single sample from another flow at the Apollo 12 site or is exotic to this site is equivocal. The question of whether the olivine and pigeonite basalt suites are co-magmatic is addressed by incompatible trace-element chemistry: the trends defined by these two suites when Co/Sm and Sm/Eu ratios are plotted against Rb/Sr ratio demonstrate that these two basaltic types cannot be co-magmatic. Crystal fractionation/accumulation paths have been calculated and show that neither the pigeonite, olivine, or ilmenite basalts are related by this process. Each suite requires a distinct and separate source region. This study also examines sample heterogeneity and the degree to which whole-rock analyses are representative, which is critical when petrogenetic interpretation is undertaken. Sample heterogeneity has been investigated petrographically (inhomogeneous mineral distribution) with consideration of duplicate analyses, and whether a specific sample (using average data) plots consistently upon a fractionation trend when a number of different compostional parameters are considered. Using these criteria, four basalts have been identified where reported analyses are not representative of the whole-rock composition: 12005, an ilmenite basalt; 12006 and 12036, olivine basalts; and 12031 previously classified as a feldspathic basalt, but reclassified as part of the pigeonite suite.

Neal, Clive R.↗

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗

Kinetics and Mechanism of the Reaction of Hydoxyl Radicals with Acetonitrile under Atmospheric Conditions

The pulsed laser photolysis-pulsed laser induced fluorescence technique has been employed to determine absolute rate coefficients for the reaction OH + CH3CN (1) and its isotopic variants, OH + CD3CN (2), OD + CH3CN (3), and OD + CD3CN (4). Reactions 1 and 2 were studied as a function of pressure and temperature in N2, N2/O2, and He buffer gases. In the absence of O2 all four reactions displayed well-behaved kinetics with exponential OH decays and pseudo-first rate constants which were proportional to substrate concentration. Data obtained in N2 over the range 50-700 Torr at 298 K are consistent with k(sub 1), showing a small pressure dependence. The Arrhenius expression obtained by averaging data at all pressures in k(sub 1)(T) = (1.1(sup +0.5)/(sub -0.3)) x 10(exp -12) exp[(-1130 +/- 90)/T] cu cm /(molecule s). The kinetics of reaction 2 are found to be pressure dependent with k(sub 2) (298 K) increasing from (1.21 +/- 0.12) x 10(exp -14) to (2.16 +/- 0.11) x 10(exp -14) cm(exp 3)/ (molecule s) over the pressure range 50-700 Torr of N2 at 298 K. Data at pressures greater than 600 Torr give k(sub 2)(T) = (9.4((sup +13.4)(sub -5.0))) x 10(exp -13) exp[(-1180 +/- 250)/T] cu cm/(molecule s). The rates of reactions 3 and 4 are found to be independent of pressure over the range 50-700 Torr of N2 with 298 K rate coefficient given by k(sub 3) =(3.18 +/- 0.40) x 10(exp -14) cu cm/(molecule s) and k(sub 4) = (2.25 +/-0.28) x 10(exp -14) cu cm/(molecule s). In the presence of O2 each reaction shows complex (non-pseudo-first-order) kinetic behavior and/or an apparent decrease in the observed rate constant with increasing [O2], indicating the presence of significant OH or OD regeneration. Observation of regeneration of OH in (2) and OD in (3) is indicative of a reaction channel which proceeds via addition followed by reaction of the adduct, or one of its decomposition products, with O2. The observed OH and OD decay profiles have been modeled by using a simple mechanistic scheme to extract kinetic information about the adduct reations with O2 and branching ratios for OH regeneration. A plausible mechanism for OH regeneration in (2) involves OH addition to the nitrogen atom followed by O2 addition to the cyano carbon atom, isomeriazation and decomposition to D2CO + DOCN + OH. Our results suggest that the OH + CH3CN reaction occurs via a complex mechanism involving both bimolecular and termolecular pathways, analogous to the mechanisms for the the important atmospheric reactions of OH with CO and HNO3.

Hynes, A. J.↗

Modal identification experiment

The Modal Identification Experiment (MIE) is a proposed on-orbit experiment being developed by NASA's Office of Aeronautics and Space Technology wherein a series of vibration measurements would be made on various configurations of Space Station Freedom (SSF) during its on-orbit assembly phase. The experiment is to be conducted in conjunction with station reboost operations and consists of measuring the dynamic responses of the spacecraft produced by station-based attitude control system and reboost thrusters, recording and transmitting the data, and processing the data on the ground to identify the natural frequencies, damping factors, and shapes of significant vibratory modes. The experiment would likely be a part of the Space Station on-orbit verification. Basic research objectives of MIE are to evaluate and improve methods for analytically modeling large space structures, to develop techniques for performing in-space modal testing, and to validate candidate techniques for in-space modal identification. From an engineering point of view, MIE will provide the first opportunity to obtain vibration data for the fully-assembled structure because SSF is too large and too flexible to be tested as a single unit on the ground. Such full-system data is essential for validating the analytical model of SSF which would be used in any engineering efforts associated with structural or control system changes that might be made to the station as missions evolve over time. Extensive analytical simulations of on-orbit tests, as well exploratory laboratory simulations using small-scale models, have been conducted in-house and under contract to develop a measurement plan and evaluate its potential performance. In particular, performance trade and parametric studies conducted as part of these simulations were used to resolve issues related to the number and location of the measurements, the type of excitation, data acquisition and data processing, effects of noise and nonlinearities, selection of target vibration modes, and the appropriate type of data analysis scheme. The purpose of this talk is to provide an executive-summary-type overview of the modal identification experiment which has emerged from the conceptual design studies conducted to-date. Emphasis throughout is on those aspects of the experiment which should be of interest to those attending the subject utilization conference. The presentation begins with some preparatory remarks to provide background and motivation for the experiment, describe the experiment in general terms, and cite the specific technical objectives. This is followed by a summary of the major results of the conceptual design studies conducted to define the baseline experiment. The baseline experiment which has resulted from the studies is then described.

Kvaternik, Raymond G.↗

Space Weathering in the Fine Size Fractions of Lunar Soils: Soil Maturity Effects

The effects of space weathering on the optical properties of lunar materials have been well documented. These effects include a reddened continuum slope, lowered albedo, and attenuated absorption features in reflectance spectra of lunar soils as compared to finely comminuted rocks from the same Apollo sites. However, the regolith processes that cause these effects are not well known, nor is the petrographic setting of the products of these processes fully understood. A Lunar Soil Characterization Consortium has been formed with the purpose of systematically integrating chemical and mineralogical data with the optical properties of lunar soils. Understanding space-weathering effects is critical in order to fully integrate the lunar sample collection with remotely-sensed data from recent robotic missions (e.g., Lunar Prospector, Clementine, and Galileo) We have shown that depositional processes (condensation of impact-derived vapors, sputter deposits, accreted impact material, e.g., splash glass, spherules, etc.) are a major factor in the modification of the optical surfaces of lunar regolith materials. In mature soils, it is the size and distribution of the nanophase metal in the soil grains that has the major effect on optical properties. In this report, we compare and contrast the space-weathering effects in an immature and a mature soil with similar elemental compositions. For this study, we analyzed <10 micron sieve fractions of two Apollo 17 soils, 79221 (mature, Is/FeO = 81) and 71061 (immature, Is/FeO = 14). Details of the sieving procedures and allocation scheme are given else where. The results of other detailed chemical, mineralogical, and spectroscopic analyses of these soil samples are reported elsewhere. A representative sample of each soil was embedded in low-viscosity epoxy, and thin sections (about 70nm thick) were obtained through ultra microtomy. The thin sections used for these analyses typically contained cross sections of up to 500 individual grains. The thin sections were studied using a JEOL 2010 transmission electron microscope (TEM) equipped with a thin window energy-dispersive X-ray (EDX) spectrometer. An individual thin section was selected from each soil, and for each grain in the section we determined (1) the elemental composition by EDX; (2) whether the grain was crystalline or glassy using electron diffraction and darkfield imaging; (3) the presence or absence of rims and accreted material; and (4) the distribution of nanophase Fe where present. Most of the categories are self-evident; however, we divide the agglutinate derived material into agglutinitic glass (glass with approximately the same composition as the bulk soil that contains nanophase Fe with or without vesicles) and agglutinate fragments, which are composed of crystalline grains and agglutinitic glass. Lithic fragments are defined as polymineralic grains with no glass. Pyroxene grains have been divided into high- and low-Ca groups. As expected, there are a number of differences in the petrography of the <10-microns fractions of 79221 and 71061 given the great difference in their respective maturities, but we focus here on two major distinctions: agglutinate content and the number of grains with micropatina. Slightly over 50% of the particles in 79221 consist of agglutinitic glass and agglutinate fragments, while the remainder are predominantly crystalline mineral grains. The agglutinic glass particles contain abundant nanophase Fe and vesicles. Angular particles are rare, with most showing smooth, rounded exteriors, Of the mineral grains analyzed thus far, over 90% of the grains have amorphous rims that contain nanophase Fe (these rims are believed to have formed by vapor deposition and irradiation effects). The nanophase Fe in these rims probably accounts for a significant fraction of the increase in Is/FeO measured in these size fractions. In addition to the rims, the majority of particles also show abundant accreted material in the form of glass splashes and spherules that also contain nanophase Fe. In stark contrast, the surfaces of the mineral grains in the 71061 sample are relatively prisitine, as only about 14% of the mineral grains in the sample exhibited amorphous rims. Furthermore, the mineral particles are more angular and show greater surface roughness than in the mature sample. Accreted material on particle surfaces is rare. Agglutinitic material is a major component of the 71061 sample; however, nanophase Fe and vesicles are not as well developed as in the 79221 sample. It is now recognized that nanophase Fe is probably the main agent in modifying the optical properties of lunar soil grains. The most important result of this study is the observation that in the fine size fractions of mature soils, nearly every grain has nanophase Fe within 100 run of the particle surface. (Additional Information contained in original)

Keller, L. P.↗

Ion Storage Tests with the High Performance Antimatter Trap (HiPAT)

The NASA/Marshall Space Flight Centers (NASA/MSFC) Propulsion Research Center (PRC) is evaluating an antiproton storage system, referred to as the High Performance Antiproton Trap (HiPAT). This interest stems from the sheer energy represented by matter/antimatter annihilation process with has an energy density approximately 10 order of magnitude above that of chemical propellants. In other terms, one gram of antiprotons contains the equivalent energy of approximately 23 space shuttle external tanks or ET's (each ET contains roughly 740,000 kgs of fuel and oxidizer). This incredible source of stored energy, if harnessed, would be an enabling technology for deep space mission where both spacecraft weight and propulsion performance are key to satisfying aggressive mission requirements. The HiPAT hardware consists of a 4 Tesla superconductor system, an ultra high vacuum test section (vacuum approaching 10(exp -12) torr), and a high voltage confinement electrode system (up to 20 kvolts operation). The current laboratory layout is illustrated. The HiPAT designed objectives included storage of up to 1 trillion antiprotons with corresponding lifetimes approaching 18 days. To date, testing has centered on the storage of positive hydrogen ions produced in situ by a stream of high-energy electrons that passes through the trapping region. However, due to space charge issues and electron beam compression as it passes through the HiPAT central field, current ion production is limited to less then 50,000 ions. Ion lifetime was determined by counting particle populations at the end of various storage time intervals. Particle detection was accomplished by destructively expelling the ions against a micro-channel plate located just outside the traps magnetic field. The effect of radio frequency (RF) stabilization on the lifetime of trapped particles was also examined. This technique, referred to as a rotating wall, made use of a segmented electrode located near the center of the trap on which various phases of a particular frequency were applied. Various experiments were performed illustrating the ability of an RF drive to both prolong and reduced the lifetimes of various ion species depending on the selected frequency. HiPAT is now being reconfigured for testing with an ion source that will provide both positive and negative hydrogen ions from an external source. This ion system shall provide higher fill capacity (order of million of ions per shot), stacking of multiple shots, and injection schemes typical of a realistic antiproton delivery system.

Martin, James J.↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Satellite Ocean-Color Validation Using Ships of Opportunity

The investigation s main objective is to collect from platforms of opportunity (merchant ships, research vessels) concomitant normalized water-leaving radiance and aerosol optical thickness data over the world s oceans. A global, long-term data set of these variables is needed to verify whether satellite retrievals of normalized water-leaving radiance are within acceptable error limits and, eventually, to adjust atmospheric correction schemes. To achieve this objective, volunteer officers, technicians, and scientists onboard the selected ships collect data from portable SIMBAD and Advanced SIMBAD (SIMBADA) radiometers. These instruments are specifically designed for evaluation of satellite-derived ocean color. They measure radiance in spectral bands typical of ocean-color sensors. The SIMBAD version measures in 5 spectral bands centered at 443, 490, 560, 670, and 870 nm, and the Advanced SIMBAD version in 11 spectral bands centered at 350, 380, 412, 443, 490, 510, 565, 620, 670, 750, and 870 nm. Aerosol optical thickness is obtained by viewing the sun disk like a classic sun photometer. Normalized water-leaving radiance, or marine reflectance, is obtained by viewing the ocean surface through a vertical polarizer in a specific geometry (nadir angle of 45o and relative azimuth angle of 135deg) to minimize direct sun glint and reflected sky radiation. The SIMBAD and SIMBADA data, after proper quality control and processing, are delivered to the SIMBIOS project office for inclusion in the SeaBASS archive. They complement data collected in a similar way by the Laboratoire d'Optique Atmospherique of the University of Lille, France. The SIMBAD and SIMBADA data are used to check the radiometric calibration of satellite ocean-color sensors after launch and to evaluate derived ocean-color variables (i.e., normalized water-leaving radiance, aerosol optical thickness, and aerosol type). Analysis of the SIMBAD and SIMBADA data provides information on the accuracy of satellite retrievals of normalized water-leaving radiance, an understanding of the discrepancies between satellite and in situ data, and algorithms that reduce the discrepancies, contributing to more accurate and consistent global ocean color data sets.

Frouin, Robert↗

Winds Measured by the Rover Environmental Monitoring Station (REMS) During the Mars Science Laboratory (MSL) Rover's Bagnold Dunes Campaign and Comparison with Numerical Modeling Using MarsWRF

A high density of REMS wind measurements were collected in three science investigations during MSL's Bagnold Dunes Campaign, which took place over approx. 80 sols around southern winter solstice (Ls approx. 90deg) and constituted the first in situ analysis of the environmental conditions, morphology, structure, and composition of an active dune field on Mars. The Wind Characterization Investigation was designed to fully characterize the near-surface wind field just outside the dunes and confirmed the primarily upslope/downslope flow expected from theory and modeling of the circulation on the slopes of Aeolis Mons in this season. The basic pattern of winds is 'upslope' (from the northwest, heading up Aeolis Mons) during the daytime (approx. 09:00-17:00 or 18:00) and 'downslope' (from the southeast, heading down Aeolis Mons) at night (approx. 20:00 to some time before 08:00). Between these times the wind rotates largely clockwise, giving generally westerly winds mid-morning and easterly winds in the early evening. The timings of these direction changes are relatively consistent from sol to sol; however, the wind direction and speed at any given time shows considerable intersol variability. This pattern and timing is similar to predictions from the MarsWRF numerical model, run at a resolution of approx. 490 m in this region, although the model predicts the upslope winds to have a stronger component from the E than the W, misses a wind speed peak at approx. 09:00, and under-predicts the strength of daytime wind speeds by approx. 2-4 m/s. The Namib Dune Lee Investigation reveals 'blocking' of northerly winds by the dune, leaving primarily a westerly component to the daytime winds, and also shows a broadening of the 1 Hz wind speed distribution likely associated with lee turbulence. The Namib Dune Side Investigation measured primarily daytime winds at the side of the same dune, in support of aeolian change detection experiments designed to put limits on the saltation threshold, and also appears to show the influence of the dune body on the local flow, though less clearly than in the lee. Using a vertical grid with lower resolution near the surface reduces the relative strength of nighttime winds predicted by MarsWRF and produces a peak in wind speed at approx. 09:00, improving the match to the observed diurnal variation of wind speed, albeit with an offset in magnitude. The annual wind field predicted using this grid also provides a far better match to observations of aeolian dune morphology and motion in the Bagnold Dunes. However, the lower overall wind speeds than observed and disagreement with the observed wind direction at approx. 09:00 suggest that the problem has not been solved and that alternative boundary layer mixing schemes should be explored which may result in more mixing of momentum down to the near-surface from higher layers. These results demonstrate a strong need for in situ wind data to constrain the setup and assumptions used in numerical models, so that they may be used with more confidence to predict the circulation at other times and locations on Mars.

Newman, Claire E.↗

An Adaptable Power System with Software Control Algorithm

A low cost, flexible and modular spacecraft power system design was developed in response to a call for an architecture that could accommodate multiple missions in the small to medium load range. Three upcoming satellites will use this design, with one launch date in 1999 and two in the year 2000. The design consists of modular hardware that can be scaled up or down, without additional cost, to suit missions in the 200 to 600 Watt orbital average load range. The design will be applied to satellite orbits that are circular, polar elliptical and a libration point orbit. Mission unique adaptations are accomplished in software and firmware. In designing this advanced, adaptable power system, the major goals were reduction in weight volume and cost. This power system design represents reductions in weight of 78 percent, volume of 86 percent and cost of 65 percent from previous comparable systems. The efforts to miniaturize the electronics without sacrificing performance has created streamlined power electronics with control functions residing in the system microprocessor. The power system design can handle any battery size up to 50 Amp-hour and any battery technology. The three current implementations will use both nickel cadmium and nickel hydrogen batteries ranging in size from 21 to 50 Amp-hours. Multiple batteries can be used by adding another battery module. Any solar cell technology can be used and various array layouts can be incorporated with no change in Power System Electronics (PSE) hardware. Other features of the design are the standardized interfaces between cards and subsystems and immunity to radiation effects up to 30 krad Total Ionizing Dose (TID) and 35 Mev/cm(exp 2)-kg for Single Event Effects (SEE). The control algorithm for the power system resides in a radiation-hardened microprocessor. A table driven software design allows for flexibility in mission specific requirements. By storing critical power system constants in memory, modifying the system code for other programs is simple. These constants can be altered also by ground command, or in response to an anomolous event. All critical power system functions have backup hardware functions to prevent a software or computer glitch from propagating. A number of battery charge control schemes can be implemented by selecting the proper control terms in the code. The architecture allows the design engineer to tune the system response to various system components and anticipated load profiles without costly alterations. A design trade was made with the size, weight and power dissipation of the electronics versus the performance of the power bus to load variations. Linear, fine control is maintained with a streamlined electronics design. This paper describes the hardware design as well as the software control algorithm. The challenges of closing the system control loop digitally is discussed. Control loop margin and power system performance is presented. Lab measurements are shown and compared to the system response of a hardware model running actual flight software.

Castell, Karen↗

A Global End-Member Approach to Derive aCDOM(440) from Near-Surface Optical Measurements

This study establishes an optical inversion scheme for deriving the absorption coefficient of colored (or chromophoric, depending on the literature) dissolved organic material (CDOM) at the 440 nm wavelength, which can be applied to global water masses with near-equal efficacy. The approach uses a ratio of diffuse attenuation coefficient spectral end members, i.e., a short and long wavelength pair. The global perspective is established by sampling "extremely" clear water plus a generalized extent in turbidity and optical properties that each span three decades of dynamic range. A unique data set was collected in oceanic, coastal, and inland waters (as shallow as 0.6 m) from the North Pacific Ocean, the Arctic Ocean, Hawaii, Japan, Puerto Rico, and the east and west coasts of the United States. The data were partitioned using subjective categorizations to define a validation quality subset of conservative water masses, i.e., the inflow and outflow of properties constrain the range in the gradient of a constituent, plus 15 subcategories of water masses that were not evolving conservatively. The dependence on subcategories was confirmed with an objective methodology based on cluster analysis techniques. The latter defined five distinct classes with validation quality data present in all classes, but which also decreased in percent composition as a function of increasing class number and optical complexity. Four different algorithms based on different validation quality end members were validated with accuracies of 1.–6.2 %, wherein pairs with the largest spectral span were most accurate. Although algorithm accuracy decreased with the inclusion of more subcategories containing non-conservative water masses, changes to the algorithm fit were small when a preponderance of subcategories were included. The high accuracy for all end-member algorithms was the result of data acquisition and data processing improvements, e.g., increased vertical sampling resolution to less than 1mm and a boundary constraint to mitigate wave focusing effects, respectively. An independent evaluation with a historical database confirmed the consistency of the algorithmic approach and its application to quality assurance, e.g., to flag data outside expected ranges, identify suspect spectra, and objectively determine the in-water extrapolation interval by converging agreement for all applicable end-member algorithms. The legacy data exhibit degraded performance (as 44 % uncertainty) due to a lack of high-quality near-surface observations, especially for clear waters wherein wave-focusing effects are problematic. The novel optical approach allows the in situ estimation of an in-water constituent in keeping with the accuracy obtained in the laboratory.

Stanford B Hooker↗

Acceleration Modes and Transitions in Pulsed Plasma Accelerators

Pulsed plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas, ionizing and accelerating it through the Lorentz body force. Two plasma accelerator types employing this general scheme have typically been studied: the gas-fed pulsed plasma thruster and the quasi-steady magnetoplasmadynamic (MPD) accelerator. The gas-fed pulsed plasma accelerator is generally represented as a completely transient device discharging in approximately 1-10 microseconds. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j (current density) by B (magnetic field) body force, entraining propellant it encounters. This process is sometimes referred to as detonation-mode acceleration because the current sheet representation approximates that of a strong shock propagating through the gas. Acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new current sheet at the breech and depriving the initial sheet of additional acceleration. In the quasi-steady MPD accelerator, the pulse is lengthened to approximately 1 millisecond or longer and maintained at an approximately constant level during discharge. The time over which the transient phenomena experienced during startup typically occur is short relative to the overall discharge time, which is now long enough for the plasma to assume a relatively steady-state configuration. The ionized gas flows through a stationary current channel in a manner that is sometimes referred to as the deflagration-mode of operation. The plasma experiences electromagnetic acceleration as it flows through the current channel towards the exit of the device. A device that had a short pulse length but appeared to operate in a plasma acceleration regime different from the gas-fed pulsed plasma accelerators was developed by Cheng, et al. The Coaxial High ENerGy (CHENG) thruster operated on the 10-microseconds timescales of pulsed plasma thrusters, but claimed high thrust density, high efficiency and low electrode erosion rates, which are more consistent with the deflagration mode of acceleration. Separate work on gas-fed pulsed plasma thrusters (PPTs) by Ziemer, et al. identified two separate regimes of performance. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency (ratio of jet kinetic energy to input electrical energy) as a function of mass bit. In the second regime at very low mass bits (termed Mode II), the efficiency increased with decreasing mass bit. Work by Poehlmann et al. and by Sitaraman and Raja sought to understand the performance of the CHENG thruster and the Mode I / Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes representing two distinct sets of solutions to the relevant conservation laws. These works studied the proposal that, depending upon the values of the various controllable parameters, the accelerator would operate in either the detonation or deflagration mode. In the present work, we propose a variation on the explanation for the differences in performance between the various pulsed plasma accelerators. Instead of treating the accelerator as if it were only operating in one mode or the other during a pulse, we model the initial stage of the discharge in all cases as an accelerating current sheet (detonation mode). If the current sheet reaches the exit of the accelerator before the discharge is completed, the acceleration mode transitions to the deflagration mode type found in the quasi-steady MPD thrusters. This modeling method is used to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗

Beneficiation of ilmenite from lumar analog

The results reported were obtained on a meteoric eucrite sample called Millbillillie Sample no. 173. Optical microscopy studies of the sample showed it to consist of ilmenite, troilite, and transparent gangue. The transparent gangue consisted of feldspar (anorthite), pyroxenes, olivines, and opaques. Troilite was present in minor quantities. Screen assay analyses of the 30, 100, 200, and 400 US mesh screen fractions showed that minor concentration of titanium occurred in the 200 x 400 and -400 mesh screen fractions. Scanning electron microscopy (SEM) studies of the bulk sample showed the presence of a variety of ilmenite grains, ranging from 50 microns down to less than 1 micron without any evidence of liberation. Electron Diffraction Scans (EDS) confirmed the ratio of Fe to Ti in the ilmenite grains. Dry magnetic separation in a Frantz Isodynamic Separator was found to be effective only at sizes finer than 150 microns (100 US mesh) and more so at 200 mesh (74 microns). In each case, dedusting of the sample to remove -400 mesh (-0.037 microns) fines was required. Liberation size was determined to be 200 mesh and finer. The highest grade concentrate assaying 3.45 percent Ti was produced by magnetic separation of the -200 + 400 mesh screen fraction assaying 0.44 Ti (from a -30 mesh sample) at a current setting of 0.35 AMP. This concentrate contained 21.2 percent of the Ti values in the screen fraction with 2.72 weight percent of feed to test. The results can be projected to a sample stage ground to -200 mesh. Magnetic separation of the 200 + 400 mesh (-0.074 + 0.037 microns) should produce a concentrate accounting for 1.41 weight percent of the feed. This concentrate will analyze 3.45 percent Ti and contain 10.3 percent of the Ti values in the feed. By changing the Frantz Magnetic Separator settings, a lower grade concentrate analyzing 0.98 percent Ti can be produced at an increased recovery of 25.4 percent. The concentrate weight will be 11.7 percent of the feed. It must be emphasized that improved grades and recoveries can be obtained with the -400 mesh fines. However, beneficiation of these extremely fine materials is not possible in a practical process scheme.

Ramadorai, G.↗

Deflections of anisotropic sandwich beams with variable face sheets and core thicknesses

A sandwich construction consists of a low-density core material with high strength face sheets bounded to the top and bottom surfaces. The construction has been widely used in the aerospace and marine industries due to its outstanding characteristics such as noise absorption, weight minimization, heat insulation, and better bending stiffness. In sandwich structures used in high-performance aircraft, the face sheets are often made of fiber-reinforced composite materials and the core is made of honeycomb. The structures may also have variable thickness so as to satisfy aerodynamic requirements. In the stress analysis, the constant-thickness face sheets are usually considered as membrane and the core is assumed to be inextensible but deformable in the thickness direction. The static behavior of variable-thickness, isotropic and homogeneous sandwich beams was successfully studied by employing a constant-thickness theory but allowing stiffnesses to vary in accordance with local thickness variations. It has been recently found in a refined theory that the analyses based on the constant thickness theory locally can lead to significant errors in structural responses if the sandwich beam is thickness-tapered and the cores are deformable in transverse shear. The errors arise mainly from two factors: (1) the transverse shear components of the membrane forces in the face sheets alter the transverse shears carried by the core; and (2) the face-sheet membrane strains arise from transverse shear deformation of the core. In practice the variable thickness may not only exist in core but also in face sheets. The thickness-variations may even be a type of step function. In this case the transverse shear stress in the face sheets and bending stress in the core should be taken into account in the refined theory mentioned. In the present study, energy principles are employed in deriving governing equations for general bending of anisotropic sandwich beams with variable thickness in both face sheets and cores. Solutions to these equations are based on a finite difference scheme. As an example in application, a simply supported thickness-tapered sandwich beam subject to a concentrated load at its center is considered. Let W' be the maximum deflection of the beam in which face sheets are considered as membrane, while W'' is that based on using the modified refined theory. It is found that W' is always larger than W'', however, the magnitude of (W'- W'') appears to be insensitive to the change of the taper of the beam.

Lu, Chu-Ho↗

Thermal Performance of ATLAS Laser Thermal Control System Demonstration Unit

The second Ice, Cloud, and Land Elevation Satellite mission currently planned by National Aeronautics and Space Administration will measure global ice topography and canopy height using the Advanced Topographic Laser Altimeter System {ATLAS). The ATLAS comprises two lasers; but only one will be used at a time. Each laser will generate between 125 watts and 250 watts of heat, and each laser has its own optimal operating temperature that must be maintained within plus or minus 1 degree Centigrade accuracy by the Laser Thermal Control System (LTCS) consisting of a constant conductance heat pipe (CCHP), a loop heat pipe (LHP) and a radiator. The heat generated by the laser is acquired by the CCHP and transferred to the LHP, which delivers the heat to the radiator for ultimate rejection. The radiator can be exposed to temperatures between minus 71 degrees Centigrade and minus 93 degrees Centigrade. The two lasers can have different operating temperatures varying between plus 15 degrees Centigrade and plus 30 degrees Centigrade, and their operating temperatures are not known while the LTCS is being designed and built. Major challenges of the LTCS include: 1) A single thermal control system must maintain the ATLAS at 15 degrees Centigrade with 250 watts heat load and minus 71 degrees Centigrade radiator sink temperature, and maintain the ATLAS at plus 30 degrees Centigrade with 125 watts heat load and minus 93 degrees Centigrade radiator sink temperature. Furthermore, the LTCS must be qualification tested to maintain the ATLAS between plus 10 degrees Centigrade and plus 35 degrees Centigrade. 2) The LTCS must be shut down to ensure that the ATLAS can be maintained above its lowest desirable temperature of minus 2 degrees Centigrade during the survival mode. No software control algorithm for LTCS can be activated during survival and only thermostats can be used. 3) The radiator must be kept above minus 65 degrees Centigrade to prevent ammonia from freezing using no more than 135 watts of heater power. 4) The LHP reservoir control heater power is limited to 15 watts with a 70 percent duty cycle. 5) The voltage of the power supply can vary between 26 volts direct current and 34 volts direct current during the spacecraft lifetime. A design analysis shows that a single LTCS can satisfy these requirements. However, shutdown of· the LHP is particularly challenging and the shutdown heater must be wired in series with two reservoir thermostats and two CCHP thermostats at different set points. An LTCS demonstration unit has been tested to verify these performance characteristics experimentally prior to proceeding to the final LTCS design and fabrication. Test results showed that the LHP shutdown scheme would be able to shut down the LHP as designed and the reservoir control heater can maintain the ATLAS mass simulator within the plus or minus 1 degrees Centigrade accuracy under various combinations of the heat load, sink temperature, and power supply voltage.

Ku, Jentung↗