SEARCH · Engineering Papers

Results for “Numerical Propulsion System Simulation”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS2023

CFD Modeling of Bi-Directional PMD inside Cryogenic Propellant Tanks Onboard Parabolic Flights

Future cryogenic propulsion systems will require efficient methods with which to transfer cryogenic propellants from a depot storage tank to a customer receiver tank to minimize cost and maximize reusability. The Reduced Gravity Cryogenic Transfer project is currently developing advanced cryogenic fluid management technology and developing and validating new numerical models for three phases of transfer: line chilldown, tank chilldown, and tank fill. Additionally, multiple liquid nitrogen (LN 2 ) parabolic flight transfer rigs are being designed by universities and NASA to investigate the gravitational sensitivities that exist in these three technologies. In order to maximize the collection of low-g data during flights, it is required to extract as much (LN 2 as possible from the supply tank, despite variable gravity levels. The purpose of this paper is to present computational fluid dynamics (CFD) volume of fluid simulations of (LN 2 behavior in the supply tank onboard parabolic flights to validate the optimal design of a bi-directional propellant management device (PMD) using the commercial software FLOW-3D. A parametric study is conducted on the effects of gravity level, fill level, pore size, open area, thickness, and type of baffle on PMD performance. Based on results, the PMD as designed exceeds the targeted expulsion efficiency.

Jason Hartwig↗

NASA NTRS1993

Research activities at the Center for Modeling of Turbulence and Transition

The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.

Shih, Tsan-Hsing↗

NASA NTRS2018

Computational Investigation of the Near-Field Plasma Plume in Ion-Ion Propulsion

A two-fluid numerical model of plasma flows was developed to investigate the plume of an ion-ion propulsion system. The densities of positive and negative ions, and the associated values of net charge, electric field, and electric potential were calculated as a function of time throughout the domain. The computational domain was chosen to be large enough (25 thruster diameters downstream of the exit plane) to allow for examining the neutralization of the plume. The resulting plasma electric potential and charge neutrality at the downstream end of the domain are shown and they indicate that it is possible to alternatively accelerate oppositely charged ions without the need for an electron-emitting neutralizer and without facing any electric potential hills that could cause stagnation. However, compared to existing literature on ion-ion plasma thrusters, the results from this simulation predict a longer length-scale for voltage decay.

Sankaran, K.↗

NASA NTRS2007

Emergency Flight Control of a Twin-Jet Commercial Aircraft using Manual Throttle Manipulation

The Department of Homeland Security (DHS) created the PCAR (Propulsion-Controlled Aircraft Recovery) project in 2005 to mitigate the ManPADS (man-portable air defense systems) threat to the commercial aircraft fleet with near-term, low-cost proven technology. Such an attack could potentially cause a major FCS (flight control system) malfunction or other critical system failure onboard the aircraft, despite the extreme reliability of current systems. For the situations in which nominal flight controls are lost or degraded, engine thrust may be the only remaining means for emergency flight control [ref 1]. A computer-controlled thrust system, known as propulsion-controlled aircraft (PCA), was developed in the mid 1990s with NASA, McDonnell Douglas and Honeywell. PCA's major accomplishment was a demonstration of an automatic landing capability using only engine thrust [ref 11. Despite these promising results, no production aircraft have been equipped with a PCA system, due primarily to the modifications required for implementation. A minimally invasive option is TOC (throttles-only control), which uses the same control principles as PCA, but requires absolutely no hardware, software or other aircraft modifications. TOC is pure piloting technique, and has historically been utilized several times by flight crews, both military and civilian, in emergency situations stemming from a loss of conventional control. Since the 1990s, engineers at NASA Dryden Flight Research Center (DFRC) have studied TOC, in both simulation and flight, for emergency flight control with test pilots in numerous configurations. In general, it was shown that TOC was effective on certain aircraft for making a survivable landing. DHS sponsored both NASA Dryden Flight Research Center (Edwards, CA) and United Airlines (Denver, Colorado) to conduct a flight and simulation study of the TOC characteristics of a twin-jet commercial transport, and assess the ability of a crew to control an aircraft down to a survivable runway landing using TOC. The PCAR project objective was a set of pilot procedures for operation of a specific aircraft without hydraulics that (a) have been validated in both simulation and flight by relevant personnel, and (b) mesh well with existing commercial operations, maintenance, and training at a minimum cost. As a result of this study, a procedure has been developed to assist a crew in making a survivable landing using TOC. In a simulation environment, line pilots with little or no previous TOC experience performed survivable runway landings after a few practice TOC approaches. In-flight evaluations put line pilots in a simulated emergency situation where TOC was used to recover the aircraft, maneuver to a landing site, and perform an approach down to 200 feet AGL. The results of this research, including pilot observations, procedure comments, recommendations, future work and lessons learned, will he discussed. Flight data and video footage of TOC approaches may also be shown.

Cole, Jennifer H.↗

NASA NTRS2021

High Contrast Imaging With MEMS Deformable Mirrors in the Decadal Survey Testbed

Deformable mirrors (DMs) are an essential part of any coronagraphic, high contrast instrument. They mitigate optical aberrations in the system and can even be used to generate contrast for the coronagraph. MEMS DMs from Boston Micromachines have been selected as the baseline for two flagship space telescopes proposed to the 2020 Decadal Survey. Although MEMS DMs have over a decade of heritage on ground-based telescopes and in in-air testbeds around the globe, they have not been tested in vacuum down to the ∼10−10 contrast level needed to image terrestrial exoplanets. In this paper, we describe vacuum tests of MEMS DMs in the Decadal Survey Testbed at the Jet Propulsion Laboratory. The first challenge was a bright, temporally incoherent signal, which was identified as electronics noise and removed with a low-pass filter. After that, the contrast has been limited in broadband light by the strong print-through on the DM surfaces. We performed numerical simulations to confirm that conclusion and to characterize the improvements needed to the MEMS DM surfaces and the testbed layout to attain our goal of 10−10 contrast.

Seo, Byoung-Joon↗

NASA NTRS2015

Investigation of Unsteady Flow Interaction Between an Ultra-Compact Inlet and a Transonic Fan

In the study presented, unsteady flow interaction between an ultra-compact inlet and a transonic fan stage is investigated. Future combat aircraft engines require ultra-compact inlet ducts as part of an integrated, advanced propulsion system to improve air vehicle capability and effectiveness to meet future mission needs. The main purpose of the current study is to advance the understanding of the flow interaction between a modern ultra-compact inlet and a transonic fan for future design applications. Many experimental/ analytical studies have been reported on the aerodynamics of compact inlets in aircraft engines. On the other hand, very few studies have been reported on the effects of flow distortion from these inlets on the performance of the following fan/compressor stages. The primary goal of the study presented is to investigate how flow interaction between an ultra-compact inlet and a transonic compressor influence the operating margin of the compressor. Both Unsteady Reynolds-averaged Navier-Stokes (URANS) and Large Eddy Simulation (LES) approaches are used to calculate the unsteady flow field, and the numerical results are used to study the flow interaction. The present study indicates that stall inception of the following compressor stage is affected directly based on how the distortion pattern evolves before it interacts with the fan/compressor face. For the present compressor, the stall initiates at the tip section with clean inlet flow and distortion pattern away from the casing itself seems to have limited impacts on the stall inception of the compressor. A counter-rotating swirl, which is generated due to flow separation inside the s-shaped compact duct, generates an increased flow angle near the blade tip. This increased flow angle near the rotor tip due to the secondary flow from the counter-rotating vortices is the primary reason for the reduced compressor stall margin.

Compressor↗

NASA NTRSDate not supplied

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown↗

NASA NTRS1995

Progress in incompressible Navier-Stokes computations for propulsion flows and its dual-use applications

Development of an incompressible Navier-Stokes solution procedure was performed for the analysis of a liquid rocket engine pump components and for the mechanical heart assist devices. The solution procedure for the propulsion systems is applicable to incompressible Navier-Stokes flows in a steadily rotating frame of reference for any general complex configurations. The computer codes were tested on different complex configurations such as liquid rocket engine inducer and impellers. As a spin-off technology from the turbopump component simulations, the flow analysis for an axial heart pump was conducted. The baseline Left Ventricular Assist Device (LVAD) design was improved by adding an inducer geometry by adapting from the liquid rocket engine pump. The time-accurate mode of the incompressible Navier-Stokes code was validated with flapping foil experiment by using different domain decomposition methods. In the flapping foil experiment, two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions for a downstream larger stationary foil. Fairly good agreement was obtained between unsteady experimental data and numerical results from two different moving boundary procedures. Incompressible Navier-Stokes code (INS3D) has been extended for heat transfer applications. The temperature equation was written for both forced and natural convection phenomena. Flow in a square duct case was used for the validation of the code in both natural and forced convection.

Kiris, Cetin↗

NASA NTRS2017

Conceptual Feasibility Study of the Hyperloop Vehicle for Next-Generation Transport

The Hyperloop concept is proposed as a faster, cheaper alternative to high-speed rail and traditional short-haul aircraft. It consists of a passenger pod traveling through a tube under light vacuum while being propelled and levitated by a combination of permanent and electro-magnets. The concept addresses NASA's research thrusts for growth in demand, sustainability, and technology convergence for high-speed transport. Hyperloop is a radical departure from other advanced aviation concepts, however it remains an aeronautics concept that tackles the same strategic goals of low-carbon propulsion and ultra-effcient vehicles. System feasibility was investigated by building a multidisciplinary vehicle sizing model that takes into account aerodynamic, thermodynamic, structures, electromagnetic, weight, and mission analyses. The sizing process emphasized the strong coupling between the two largest systems: the tube and the passenger pod. The model was then exercised to examine Hyperloop from a technical and cost perspective. The structural sizing analysis of the travel tube demonstrates potential for signi cant capital cost reductions by considering an underwater route. Examination of varying passenger capacity indicates that the system can be operated with a wide range of passenger loads without significant change in operating expenses. Lastly, a high-level sizing study simulated variations in tube area, pressure, pod speed, and passenger capacity showing that there is a tube pressure that minimizes operating energy usage. The value of this optimal tube pressure is highly sensitive to numerous design details. These combined estimates of energy consumption, passenger throughput, and mission analyses all support Hyperloop as a faster and cheaper alternative to short-haul flights. The tools and expertise used to quantify these results also demonstrate how traditional aerospace design methods can be leveraged to handle the complex and coupled design process. Much of the technology development required for the Hyperloop is shared with next-generation aircraft. Furthermore, the substantial public interest and active commercial development make it an ideal candidate as an aircraft technology driver and test bed.

Hyperloop Vehicle↗

NASA NTRS2003

Large-Eddy Simulation Code Developed for Propulsion Applications

A large-eddy simulation (LES) code was developed at the NASA Glenn Research Center to provide more accurate and detailed computational analyses of propulsion flow fields. The accuracy of current computational fluid dynamics (CFD) methods is limited primarily by their inability to properly account for the turbulent motion present in virtually all propulsion flows. Because the efficiency and performance of a propulsion system are highly dependent on the details of this turbulent motion, it is critical for CFD to accurately model it. The LES code promises to give new CFD simulations an advantage over older methods by directly computing the large turbulent eddies, to correctly predict their effect on a propulsion system. Turbulent motion is a random, unsteady process whose behavior is difficult to predict through computer simulations. Current methods are based on Reynolds-Averaged Navier- Stokes (RANS) analyses that rely on models to represent the effect of turbulence within a flow field. The quality of the results depends on the quality of the model and its applicability to the type of flow field being studied. LES promises to be more accurate because it drastically reduces the amount of modeling necessary. It is the logical step toward improving turbulent flow predictions. In LES, the large-scale dominant turbulent motion is computed directly, leaving only the less significant small turbulent scales to be modeled. As part of the prediction, the LES method generates detailed information on the turbulence itself, providing important information for other applications, such as aeroacoustics. The LES code developed at Glenn for propulsion flow fields is being used to both analyze propulsion system components and test improved LES algorithms (subgrid-scale models, filters, and numerical schemes). The code solves the compressible Favre-filtered Navier- Stokes equations using an explicit fourth-order accurate numerical scheme, it incorporates a compressible form of Smagorinsky s model for the subgrid-scale turbulence, and it uses generalized curvilinear coordinates to allow analysis of a wide range of geometries. The code runs in parallel on shared memory multiprocessor computers and is written in Fortran 90 with dynamic memory allocation. A sample result for a Mach-1.4 round jet is presented in the figure. Instantaneous Mach number contours in several cross-planes downstream of the nozzle exit are shown, illustrating how an LES captures the large unsteady three-dimensional turbulent structures present in the jet.

DeBonis, James R.↗

NASA NTRS2020

Multibody Simulation of REMORA CubeSat Docking to and Pushing a Spent Rocket Booster

This paper details the multibody simulation of three phases of critical importance to the feasibility of the REMORA CubeSat space debris mitigation concept: the final approach of the CubeSat to a spent rocket booster; the grappling of the spent rocket booster using a robotic arm; and the pushing of the spent rocket booster to divert its course from another on-orbit asset. The extension of a robotic mobility and manipulation modeling toolkit (M3TK) from multibody dynamics simulation of manipulators and ground vehicles to simulation of orbital robotics is outlined. This includes the identification of the appropriate parameters required to concisely and generically describe thruster loads, thruster mixing, spacecraft control, and spacecraft navigation for the purpose of on-orbit robotics simulation. A high-level spacecraft navigator commands maneuvers to target spacecraft states. A PID spacecraft controller takes the target states and calculates desired forces and torques. A thruster mixer solves a quadratic program to determine the optimal thruster firing times for the propulsion system. Pulse width modulated actuation of eight canted cold gas thrusters is used in the simulated approach to a rocket nozzle from a distance of 200 m. A five degree of freedom robotic arm is controlled to position a pair of pincers to grasp the rocket nozzle. Contact dynamics are used to accurately simulate the grasping of the rocket nozzle by the pincers. A similar, but separate simulation is performed to assess the ability of the REMORA CubeSat to push the large spent rocket booster. This diversion maneuver makes use of an additional, larger thruster, and pushes the rocket booster in excess of 400 m. Appropriate motor control gains on the robotic arm are found to be higher during the pushing phase than those which are appropriate during free motion; this increase promotes rigidity of the arm and allows it to properly direct the pushing force. Challenges encountered in time step selection for numerical stability of the simulation are also discussed.

Mukherjee, Rudranarayan↗

NASA NTRS2001

Filtered Mass Density Function for Design Simulation of High Speed Airbreathing Propulsion Systems

The objective of this research is to develop and implement new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. We have just completed two (2) years of Phase I of this research. This annual report provides a brief and up-to-date summary of our activities during the period: September 1, 2000 through August 31, 2001. In the work within the past year, a methodology termed "velocity-scalar filtered density function" (VSFDF) is developed and implemented for large eddy simulation (LES) of turbulent flows. In this methodology the effects of the unresolved subgrid scales (SGS) are taken into account by considering the joint probability density function (PDF) of all of the components of the velocity and scalar vectors. An exact transport equation is derived for the VSFDF in which the effects of the unresolved SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source terms appear in closed form. The remaining unclosed terms in this equation are modeled. A system of stochastic differential equations (SDEs) which yields statistically equivalent results to the modeled VSFDF transport equation is constructed. These SDEs are solved numerically by a Lagrangian Monte Carlo procedure. The consistency of the proposed SDEs and the convergence of the Monte Carlo solution are assessed by comparison with results obtained by an Eulerian LES procedure in which the corresponding transport equations for the first two SGS moments are solved. The unclosed SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source in the Eulerian LES are replaced by corresponding terms from VSFDF equation. The consistency of the results is then analyzed for a case of two dimensional mixing layer.

Drozda, T. G.↗

NASA NTRS2015

Plume Characterization of a Laboratory Model 22 N GPIM Thruster via High-Frequency Raman Spectroscopy

The Green Propellant Infusion Mission (GPIM) will demonstrate the capability of a green propulsion system, specifically, one using the monopropellant, AF-M315E. One of the risks identified for GPIM is potential contamination of sensitive areas of the spacecraft from the effluents in the plumes of AF-M315E thrusters. Plume characterization of a laboratory-model 22 N thruster via optical diagnostics was conducted at NASA GRC in a space-simulated environment. A high-frequency pulsed laser was coupled with an electron-multiplied ICCD camera to perform Raman spectroscopy in the near-field, low-pressure plume. The Raman data yielded plume constituents and temperatures over a range of thruster chamber pressures and as a function of thruster (catalyst) operating time. Schlieren images of the near-field plume enabled calculation of plume velocities and revealed general plume structure of the otherwise invisible plume. The measured velocities are compared to those predicted by a two-dimensional, kinetic model. Trends in data and numerical results are presented from catalyst mid-life to end-of-life. The results of this investigation were coupled with the Raman and Schlieren data to provide an anchor for plume impingement analysis presented in a companion paper. The results of both analyses will be used to improve understanding of the nature of AF-M315E plumes and their impacts to GPIM and other future missions.

Raman Spectroscopy↗

NASA NTRS2006

Evaluation of Open Cell Foam Heat Transfer Enhancement for Liquid Rocket Engine

As NASA pursues the exploration mission, advanced propulsion for the next generation of spacecraft will be needed. These new propulsion systems will require higher performance and increased durability, despite current limitations on materials. A break-through technology is needed in the thrust chamber. In this paper the idea of using a porous metallic foam is examined for its potential cooling enhancement capabilities. The goal is to increase the chamber wall cooling without creating an additional pressure drop penalty. A feasibility study based on experiments at laboratory-scale conditions was performed and analysis at rocket conditions is underway. In the experiment, heat transfer and pressure drop data were collected using air as the coolant in a copper or nickel foam filled annular channel. The foam-channel performance was evaluated based on comparison with conventional microchannel cooling passages under equal pressure drop conditions. The heat transfer enhancement of the foam channel over the microchannel ranges from 130% to 172%. The enhancement is relatively independent of the pressure drop and increases with decreasing pore size. A direct numerical simulation model of the foam heat exchange has been built. The model is based on the actual metal foam microstructure of thin ligaments (0.2- 0.3 mm in diameter) that form a network of interconnected open-cells. The cell dimension is around 2 mm. The numerical model was built using the FLUENT CFD code. Comparison of the pressure drop results predicted by the current model with those experimental data of Leong and Jin [8] shows favorable comparisons. Pressure drop predictions have been made using hydrogen as a coolant at typical rocket conditions. Conjugate heat transfer analysis using the foam filled channel is planned for the future.

Chung, J. N.↗

NASA NTRS1997

The Effects of Magnetic Nozzle Configurations on Plasma Thrusters

Over the course of eight years, the Ohio State University has performed research in support of electric propulsion development efforts at the NASA Lewis Research Center, Cleveland, OH. This research has been largely devoted to plasma propulsion systems including MagnetoPlasmaDynamic (MPD) thrusters with externally-applied, solenoidal magnetic fields, hollow cathodes, and Pulsed Plasma Microthrusters (PPT's). Both experimental and theoretical work has been performed, as documented in four master's theses, two doctoral dissertations, and numerous technical papers. The present document is the final report for the grant period 5 December 1987 to 31 December 1995, and summarizes all activities. Detailed discussions of each area of activity are provided in appendices: Appendix 1 - Experimental studies of magnetic nozzle effects on plasma thrusters; Appendix 2 - Numerical modeling of applied-field MPD thrusters; Appendix 3 - Theoretical and experimental studies of hollow cathodes; and Appendix 4 -Theoretical, numerical and experimental studies of pulsed plasma thrusters. Especially notable results include the efficacy of using a solenoidal magnetic field downstream of a plasma thruster to collimate the exhaust flow, the development of a new understanding of applied-field MPD thrusters (based on experimentally-validated results from state-of-the art, numerical simulation) leading to predictions of improved performance, an experimentally-validated, first-principles model for orificed, hollow-cathode behavior, and the first time-dependent, two-dimensional calculations of ablation-fed, pulsed plasma thrusters.

Turchi, P. J.↗

NASA NTRS1996

Optimum combining of residual carrier array signals in correlated noises

An array feed combining system for the recovery of signal-to-noise ratio (SNR) loss due to antenna reflector deformation has been implemented and is currently being evaluated on the Jet Propulsion Laboratory 34-m DSS-13 antenna. The current signal-combining system operates under the assumption that the white Gaussian noise processes in the received signals from different array elements are mutually uncorrelated. However, experimental data at DSS 13 indicate that these noise processes are indeed mutually correlated. The objective of this work is to develop a signal-combining system optimized to account for the mutual correlations between these noise processes. The set of optimum combining weight coefficients that maximizes the combined signal SNR in the correlated noises environment is determined. These optimum weights depend on unknown signal and noise covariance parameters. A maximum-likelihood approach is developed to estimate these unknown parameters to obtain estimates of the optimum weight coefficients based on residual carrier signal samples. The actual combined signal SNR using the estimated weight coefficients is derived and shown to converge to the maximum achievable SNR as the number of signal samples increases. These results are also verified by simulation. A numerical example shows a significant improvement in SNR performance can be obtained, especially when the amount of correlation increases.

Liang, R.↗

NASA NTRS1996

Optimum Combining of Residual Carrier Array Signals in Correlated Noises

Tan, H. H.↗

NASA NTRS2015

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

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

Majumdar, A. K.↗