SEARCH · Engineering Papers

Results for “Optimal Power Flow”

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.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Development of a High-Order Space-Time Matrix-Free Adjoint Solver

The growth in computational power and algorithm development in the past few decades has granted the science and engineering community the ability to simulate flows over complex geometries, thus making Computational Fluid Dynamics (CFD) tools indispensable in analysis and design. Currently, one of the pacing items limiting the utility of CFD for general problems is the prediction of unsteady turbulent ows.1{3 Reynolds-averaged Navier-Stokes (RANS) methods, which predict a time-invariant mean flowfield, struggle to provide consistent predictions when encountering even mild separation, such as the side-of-body separation at a wing-body junction. NASA's Transformative Tools and Technologies project is developing both numerical methods and physical modeling approaches to improve the prediction of separated flows. A major focus of this e ort is efficient methods for resolving the unsteady fluctuations occurring in these flows to provide valuable engineering data of the time-accurate flow field for buffet analysis, vortex shedding, etc. This approach encompasses unsteady RANS (URANS), large-eddy simulations (LES), and hybrid LES-RANS approaches such as Detached Eddy Simulations (DES). These unsteady approaches are inherently more expensive than traditional engineering RANS approaches, hence every e ort to mitigate this cost must be leveraged. Arguably, the most cost-effective approach to improve the efficiency of unsteady methods is the optimal placement of the spatial and temporal degrees of freedom (DOF) using solution-adaptive methods.

Adjoint↗

NASA NTRS2011

Feasibility of Ultraviolet Light Emitting Diodes as an Alternative Light Source for Photocatalysis

The objective of this study was to determine whether ultraviolet light emitting diodes (UV-LEDs) could serve as an alternative photon source efficiently for heterogeneous photocatalytic oxidation (PCO). An LED module consisting of 12 high-power UV-A LEDs was designed to be interchangeable with a UV-A fluorescent black light blue (BLB) lamp in a Silica-Titania Composite (STC) packed bed annular reactor. Lighting and thermal properties were characterized to assess the uniformity and total irradiant output. A forward current of (I(sub F)) 100 mA delivered an average irradiance of 4.0 m W cm(exp -2), which is equivalent to the maximum output of the BLB, but the irradiance of the LED module was less uniform than that of the BLB. The LED- and BLB-reactors were tested for the oxidization of 50 ppmv ethanol in a continuous flow-through mode with 0.94 sec space time. At the same irradiance, the UV-A LED reactor resulted in a lower PCO rate constant than the UV-A BLB reactor (19.8 vs. 28.6 nM CO2 sec-I), and consequently lower ethanol removal (80% vs. 91%) and mineralization efficiency (28% vs. 44%). Ethanol mineralization increased in direct proportion to the irradiance at the catalyst surface. This result suggests that reduced ethanol mineralization in the LED- reactor could be traced to uneven irradiance over the photocatalyst, leaving a portion of the catalyst was under-irradiated. The potential of UV-A LEDs may be fully realized by optimizing the light distribution over the catalyst and utilizing their instantaneous "on" and "off' feature for periodic irradiation. Nevertheless, the current UV-A LED module had the same wall plug efficiency (WPE) of 13% as that of the UV-A BLB. These results demonstrated that UV-A LEDs are a viable photon source both in terms of WPE and PCO efficiency.

Levine, Langanf H.↗

NASA NTRS2016

Computational Study of Flow Interactions in Coaxial Rotors

Although the first idea of coaxial rotors appeared more than 150 years ago, most helicopters have used single main-rotor/tail-rotor combination. Since reactive moments of coaxial rotors are canceled by contra-rotation, no tail rotor is required to counter the torque generated by the main rotor. Unlike the single main rotor design that distributes power to both main and tail rotors, all of the power for coaxial rotors is used for vertical thrust. Thus, no power is wasted for anti-torque or directional control. The saved power helps coaxial rotors reach a higher hover ceiling than single rotor helicopters. Another advantage of coaxial rotors is that the overall rotor diameter can be reduced for a given vehicle gross weight because each rotor provides a maximum contribution to vertical thrust to overcome vehicle weight. However, increased mechanical complexity of the hub has been one of the challenges for manufacturing coaxial rotorcraft. Only the Kamov Design Bureau of Russia had been notably successful in production of coaxial helicopters until Sikorsky built X2, an experimental compound helicopter. Recent developments in unmanned aircraft systems and high-speed rotorcraft have renewed interest in the coaxial configuration. Multi-rotors are frequently used for small electric unmanned rotorcraft partly due to mechanical simplicity. The use of multiple motors provides redundancy as well as cost-efficiency. The multi-rotor concept has rarely been used until recently because of its inherent stability and control problems. However, advances in inexpensive electronic flight control systems have opened the floodgates for small drones using multirotors. Coaxial rotors have started to appear in some multi-rotor configurations. Small coaxial rotors have often been designed using a hundred year old approach that is "sketch, build, fly, and iterate." In that approach, there is no systematic way to explore trade-offs or determine logical next steps. It is neither possible to account for multiple real-world constraints up front in design nor possible to know what performance is possible with a given design. Since unmanned vehicles are sized and optimized for the particular mission, a modern low-fidelity conceptual design and sizing tool that has been used for the design of large helicopters can be used for design of small coaxial rotorcraft. However, unlike most helicopters with single main rotor, the interactions between the upper and lower rotors emerge as an important factor to consider in design because an increase in performance of a multi-rotor system is not proportional to the number of rotors. Interference losses and differences in thrusts between the upper and lower rotors were investigated by theoretical methods as well as a computational fluid dynamics (CFD) method using the Reynolds-Averaged Navier-Stokes (RANS) equations. In this work, hybrid turbulence models are used to investigate the physics of interactions between coaxial rotors and a fuselage that are not well understood. Present study covers not only small-scale drones but also large-scale coaxial rotors for heavy-lifting missions. Considering the recently proposed FAA drone rules that require the flight only in visual line-of-sight, a large multirotor might be used as an airborne carrier for launch and recovery of unmanned aircraft systems with a human operator onboard. For applications to civil operations, their aerodynamic performance and noise levels need to be assessed. Noise is one of the largest limiting factors to rotorcraft operations in urban area. Since the high-frequency noise of multi-rotors may increase the annoyance, noise may turn out to be a key issue that must be addressed for market acceptability. One of the objectives of the present work is to study the effects of inter-rotor spacing and collectives on the performance, efficiency, and acoustics of coaxial rotor systems.

Study↗

NASA NTRS2010

Optimal and Local Connectivity Between Neuron and Synapse Array in the Quantum Dot/Silicon Brain

This innovation is used to connect between synapse and neuron arrays using nanowire in quantum dot and metal in CMOS (complementary metal oxide semiconductor) technology to enable the density of a brain-like connection in hardware. The hardware implementation combines three technologies: 1. Quantum dot and nanowire-based compact synaptic cell (50x50 sq nm) with inherently low parasitic capacitance (hence, low dynamic power approx.l0(exp -11) watts/synapse), 2. Neuron and learning circuits implemented in 50-nm CMOS technology, to be integrated with quantum dot and nanowire synapse, and 3. 3D stacking approach to achieve the overall numbers of high density O(10(exp 12)) synapses and O(10(exp 8)) neurons in the overall system. In a 1-sq cm of quantum dot layer sitting on a 50-nm CMOS layer, innovators were able to pack a 10(exp 6)-neuron and 10(exp 10)-synapse array; however, the constraint for the connection scheme is that each neuron will receive a non-identical 10(exp 4)-synapse set, including itself, via its efficacy of the connection. This is not a fully connected system where the 100x100 synapse array only has a 100-input data bus and 100-output data bus. Due to the data bus sharing, it poses a great challenge to have a complete connected system, and its constraint within the quantum dot and silicon wafer layer. For an effective connection scheme, there are three conditions to be met: 1. Local connection. 2. The nanowire should be connected locally, not globally from which it helps to maximize the data flow by sharing the same wire space location. 3. Each synapse can have an alternate summation line if needed (this option is doable based on the simple mask creation). The 10(exp 3)x10(exp 3)-neuron array was partitioned into a 10-block, 10(exp 2)x10(exp 3)-neuron array. This building block can be completely mapped within itself (10,000 synapses to a neuron).

Duong, Tuan A.↗

NASA NTRS2007

NASA Tech Briefs, January 2007

Topics covered include: Flexible Skins Containing Integrated Sensors and Circuitry; Artificial Hair Cells for Sensing Flows; Video Guidance Sensor and Time-of-Flight Rangefinder; Optical Beam-Shear Sensors; Multiple-Agent Air/Ground Autonomous Exploration Systems; A 640 512-Pixel Portable Long-Wavelength Infrared Camera; An Array of Optical Receivers for Deep-Space Communications; Microstrip Antenna Arrays on Multilayer LCP Substrates; Applications for Subvocal Speech; Multiloop Rapid-Rise/Rapid Fall High-Voltage Power Supply; The PICWidget; Fusing Symbolic and Numerical Diagnostic Computations; Probabilistic Reasoning for Robustness in Automated Planning; Short-Term Forecasting of Radiation Belt and Ring Current; JMS Proxy and C/C++ Client SDK; XML Flight/Ground Data Dictionary Management; Cross-Compiler for Modeling Space-Flight Systems; Composite Elastic Skins for Shape-Changing Structures; Glass/Ceramic Composites for Sealing Solid Oxide Fuel Cells; Aligning Optical Fibers by Means of Actuated MEMS Wedges; Manufacturing Large Membrane Mirrors at Low Cost; Double-Vacuum-Bag Process for Making Resin- Matrix Composites; Surface Bacterial-Spore Assay Using Tb3+/DPA Luminescence; Simplified Microarray Technique for Identifying mRNA in Rare Samples; High-Resolution, Wide-Field-of-View Scanning Telescope; Multispectral Imager With Improved Filter Wheel and Optics; Integral Radiator and Storage Tank; Compensation for Phase Anisotropy of a Metal Reflector; Optical Characterization of Molecular Contaminant Films; Integrated Hardware and Software for No-Loss Computing; Decision-Tree Formulation With Order-1 Lateral Execution; GIS Methodology for Planning Planetary-Rover Operations; Optimal Calibration of the Spitzer Space Telescope; Automated Detection of Events of Scientific Interest; Representation-Independent Iteration of Sparse Data Arrays; Mission Operations of the Mars Exploration Rovers; and More About Software for No-Loss Computing.

Source record↗

NASA NTRSDate not supplied

Comprehensive 3D Multiphysics Model on Electrochemical Recovery of O 2 from Metabolic CO 2 at the International Space Station (ISS)

The International Space Station (ISS) is presently equipped with an elaborate, heavy, and high-power consuming system that recovers approximately 50% of O 2 from metabolic CO 2 as part of the atmospheric revitalization (AR) at the ISS habitat. Future long-duration missions will require a sustainable and efficient system capable of yielding a minimum of 75% O 2 recovery to reach the self-sufficiency required for long space missions beyond earth’s low orbit. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to not only increase significantly current O 2 recovery efficiency, improving self-sufficiency on AR at the ISS habitat and future long missions, but also reduce the complexity of the system. The authors have developed and deployed a comprehensive 3D multiphysics model that thoroughly replicates the actual configuration and fluid/material domains of the MFECR. The coupled physics in this multiphysics model include multicomponent-multiphase electrochemical-driven reactions, non-ideal mass transport mechanism, free and porous flow, heat transfer, CO 2 solubility on alkaline electrolyte, water condensation on porous medium, and DC electrical current generation along with Joule heating effect. This model is aimed to conduct quantitive benchmark on three different MFECR’s layouts, one without serpentine paths (plain) and two with serpentines leading to four and twelve paths respectively. Once experimental data is generated via a test matrix of 200 tests, the model will be validated to conduct MFECR’s process optimization and revalidate the quantitive benchmark on three different MFECR’s layouts.

Jesus A. Dominguez↗

NASA NTRS2008

Investigation of Condensing Ice Heat Exchangers for MTSA Technology Development

Metabolic heat regenerated Temperature Swing Adsorption (MTSA) technology is being developed for thermal, carbon dioxide (CO2) and humidity control for a Portable Life Support Subsystem (PLSS). Metabolically-produced CO2 present in the ventilation gas of a PLSS is collected using a CO2selective adsorbent via temperature swing adsorption. The temperature swing is initiated through cooling to well below metabolic temperatures. Cooling is achieved with a sublimation heat exchanger using water or liquid carbon dioxide (LCO2) expanded below sublimation temperature when exposed to low pressure or vacuum. Subsequent super heated vapor, as well as additional coolant, is used to further cool the astronaut. The temperature swing on the adsorbent is then completed by warming the adsorbent with a separate condensing ice heat exchanger (CIHX) using metabolic heat from moist ventilation gas. The condensed humidity in the ventilation gas is recycled at the habitat. The water condensation from the ventilation gas is a significant heat transfer mechanism for the warming of the adsorbent bed because it represents as much as half of the energy potential in the moist ventilation gas. Designing a heat exchanger to efficiently transfer this energy to the adsorbent bed and allow the collection of the water is a challenge since the CIHX will operate in a temperature range from 210K to 280K. The ventilation gas moisture will first freeze and then thaw, sometimes existing in three phases simultaneously. A NASA Small Business Innovative Research (SBIR) Phase 1 contract was performed to investigate condensing and icing as applied to MTSA to enable higher fidelity modeling and assess the impact of geometry variables on CIHX performance for future CIHX design optimization. Specifically, a design tool was created using analytical relations to explore the complex, interdependent design space of a condensing ice heat exchanger. Numerous variables were identified as having nontrivial contributions to performance such as hydraulic diameter, heat exchanger effectiveness, ventilation gas mass flow rate and surface roughness. Using this tool, four test articles were designed and manufactured to map to a full MTSA subassembly (the adsorbent bed, the sublimation heat exchanger for cooling and the condensing ice heat exchanger for warming). The design mapping considered impacts due to CIHX geometry as well as subassembly impacts such as thermal mass and thermal resistance through the adsorbent bed. The test articles were tested at simulated PLSS ventilation loop temperature, moisture content and subambient pressure. Ice accumulation and melting were observed. Data and test observations were analyzed to identify drivers of the condensing ice heat exchanger performance. This paper will discuss the analytical models, the test article designs, and testing procedures. Testing issues will be discussed to better describe data and share lessons learned. Data analysis and subsequent conclusions will be presented.

Padilla, Sebastian↗

NASA NTRS1999

The Effect of Penetration Depth on Thermal Contrast of NDT by Thermography

Nondestructive evaluation by Thermography (TNDE) is generally classified into two categories, the passive approach and the active approach. The passive approach is usually performed by measuring the natural temperature difference between the ambient and the material or structure to be tested. The active approach, on the other hand, requires the application of an external energy source to stimulate the material for inspection. A laser, a heater, a hot air blower, a high power thermal pulse, mechanical, or electromagnetic energy may provide the energy sources. For the external heating method to inspect materials for defects and imperfection at ambient temperature, a very short burst of heat can be introduced to one of the surfaces or slow heating of the side opposite to the side being observed. Due to the interruption of the heat flow through the defects, the thermal images will reveal the defective area by contrasting against the surrounding good materials. This technique is called transient Thermography, pulse video Thermography, or thermal wave imaging. As an empirical rule, the radius of the smallest defect should be at least one to two times larger than its depth under the surface. Thermography is being used to inspect void, debond, impact damage, and porosity in composite materials. It has been shown that most of the defects and imperfection can be detected. However, the current method of inspection using thermographic technique is more of an art than a practical scientific and engineering approach. The success rate of determining the defect location and defect type is largely depend on the experience of the person who operates thermography system and performs the inspection. The operator has to try different type of heat source, different duration of its application time, as well as experimenting with the thermal image acquisition time and interval during the inspection process. Further-more, the complexity of the lay-up and structure of composites makes it more difficult to determine the optimal operating condition for revealing the defects. In order to develop an optimal thermography inspection procedure, we must understand the thermal behavior inside the material subjected to transient heat in order to interpret the thermal images correctly. Fabrication of finite element models of characteristic defects in composite materials subjected to transient heat will enable the development of appropriate procedure for thermography inspection. Design of phantom defects could be modeled and behavior characterized prior to physically building these test parts. Since production of phantom test parts can be very time consuming and laborious, it is important to design good representative defects.

Chu, Tsuchin Philip↗

NASA NTRS2020

PACE Propulsion Subsystem Surge Analysis and Testing

Propulsion subsystem surge analysis and testing was performed at NASA Goddard Space Flight Center (GSFC) in order to analytically and experimentally determine the potential for damaging dynamic pressure events in the Plankton Aerosol Cloud and ocean Ecosystem (PACE) propulsion subsystem. Surge pressures occur when propellant flow is initiated following the actuation of an isolation valve separating an upstream high pressure fluid from a downstream low pressure manifold. By design, the PACE propulsion subsystem mitigates surge pressures using an orifice, or a cavitating venturi, placed upstream of each thruster bank isolation latch valve. Venturi sizing is a balancing act due to the competing performance goals of staying within the margins of an acceptable pressure loss and maintaining surge pressures below the proof pressure rating of the subsystem components. Subsystem pressure loss is dependent upon the venturi throat diameter and is higher for smaller venturi sizes. A venturi with a smaller throat diameter incurs an increased pressure loss through the orifice resulting in a reduced thruster inlet pressure and diminished thruster performance. The PACE propulsion subsystem surge analysis and testing supports system optimization, i.e. pressure drop minimization, while still meeting the surge pressure requirements to maximize overall system performance. In this paper a new approach to surge pressure analytical gauging is discussed, as well as its accuracy when compared to traditional surge testing campaign.

Space↗

NASA NTRS1959

Analysis and Evaluation of Supersonic Underwing Heat Addition

The linearized theory for heat addition under a wing has been developed to optimize wing geometry, heat addition, and angle of attack. The optimum wing has all of the thickness on the underside of the airfoil, with maximum-thickness point well downstream, has a moderate thickness ratio, and operates at an optimum angle of attack. The heat addition is confined between the fore Mach waves from under the trailing surface of the wing. By linearized theory, a wing at optimum angle of attack may have a range efficiency about twice that of a wing at zero angle of attack. More rigorous calculations using the method of characteristics for particular flow models were made for heating under a flat-plate wing and for several wings with thickness, both with heat additions concentrated near the wing. The more rigorous calculations yield in practical cases efficiencies about half those estimated by linear theory. An analysis indicates that distributing the heat addition between the fore waves from the undertrailing portion of the wing is a way of improving the performance, and further calculations appear desirable. A comparison of the conventional ramjet-plus wing with underwing heat addition when the heat addition is concentrated near the wing shows the ramjet to be superior on a range basis up to Mach number of about B. The heat distribution under the wing and the assumed ramjet and airframe performance may have a marked effect on this conclusion. Underwing heat addition can be useful in providing high-altitude maneuver capability at high flight Mach numbers for an airplane powered by conventional ramjets during cruise.

wings↗

NASA NTRS2011

NASA Tech Briefs, May 2011

Topics covered include: 1) Method to Estimate the Dissolved Air Content in Hydraulic Fluid; 2) Method for Measuring Collimator-Pointing Sensitivity to Temperature Changes; 3) High-Temperature Thermometer Using Cr-Doped GdAlO3 Broadband Luminescence; 4)Metrology Arrangement for Measuring the Positions of Mirrors of a Submillimeter Telescope; 5) On-Wafer S-Parameter Measurements in the 325-508-GHz Band; 6) Reconfigurable Microwave Phase Delay Element for Frequency Reference and Phase-Shifter Applications; 7) High-Speed Isolation Board for Flight Hardware Testing; 8) High-Throughput, Adaptive FFT Architecture for FPGA-Based Spaceborne Data Processors; 9) 3D Orbit Visualization for Earth-Observing Missions; 10) MaROS: Web Visualization of Mars Orbiting and Landed Assets; 11) RAPID: Collaborative Commanding and Monitoring of Lunar Assets; 12) Image Segmentation, Registration, Compression, and Matching; 13) Image Calibration; 14) Rapid ISS Power Availability Simulator; 15) A Method of Strengthening Composite/Metal Joints; 16) Pre-Finishing of SiC for Optical Applications; 17) Optimization of Indium Bump Morphology for Improved Flip Chip Devices; 18) Measuring Moisture Levels in Graphite Epoxy Composite Sandwich Structures; 19) Marshall Convergent Spray Formulation Improvement for High Temperatures; 20) Real-Time Deposition Monitor for Ultrathin Conductive Films; 21) Optimized Li-Ion Electrolytes Containing Triphenyl Phosphate as a Flame-Retardant Additive; 22) Radiation-Resistant Hybrid Lotus Effect for Achieving Photoelectrocatalytic Self-Cleaning Anticontamination Coatings; 23) Improved, Low-Stress Economical Submerged Pipeline; 24) Optical Fiber Array Assemblies for Space Flight on the Lunar Reconnaissance Orbiter; 25) Local Leak Detection and Health Monitoring of Pressurized Tanks; 26) Dielectric Covered Planar Antennas at Submillimeter Wavelengths for Terahertz Imaging; 27) Automated Cryocooler Monitor and Control System; 28) Broadband Achromatic Phase Shifter for a Nulling Interferometer; 29) Super Dwarf Wheat for Growth in Confined Spaces; 30) Fine Guidance Sensing for Coronagraphic Observatories; 31) Single-Antenna Temperature- and Humidity-Sounding Microwave Receiver; 32) Multi-Wavelength, Multi-Beam, and Polarization-Sensitive Laser Transmitter for Surface Mapping; 33) Optical Communications Link to Airborne Transceiver; 34) Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings; 35) Entanglement in Self-Supervised Dynamics; 36) Prioritized LT Codes; 37) Fast Image Texture Classification Using Decision Trees; 38) Constraint Embedding Technique for Multibody System Dynamics; 39) Improved Systematic Pointing Error Model for the DSN Antennas; 40) Observability and Estimation of Distributed Space Systems via Local Information-Exchange Networks; 41) More-Accurate Model of Flows in Rocket Injectors; 42) In-Orbit Instrument-Pointing Calibration Using the Moon as a Target; 43) Reliability of Ceramic Column Grid Array Interconnect Packages Under Extreme Temperatures; 44) Six Degrees-of-Freedom Ascent Control for Small-Body Touch and Go; and 45) Optical-Path-Difference Linear Mechanism for the Panchromatic Fourier Transform Spectrometer.

Source record↗

NASA NTRSDate not supplied

Towards Sheaf Theoretic Analyses for Delay Tolerant Networking

The goal of Delay Tolerant Networking (DTN) is to take a collection of heterogeneous, disparate connections between satellites, space assets, ground stations, and ground infrastructure and bring it together into a cohesive, functioning overlay network. Depending on the systems being considered, one can find links with a one-way light time exceeding minutes (and hours),periodic links which can sometimes be predicted by orbital mechanics, and restrictions based on the variety of capabilities built into these systems. These characteristics preclude traditional network models and routing techniques and have classically led to either rigid routing tables or purely probabilistic models. As the deeper underlying structures remain unknown, development of more DTN-optimized algorithms has lacked the necessary foundation. In a continuation of previous work, the goal of this paper is to identify and study these fundamental structures that exist in delay tolerant networks (DTN), with a focus on space networks. The current routing methodology has been to use contact graph routing (CGR) algorithms. CGR models a series of known contacts as a static graph. For CGR to work, this graph must be globally consistent and must have an accurate picture of the network. Because this is a globally controlled structure, there is little room for flexibility in the event of changes to the network which would naturally occur as the network grows. As a response to the desire for flexibility as the network changes, we introduced the mathematical structure known as sheaves to DTNs last year. The tag-line for sheaves is that they are a mathematically precise way of gluing local data together into unique global data. Thus, sheaves lend extra power to traditional models(and routing algorithms) by taking additional information and merging it, in as consistent a manner as possible, with the representation itself. The clearest example of how Earth-bound networks exhibit behavior that is “sheafy” is link state routers, which build a local-to-global picture of their network by gluing local information together into a global network, exactly as a sheaf would do. For routing within delay tolerant networks to truly exploit this structure, a deeper structure than a graph is required. In this paper, we develop sheaves that can work over directed graphs such as temporal flow networks, we construct a sheaf representation for Dijkstra’s algorithm, and we outline a construction for routing sheaves capable of modeling multicast scenarios. Finally, there is a section of future work suggesting follow-on research.

Robert Short↗

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 NTRS2015

Development of Long-Lifetime Pulsed Gas Valves for Pulsed Electric Thrusters

It is advantageous for gas-fed pulsed electric thrusters to employ pulsed valves so propellant is only flowing to the device during operation. The propellant utilization of the thruster will be maximized when all the gas injected into the thruster is acted upon by the fields produced by the electrical pulse. Gas that is injected too early will diffuse away from the thruster before the electrical pulse can act to accelerate the propellant. Gas that is injected too late will miss being accelerated by the already-completed electrical pulse. As a consequence, the valve must open quickly and close equally quickly, only remaining open for a short duration. In addition, the valve must have only a small amount of volume between the sealing body and the thruster so the front and back ends of the pulse are as coincident as possible with the valve cycling, with very little latent propellant remaining in the feed lines after the valve is closed. For a real mission of interest, a pulsed thruster can be expected to pulse at least 10(exp 10) - 10(exp 11) times, setting the range for the number of times a valve must open and close. The valves described in this paper have been fabricated and tested for operation in an inductive pulsed plasma thruster (IPPT) for in-space propulsion. In general, an IPPT is an electrodeless space propulsion device where a capacitor is charged to an initial voltage and then discharged, producing a high-current pulse through a coil. The field produced by this pulse ionizes propellant, inductively driving current in a plasma located near the face of the coil. Once the plasma is formed, it can be accelerated and expelled at a high exhaust velocity by the electromagnetic Lorentz body force arising from the interaction of the induced plasma current and the magnetic field produced by the current in the coil. The valve characteristics needed for the IPPT application require a fast-acting valve capable of a minimum of 10(exp 10) valve actuation cycles. Since even 10(exp 9) cycles is well above anything demonstrated, this lower value was selected as the design point for the present work. The valve seal must remain leak-tight throughout operation, and the body must maintain a low internal leakage at relatively high operating temperatures. The full set of design requirements used for this program are summarized in Table 1. In this work, we describe two pulsed gas valves that have been fabricated to have long lifetime and demonstrate the characteristics listed above. The first is a miniaturized, conventional electromagnet-based valve while the second is a piezoelectric-based valve design. The conventional valve, shown in Fig. 1, is opened by use of a solenoid electromagnetic actuator. When current is applied to the solenoid coil, magnetic forces pull the plunger away from the valve seat, allowing fluid to flow through the valve. Removal of electrical current permits the spring and fluid pressure to seat the plunger, halting the flow of fluid. The valve body is fabricated from 304L corrosion resistant steel (CRES) and while the parts that form the magnetic circuit are fabricated from 430 CRES. This material does not have optimum magnetic properties, but its corrosion resistance permits incorporation into a design without requiring an additional plating process. A viton O-ring compound (Parker V0884-75), selected for its mechanical strength at elevated temperatures, was used for the valve seat seal. The design was based solely on the use of analytical sizing calculations, as opposed to a more rigorous finite element analysis. While this valve is small and relatively lightweight, it does not represent a design that is optimized for mass and/or a given volume envelope. The piezoelectric valve is a "puller" valve design. Applying a voltage to the piezo crystal causes it to elongate and pull a pintle off the seat, opening the valve. The valve seal consists of the pintle with an external, spherically-formed tip fabricated from Inconel 625 and the seat fabricated from Inconel 718 and possessing a spherical surface that permits the pintle to align with the valve body centerline after each cycle, maintaining the circumferential metal-on-metal seal. The contact stress on the valve pintle is expected in the worst case to be 96.5 MPa (14 ksi), which is well below the work hardening value of 413.7 MPa (60 ksi) for Inconel 615. Results presented in this work will compare the performance of each valve, both to each other and to the valve design requirements given in Table 1. In general, both valves meet all the performance envelope requirements, but the piezoelectric valve consumes considerably less power and has a faster response time over the conventional solenoid-driven valves.

Burkhardt, Wendel M.↗

NASA NTRS1989

The transient performance of a two-phase fluid reservoir

Thermal control of future large, high power spacecraft will require a two-phase fluid central bus. The two-phase fluid reservoir is a critical component in the two-phase fluid bus. It both controls the saturation temperature and provides a space for volumetric changes. A dynamic reservoir simulation model does not currently exist, but it is needed to expedite efforts and reduce risk. During 1989 an effort was made to develop a simulation model of the transient performance of a two-phase fluid reservoir. As a beginning, a preliminary model was developed. It is based upon component mathematical models in lumped parametric form and build upon five component mathematical models for calculating dynamic responses of two-phase fluid reservoirs, primary feedback elements, controller commands, heater actuators, and reservoir heaters. As much as possible, the model took advantage of the available SINDA'85/FLUINT thermal/fluid integrator. Additional calculation logic and computer subroutines were developed to complete implementation of the model. The model is capable of simulating dynamic response of an equilibrium two-phase fluid reservoir. Modification of the model to include the liquid/vapor nonequilibrium is required for applications of the model to simulate performance of reservoir in which the liquid and vapor phases of the reservoir fluid are not in equilibrium. In addition, the model in its present form, needs to be refined in several respects. More empirical data are needed to guide the model development. The model may then be used to conduct a full parametric study of two-phase fluid reservoirs. More complexities in two-phaes flow regions in laboratory and flight conditions may have to be considered eventually if empirical data cannot be simulated satisfactorily. System with other components arrangement also need to be simulated if optimization is ever to be attained. The present model does, however, preliminarily demonstrates that such analyses are quite possible and offers a far less expensive method to understand the transient of a two-phase fluid reservoir system than a totally headware approach.

Chi, Joseph↗

NASA NTRSDate not supplied

Quantum-Assisted Variational Segmentation for Image-to-Image Wildfire Detection Using Satellite Data

The quantum computing community has been searching for suitable applications to demonstrate the potential of near-term quantum devices. Quantum machine learning is a potential candidate, particularly using models that cannot be efficiently simulated with classical computers [1, 2]. This work focuses on a transition phase of quantum computers where the quantum machine learning model is still simulable classically but projected not to be simulable as the size of the model grows. Ultimately quantum computers may have advantages for high-dimensional real-world problems. Due to the limited number of qubits in current noisy intermediate-scale quantum (NISQ) devices, the direct application of quantum computers in high dimensional data is not feasible. To remedy this problem, an encoder-decoder architecture can be utilized. The encoder model would transform the high-dimensional data into a compact representation, to a level that small quantum computers can be used today (or in the near future), and the decoder would take the quantum processed outputs back to the high-dimensional space. Addressing the two challenges of quantum machine learning, this work investigates a hybrid supervised generative model with a quantum Ising Born machine embedded as the latent distribution. The model contains four main parts (Figure 1.a.): (1) a U-NET architecture responsible for learning segmentation flow, (2) a Prior network responsible for learning an encoded latent distribution of the input data, (3) a Born machine which represents the latent distribution, and (4) a Posterior network in charge of learning the joint encoded latent distribution of inputs and target data. The initial model, proposed by [3], is optimized by (1) maximizing the overlap of the prior and posterior latent distributions, and (2) minimizing the segmentation loss. The proposed model is designed to be investigated in a simulation environment applied to the real-world application of wildfire segmentation. Specifically, the model is designed to solve the patchy wildfire segmentations of Moderate Resolution Imaging Spectroradiometer (MODIS) by taking the MODIS observations and using Visible Infrared Imaging Radiometer Suite’s (VIIRS) consistent wildfire product as the target. The model solves patchy wildfire segmentations and provides insight into the epistemic errors sourced from model variation. The model utilizes the Born machine as a QUBO solver to represent the latent space as a Bernoulli distribution. The proposed configuration allows the variational segmentation model to leverage the true quantum probabilistic nature and derive a more expressive latent configuration, increasing the model performance in describing wildfire segmentations. The quantum probabilistic information of the Born machine is directly incorporated in the Kullback-Leibler divergence loss in the prior and posterior distributions, forcing the Bernoulli latent distribution to maximize the overlap of input and joint input-target distributions. The proposed model is then trained and compared with a baseline only consisting of direct Bernoulli latent distribution with no Born machine representing the latent space. The models are evaluated based on the segmentation metrics, such as precision, recall, intersect of union, with uncertainty boundaries accounting for the stochastic nature of the model. Our findings show that even in low latent-dimensional space (due to the limit in computational power of the classical quantum simulator), we are able to effectively capture the latent representation and hence the model performs better than the baseline. The findings are a projection for scaling the model into higher dimensional latent space with the Born machine surpassing the baseline performance. Figure 1. Sub-figure (a) demonstrates the architecture for the training phase. The model consists of a Prior and Posterior network that encode inputs and joint input-target data into compact representations, respectively. The Born machine represents the latent distribution, and the U-NET branch learns the segmentation patterns of the data. The stochasticity is introduced to the U-NET through its last layer to create meaningful but stochastic segmentations. Sub-figure (b) represents the inference phase where the model takes the stochastic behavior from the prior network and injects that into the U-NET. Each attempt of inference will generate different but similar segmentations from the same distribution of the wildfire event. REFERENCES [1] Coyle, B., Mills, D., Danos, V., & Kashefi, E. (2020). The Born supremacy: quantum advantage and training of an Ising Born machine. npj Quantum Information, 6(1), 1-11. [2] Liu, J. G., & Wang, L. (2018). Differentiable learning of quantum circuit born machines. Physical Review A, 98(6), 062324. [3] Kohl, S., Romera-Paredes, B., Meyer, C., De Fauw, J., Ledsam, J. R., Maier-Hein, K., ... & Ronneberger, O. (2018). A probabilistic u-net for segmentation of ambiguous images. Advances in neural information processing systems, 31.

quantum machine learning↗

NASA NTRS2007

Biomedical Support of U.S. Extravehicular Activity

The world's first extravehicular activity (EVA) was performed by A. A. Leonov on March 18, 1965 during the Russian Voskhod-2 mission. The first US EVA was executed by Gemini IV astronaut Ed White on June 3, 1965, with an umbilical tether that included communications and an oxygen supply. A hand-held maneuvering unit (HHMU) also was used to test maneuverability during the brief EVA; however the somewhat stiff umbilical limited controlled movement. That constraint, plus difficulty returning through the vehicle hatch, highlighted the need for increased thermal control and improved EVA ergonomics. Clearly, requirements for a useful EVA were interrelated with the vehicle design. The early Gemini EVAs generated requirements for suits providing micro-meteor protection, adequate visual field and eye protection from solar visual and infrared radiation, gloves optimized for dexterity while pressurized, and thermal systems capable of protecting the astronaut while rejecting metabolic heat during high workloads. Subsequent Gemini EVAs built upon this early experience and included development of a portable environmental control and life support systems (ECLSS) and an astronaut maneuvering unit. The ECLSS provided a pressure vessel and controller with functional control over suit pressure, oxygen flow, carbon dioxide removal, humidity, and temperature control. Gemini EVA experience also identified the usefulness of underwater neutral buoyancy and altitude chamber task training, and the importance of developing reliable task timelines. Improved thermal management and carbon dioxide control also were required for high workload tasks. With the Apollo project, EVA activity was primarily on the lunar surface; and suit durability, integrated liquid cooling garments, and low suit operating pressures (3.75 pounds per square inch absolute [psia] or 25.8 kilopascal [kPa],) were required to facilitate longer EVAs with ambulation and significant physical workloads with average metabolic rates of 1000 BTU/hr and peaks of up to 2200 BTU/hr. Mobility was further augmented with the Lunar Roving Vehicle. The Apollo extravehicular mobility unit (EMU) was made up of over 15 components, ranging from a biomedical belt for capturing and transmitting biomedical data, urine and fecal containment systems, a liquid cooling garment, communications cap, a modular portable life support system (PLSS), a boot system, thermal overgloves, and a bubble helmet with eye protection. Apollo lunar astronauts performed successful EVAs on the lunar surface from a 5 psia (34.4 kPa) 100% oxygen environment in the Lunar Lander. A maximum of three EVAs were performed on any mission. For Skylab a modified A7LB suit, used for Apollo 15, was selected. The Skylab astronaut life support assembly (ALSA) provided umbilical support through the life support umbilical (LSU) and used open loop oxygen flow, rather than closed-loop as in Apollo missions. Thermal control was provided by liquid water circulated by spacecraft pumps and electrical power also was provided from the spacecraft via the umbilical. The cabin atmosphere of 5 psia (34.4 kPa), 70% oxygen, provided a normoxic atmosphere and because of the very low nitrogen partial pressures, no special protocols were required to protect against decompression sickness (DCS) as was the case with the Apollo spacecraft with a 5 psi, 100% oxygen environment.

Gernhardt, Michael L.↗

NASA NTRS2022

In Situ Measurements of Surface Texture with Virtual Environments Support Science-Driven Human Surface Operations on the Moon and Beyond

Visualization tools enabling real-time scientific analysis are important for supporting future astronaut operations on the lunar surface. Such tools can be built into virtual environments to support scientific investigations, as well as situational awareness, real-time decision making, and efficient communication between astronauts and ground and support systems. Understanding how these tools can be optimized for science is essential for upcoming Artemis missions. In this contribution, we discuss how measurements of surface texture at multiple length scales can greatly enhance in situ science on/of the Moon, and eventually Mars, asteroids, and beyond. Roughness measurements at various wavelengths directly support objectives defined in the Artemis Science Plan, including (O1) “understanding planetary processes,” (O2) “understanding volatile cycles,” and (O3) “interpreting the impact history of the Earth-Moon system” . Key scientific analyses enabled by texture measurements at different length scales include: ● Sub-centimeter scales: Texture measurements can help constrain lava flow crystallinity, lava rheology, emplacement flow dynamics, and cooling histories (O1). Measurements of lacunarity (voids in fractal fill space) can shed light on eruptive volatile content, residence time of migrating volatiles, and near-surface volume available for micro-cold trapping of volatiles (O1, O2). ● Centimeter–meter scales: Texture measurements can be used for the differentiation of individual lava flows, the reconstruction of local stratigraphies and emplacement sequences, characterization of post-emplacement surface modification processes (O1, O3). Derived roughness (polarization) metrics can be used in the detection of water ice and characterization of ice properties (e.g., purity, grade, depth, abundance). ● Hectometer–Kilometer scales: Texture measurements can be used to differentiate major geologic surface units and surface structures (O1), constrain the presence of abundant ground ices (O2), and analyze surface modification and estimate surface age (O3). Real-time measurements of surface texture across these multiple length scales will enable efficient sample identification and scientific investigations by future astronauts. To support these investigations and the objective classification of surface texture, virtual environments employed by astronauts should be able to instantaneously convert raw data into processed data (e.g., digital terrain and elevation models) and derived metrics (e.g., RMS, std, Hurst, CPR) and perform statistical analyses (e.g., PCA, outliers, correlation matrices). Such tools are being developed and tested by the Resource Exploration and Science of our Cosmic Environment (RESOURCE) team, a node of NASA’s Solar System Exploration Research Virtual Institute (SSERVI), and are an excellent example of the powerful synergies of human and robotic ground assets critical in the return of humans to the Moon.

Ariel N. Deutsch↗