Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermal Energy Distribution System”

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.

299 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Grant Proposal for the Continuation of the Voyager Interstellar Mission: LECP Investigation

This proposal documents the plans of the Low Energy Charged Particle (LECP) investigation team for participation in NASA's Voyager Interstellar Mission (VIM) as the Voyager 1 and 2 spacecraft explore the outer reaches of the heliosphere and search for the termination shock and the heliopause. The proposal covers the four year period from 1 January 1997 to 31 December 2000. The LECP instruments on Voyager 1 and 2 measure in situ intensities of charged particles with energies from about 30 keV to 100 MeV for ions, and about 20 keV to greater than 10 MeV for electrons. The instruments provide detailed spectral, angular, and compositional information about the particles. Composition is available for greater than 200 keV/nuc using multi-parameter measurements. Angular information is obtained by a mechanically scanned platform that rotates at various commanded rates. Measurements of low energy ion and electron intensities versus time and spatial location within the heliosphere contain an abundance of information regarding various transport and acceleration processes on both local (approx. 1 hr, approx. 0.01 AU) and global (approx. 11 yrs, approx. 100 AU) scales. The LECP instruments provide unique observations of such dynamical processes, and we anticipate that it will return critical information regarding the boundaries of the heliosphere. Several recent and exciting discoveries based on LECP measurements emphasize the important role that low energy charged particle distributions play in physical processes in the interplanetary medium. Yet, at the same time, these discoveries also underscore the fact that our understanding of processes in the outer heliosphere is, in most cases, incomplete, and in others, only rudimentary at best. Among the discoveries referred to above are the following: (1) Shocks: Examination of greater than 30 keV ion intensities have revealed: (a) a total absence of acceleration beyond only -100-200 keV at a strong transient shock in May 1991 at 35 AU, despite an enhanced level of seed particles; (b) a large transient shock in September 1991 of global scale, with intensities of shock-accelerated ions greater than or equal to 30 keV to approx. 30 MeV showing complex, highly energy-dependent spatial evolution, and small-scale (approx. few gyroradii), often anisotropic, micro-structures; (c) recurrent intensity increases in greater than or equal to 30 keV to -few MeV ions, with structures that, in some cases, show no correlation with the associated corotating shock. (2) Superthermal ion pressure: A global merged interaction region with a leading shock, downstream of which the superthermal ion (greater than or equal to 30 keV to approx. 4 MeV) pressure is comparable to that of the thermal plasma, and the total particle pressure yields a plasma beta of order unity. (3) Pickup ions: Measurements of the C/O ratio within transient structures at 35-45 AU showing the first clear evidence that transient shocks can pre-accelerate interstellar pickup ions from approx. 1 keV/nuc to at least 1 MeV/nuc. (4) Seed particles: Injection of ions for acceleration to high energies at the termination shock is unlikely to be a problem, since interplanetary transient and recurrent shocks are continually accelerating ions, of solar wind or interstellar origin, to highly superthermal energies. (5) Precursor electrons: Ambient solar electrons (greater than or equal to few tens of keV) that exist in the outer heliosphere ca form a broad precursor, several days wide, that is upstream of the termination shock and potentially observable a few months prior to the shock crossing. (6) Solar wind velocity at Voyager 1: We can use LECP ion data to obtain the solar wind velocity at Voyager 1, enabling us to provide critical measurement of the plasma flow as we approach and encounter the termination shock and other regions (necessary due to the partial failure of the Voyager 1 PLS experiment). The work of the LECP investigator team during the VIM will include: (1) Continuing operations with regard to the receipt, processing, verification, cataloging, display, and distribution of the data from the LECP instruments on Voyager 1 and 2, (2) Monitoring the health and performance of the LECP instruments, and evaluating and characterizing the response of the LECP instruments to various energetic particle and plasma environments, (3) Participating in, and supporting Voyager Project planning exercises and other coordinated activities relevant to exploration of the outer heliosphere, (4) Developing analysis techniques and operational procedures suitable for searching for and characterizing the boundaries and unique regions of the outher heliosphere, (5) Continuing the preparation of data sets appropriate for submission to the National Space Sciences Data Center (NSSDC) and, where appropriate, the Planetary Data System (PDS), (6) Maintaining direct Web access to online LECP data through the JHU/APL Voyager LECP home page, (7) Performing scientific evaluations of the Voyager 1 and 2 LECP data sets in conjunction with other data sets and other investigators, with particular focus on the outer regions of the heliosphere, and (8) Publishing the results of these evaluations in the scientific literature and presenting the results in scientific conferences.

Krimigis, Stamatios M.↗
NASA NTRS1994

Relativistic particle acceleration in plerions

We discuss recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is either pure electrons and positrons or primarily electrons and positrons with an admixture of heavy ions. Particle-in-cell simulation techniques as well as analytic theory have been used to show that such shocks in pure pair plasmas are fully thermalized -- the downstream particle spectra are relativistic Maxwellians at the temperature expected from the jump conditions. On the other hand, shocks containing heavy ions which are a minority constituent by number but which carry most of the energy density in the upstream medium do put approximately 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) varies as E(exp -2), where N(E)dE is the number of particles with energy between E and E+dE. The mechanism of thermalization and particle acceleration is found to be synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front. The synchrotron maser modes radiated by the heavy ions are absorbed by the pairs at their (relativistic) cyclotron frequencies, allowing the maximum energy achievable by the pairs to be gamma(sub +/-)m(sub +/-)c squared = m(sub i)c squared gamma(sub 1)/Z(sub i), where gamma(sub 1) is the Lorentz factor of the upstream flow and Z(sub i) is the atomic number of the ions. The shock's spatial structure is shown to contain a series of 'overshoots' in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. This shock model is applied to an interpretation of the structure of the inner regions of the Crab Nebula, in particular to the 'wisps,' surface brightness enhancements near the pulsar. We argue that these surface brightness enhancements are the regions of magnetic overshoot, which appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure gamma(sub 1) and a number of other properties of the pulsar's wind. We also discuss applications of the shock theory to the termination shocks of the winds from rotation-powered pulsars embedded in compact binaries. We show that this model adequately accounts for (and indeed predicted) the recently discovered X-ray flux from PSR 1957+20, and we discuss several other applications to other examples of these systems.

Arons, Jonathan↗
NASA NTRS2019

Electron Beam Welding of Pure Tungsten Hex Cans for Nuclear Thermal Propulsion Engines

Nuclear thermal propulsion (NTP) is an in-space propulsion method currently being developed at the NASA Marshall Space Flight Center (MSFC). NTP systems are a high specific impulse (750–1,100 s), high thrust (15,000–250,000 lbf ) method of propulsion which have the potential to allow for faster transit times when optimizing for high ΔV. In the nuclear rocket engine, the heat from the nuclear fission reaction is transferred to a low molecular mass propellant (such as hydrogen). Hot propellant is expanded through a nozzle to generate thrust. Development of ceramic metal (cermet) fuel systems for NTP applications is currently ongoing at MSFC. In cermet fuel systems, ceramic fissile fuel particles such as uranium nitride or uranium dioxide are dispersed within a net-shaped, high-density structural matrix. The composite material is cladded by a protective metal structure to make up an NTP fuel element. Cladding materials must be able to withstand the demanding operating conditions required of the engine as well as retain a hermetic seal to allow for retention of fuel element structural integrity, prevent hydrogen attack or migration of the ceramic fuel, and limit release of fission products during operation. For NTP applications, tungsten is a prime material for both the metal matrix and cladding in cermet fuel systems because of its high melting point, high temperature strength, and compatibility with hot hydrogen. If a weld in tungsten with the capability of holding a hermetic seal is achievable, tungsten becomes a strong candidate for NTP applications. This Technical Memorandum focuses on determining the weldability of pure tungsten using electron beam welding (EBW). Tungsten appears well suited for NTP applications, but it has a high ductile to brittle transition temperature (DBTT) dependent upon chemical composition, structure/stress distribution, and mechanical conditions. Therefore, it is highly subject to brittle fracture. Because of its high susceptibility to brittle fracture, it is very difficult to weld. EBW was chosen for joining pure tungsten because of its low heat input compared to gas tungsten arc welding. Reduced heat input can be directly correlated with an increase in ductility of a tungsten weld. EBW is a high energy density welding process in which a stream of electrons penetrates a weld joint in a deep, narrow spike in contrast to a broad gas tungsten arc weld pool. The investigation initially focused on EBW of tungsten plates of both 0.01 in and 0.03 in thickness to determine if EBW could weld pure tungsten without the presence of visual defects—particularly cracking—in the welds. Variation in the weld procedure and post-weld heat treatment (PWHT) was used to improve the surface appearance of flat EBWs on a pure tungsten sheet. The investigation moved on to weld 0.05-in-thick hexagonal tungsten cans with a weld joint thickness of 0.025 in. The goal for welding the pure tungsten hex cans was to avoid any visual surface defects and generate a weld capable of a hermetic seal. This proved difficult. Cold welds commonly exhibited porosity that leaked air. Hot welds exhibited cracks, typically observed immediately after welding. Later welds were preheated to increase ductility and decrease the likelihood of through-thickness cracking. PWHT was used to arrest microcrack growth both in the flat weld samples and hexagonal weld samples.

Courtright, Z. S.↗
NASA NTRSDate not supplied

Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models

The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.

Automation↗
NASA NTRS2001

Solar-Thermal Engine Testing

A solar-thermal engine serves as a high-temperature solar-radiation absorber, heat exchanger, and rocket nozzle. collecting concentrated solar radiation into an absorber cavity and transferring this energy to a propellant as heat. Propellant gas can be heated to temperatures approaching 4,500 F and expanded in a rocket nozzle, creating low thrust with a high specific impulse (I(sub sp)). The Shooting Star Experiment (SSE) solar-thermal engine is made of 100 percent chemical vapor deposited (CVD) rhenium. The engine 'module' consists of an engine assembly, propellant feedline, engine support structure, thermal insulation, and instrumentation. Engine thermal performance tests consist of a series of high-temperature thermal cycles intended to characterize the propulsive performance of the engines and the thermal effectiveness of the engine support structure and insulation system. A silicone-carbide electrical resistance heater, placed inside the inner shell, substitutes for solar radiation and heats the engine. Although the preferred propellant is hydrogen, the propellant used in these tests is gaseous nitrogen. Because rhenium oxidizes at elevated temperatures, the tests are performed in a vacuum chamber. Test data will include transient and steady state temperatures on selected engine surfaces, propellant pressures and flow rates, and engine thrust levels. The engine propellant-feed system is designed to Supply GN2 to the engine at a constant inlet pressure of 60 psia, producing a near-constant thrust of 1.0 lb. Gaseous hydrogen will be used in subsequent tests. The propellant flow rate decreases with increasing propellant temperature, while maintaining constant thrust, increasing engine I(sub sp). In conjunction with analytical models of the heat exchanger, the temperature data will provide insight into the effectiveness of the insulation system, the structural support system, and the overall engine performance. These tests also provide experience on operational aspects of the engine and associated subsystems, and will include independent variation of both steady slate heat-exchanger temperature prior to thrust operation and nitrogen inlet pressure (flow rate) during thrust operation. Although the Shooting Star engines were designed as thermal-storage engines to accommodate mission parameters, they are fully capable of operating as scalable, direct-gain engines. Tests are conducted in both operational modes. Engine thrust and propellant flow rate will be measured and thereby I(sub sp). The objective of these tests is to investigate the effectiveness of the solar engine as a heat exchanger and a rocket. Of particular interest is the effectiveness of the support structure as a thermal insulator, the integrity of both the insulation system and the insulation containment system, the overall temperature distribution throughout the engine module, and the thermal power required to sustain steady state fluid temperatures at various flow rates.

Tucker, Stephen↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
NASA NTRS2011

Catalytic Microtube Rocket Igniter

Devices that generate both high energy and high temperature are required to ignite reliably the propellant mixtures in combustion chambers like those present in rockets and other combustion systems. This catalytic microtube rocket igniter generates these conditions with a small, catalysis-based torch. While traditional spark plug systems can require anywhere from 50 W to multiple kW of power in different applications, this system has demonstrated ignition at less than 25 W. Reactants are fed to the igniter from the same tanks that feed the reactants to the rest of the rocket or combustion system. While this specific igniter was originally designed for liquid methane and liquid oxygen rockets, it can be easily operated with gaseous propellants or modified for hydrogen use in commercial combustion devices. For the present cryogenic propellant rocket case, the main propellant tanks liquid oxygen and liquid methane, respectively are regulated and split into different systems for the individual stages of the rocket and igniter. As the catalyst requires a gas phase for reaction, either the stored boil-off of the tanks can be used directly or one stream each of fuel and oxidizer can go through a heat exchanger/vaporizer that turns the liquid propellants into a gaseous form. For commercial applications, where the reactants are stored as gases, the system is simplified. The resulting gas-phase streams of fuel and oxidizer are then further divided for the individual components of the igniter. One stream each of the fuel and oxidizer is introduced to a mixing bottle/apparatus where they are mixed to a fuel-rich composition with an O/F mass-based mixture ratio of under 1.0. This premixed flow then feeds into the catalytic microtube device. The total flow is on the order of 0.01 g/s. The microtube device is composed of a pair of sub-millimeter diameter platinum tubes connected only at the outlet so that the two outlet flows are parallel to each other. The tubes are each approximately 10 cm long and are heated via direct electric resistive heating. This heating brings the gasses to their minimum required ignition temperature, which is lower than the auto-thermal ignition temperature, and causes the onset of both surface and gas phase ignition producing hot temperatures and a highly reacting flame. The combustion products from the catalytic tubes, which are below the melting point of platinum, are injected into the center of another combustion stage, called the primary augmenter. The reactants for this combustion stage come from the same source but the flows of non-premixed methane and oxygen gas are split off to a secondary mixing apparatus and can be mixed in a near-stoichiometric to highly lean mixture ratio. The primary augmenter is a component that has channels venting this mixed gas to impinge on each other in the center of the augmenter, perpendicular to the flow from the catalyst. The total crosssectional area of these channels is on a similar order as that of the catalyst. The augmenter has internal channels that act as a manifold to distribute equally the gas to the inward-venting channels. This stage creates a stable flame kernel as its flows, which are on the order of 0.01 g/s, are ignited by the combustion products of the catalyst. This stage is designed to produce combustion products in the flame kernel that exceed the autothermal ignition temperature of oxygen and methane.

Schneider, Steven J.↗
NASA NTRS1995

Titan's stratospheric temperature asymmetry: a radiative origin?

During the 1981 Voyager encounter, Titan's stratosphere exhibited a large thermal asymmetry, with high northern latitudes being colder than comparable southern latitudes. Given the short radiative time constant, this asymmetry would not be expected at the season of the Voyager observations (spring equinox), if the infrared and solar opacity sources were distributed symmetrically. We have investigated the radiative budget of Titan's stratosphere, using two selections of Voyager IRIS spectra recorded at symmetric northern and southern latitudes. In the region 0.1-1 mbar, temperatures are 7 K colder at 50 degrees N than at 53 degrees S and the difference reaches approximately 13 K at 5 mbar. On the other hand, the northern region is strongly enriched in nitriles and hydrocarbons, and the haze optical depth derived from the continuum emission between 8 and 15 micrometers is twice as large as in the south. Cooling rate profiles have been computed at the two locations, using the gas and haze abundances derived from the IRIS measurements. We find that, despite lower temperatures, the cooling rate profiles in the pressure range 0.15-5 mbar are 20 to 40% larger in the north than in the south, because of the enhanced concentrations of infrared radiators. Because the northern hemisphere appears darker than the southern one in the Voyager images, enhanced solar heating is also expected to take place at 50 degrees N. Solar heating rate profiles have been calculated, with two different assumptions on the origin of the hemispheric asymmetry. In the most likely case where it results from a variation in the absorbance of the haze material, the heating rates are found to be 12-15% larger at the northern location than at the southern one, a smaller increase than that in the cooling rates. If the lower albedo in the north results from an increase in the particle number density, a 55 to 75% difference is found for the pressure range 0.15-5 mbar, thus larger than that calculated for the cooling rates. Considering the uncertainties in the haze model, dynamical heat transport may significantly contribute to the meridional temperature gradients observed in the stratosphere. On the other hand, the latitudinal variation in gas and haze composition may be sufficient to explain the entire temperature asymmetry observed, without invoking a lag in the thermal response of the atmosphere due to dynamical inertia.

NASA Program Exobiology↗
NASA NTRSDate not supplied

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗
NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗
NASA NTRS1997

Cratering on Titan: A Pre-Cassini Perspective

The NASA-ESA Cassini mission, comprising a formidably instrumented orbiter and parachute-borne probe to be launched this October, promises to reveal a crater population on Titan that has been heretofore hidden by atmospheric haze. This population on the largest remaining unexplored surface in the solar system will be invaluable in comparative planetological studies, since it introduces evidence of the atmospheric effects of cratering on an icy satellite. Here, I highlight some impact features we may hope to find and could devote some modeling effort toward. Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s Atmosphere -94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface.Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s; Atmosphere about 94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface. Titan is comparable to Callisto and Ganymede for strength/gravity, Mars/Earth/Venus for atmospheric interaction, and Hyperion, Rhea, and Iapetus for impactor distribution. The leading/trailing asymmetry of crater density from heliocentric impactors is expected to be about 5-6, in the absence of resurfacing. Any Saturnocentric impactor population is likely to alter this. In particular the impact disruption of Hyperion is noted; because of the 3:4 orbital resonance with Titan, fragments from the proto-Hyperion breakup would have rapidly accreted onto Titan. Titan's resurfacing history is of course unknown. The disruption of impactors into fragments that individually create small craters is expected to occur. A crude estimate suggests a maximum separation of about 2 km (compared with 4 km on Venus, or 0.5 km on Earth). Crater chains are unlikely on Titan, since impactors must pass close enough to Saturn to be tidally disrupted; as a result, they would suffer aerodynamic disruption. Crater counting on adjacent satellites gives densities of about 200 per 10 (exp 6) square km for 20-km-diameter craters. However, the presence of a thick atmosphere leads to atmospheric shielding, depleting the relative abundance of small craters. This has been evaluated by models, and the relative abundance of small craters may be due to a diagnostic atmospheric collapse. A number of radar-dark "splotches" have been detected on Venus; these have been attributed to the interaction of the surface with the atmospheric shockwave produced by the Tunguska-like explosion of a bolide in the atmosphere. Simple analogy suggests that similar features might occur on Titan, but the shocked mass density (which controls the momentum coupling between the surface and the shockwave) of Titan's cold N2 atmosphere is about 20x smaller than that of Venus's hot CO2 atmosphere. Unless ice is much more easily turned to rubble than is rock, such features seem less probable on Titan. When the energy deposited by an impact forms a fireball with an equilibrate greater than one scale height, the fireball expands upward and can distribute ejecta. on ballistic exoatmospheric trajectories. On Venus this process is believed to be responsible for the parabolic features; the interaction of various-sized particles falling through the atmosphere with the zonal wind field winnows the particles to form a parabolic deposit. Although such a process is possible on Titan, the large scale height at higher altitudes would make it more difficult. Comparison with craters on other icy satellites suggests that craters on Titan will be fairly shallow (depth/diameter about 0.1) and craters greater than 10 km in diameter will have central peaks or domed bases, perhaps with central pits. The formation of ejecta. blankets may involve the atmosphere in a significant way, both by restraining the expansion of the ejecta cloud and by influencing the thermal history of the ejecta. Compared with Venus, Titan's atmosphere will chill an impact melt somewhat quickly, so the long ejecta flows seen on Venus seem less likely, detailed modeling needs to be performed to determine the impact melt production. Crater topography on Titan may be highlighted by the influence of liquids forming crater lakes. Craters with central peaks will typically form ring-shaped lakes, although horseshoe-shaped takes may be common; domed craters with central pits may even form bullseye lakes with islands with central ponds. If liquids have covered a substantial part of Titan's surface for a substantial period, hydroblemes and tsunami deposits may be common.

Lorenz, R. D.↗