Search NASA⌕ Search

SEARCH · Search NASA

Results for “fuel performance codes”

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.

267 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

NDARC - NASA Design and Analysis of Rotorcraft: Input - Appendix 3

The NDARC code performs design and analysis tasks. The design task involves sizing the rotorcraft to satisfy specified design conditions and missions. The analysis tasks can include off-design mission performance analysis, flight performance calculation for point operating conditions, and generation of subsystem or component performance maps. Figure 1-1 illustrates the tasks. The principal tasks (sizing, mission analysis, flight performance analysis) are shown in the figure as boxes with heavy borders. Heavy arrows show control of subordinate tasks. The aircraft description (figure 1-1) consists of all the information, input and derived, that defines the aircraft. The aircraft consists of a set of components, including fuselage, rotors, wings, tails, and propulsion. This information can be the result of the sizing task; can come entirely from input, for a fixed model; or can come from the sizing task in a previous case or previous job. The aircraft description information is available to all tasks and all solutions (indicated by light arrows). The sizing task determines the dimensions, power, and weight of a rotorcraft that can perform a specified set of design conditions and missions. The aircraft size is characterized by parameters such as design gross weight, weight empty, rotor radius, and engine power available. The relations between dimensions, power, and weight generally require an iterative solution. From the design flight conditions and missions, the task can determine the total engine power or the rotor radius (or both power and radius can be fixed), as well as the design gross weight, maximum takeoff weight, drive system torque limit, and fuel tank capacity. For each propulsion group, the engine power or the rotor radius can be sized. Missions are defined for the sizing task, and for the mission performance analysis. A mission consists of a number of mission segments, for which time, distance, and fuel burn are evaluated. For the sizing task, certain missions are designated to be used for design gross weight calculations; for transmission sizing; and for fuel tank sizing. The mission parameters include mission takeoff gross weight and useful load. For specified takeoff fuel weight with adjustable segments, the mission time or distance is adjusted so the fuel required for the mission (burned plus reserve) equals the takeoff fuel weight. The mission iteration is on fuel weight or energy. Flight conditions are specified for the sizing task, and for the flight performance analysis. For the sizing task, certain flight conditions are designated to be used for design gross weight calculations; for transmission sizing; for maximum takeoff weight calculations; and for antitorque or auxiliary thrust rotor sizing. The flight condition parameters include gross weight and useful load. For flight conditions and mission takeoff, the gross weight can be maximized, such that the power required equals the power available. A flight state is defined for each mission segment and each flight condition. The aircraft performance can be analyzed for the specified state, or a maximum effort performance can be identified. The maximum effort is specified in terms of a quantity such as best endurance or best range, and a variable such as speed, rate of climb, or altitude. The aircraft must be trimmed, by solving for the controls and motion that produce equilibrium in the specified flight state. Different trim solution definitions are required for various flight states. Evaluating the rotor hub forces may require solution of the blade flap equations of motion.

NDARC↗
NASA NTRS1991

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗
NASA NTRS2010

Optical Characterization of a Multipoint Lean Direct Injector for Gas Turbine Combustors: Velocity and Fuel Drop Size Measurements

Performance of a multipoint, lean direct injection (MP-LDI) strategy for low emission aero-propulsion systems has been tested in a Jet-A fueled, lean flame tube combustion rig. Operating conditions for the series of tests included inlet air temperatures between 672 and 828 K, pressures between 1034 and 1379 kPa and total equivalence ratios between 0.41 and 0.45, resulting in equilibrium flame temperatures approaching 1800 K. Ranges of operation were selected to represent the spectrum of subsonic and supersonic flight conditions projected for the next-generation of commercial aircraft. This document reports laser-based measurements of in situ fuel velocities and fuel drop sizes for the NASA 9-point LDI hardware arranged in a 3 3 square grid configuration. Data obtained represent a region of the flame tube combustor with optical access that extends 38.1-mm downstream of the fuel injection site. All data were obtained within reacting flows, without particle seeding. Two diagnostic methods were employed to evaluate the resulting flow path. Three-component velocity fields have been captured using phase Doppler interferometry (PDI), and two-component velocity distributions using planar particle image velocimetry (PIV). Data from these techniques have also offered insight into fuel drop size and distribution, fuel injector spray angle and pattern, turbulence intensity, degree of vaporization and extent of reaction. This research serves to characterize operation of the baseline NASA 9- point LDI strategy for potential use in future gas-turbine combustor applications. An additional motive is the compilation of a comprehensive database to facilitate understanding of combustor fuel injector aerodynamics and fuel vaporization processes, which in turn may be used to validate computational fluid dynamics codes, such as the National Combustor Code (NCC), among others.

Heath, Christopher M.↗
NASA NTRS2001

Advanced Methods for Aircraft Engine Thrust and Noise Benefits: Nozzle-Inlet Flow Analysis

The Fluid Mechanics and Acoustics Laboratory at Hampton University (HU/FM&AL) jointly with the NASA Glenn Research Center has conducted four connected subprojects under the reporting project. Basically, the HU/FM&AL Team has been involved in joint research with the purpose of theoretical explanation of experimental facts and creation of accurate numerical simulation techniques and prediction theory for solution of current problems in propulsion systems of interest to the NAVY and NASA agencies. This work is also supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations and possible experimental tests at the Hampton University campus. The fundamental idea uniting these subprojects is to use nontraditional 3D corrugated and composite nozzle and inlet designs and additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. These subprojects are: (1) Aeroperformance and acoustics of Bluebell-shaped and Telescope-shaped designs; (2) An analysis of sharp-edged nozzle exit designs for effective fuel injection into the flow stream in air-breathing engines: triangular-round, diamond-round and other nozzles; (3) Measurement technique improvement for the HU Low Speed Wind Tunnel; a new course in the field of aerodynamics, teaching and training of HU students; experimental tests of Mobius-shaped screws: research and training; (4) Supersonic inlet shape optimization. The main outcomes during this reporting period are: (l) Publications: The AIAA Paper #00-3170 was presented at the 36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, 17-19 June, 2000, Huntsville, AL. The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01 -3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, being held on 08-11 July, in Salt Lake City, UT; (2) A U.S. patent #6,082,635 was granted on July 4, 2000; (3) Grants and proposals: The H U/ FM&AL was awarded the NASA grant NAG-3-2495 in October 2000 and the laboratory is a primary U.S. research team in a joint project under the CRDF award granted to the NASA GRC and IM/MSU (Russia) in July 2000; (4) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc., have been conducted for several complicated 2D and 3D nozzle and inlet designs using NASA, ICASE, and IM/MSU codes based on full Euler and Navier-Stokes solvers: CFL3D, FLUENT, and GODUNOV, and others; (5) Experimental Tests: (a) A new course: "Advanced Aerodynamics and Aircraft Performance" presented in spring semester, 2001; training and experimental test research using the HU LSWT. (b) Small-scale M6bius-shaped screws were tested in different conditions and their application has shown essential benefits by comparison with traditional designs; (6) Installation in the FM&AL computer system: second software TECPLOT 8.0 for the UNIX SGI workstation and free TECPLOT 7.5 for the PC Dell computer, and 2D and 3D GRIDGEN (version 9) for the UNIX SGI as well as installation of two free NASA codes, 3D MAG and VULCAN; (7) Student Research Activity: Involvement of two undergraduate students as research assistants in the current research project.

Gilinsky, Mikhail↗
NASA NTRS2001

CFD Simulation of Liquid Rocket Engine Injectors

Detailed design issues associated with liquid rocket engine injectors and combustion chamber operation require CFD methodology which simulates highly three-dimensional, turbulent, vaporizing, and combusting flows. The primary utility of such simulations involves predicting multi-dimensional effects caused by specific injector configurations. SECA, Inc. and Engineering Sciences, Inc. have been developing appropriate computational methodology for NASA/MSFC for the past decade. CFD tools and computers have improved dramatically during this time period; however, the physical submodels used in these analyses must still remain relatively simple in order to produce useful results. Simulations of clustered coaxial and impinger injector elements for hydrogen and hydrocarbon fuels, which account for real fluid properties, is the immediate goal of this research. The spray combustion codes are based on the FDNS CFD code' and are structured to represent homogeneous and heterogeneous spray combustion. The homogeneous spray model treats the flow as a continuum of multi-phase, multicomponent fluids which move without thermal or velocity lags between the phases. Two heterogeneous models were developed: (1) a volume-of-fluid (VOF) model which represents the liquid core of coaxial or impinger jets and their atomization and vaporization, and (2) a Blob model which represents the injected streams as a cloud of droplets the size of the injector orifice which subsequently exhibit particle interaction, vaporization, and combustion. All of these spray models are computationally intensive, but this is unavoidable to accurately account for the complex physics and combustion which is to be predicted, Work is currently in progress to parallelize these codes to improve their computational efficiency. These spray combustion codes were used to simulate the three test cases which are the subject of the 2nd International Workshop on-Rocket Combustion Modeling. Such test cases are considered by these investigators to be very valuable for code validation because combustion kinetics, turbulence models and atomization models based on low pressure experiments of hydrogen air combustion do not adequately verify analytical or CFD submodels which are necessary to simulate rocket engine combustion. We wish to emphasize that the simulations which we prepared for this meeting are meant to test the accuracy of the approximations used in our general purpose spray combustion models, rather than represent a definitive analysis of each of the experiments which were conducted. Our goal is to accurately predict local temperatures and mixture ratios in rocket engines; hence predicting individual experiments is used only for code validation. To replace the conventional JANNAF standard axisymmetric finite-rate (TDK) computer code 2 for performance prediction with CFD cases, such codes must posses two features. Firstly, they must be as easy to use and of comparable run times for conventional performance predictions. Secondly, they must provide more detailed predictions of the flowfields near the injector face. Specifically, they must accurately predict the convective mixing of injected liquid propellants in terms of the injector element configurations.

Farmer, Richard↗
NASA NTRS2000

A Collaborative Analysis Tool for Integrated Hypersonic Aerodynamics, Thermal Protection Systems, and RBCC Engine Performance for Single Stage to Orbit Vehicles

Presented is a computer-based tool that connects several disciplines that are needed in the complex and integrated design of high performance reusable single stage to orbit (SSTO) vehicles. Every system is linked to every other system, as is the case of SSTO vehicles with air breathing propulsion, which is currently being studied by NASA. An RBCC propulsion system integrates airbreathing and rocket propulsion into a single engine assembly enclosed within a cowl or duct. A typical RBCC propulsion system operates as a ducted rocket up to approximately Mach 3. Then there is a transition to a ramjet mode for supersonic-to-hypersonic acceleration. Around Mach 8 the engine transitions to a scramjet mode. During the ramjet and scramjet modes, the integral rockets operate as fuel injectors. Around Mach 10-12 (the actual value depends on vehicle and mission requirements), the inlet is physically closed and the engine transitions to an integral rocket mode for orbit insertion. A common feature of RBCC propelled vehicles is the high degree of integration between the propulsion system and airframe. At high speeds the vehicle forebody is fundamentally part of the engine inlet, providing a compression surface for air flowing into the engine. The compressed air is mixed with fuel and burned. The combusted mixture must be expanded to an area larger than the incoming stream to provide thrust. Since a conventional nozzle would be too large, the entire lower after body of the vehicle is used as an expansion surface. Because of the high external temperatures seen during atmospheric flight, the design of an airbreathing SSTO vehicle requires delicate tradeoffs between engine design, vehicle shape, and thermal protection system (TPS) sizing in order to produce an optimum system in terms of weight (and cost) and maximum performance. To adequately determine the performance of the engine/vehicle, the Hypersonic Flight Inlet Model (HYFIM) module was designed to interface with the RBCC engine model. HYFIM performs the aerodynamic analysis of forebodies and inlet characteristics of RBCC powered SSTO launch vehicles. HYFIM is applicable to the analysis of the ramjet/scramjet engine operations modes (Mach 3-12), and provides estimates of parameters such as air capture area, shock-on-lip Mach number, design Mach number, compression ratio, etc., based on a basic geometry routine for modeling axisymmetric cones, 2-D wedge geometries. HYFIM also estimates the variation of shock layer properties normal to the forebody surface. The thermal protection system (TPS) is directly linked to determination of the vehicle moldline and the shaping of the trajectory. Thermal protection systems to maintain the structural integrity of the vehicle must be able to mitigate the heat transfer to the structure and be lightweight. Herein lies the interdependency, in that as the vehicle's speed increases, the TPS requirements are increased. And as TPS masses increase the effect on the propulsion system and all other systems is compounded. The need to analyze vehicle forebody and engine inlet is critical to be able to design the RBCC vehicle. To adequately determine insulation masses for an RBCC vehicle, the hypersonic aerodynamic environment and aeroheating loads must be calculated and the TPS thicknesses must be calculated for the entire vehicle. To accomplish this an ascent or reentry trajectory is obtained using the computer code Program to Optimize Simulated Trajectories (POST). The trajectory is then used to calculate the convective heat rates on several locations on the vehicles using the Miniature Version of the JA70 Aerodynamic Heating Computer Program (MINIVER). Once the heat rates are defined for each body point on the vehicle, then insulation thicknesses that are required to maintain the vehicle within structural limits are calculated using Systems Improved Numerical Differencing Analyzer (SINDA) models. If the TPS masses are too heavy for the performance of the vehicle the process may be repeated altering the trajectory or some other input to reduce the TPS mass. E-PSURBCC is an "engine performance" model and requires the specification of inlet air static temperature and pressure as well as Mach number (which it pulls from the HYFIM and POST trajectory files), and calculates the corresponding stagnation properties. The engine air flow path geometry includes inlet, a constant area section where the rocket is positioned, a subsonic diffuser, a constant area afterburner, and either a converging nozzle or a converging-diverging nozzle. The current capabilities of E-PSURBCC ejector and ramjet mode treatment indicated that various complex flow phenomena including multiple choking and internal shocks can occur for combinations of geometry/flow conditions. For a given input deck defining geometry/flow conditions, the program first goes through a series of checks to establish whether the input parameters are sound in terms of a solution path. If the vehicle/engine performance fails mission goals, the engineer is able to collaboratively alter the vehicle moldline to change aerodynamics, or trajectory, or some other input to achieve orbit. The problem described is an example of the need for collaborative design and analysis. RECIPE is a cross-platform application capable of hosting a number of engineers and designers across the Internet for distributed and collaborative engineering environments. Such integrated system design environments allow for collaborative team design analysis for performing individual or reduced team studies. To facilitate the larger number of potential runs that may need to be made, RECIPE connects the computer codes that calculate the trajectory data, aerodynamic data based on vehicle geometry, heat rate data, TPS masses, and vehicle and engine performance, so that the output from each tool is easily transferred to the model input files that need it.

Stanley, Thomas Troy↗
NASA NTRS1996

Tilt Rotor Aircraft Aeroacoustics

A fleet of civil tilt rotor transports offers a means of reducing airport congestion and point-to-point travel time. The speed, range, and fuel economy of these aircraft, along with their efficient use of vertiport area, make them good candidates for short-to-medium range civil transport. However, to be successfully integrated into the civilian community, the tilt rotor must be perceived as a quiet, safe, and economical mode of transportation that does not harm the environment. In particular, noise impact has been identified as a possible barrier to the civil tilt rotor. Along with rotor conversion-mode flight, and blade-vortex interaction noise during descent, hover mode is a noise problem for tilt rotor operations. In the present research, tilt rotor hover aeroacoustics have been studied analytically, experimentally, and computationally. Various papers on the subject were published as noted in the list of publications. More recently, experimental measurements were made on a 1/12.5 scale model of the XV-15 in hover and analyses of this data and extrapolations to full scale were also carried out. A dimensional analysis showed that the model was a good aeroacoustic approximation to the full-scale aircraft, and scale factors were derived to extrapolate the model measurements to the full-scale XV-15. The experimental measurements included helium bubble flow visualization, silk tuft flow visualization, 2-component hot wire anemometry, 7-hole pressure probe measurements, vorticity measurements, and outdoor far field acoustic measurements. The hot wire measurements were used to estimate the turbulence statistics of the flow field into the rotors, such as length scales, velocity scales, dissipation, and turbulence intermittency. Several different configurations of the model were tested: (1) standard configurations (single isolated rotor, two rotors without the aircraft, standard tilt rotor configuration); (2) flow control devices (the 'plate', the 'diagonal fences'); (3) basic configuration changes (increasing the rotor/rotor spacing, reducing the rotor plane/wing clearance. operating the rotors out of phase). Also, an approximation to Sikorsky's Variable Diameter Tilt Rotor (VDTR) configuration was tested, and some flow measurements were made on a semi-span configuration of the model. Acoustic predictions were made using LOWSON.M, a Mathematica code Mean aerodynamic models were developed based on hover performance predictions from HOVER.FOR. This hover prediction code used blade element theory for the aerodynamics, and Prandtl's Vortex theory to model the wake, along with empirical formulas for the effects of Reynolds number, Mach number, and stall. Aerodynamic models were developed from 7-hole pressure probe measurements of the mean velocity into the model rotors.

George, Albert R.↗
NASA NTRS2013

NASA Tech Briefs, January 2013

Topics include: Single-Photon-Sensitive HgCdTe Avalanche Photodiode Detector; Surface-Enhanced Raman Scattering Using Silica Whispering-Gallery Mode Resonators; 3D Hail Size Distribution Interpolation/Extrapolation Algorithm; Color-Changing Sensors for Detecting the Presence of Hypergolic Fuels; Artificial Intelligence Software for Assessing Postural Stability; Transformers: Shape-Changing Space Systems Built with Robotic Textiles; Fibrillar Adhesive for Climbing Robots; Using Pre-Melted Phase Change Material to Keep Payloads in Space Warm for Hours without Power; Development of a Centrifugal Technique for the Microbial Bioburden Analysis of Freon (CFC-11); Microwave Sinterator Freeform Additive Construction System (MS-FACS); DSP/FPGA Design for a High-Speed Programmable S-Band Space Transceiver; On-Chip Power-Combining for High-Power Schottky Diode-Based Frequency Multipliers; FPGA Vision Data Architecture; Memory Circuit Fault Simulator; Ultra-Compact Transputer-Based Controller for High-Level, Multi-Axis Coordination; Regolith Advanced Surface Systems Operations Robot Excavator; Magnetically Actuated Seal; Hybrid Electrostatic/Flextensional Mirror for Lightweight, Large-Aperture, and Cryogenic Space Telescopes; System for Contributing and Discovering Derived Mission and Science Data; Remote Viewer for Maritime Robotics Software; Stackfile Database; Reachability Maps for In Situ Operations; JPL Space Telecommunications Radio System Operating Environment; RFI-SIM: RFI Simulation Package; ION Configuration Editor; Dtest Testing Software; IMPaCT - Integration of Missions, Programs, and Core Technologies; Integrated Systems Health Management (ISHM) Toolkit; Wind-Driven Wireless Networked System of Mobile Sensors for Mars Exploration; In Situ Solid Particle Generator; Analysis of the Effects of Streamwise Lift Distribution on Sonic Boom Signature; Rad-Tolerant, Thermally Stable, High-Speed Fiber-Optic Network for Harsh Environments; Towed Subsurface Optical Communications Buoy; High-Collection-Efficiency Fluorescence Detection Cell; Ultra-Compact, Superconducting Spectrometer-on-a-Chip at Submillimeter Wavelengths; UV Resonant Raman Spectrometer with Multi-Line Laser Excitation; Medicine Delivery Device with Integrated Sterilization and Detection; Ionospheric Simulation System for Satellite Observations and Global Assimilative Model Experiments - ISOGAME; Airborne Tomographic Swath Ice Sounding Processing System; flexplan: Mission Planning System for the Lunar Reconnaissance Orbiter; Estimating Torque Imparted on Spacecraft Using Telemetry; PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils; Multiple-Frame Detection of Subpixel Targets in Thermal Image Sequences; Metric Learning to Enhance Hyperspectral Image Segmentation; Basic Operational Robotics Instructional System; Sheet Membrane Spacesuit Water Membrane Evaporator; Advanced Materials and Manufacturing for Low-Cost, High-Performance Liquid Rocket Combustion Chambers; Motor Qualification for Long-Duration Mars Missions.

Source record↗
NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗
NASA NTRS2012

Analysis of Roll Steering for Solar Electric Propulsion Missions

Nothing is more vital to a spacecraft than power. Solar Electric Propulsion (SEP) uses that power to provide a safe, reliable, and, most importantly, fuel efficient means to propel a spacecraft to its destination. The power performance of an SEP vehicle s solar arrays and electrical power system (EPS) is largely influenced by the environment in which the spacecraft is operating. One of the most important factors that determines solar array power performance is how directly the arrays are pointed to the sun. To get the most power from the solar arrays, the obvious solution is to point them directly at the sun at all times. Doing so is not a problem in deep space, as the environment and pointing conditions that a spacecraft faces are fairly constant and are easy to accommodate, if necessary. However, large and sometimes rapid variations in environmental and pointing conditions are experienced by Earth orbiting spacecraft. SEP spacecraft also have the additional constraint of needing to keep the thrust vector aligned with the velocity vector. Thus, it is important to analyze solar array power performance for any vehicle that spends an extended amount of time orbiting the Earth, and to determine how much off-pointing can be tolerated to produce the required power for a given spacecraft. This paper documents the benefits and drawbacks of perfectly pointing the solar arrays of an SEP spacecraft spiraling from Earth orbit, and how this might be accomplished. Benefits and drawbacks are defined in terms of vehicle mass, power, volume, complexity, and cost. This paper will also look at the application of various solar array pointing methods to future missions. One such pointing method of interest is called roll steering . Roll steering involves rolling the entire vehicle twice each orbit. Roll steering, combined with solar array gimbal tracking, is used to point the solar arrays perfectly towards the sun at all points in the orbit, while keeping the vehicle thrusters aligned in the velocity direction. Roll steering is particularly attractive for a recently proposed mission that involves a spiral trajectory from low Earth orbit (LEO) to the Earth-Moon Lagrange Point 1 (E-M L1). During the spiral, the spacecraft will spend over 300 days experiencing the full spectrum of near-earth environments and solar array pointing conditions. An extensive study of the application of SEP (and roll steering) to this spiral mission is included, highlighting the ultimate goal of reduced vehicle cost and mass. Tools used for this analysis include the Systems Power Analysis for Capability Evaluation (Refs. 1 and 2) (SPACE) electrical power systems code, and SEP trajectory simulation tools developed at NASA Glenn Research Center.

Pederson, Dylan, M.↗
NASA NTRSDate not supplied

Initial Development of A Digital Twin Model for an Electrified Aircraft Propulsion Emulation Rig

In support of aviation fuel burn and emission reduction goals, NASA is pursing high-payoff research investments that promise to transform aviation. This includes investments in Electrified Aircraft Propulsion (EAP), which relies on the generation, storage, transmission, and use of electrical power for producing thrust and optimizing propulsion system efficiency. Multiple technology challenges must be addressed to unlock the full potential of EAP. This includes advances in propulsion controls, which will be vital for ensuring coordinated efficient operation of the complex integrated subsystems that comprise EAP architectures. To support EAP controls research, the NASA Glenn Research Center has developed the Hybrid Propulsion Emulation Rig (HyPER). The HyPER laboratory hardware includes shaft-mounted electric machines, power converters, power supplies, power distribution cables, and an energy storage device that can be reconfigured to represent a variety of EAP architectures. It also includes an integrated real-time computer system that hosts developed EAP control software and turbomachinery simulations. This enables the electrical system and rotating shafts of EAP designs to be implemented in actual hardware and integrated with turbomachinery simulations and system-level EAP control logic implemented in software. In this form, the HyPER laboratory provides a partially simulated, partially hardware-in-the-loop test environment enabling the initial development and evaluation of EAP control technology. A prerequisite for the development of EAP control designs is the availability of a system model that accurately reflects the operation of the electrical system hardware. To support this need, a digital twin model of the HyPER electrical system hardware is under development. This model is being coded in the MATLAB Simulink environment and uses the NASA-developed Electrical Modeling and Thermal Analysis Toolbox (EMTAT) to construct a digital twin framework. EMTAT contains generic electrical component building blocks that are simulated at turbomachinery timescales. Associated inputs and outputs allow the blocks to be combined to model complete electrical systems. The EMTAT blocks also contain adjustable internal maps and parameters that can be set to reflect the operation of a specific electrical component. For the HyPER digital twin, the settings of these EMTAT block internal maps and parameters is determined through machine learning approaches applied to characterization run data collected from the laboratory. During characterization runs the laboratory electrical system hardware is subjected to a full range of torque, speed, and power settings. Acquired data is then used to estimate EMTAT block parameters using a variety of machine learning techniques. The resulting digital twin model is found to match the operation of actual HyPER hardware with an accuracy suitable for control development purposes. It also holds promise for other applications including modeling the performance of HyPER laboratory reconfigurations and model-based anomaly detection. Planned follow-on work to automate post-processing of acquired laboratory data to update the HyPER digital twin model will also be presented and discussed.

Electrified Aircraft Propulsion↗
NASA NTRS1998

The NASA Low-Pressure Turbine Flow Physics Program

An overview of the NASA Lewis Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The program was established in response to the aero-engine industry's need for improved LPT efficiency and designs. Modern jet engines have four to seven LPT stages, significantly contributing to engine weight. In addition, there is a significant efficiency degradation between takeoff and cruise conditions, of up to 2 points. Reducing the weight and part count of the LPT and minimizing the efficiency degradation will translate into fuel savings. Accurate prediction methods of LPT flows and losses are needed to accomplish those improvements. The flow in LPT passages is at low Reynolds number, and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The affecting parameters traditionally considered are Reynolds number, freestream turbulence intensity, wake frequency parameter, and the pressure distribution (loading). Three-dimensional effects and additional parameters, particularly turbulence characteristics like length scales, spectra and other statistics, as well as wake turbulence intensity and properties also play a role. The flow of most interest is on the suction surface, where large losses are generated as the flow tends to separate at the low Reynolds numbers. Ignoring wakes, a common flow scenario, there is laminar separation, followed by transition on the separation bubble and turbulent reattachment. If transition starts earlier the separation will be eliminated and the boundary layer will be attached leading to the well known bypass transition issues. In contrast, transition over a separation bubble is closer to free shear layer transition and was not investigated as well, particularly in the turbine environment. Unsteadiness created by wakes complicates the picture. Wakes induce earlier transition, and the calmed regions trailing the induced turbulent spots can delay or eliminate separation via shear stress modification. Three-dimensional flow physics and geometry will have strong effects. Altogether a very complex and challenging problem emerges. The objective of the program is to provide improved models and physical understanding of the complex flow, which are essential for accurate prediction of flow and losses in the LPT. Experimental, computational and analytical work as complementing and augmenting approaches are used. The program involves industry, universities and research institutes, and other government laboratories. It is characterized by strong interaction among participants, quick dissemination of results, and responsiveness to industry's needs. The presentation will describe the work elements. Highlighting some activities in progress are experiments on simulated blade suction surface in low-speed wind tunnels, on curved wall, and on a flat-plate, both with pressure gradient. In the area of computation, assessment of existing models is performed using RANS (Reynolds Averaged Navier Stokes) simulations. Laminar flow DNS was completed. Analytical studies of instability and receptivity in attached and separated flows were started. In the near future the program is moving to include wake effects and development of improved modeling. Experimental work in preparation stages are: (1) Addition of wakes to the curved tunnel experiment; (2) Low-speed rotating rig experiment on GE90 engine LPT; and (3) Transonic cascade. In the area of computation, it is expected to move from model assessment towards development of improved models. In addition, a new project of Large Eddy Simulation (LES) of LPT is to begin and will provide numerical data bases. It is planned to implement the emerging improved models in a multistage turbomachinery code and to validate against the GE90 engine LPT.

Ashpis, David E.↗
NASA NTRS2018

The Fusion Driven Rocket: Nuclear Propulsion through Direct Conversion of Fusion Energy

The future of manned space exploration and development of space depends critically on the creation of a dramatically more efficient propulsion architecture for in-space transportation. A very persuasive reason for investigating the applicability of nuclear power in rockets is the vast energy density gain of nuclear fuel when compared to chemical combustion energy. The Fusion Driven rocket (FDR) represents a revolutionary approach to fusion propulsion where the power source releases its energy directly into the propellant, not requiring conversion to electricity. It employs a solid lithium propellant that requires no significant tankage mass. The propellant is rapidly heated and accelerated to high exhaust velocity (> 30 km/s), while having no substantial physical interaction with the spacecraft thereby avoiding damage to the rocket and limiting both the thermal heat load and radiator mass. The key to achieving this stems from research at MSNW and the UW on the magnetically driven implosion of metal foils onto a magnetized plasma target to obtain fusion conditions. A logical extension of this work leads to a method that utilizes these metal shells (or liners) to not only achieve fusion conditions, but to serve as the propellant as well. Several low-mass, magnetically driven metal liners are inductively driven to converge radially and axially and form a thick blanket surrounding the target plasmoid and compress the plasmoid to fusion conditions. Virtually all of the radiant, neutron and particle energy from the plasma is absorbed by the encapsulating, thick metal blanket thereby isolating the spacecraft from the fusion. This energy, in addition to the intense Ohmic heating at peak magnetic field compression, is adequate to vaporize and ionize the metal blanket. The expansion of this hot, ionized metal propellant through a magnetically insulated nozzle produces high thrust at the optimal Isp. The energy from the fusion process, along with the waste heat, is thus utilized at very high efficiency. The basic scheme for FDR is illustrated and described in the report (see Fig. 2) The two most critical issues in meeting challenges introduced employing magneto-inertial fusion as the power source is driver efficiency and “stand-off” – the ability to isolate and protect fusion and thruster from the resultant fusion energy. By employing metal shells for compression, it is possible to produce the desired convergent motion inductively by inserting the metal sheets along the inner surface of cylindrical or conically tapered coils. Both stand-off and energy efficiency issues are solved by this arrangement. 3 This two year effort focused on achieving three key criteria for the Fusion Driven Rocket to move forward for technological development: (1) the physics of the FDR must be fully understood and validated, (2) the design and technology development for the FDR required for its implementation in space must be fully characterized, and (3) an in-depth analysis of the rocket design and spacecraft integration as well as mission architectures enabled by the FDR need to be performed. A subscale, laboratory liner compression test facility was assembled at the University of Washington Plasma Dynamics Laboratory with sufficient liner kinetic energy (~ 0.5 MJ) to reach conditions required for fusion breakeven conditions. Detailed experimental studies of the dynamic behavior of the driven liners as well as liner convergence and magnetic compression were performed. The development of both the 1D liner dynamics code and the full 3D ANSYS® liner calculations was achieved. The characterization of both the FDR and spacecraft as well as a design architecture analysis was conducted that included an examination of a wide range of mission architectures and destinations for which this fusion propulsion system would be enabling or critical. In particular a rapid, single launch manned Mars mission was developed.

Energy↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗
NASA NTRS1997

Comparative Soot Diagnostics: Preliminary Results

The motivation for the Comparative Soot Diagnostics (CSD) experiment lies in the broad practical importance of understanding combustion generated particulate. Depending upon the circumstances, particulate matter can affect the durability and performance of combustion equipment, can be a pollutant, can be used to detect fires and, in the form of soot, can be the dominant source of radiant energy from flames. The nonbuoyant structure of most flames of practical interest makes understanding of soot processes in low gravity flames important to our ability to predict fire behavior on earth. These studies also have direct applications to fire safety in human-crew spacecraft, since smoke is the indicator used for automated detection in current spacecraft. In the earliest missions (Mercury, Gemini and Apollo), the crew quarters were so cramped that it was considered reasonable that the astronauts would rapidly detect any fire. The Skylab module, however, included approximately 20 UV-sensing fire detectors. The Space Shuttle has 9 particle-ionization smoke detectors in the mid-deck and flight deck and Spacelab has six additional particle-ionization smoke detectors. The designated detectors for the ISS are laser-diode, forward-scattering, smoke or particulate detectors. Current plans for the ISS call for two detectors in the open area of the module, and detectors in racks that have both cooling air flow and electrical power. Due to the complete absence of data concerning the nature of particulate and radiant emission from incipient and fully developed low-g fires, all three of these detector systems were designed based upon l-g test data and experience. As planned mission durations and complexity increase and the volume of spacecraft increases, the need for and importance of effective, crew-independent, fire detection grows significantly. To provide this level of protection, more knowledge is needed concerning low-gravity fire phenomena and, in particular, how they might be detected and suppressed. Prior to CSD, no combustion-generated particulate samples had been collected near the flame zone for well-developed microgravity flames. All of the extant data either came from drop tower tests and therefore only corresponded to the early stages of a fire or were collected far from the flame zone. The fuel sources in the drop tower tests were restricted to laminar gas-jet diffusion flames and very rapidly overheated wire insulation. The gas-jet tests indicated, through thermophoretic sampling, (2) that soot primaries and aggregates (groups of primary particles) in low-gravity may be significantly larger than those in normal gravity (1-g). This raises new scientific questions about soot processes as well as practical issues for particulate size sensitivity and detection alarm threshold levels used in on-orbit smoke detectors. Preliminary tests in the 2.2 second drop tower suggest that particulate generated by overheated wire insulation may be larger in low-g than in 1-g. Transmission Electron Microscope (TEM) grids downstream of the fire region in the Wire Insulation Flammability experiment as well as visual observation of long string-like aggregates, further confirm this suggestion. The combined impact of these limited results and theoretical predictions is that, as opposed to extrapolation from l-g data, direct knowledge of low-g combustion particulate is needed for more confident design of smoke detectors for spacecraft. This paper describes the operation and preliminary results of the CSD, a project conceived and developed at NASA Lewis Research Center. The CSD flight experiment was conducted in the Middeck Glovebox Facility (MGBX) on USMP-3. The project is support by NASA Headquarters Microgravity Science and Applications Division and Code Q. The results presented here are from the microgravity portion of the experiment, including the temporal response of the detectors and average sizes of the primary and aggregate particles captured on the thermophoretic probes.

Urban, David L.↗