SEARCH · Engineering Papers
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.
Boiling heat transfer coefficients, interface behavior, and vapor quality in heated wall rotating boilers with continuous through-flow of water as test fluid
Fluid flow problems during orbital refueling vehicle maneuvers and system operations in propellant transfer, including disturbances, liquid-vapor interface instability, vapor ingestion, etc
A long-life hydraulic converter with unique features was conceptually designed to interface with a specified 15 kW(e) free-piston Stirling engine in a solar thermal dish application. Hydraulic fluid at 34.5 MPa (5000 psi) is produced to drive a conventional hydraulic motor and rotary alternator. Efficiency of the low-maintenance converter design was calculated at 93.5% for a counterbalanced version and 97.0% without the counterbalance feature. If the converter were coupled to a Stirling engine with design parameters more typcial of high-technology Stirling engines, counterbalanced converter efficiency could be increased to 99.6%. Dynamic computer simulation studies were conducted to evaluate performance and system sensitivities. Production costs of the complete Stirling hydraulic/electric power system were evaluated at $6506 which compared with $8746 for an alternative Stirling engine/linear alternator system.
Many previous researchers have modeled sheet cavitation by means of a constant pressure solution in the cavity region coupled with a velocity potential formulation for the outer flow. The present paper discusses the issues involved in extending these cavitation models to Euler or Navier-Stokes codes. The approach taken is to start from a velocity potential model to ensure our results are compatible with those of previous researchers and available experimental data, and then to implement this model in both Euler and Navier-Stokes codes. The model is then augmented in the Navier-Stokes code by the inclusion of the energy equation which allows the effect of subcooling in the vicinity of the cavity interface to be modeled to take into account the experimentally observed reduction in cavity pressures that occurs in cryogenic fluids such as liquid hydrogen. Although our goal is to assess the practicality of implementing these cavitation models in existing three-dimensional, turbomachinery codes, the emphasis in the present paper will center on two-dimensional computations, most specifically isolated airfoils and cascades. Comparisons between velocity potential, Euler and Navier-Stokes implementations indicate they all produce consistent predictions. Comparisons with experimental results also indicate that the predictions are qualitatively correct and give a reasonable first estimate of sheet cavitation effects in both cryogenic and non-cryogenic fluids. The impact on CPU time and the code modifications required suggests that these models are appropriate for incorporation in current generation turbomachinery codes.
This work describes the application of a control theory-based aerodynamic shape optimization method to the problem of supersonic aircraft design. The design process is greatly accelerated through the use of both control theory and a parallel implementation on distributed memory computers. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is then implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) Standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on higher order computational fluid dynamics methods (CFD). In our earlier studies, the serial implementation of this design method was shown to be effective for the optimization of airfoils, wings, wing-bodies, and complex aircraft configurations using both the potential equation and the Euler equations. In our most recent paper, the Euler method was extended to treat complete aircraft configurations via a new multiblock implementation. Furthermore, during the same conference, we also presented preliminary results demonstrating that this basic methodology could be ported to distributed memory parallel computing architectures. In this paper, our concern will be to demonstrate that the combined power of these new technologies can be used routinely in an industrial design environment by applying it to the case study of the design of typical supersonic transport configurations. A particular difficulty of this test case is posed by the propulsion/airframe integration.
This work describes the application of a control theory-based aerodynamic shape optimization method to the problem of supersonic aircraft design. The design process is greatly accelerated through the use of both control theory and a parallel implementation on distributed memory computers. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods (13, 12, 44, 38). The resulting problem is then implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) Standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on higher order computational fluid dynamics methods (CFD). In our earlier studies, the serial implementation of this design method (19, 20, 21, 23, 39, 25, 40, 41, 42, 43, 9) was shown to be effective for the optimization of airfoils, wings, wing-bodies, and complex aircraft configurations using both the potential equation and the Euler equations (39, 25). In our most recent paper, the Euler method was extended to treat complete aircraft configurations via a new multiblock implementation. Furthermore, during the same conference, we also presented preliminary results demonstrating that the basic methodology could be ported to distributed memory parallel computing architectures [241. In this paper, our concem will be to demonstrate that the combined power of these new technologies can be used routinely in an industrial design environment by applying it to the case study of the design of typical supersonic transport configurations. A particular difficulty of this test case is posed by the propulsion/airframe integration.
Flow, heat and mass transfer in a two-layer liquid system is investigated when the layers are subjected to a static magnetic field. The physical system consists of two immiscible liquids in a layer configuration with a free surface and an interfacial layer simulating two crystal growth configurations, namely, floating zone and horizontal zone melting. Mass transfer effects due to crystallization of one of the layers is also considered. The resulting flows are categorized into three types depending on: (a) the dominance of the free surface effect top layer, (b) a balance between the thermocapillary effect on the free surface and the interface effect between the two layers, or (c) the dominance of interfacial effects, in determining the system thermo-fluid characteristics. The use of an additional liquid layer atop the crystallizing layer leads to a reduction of the radial dopant inhomogeneity in comparison to the one-layer case. The effect an externally imposed magnetic field on dopant distribution in this two-layer configuration is rather ambiguous and depends on its direction and intensity resulting in the necessity of optimizing the applied magnetic field for a specific situation. In case of electrically non-conducting melts, the calculations show that some measure of flow control can be achieved by using an electrically conducting encapsulant layer of the desired thickness.
Many thermoplastic materials in common use for a wide range of applications, including spacecraft, develop bubbles internally as they burn due to chemical reactions taking place within the bulk. These bubbles grow and migrate until they burst at the surface, forceably ejecting volatile gases and, occasionally, molten fuel. In experiments in normal gravity, Kashiwagi and Ohlemiller observed vapor jets extending a few centimeters from the surface of a radiatively heated polymethylmethacrylate (PMMA) sample, with some molten material ejected into the gas phase. These physical phenomena complicated the combustion process considerably. In addition to the non-steady release of volatiles, the depth of the surface layer affected by oxygen was increased, attributed to the roughening of the surface by bursting events. The ejection of burning droplets in random directions presents a potential fire hazard unique to microgravity. In microgravity combustion experiments on nylon Velcro fasteners and on polyethylene wire insulation, the presence of bursting fuel vapor bubbles was associated with the ejection of small particles of molten fuel as well as pulsations of the flame. For the nylon fasteners, particle velocities were higher than 30 cm/sec. The droplets burned robustly until all fuel was consumed, demonstrating the potential for the spread of fire in random directions over an extended distance. The sequence of events for a bursting bubble has been photographed by Newitt et al.. As the bubble reaches the fluid surface, the outer surface forms a dome while the internal bubble pressure maintains a depression at the inner interface. Liquid drains from the dome until it breaks into a cloud of droplets on the order of a few microns in size. The bubble gases are released rapidly, generating vortices in the quiescent surroundings and transporting the tiny droplets. The depression left by the escaping gases collapses into a central jet, which rises with a high velocity and may break up, releasing one or more relatively large drops (on the order of a millimeter in these experiments). A better understanding of bubble development and bursting processes, the effects of bursting behavior on burning rate of the bulk material, and the circumstances under which large droplets are expelled, as well as their trajectories, sizes, and burning rates, is sought through computer modeling compared with experiment.
The KSC Autonomous Test Engineer (KATE) program has a long history at KSC. Now a part of the Autonomous Cryogenic Load Operations (ACLO) mission, this software system has been sporadically developed over the past 20+ years. Originally designed to provide health and status monitoring for a simple water-based fluid system, it was proven to be a capable autonomous test engineer for determining sources of failure in. the system, As part.of a new goal to provide this same anomaly-detection capability for a complicated cryogenic fluid system, software engineers, physicists, interns and KATE experts are working to upgrade the software capabilities and graphical user interface. Much progress was made during this effort to improve KATE. A display ofthe entire cryogenic system's graph, with nodes for components and edges for their connections, was added to the KATE software. A searching functionality was added to the new graph display, so that users could easily center their screen on specific components. The GUI was also modified so that it displayed information relevant to the new project goals. In addition, work began on adding new pneumatic and electronic subsystems into the KATE knowledgebase, so that it could provide health and status monitoring for those systems. Finally, many fixes for bugs, memory leaks, and memory errors were implemented and the system was moved into a state in which it could be presented to stakeholders. Overall, the KATE system was improved and necessary additional features were added so that a presentation of the program and its functionality in the next few months would be a success.
The Knowledge-Based Autonomous Test Engineer (KATE) program has a long history at KSC. Now a part of the Autonomous Cryogenic Load Operations (ACLO) mission, this software system has been sporadically developed over the past 20 years. Originally designed to provide health and status monitoring for a simple water-based fluid system, it was proven to be a capable autonomous test engineer for determining sources of failure in the system. As part of a new goal to provide this same anomaly-detection capability for a complicated cryogenic fluid system, software engineers, physicists, interns and KATE experts are working to upgrade the software capabilities and graphical user interface. Much progress was made during this effort to improve KATE. A display of the entire cryogenic system's graph, with nodes for components and edges for their connections, was added to the KATE software. A searching functionality was added to the new graph display, so that users could easily center their screen on specific components. The GUI was also modified so that it displayed information relevant to the new project goals. In addition, work began on adding new pneumatic and electronic subsystems into the KATE knowledge base, so that it could provide health and status monitoring for those systems. Finally, many fixes for bugs, memory leaks, and memory errors were implemented and the system was moved into a state in which it could be presented to stakeholders. Overall, the KATE system was improved and necessary additional features were added so that a presentation of the program and its functionality in the next few months would be a success.
The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.
The implementation of a general three-dimensional hydroelastic capability in the NASTRAN computer program is described. A method for analyzing the combined mode shapes of arbitrary fluid and structure finite element models was developed. The fluid was modeled with three-dimensional solid elements with options for tetrahedron, wedge, and hexahedron shapes. The elements were then connected to fluid grid points which defined the pressure in the fluid at the specified location. Matrices were formulated to provide efficient solutions for large-order problems. The structure matrices were processed separately and could be reduced using matrix condensation procedures or through a modal formulation using the normal modes of the empty structure as solution coordinates. The fluid/structure interface and the free surface were defined by the user with special NASTRAN boundary elements. A special purpose mesh generator program was provided to generate the actual NASTRAN data cards for the fluid, the structure, and the boundary elements for typical tank-type models.
A method has been developed which integrates a fluid flow analyzer and a thermal analyzer to produce both steady state and transient results of 1-D, 2-D, and 3-D analysis models. The Generalized Fluid System Simulation Program (GFSSP) is a one dimensional, general purpose fluid analysis code which computes pressures and flow distributions in complex fluid networks. The MSC Systems Improved Numerical Differencing Analyzer (MSC.SINDA) is a one dimensional general purpose thermal analyzer that solves network representations of thermal systems. Both GFSSP and MSC.SINDA have graphical user interfaces which are used to build the respective model and prepare it for analysis. The SINDA/GFSSP Conjugate Integrator (SGCI) is a formbase graphical integration program used to set input parameters for the conjugate analyses and run the models. The contents of this paper describes SGCI and its thermo-fluids conjugate analysis techniques and capabilities by presenting results from some example models including the cryogenic chill down of a copper pipe, a bar between two walls in a fluid stream, and a solid plate creating a phase change in a flowing fluid.
Cryogenic tank chill and fill is an important cryogenic fluid management (CFM) technology that supports and enables many of NASA’s long-duration space missions. For no-vent fill, the receiver tank pressure remains below the supply tank pressure during the entire duration of the fill so that the tank does not require venting. This is especially advantageous for tank fill operations in low gravity where the position of the liquid is not always known and venting the tank may cause loss of propellant by venting liquid. In lieu of expensive tests conducted on-orbit, accurate computational models capable of predicting receiver tank pressure during cryogenic propellant tank fill may be used to reduce system and propellant mass as well as mission risk. However, these numerical models must be validated or anchored to test data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict a liquid hydrogen tank chill and no-vent fill ground test conducted at Lewis Research Center (now Glenn Research Center) in 1991. The specific test case chosen for model validation implemented an upward-facing jet near the bottom of a cylindrical 34 liter tank. CFD predictions are compared to experimental measurements of tank pressure, fill level, fluid temperatures, and wall temperatures. The CFD results show reasonable agreement to the test data but overpredict the pressure collapse near the end of the fill despite the liquid jet penetrating the liquid-vapor interface. Several sensitivity studies are considered due to notable uncertainties in the experiment.
During my Spring 2016 internship, my two major contributions were the design of the Commodities Exchange Pallet and the design of a test stand for a 100 pounds-thrust sounding rocket. The Commodities Exchange Pallet is a prototype developed for the Human Exploration Spacecraft Testbed for Integration and Advancement (HESTIA) program. Under the HESTIA initiative the Commodities Exchange Pallet was developed as a method for demonstrating multi-system integration thru the transportation of In-Situ Resource Utilization produced oxygen and water to a human habitat. Ultimately, this prototype's performance will allow for future evaluation of integration, which may lead to the development of a flight capable pallet for future deep-space exploration missions. For HESTIA, my main task was to design the Commodities Exchange Pallet system to be used for completing an integration demonstration. Under the guidance of my mentor, I designed, both, the structural frame and fluid delivery system for the commodities pallet. The fluid delivery system includes a liquid-oxygen to gaseous-oxygen system, a water delivery system, and a carbon-dioxide compressors system. The structural frame is designed to meet safety and transportation requirements, as well as the ability to interface with the ER division's Portable Utility Pallet. The commodities pallet structure also includes independent instrumentation oxygen/water panels for operation and system monitoring. My major accomplishments for the commodities exchange pallet were the completion of the fluid delivery systems and the structural frame designs. In addition, parts selection was completed in order to expedite construction of the prototype, scheduled to begin in May of 2016. Once the commodities pallet is assembled and tested it is expected to complete a fully integrated transfer demonstration with the ISRU unit and the Environmental Control and Life Support System test chamber in September of 2016. In addition to the development of the Commodities Exchange Pallet, I also assisted in preparation for testing the upper stage of a sounding rocket developed as a Center Innovation Fund project. The main objective of this project is to demonstrate the integration between a propulsion system and a solid oxide fuel cell (SOFC). The upper stage and SOFC are scheduled to complete an integrated test in August of 2016. As part of preparation for scheduled testing, I was responsible for designing the upper stage's test stand/support structure and main engine plume deflector to be used during hot-fire testing (fig. 3). The structural components of the test stand need to meet safety requirements for operation of the propulsion system, which consist of a 100 pounds-thrust main engine and two 15 pounds-thrust reaction control thrusters. My main accomplishment for this project was the completion of the design and the parts selection for construction of the structure, scheduled to begin late April of 2016.
NASA's effort for the thermal environmental control of the Space Station Freedom is directed towards the design, analysis, and development of an Active Thermal Control System (ATCS). A two phase, flow through condenser/radiator concept was baselined, as a part of the ATCS, for the radiation of space station thermal load into space. The proposed condenser rejects heat through direct condensation of ATCS working fluid (ammonia) in the small diameter radiator tubes. Analysis of the condensation process and design of condenser tubes are based on the available two phase flow models for the prediction of flow regimes, heat transfer, and pressure drops. The prediction formulas use the existing empirical relationships of friction factor at gas-liquid interface. An attempt is made to study the stability of interfacial waves in two phase annular flow. The formulation is presented of a stability problem in cylindrical coordinates. The contribution of fluid viscosity, surface tension, and transverse radius of curvature to the interfacial surface is included. A solution is obtained for Kelvin-Helmholtz instability problem which can be used to determine the critical and most dangerous wavelengths for interfacial waves.
FMG3D (full multigrid 3 dimensions) is a pilot computer program that solves equations of fluid flow using a finite difference representation on a structured grid. Infrastructure exists for three dimensions but the current implementation treats only two dimensions. Written in Fortran 90, FMG3D takes advantage of the recursive subroutine feature, dynamic memory allocation, and structured-programming constructs of that language. FMG3D supports multi-block grids with three types of block-to-block interfaces: periodic, C-zero, and C-infinity. For all three types, grid points must match at interfaces. For periodic and C-infinity types, derivatives of grid metrics must be continuous at interfaces. The available equation sets are as follows: scalar elliptic equations, scalar convection equations, and the pressure-Poisson formulation of the Navier-Stokes equations for an incompressible fluid. All the equation sets are implemented with nonzero forcing functions to enable the use of user-specified solutions to assist in verification and validation. The equations are solved with a full multigrid scheme using a full approximation scheme to converge the solution on each succeeding grid level. Restriction to the next coarser mesh uses direct injection for variables and full weighting for residual quantities; prolongation of the coarse grid correction from the coarse mesh to the fine mesh uses bilinear interpolation; and prolongation of the coarse grid solution uses bicubic interpolation.
A low noise electrode suited for sensing electrocardiograms when chronically and subcutaneously implanted in a free ranging subject is described. The electrode comprises a pocket shaped electrically conductive member with a single entrance adapted to receive body fluids. The exterior of the member and the entrance region is coated with electrical insulation so that the only electrolyte/electrode interface is within the member, remote from artifact-generating tissue. Cloth straps are bonded to the member to permit the electrode to be sutured to tissue and to provide electrical lead flexure relief.