SEARCH · Search NASA

Results for “electrical distribution grids”

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.

SourceResearch connectionFrom yearThrough year

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

69 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Preliminary designs for 25 kWe advanced Stirling conversion systems for dish electric applications

Under the Department of Energy's (DOE) Solar Thermal Technology Program, Sandia National Laboratories is evaluating heat engines for terrestrial Solar Distributed Heat Receivers. The Stirling engine has been identified by Sandia as one of the most promising engines for terrestrial applications. The Stirling engine also has the potential to meet DOE's performance and cost goals. The NASA Lewis Research Center is conducting Stirling engine technology development activities directed toward a dynamic power source for space applications. Space power systems requirements include high reliability, very long life, low vibration and high efficiency. The free-piston Stirling engine has the potential for future high power space conversion systems, either nuclear or solar powered. Although both applications appear to be quite different, their requirements complement each other. Preliminary designs feature a free-piston Stirling engine, a liquid metal heat transport system, and a means to provide nominally 25 kW electric power to a utility grid while meeting DOE's performance and long term cost goals. The Cummins design incorporates a linear alternator to provide the electrical output, while the STC design generates electrical power indirectly through a hydraulic pump/motor coupled to an induction generator. Both designs for the ASCS's will use technology which can reasonably be expected to be available in the early 1990's

Shaltens, Richard K.↗

NASA NTRS1990

Preliminary designs for 25 kWe advanced Stirling conversion systems for dish electric applications

Shaltens, Richard K.↗

NASA NTRS1979

Electromagnetic heating of minor planets in the early solar system

Electromagnetic processes occurring in the primordial solar system are likely to have significantly affected planetary evolution. In particular, electrical coupling of the kinetic energy of a dense T-Tauri-like solar wind into the interior of the smaller planets could have been a major driver of thermal metamorphism. Accordingly a grid of asteroid models of various sizes and solar distances was constructed using dc transverse magnetic induction theory. Plausible parameterizations with no requirement for a high environmental temperature led to complete melting for Vesta with no melting for Pallas and Ceres. High temperatures were reached in the Pallas model, perhaps implying nonmelting thermal metamorphosis as a cause of its anomalous spectrum. A reversal of this temperature sequence seems implausible, suggesting that the Ceres-Pallas-Vesta dichotomy is a natural outcome of the induction mechanism. Highly localized heating is expected to arise due to an instability in the temperature-controlled current distribution. Localized metamorphosis resulting from this effect may be relevant to the production and evolution of pallasites, the large presumed metal component of S object spectra, and the formation of the lunar magma ocean.

Herbert, F.↗

NASA NTRS2017

Comparison of Aero-Propulsive Performance Predictions for Distributed Propulsion Configurations

NASA's X-57 "Maxwell" flight demonstrator incorporates distributed electric propulsion technologies in a design that will achieve a significant reduction in energy used in cruise flight. A substantial portion of these energy savings come from beneficial aerodynamic-propulsion interaction. Previous research has shown the benefits of particular instantiations of distributed propulsion, such as the use of wingtip-mounted cruise propellers and leading edge high-lift propellers. However, these benefits have not been reduced to a generalized design or analysis approach suitable for large-scale design exploration. This paper discusses the rapid, "design-order" toolchains developed to investigate the large, complex tradespace of candidate geometries for the X-57. Due to the lack of an appropriate, rigorous set of validation data, the results of these tools were compared to three different computational flow solvers for selected wing and propulsion geometries. The comparisons were conducted using a common input geometry, but otherwise different input grids and, when appropriate, different flow assumptions to bound the comparisons. The results of these studies showed that the X-57 distributed propulsion wing should be able to meet the as-designed performance in cruise flight, while also meeting or exceeding targets for high-lift generation in low-speed flight.

Borer, Nicholas K.↗

NASA NTRS1994

A PC based time domain reflectometer for space station cable fault isolation

Significant problems are faced by astronauts on orbit in the Space Station when trying to locate electrical faults in multi-segment avionics and communication cables. These problems necessitate the development of an automated portable device that will detect and locate cable faults using the pulse-echo technique known as Time Domain Reflectometry. A breadboard time domain reflectometer (TDR) circuit board was designed and developed at the NASA-JSC. The TDR board works in conjunction with a GRiD lap-top computer to automate the fault detection and isolation process. A software program was written to automatically display the nature and location of any possible faults. The breadboard system can isolate open circuit and short circuit faults within two feet in a typical space station cable configuration. Follow-on efforts planned for 1994 will produce a compact, portable prototype Space Station TDR capable of automated switching in multi-conductor cables for high fidelity evaluation. This device has many possible commercial applications, including commercial and military aircraft avionics, cable TV, telephone, communication, information and computer network systems. This paper describes the principle of time domain reflectometry and the methodology for on-orbit avionics utility distribution system repair, utilizing the newly developed device called the Space Station Time Domain Reflectometer (SSTDR).

Pham, Michael↗

NASA NTRS1991

Studies of the charging of a thin dust layer in a plasma

The present study is designed to extend Wilson's model to include refinements in the actual operational design of the simulation code itself and to include additional physics in order to assess the resulting electric field structure and grain charges. Most of the initial phases of the work have focussed on the fact that the code, though in principle reliable, suffers from severe time constraints, requiring as must as 10-50 hours of CPU time for typical runs on a VAX computer. It is likely that by staging the code into steps in which early steps are carried out crudely and later steps successively fine tuned the total CPU time can be cut significantly. Preliminary results indicate the possibility to cut CPU time to as much as 1/3 or 1/5 the former time. The original Wilson code has been carefully examined, tested and documented. Shortcuts and staging systems have been introduced and tested which should make the code produce more accurate results with a finer grid system than would have been feasible heretofor. Continued testing and implementation of these refinements will be carried out at the author's home institution in the near future. In addition to the effort at reducing CPU time to a manageable amount, work is in progress to include a dust-size distribution instead of assuming all dust grains have the same size.

Peterson, Lennart R.↗

NASA NTRS1999

A Multipole Expansion Method for Analyzing Lightning Field Changes

Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q model) or a point dipole (the P model). The Q model has four unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P model has six unknown parameters and describes many intracloud (IC) discharges. In this paper we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a dme-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer-simulated sources and on natural lightning at the NASA Kennedy Space Center and U.S. Air Force Eastern Range.

Koshak, William J.↗

NASA NTRS1998

A Multipole Expansion Method for Analyzing Lightning Field Changes

Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q-model) or a point dipole (the P-model). The Q-model has 4 unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P-model has 6 unknown parameters and describes many intracloud (IC) discharges. In this paper, we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a 3-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer simulated sources and on natural lightning at the NASA Kennedy Space Center and USAF Eastern Range.

Koshak, William J.↗

NASA NTRS2013

Analyzing Carbon Dioxide and Methane Emissions in California Using Airborne Measurements and Model Simulations

Greenhouse gas (GHG) concentrations have increased over the past decades and are linked to global temperature increases and climate change. These changes in climate have been suggested to have varying effects, and uncertain consequences, on agriculture, water supply, weather, sea-level rise, the economy, and energy. To counteract the trend of increasing atmospheric concentrations of GHGs, the state of California has passed the California Global Warming Act of 2006 (AB-32). This requires that by the year 2020, GHG (e.g., carbon dioxide (CO2) and methane (CH4)) emissions will be reduced to 1990 levels. To quantify GHG fluxes, emission inventories are routinely compiled for the State of California (e.g., CH4 emissions from the California Greenhouse Gas Emissions Measurement (CALGEM) Project). The major sources of CO2 and CH4 in the state of California are: transportation, electricity production, oil and gas extraction, cement plants, agriculture, landfills/waste, livestock, and wetlands. However, uncertainties remain in these emission inventories because many factors contributing to these processes are poorly quantified. To alleviate these uncertainties, a synergistic approach of applying air-borne measurements and chemical transport modeling (CTM) efforts to provide a method of quantifying local and regional GHG emissions will be performed during this study. Additionally, in order to further understand the temporal and spatial distributions of GHG fluxes in California and the impact these species have on regional climate, CTM simulations of daily variations and seasonality of total column CO2 and CH4 will be analyzed. To assess the magnitude and spatial variation of GHG emissions and to identify local “hot spots”, airborne measurements of CH4 and CO2 were made by the Alpha Jet Atmospheric eXperiment (AJAX) over the San Francisco Bay Area (SFBA) and San Joaquin Valley (SJV) in January and February 2013 during the Discover-AQ-CA study. High mixing ratios of GHGs were observed in-flight with a high degree of spatial variability. To provide an additional method to quantify GHG emissions, and analyze AJAX measurement data, the GEOS-Chem CTM is used to simulate SFBA/SJV GHG measurements. A nested-grid version of GEOS-Chem will be applied and utilizes varying emission inventories and model parameterizations to simulate GHG fluxes/emissions. The model considers CO2 fluxes from fossil fuel use, biomass/biofuel burning, terrestrial and oceanic biosphere exchanges, shipping and aviation, and production from the oxidation of carbon monoxide, CH4, and non-methane volatile organic carbons. The major sources of CH4 simulated in GEOS-Chem are domesticated animals, rice fields, natural gas leakage, natural gas venting/flaring (oil production), coal mining, wetlands, and biomass burning. Preliminary results from the comparison between available observations (e.g., AJAX and CALGEM CH4 emission maps) and GEOS-Chem results will be presented, along with a discussion of CO2 and CH4 source apportionment and the use of the GEOS-Chem-adjoint to perform inverse GHG modeling.

Johnson, Matthew S.↗

NASA NTRS2017

High-Lift Propeller Noise Prediction for a Distributed Electric Propulsion Flight Demonstrator

Over the past several years, the use of electric propulsion technologies within aircraft design has received increased attention. The characteristics of electric propulsion systems open up new areas of the aircraft design space, such as the use of distributed electric propulsion (DEP). In this approach, electric motors are placed in many different locations to achieve increased efficiency through integration of the propulsion system with the airframe. Under a project called Scalable Convergent Electric Propulsion Technology Operations Research (SCEPTOR), NASA is designing a flight demonstrator aircraft that employs many "high-lift propellers" distributed upstream of the wing leading edge and two cruise propellers (one at each wingtip). As the high-lift propellers are operational at low flight speeds (take-off/approach flight conditions), the impact of the DEP configuration on the aircraft noise signature is also an important design consideration. This paper describes efforts toward the development of a mulit-fidelity aerodynamic and acoustic methodology for DEP high-lift propeller aeroacoustic modeling. Specifically, the PAS, OVERFLOW 2, and FUN3D codes are used to predict the aerodynamic performance of a baseline high-lift propeller blade set. Blade surface pressure results from the aerodynamic predictions are then used with PSU-WOPWOP and the F1A module of the NASA second generation Aircraft NOise Prediction Program to predict the isolated high-lift propeller noise source. Comparisons of predictions indicate that general trends related to angle of attack effects at the blade passage frequency are captured well with the various codes. Results for higher harmonics of the blade passage frequency appear consistent for the CFD based methods. Conversely, evidence of the need for a study of the effects of increased azimuthal grid resolution on the PAS based results is indicated and will be pursued in future work. Overall, the results indicate that the computational approach is acceptable for fundamental assessment of low-noise high-lift propeller designs. The extent to which the various approaches may be used in a complementary manner will be further established as measured data becomes available for validation. Ultimately, it is anticipated that this combined approach may be used to provide realistic incident source fields for acoustic shielding/scattering studies on various aircraft configurations.

Nark, Douglas M.↗

NASA NTRS2005

NSTAR Extended Life Test Discharge Chamber Flake Analysis

The Extended Life Test (ELT) of the NASA Solar Electric Propulsion Technology Readiness (NSTAR) ion thruster was concluded after 30,352 hours of operation. The ELT was conducted using the Deep Space 1 (DS1) back-up flight engine, a 30 cm diameter xenon ion thruster. Post-test inspection of the ELT engine revealed numerous contaminant flakes distributed over the bottom of the cylindrical section of the anode within the discharge chamber (DC). Extensive analyses were conducted to determine the source of the particles, which is critical to the understanding of degradation mechanisms of long life ion thruster operation. Analyses included: optical microscopy (OM) and particle length histograms, field emission scanning electron microscopy (FESEM) combined with energy dispersive spectroscopy (EDS), and atomic oxygen plasma exposure tests. Analyses of the particles indicate that the majority of the DC flakes consist of a layered structure, typically with either two or three layers. The flakes comprising two layers were typically found to have a molybdenum-rich (Mo-rich) layer on one side and a carbon-rich (C-rich) layer on the other side. The flakes comprising three layers were found to be sandwich-like structures with Mo-rich exterior layers and a C-rich interior layer. The presence of the C-rich layers indicates that these particles were produced by sputter deposition build-up on a surface external to the discharge chamber from ion sputter erosion of the graphite target in the test chamber. This contaminant layer became thick enough that particles spalled off, and then were electro-statically attracted into the ion thruster interior, where they were coated with Mo from internal sputter erosion of the screen grid and cathode components. Atomic oxygen tests provided evidence that the DC chamber flakes are composed of a significant fraction of carbon. Particle size histograms further indicated that the source of the particles was spalling of carbon flakes from downstream surfaces. Analyses of flakes taken from the downstream surface of the accelerator grid provided additional supportive information. The production of the downstream carbon flakes, and hence the potential problems associated with the flake particles in the ELT ion thruster engine is a facility induced effect and would not occur in the space environment.

deGroh, Kim K.↗

NASA NTRS2004

The Investigation Of Carbon Contamination And Sputtering Effects Of Xenon Ion Thrusters

The Electro-Physics Branch of the NASA Glenn Research Center investigates the effect of atomic oxygen, environmental durability of high performance power materials and surfaces, and low earth orbit. One of its current projects involves the analysis of ion thrusters. Ion thrusters are devices that initiate a beam of ions to a target area. The type of ion thruster that I have been working with this Summer of 2004 emits positively charged Xenon (Xe(+)) atoms through two grids, the screen grid and the accelerator grid, after it enters an ionization chamber. Insulators are used to mechanically hold and separate these two grids. A propellant isolator, an instrument that closely resembles insulators, is placed in front of the ionization chamber. Both the insulator and isolator are made with a ceramic compound and filled with insulating beads. The main difference between the two devices is that the propellant isolator allows gas to flow through, in this case, the gas is Xe(+) and the insulators do not. In order to avoid carbon deposits and other contaminating chemicals to settle on the insulators and propellant isolator, a metal shadow shield is placed around them. These shadow shields function as a protectant and can be shaped in numerous configurations. Part of my job responsibility this summer is to investigate the effectiveness of different shadow shields that are utilized on three different ion engines: the NSTAR (NASA Solar Electric Propulsion Technology Application Readiness), JIMO (Jupiter Icy Moons Orbiter), and NEXIS (Nuclear Electric Xenon Ion System). Using calculus and other mathematical tactics, I was asked to find the total flux of carbon contamination that was able to pass the protectant shadow shield. I familiarized myself with the software program, MathCad2004, to help perform some mathematical computations such as complex integration. Another method of studying the probability of contamination is by experimental simulation. After attaining the precise parameters of the actual shadow shields, I created replicas of three types of shadow shielding to be used to undergo testing. It will be placed in a machine that produces carbon atoms at a high temperature of 200 C. or beam is aimed at a targeted material. As a result of this collision, atoms and other particles are ejected out of the target surface. Another part of my internship consisted of research on sputter ejection, or the angle distribution of sputtered material. This research entailed finding the past results of sputter ejection investigation as well as creating another type of mock simulation. Other minor projects include calculating the path of Xe(+) gas through the insulating beads of the isolators and assisting my mentor in collecting data for his paper for the Joint Propulsion Conference & Exhibit to be held July 11-14,2004 in Fort Lauderdale, Florida.

Prak, Moline K.↗

NASA NTRS2005

Research on Spectroscopy, Opacity, and Atmospheres

I propose to continue providing observers with basic data for interpreting spectra from stars, novas, supernovas, clusters, and galaxies. These data will include allowed and forbidden line lists, both laboratory and computed, for the first five to ten ions of all atoms and for all relevant diatomic molecules. I will eventually expand to all ions of the first thirty elements to treat far UV and X-ray spectra, and for envelope opacities. I also include triatomic molecules provided by other researchers. I have also made CDs with Partridge and Schwenke's water data for work on UV stars. The line data also serve as input to my model atmosphere and synthesis programs that generate energy distributions, photometry, limb darkening, and spectra that can be used for planning observations and for fitting observed spectra. The spectrum synthesis programs produce detailed plots with the lines identified. Grids of stellar spectra can be used for radial velocity-, rotation-, or abundance templates and for population synthesis. I am fitting spectra of bright stars to test the data and to produce atlases to guide observers. For each star the whole spectrum is computed from the UV to the far IR. The line data, opacities, models, spectra, and programs are freely distributed on CDs and on my Web site and represent a unique resource for many NASA programs. I am now in full production of new line lists for atoms. I am computing all ions of all elements from H to Zn and the first 5 ions of all the heavier elements, about 800 ions. For each ion I treat as many as 61 even and 61 odd configurations, computing all energy levels and eigenvectors. The Hamiltonian is determined from a scaled-Hartree-Fock starting guess by least squares fitting the observed energy levels. The average energy of each configuration is used in computing scaled-Thomas-Fermi-Dirac wavefunctions for each configuration which in turn are used to compute allowed and forbidden transition integrals. These are multiplied into the LS allowed and forbidden transition arrays. The transition arrays are transformed to the observed coupling to yield the allowed and forbidden line lists. Results are put on the web as they are finished. Provided I get funding,there will be more than 500 million lines. I will then compare ion by ion, to all the laboratory and computed data in the literature and make up a working line list for spectrum synthesis and opacity calculations with the best available data. As the laboratory spectrum analyses are improved, I will redo the calculations with the new energy levels. My original plan when I started the new calculations was to run through all the atoms using my old Cray programs from the 1980's that were limited to 1100 x 1100 arrays in the Hamiltonian for each J. Then I would go back and rerun the more complicated cases with 3000 x 3000 arrays so that I could include many more configurations and more configuration interactions. At present I am limited to 61 even and 61 odd configurations and I try to include everything up through n = 9. The current program runs on Alpha workstations. I decided to test the big program on Fe I and Fe II to see whether there was any great difference in the low configurations compared to those from the Cray program. Besides increasing the number of E1 lines by a factor of 6 to 7.7 million, there was an unexpected result: the electric quadrupole transitions were 10 times stronger than before because the transition integrals are weighted by r(exp 2) ---they become very large for high n, and because there are numerous configuration interactions that mix the low and high configurations. As a check I was able to reproduce Carstang's (1962) lower results by running his three configurations with my program. Since my model atom is still only a subset of a real Fe II ion, the true quadrupole A values are probably larger than mine. The magnetic dipole lines are affected by the mixing but the overall scale does not change. Because of this scovery I decided that there was no point in computing the small array cases. I have been running with as many configurations as I can and with thousands of parameters in the Hamiltonian. The computer runs take much longer to set up and produce than I had expected. I have concentrated on redoing the low iron group spectra, especially to get data for supernova modelers. I have done only Ca I -- Zn I, Ca II -- Zn II, CU I -- Cu XXIX, Zn I - Zn XXX, for practice at high stages of ionization, C I, C II, S I, and CL I and Ag I for people who were working on the laboratory spectra. Check my web site kurucz.harvard.edu for current additions. My latest calculations have been for carbon I and sulphur I, and silicon I is under way using the same elaborate approach as for C I, which took many months to do. These line lists greatly increase the number of lines in the ultraviolet, in the visible, and especially in the infrared. They will increase the opacity in A, F, and G stars. They will account for many unidentified lines in the sun.

Kurucz, Robert L.↗

NASA NTRS2017

Design of the Cruise and Flap Airfoil for the X-57 Maxwell Distributed Electric Propulsion Aircraft

A computational and design study on an airfoil and high-lift flap for the X-57 Maxwell Distributed Electric Propulsion (DEP) testbed aircraft was conducted. The aircraft wing sizing study resulted in a wing area of 66.67 sq ft and aspect ratio of 15 with a design requirement of V(stall) = 58 KEAS, at a gross weight of 3,000 lb. To meet this goal an aircraft C(L,max) of 4.0 was required. The design cruise condition is 150 KTAS at 8,000 ft. This resulted in airfoil requirements of c(l) is approximately 0.90 for the cruise condition at Re = 2.35 x 10 (exp 6). A flapped airfoil with a c(l,max) of approximately 2.5 or greater, at Re = 1.0 x 10 (exp 6), was needed to have enough lift to meet the stall requirement with the DEP system. MSES computational analyses were conducted on the GAW-1, GAW-2, and the NACA 5415 airfoil sections, however they had limitations in either high drag or low c(l,max) on the cruise airfoil, which was the impetus for a new design. A design was conducted to develop a low drag airfoil for the X-57 cruise conditions with high c(l,max). The final design was the GNEW5BP93B airfoil with a minimum drag coefficient of c(d) = 0.0053 at c(l) = 0.90 and achieved laminar flow back to 69% chord on the upper surface and 62% chord on the lower surface. With fully turbulent flow, the drag increases to c(d) = 0.0120. The predicted maximum lift with turbulent flow is a c(l,max) of 1.95 at alpha = 19 deg. The airfoil is characterized by relatively flat pressure gradient regions on both surfaces at alpha = 0 deg, and aft camber to get extra lift out of the lower surface concave region. A 25% chord slotted flap was designed and analyzed with MSES for a 30 deg flap deflection. Additional 30 deg and 40 deg flap deflection analyses for two flap positions were conducted with USM3D using several turbulence models, for two angles of attack, to assess near c(l,max) with varied flap position. The maximum c(l) varied between 2.41 and 3.35. An infinite-span powered high-lift study was conducted on a GAW-1 constant chord 40 deg flapped airfoil section with FUN3D to quantify the airfoil lift increment that can be expected from a DEP system. The 16.7 hp/propeller blown wing increases the maximum C(L) from 3.45 to C(L) = 6.43, which is an effective q ratio of 1.86. This indicates that if the unblown high-lift flapped airfoil of the X-57 airplane achieves a c(l,max) of 2.78, then the high-lift augmentation blowing could yield a sectional lift coefficient of approximately 4.95 at c(l,max). Finally, a computational study was conducted with FUN3D on an infinite-span constant chord GAW-1 cruise airfoil to determine the impact of high-lift propeller diameter to wing chord ratio on the lift increment of the DEP system. A constant diameter propeller and nacelle size were used in the study. Three computational grids were made with airfoil chords of 0.5*chord, 1.0*chord, and 2.0*chord. Results of the propeller diameter to wing chord ratio study indicated that the blown to unblown C(L) ratio increased as the chord was decreased. However, because of the increase in relative size of the high-lift nacelle to the wing, which impacted wing lift performance, the study indicated that a propeller diameter to wing chord ratio of 1.0 gives the overall best maximum lift on the wing with the DEP system.

Viken, Jeffrey K.↗

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↗