SEARCH · Search NASA
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.
The focus of research in the computational fluid dynamics (CFD) area is two fold: (1) to develop new approaches for turbulence modeling so that high speed compressible flows can be studied for applications to entry and re-entry flows; and (2) to perform research to improve CFD algorithm accuracy and efficiency for high speed flows. Research activities, faculty and student participation, publications, and financial information are outlined.
Results are presented from a simulation modeling the formation of a group of galaxies in a 'standard' cold, dark matter universe with delta = 1, h sub 0 = 50 km/(s(Mpc)), baryon fraction omega sub b = 0.1 and spectrum normalization sigma sub 8 = 0.6 (bias parameter b = 1.7). Initial conditions are generated within a periodic box with comoving length 16 Mpc in a manner constrained to produce a small cluster of total mass approximately 10 exp 14 solar mass. Two sets of 643 particles are used to model the dark matter and baryon fluids. Each gas particle represents 1.08 x 10 exp -8 solar mass, implying an L* galaxy is resolved by approximately 1000 particles. The system is evolved self-consistently in three dimensions using the combined N-body/hydrodynamic scheme P3MSPH up to a final redshift z = 1. Evolving to the present is prohibited by the fact that the mean density in the simulated volume is above critical and the entire volume would be going nonlinear beyond this point, We are currently analyzing another run with somewhat poorer mass resolution which was evolved to the present.
The Variable Thrust Engine (VTE) of the Orbital Maneuvering Vehicle (OMV) uses a hypergolic propellant combination of Monomethyl Hydrazine (MMH) and Nitrogen Tetroxide (NTO) as fuel and oxidizer, respectively. The performance of the VTE depends on a number of complex interacting phenomena such as atomization, spray dynamics, vaporization, turbulent mixing, convective/radiative heat transfer, and hypergolic combustion. This study involved the development of a comprehensive numerical methodology to facilitate detailed analysis of the VTE. An existing Computational Fluid Dynamics (CFD) code was extensively modified to include the following models: a two-liquid, two-phase Eulerian-Lagrangian spray model; a chemical equilibrium model; and a discrete ordinate radiation heat transfer model. The modified code was used to conduct a series of simulations to assess the effects of various physical phenomena and boundary conditions on the VTE performance. The details of the models and the results of the simulations are presented.
A computational and experimental study was conducted to investigate the effects of multiple injection ports in a two-dimensional, convergent-divergent nozzle, for fluidic thrust vectoring. The concept of multiple injection ports was conceived to enhance the thrust vectoring capability of a convergent-divergent nozzle over that of a single injection port without increasing the secondary mass flow rate requirements. The experimental study was conducted at static conditions in the Jet Exit Test Facility of the 16-Foot Transonic Tunnel Complex at NASA Langley Research Center. Internal nozzle performance was obtained at nozzle pressure ratios up to 10 with secondary nozzle pressure ratios up to 1 for five configurations. The computational study was conducted using the Reynolds Averaged Navier-Stokes computational fluid dynamics code PAB3D with two-equation turbulence closure and linear Reynolds stress modeling. Internal nozzle performance was predicted for nozzle pressure ratios up to 10 with a secondary nozzle pressure ratio of 0.7 for two configurations. Results from the experimental study indicate a benefit to multiple injection ports in a convergent-divergent nozzle. In general, increasing the number of injection ports from one to two increased the pitch thrust vectoring capability without any thrust performance penalties at nozzle pressure ratios less than 4 with high secondary pressure ratios. Results from the computational study are in excellent agreement with experimental results and validates PAB3D as a tool for predicting internal nozzle performance of a two dimensional, convergent-divergent nozzle with multiple injection ports.
A computational and experimental study was conducted to investigate the effects of multiple injection ports in a two-dimensional, convergent-divergent nozzle, for fluidic thrust vectoring. The concept of multiple injection ports was conceived to enhance the thrust vectoring capability of a convergent-divergent nozzle over that of a single injection port without increasing the secondary mass flow rate requirements. The experimental study was conducted at static conditions in the Jet Exit Test Facility of the 16-Foot Transonic Tunnel Complex at NASA Langley Research Center. Internal nozzle performance was obtained at nozzle pressure ratios up to 10 with secondary nozzle pressure ratios up to 1 for five configurations. The computational study was conducted using the Reynolds Averaged Navier-Stokes computational fluid dynamics code PAB3D with two-equation turbulence closure and linear Reynolds stress modeling. Internal nozzle performance was predicted for nozzle pressure ratios up to 10 with a secondary nozzle pressure ratio of 0.7 for two configurations. Results from the experimental study indicate a benefit to multiple injection ports in a convergent-divergent nozzle. In general, increasing the number of injection ports from one to two increased the pitch thrust vectoring capability without any thrust performance penalties at nozzle pressure ratios less than 4 with high secondary pressure ratios. Results from the computational study are in excellent agreement with experimental results and validates PAB3D as a tool for predicting internal nozzle performance of a two dimensional, convergent-divergent nozzle with multiple injection ports.
A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.
A two-stage procedure has been devised for estimating the spatial resolution achievable in the simulation of a given flow on a given computational grid by a computational fluid dynamics (CFD) code that incorporates a hybrid model of turbulence. The hybrid models to which this procedure is especially relevant are those of the Reynolds-averaged Navier-Stokes (RANS) and the partial-averaged Navier-Stokes (PANS) approaches. This procedure represents the first step toward adding variable-resolution turbulence-modeling capabilities to CFD codes as part of a continuing effort to increase the accuracy and robustness of CFD simulations of unsteady flows. Some background information is prerequisite to a meaningful summary of the procedure. Among experts in CFD, it is well known that combination of the Reynolds-averaged Navier-Stokes (RANS) approach and eddy-viscosity turbulence models offers limited capability for simulating unsteady and complex flows. The RANS approach includes an assumption that most of the energy in a given flow is modeled through turbulence-transport equations and is resolved in a computational grid used to simulate the flow. RANS also overpredicts eddy viscosity, thereby yielding excessive damping of unsteady motion. The eddy viscosity attains an unphysically large value because of unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved flow field. One approach used to overcome this deficiency is to provide a mechanism for the RANS equations to resolve motion only on the largest scales and to use a hybrid model to represent effects at smaller scales. The RANS approach involves the use of a standard two-equation turbulence model in which the effect of turbulence is summarized by a viscosity that is a function of (1) the time-averaged kinetic- energy density (k) associated with the local fluctuating (turbulent) component of flow and (2) the time-averaged rate of dissipation of the turbulent-kinetic- energy density ( ). In PANS, which was developed to overcome the grid dependency associated with other hybrid turbulence models (including that of RANS), the standard two-equation turbulence model is replaced by another two-equation model in which one solves for the previously unresolved k and , which are now allowed to vary in space and time.
A semi-empirical theoretical model and a C++ computer program that implements the model have been developed for use in predicting the noise generated by a rocket or jet engine. The computer program, entitled the Realtime Rocket and Jet Engine Noise Analysis and Prediction Software, is one of two main subsystems of the Acoustic Prediction/Measurement Tool, which comprises software, acoustic instrumentation, and electronic hardware combined to afford integrated capabilities for real-time prediction and measurement of noise emitted by rocket and jet engines. [The other main subsystem, consisting largely of acoustic instrumentation and electronic hardware, is described in Wireless Acoustic Measurement System, which appears elsewhere in this section.] The theoretical model was derived from the fundamental laws of fluid mechanics, as first was done by M. J. Lighthill in his now famous theory of aerodynamically generated sound. The far-field approximation of the Lighthill theory is incorporated into this model. Many other contributions from various researchers have also been introduced into the model. The model accounts for two noise components: shear noise and self noise. The final result of the model is expressed in terms of a volume integral of the acoustic intensities attributable to these two components, subject to various directivity coefficients. The computer program was written to solve the volume integral. The inputs required by the program are two data files from a computational fluid dynamics (CFD) simulation of the flow of interest: the computational-grid file and the solution file. The CFD solution should be one that has been obtained for conditions that closely approximate those of an experimental test that is yet to be performed. In the current state of development of the model and software, it is recommended that the observation points lie along a radius at an angle >60 from the jet axis. The software provides, and is driven via, a graphical user interface, which facilitates its use. Optionally, the program accepts additional input in the form of data on the measured sound pressure level as a function of frequency at a given far-field location, preferably at an angle of 90 from the jet axis. The user is prompted to use default empirical constants or to choose constants based the measurement data. The user can view the results and compare them with other computational or experimental data. Once satisfied with the results, the user can save a graph of the results in a file that can be imported into documents.
Full seven-equation Reynolds stress turbulence models are a relatively new and promising tool for todays aerospace technology challenges. This paper uses two stress-omega full Reynolds stress models to evaluate challenging flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSGLRR full second-moment Reynolds stress models are evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST) two-equation models.
Demanding aerodynamic modeling requirements for military and civilian aircraft have provided impetus for researchers to improve computational and experimental techniques. Model validation is a key component for these research endeavors so this study is an initial effort to extend conventional time history comparisons by comparing model parameter estimates and their standard errors using system identification methods. An aerodynamic model of an aircraft performing one-degree-of-freedom roll oscillatory motion about its body axes is developed. The model includes linear aerodynamics and deficiency function parameters characterizing an unsteady effect. For estimation of unknown parameters two techniques, harmonic analysis and two-step linear regression, were applied to roll-oscillatory wind tunnel data and to computational fluid dynamics (CFD) simulated data. The model used for this study is a highly swept wing unmanned aerial combat vehicle. Differences in response prediction, parameters estimates, and standard errors are compared and discussed
Thermal modeling of a cryogenic pump and a hot-gas turbine in a turbopump assembly proposed for the Space Shuttle is described in this paper. A model, developed by identifying the heat-transfer regimes and incorporating their dependencies into a turbopump system model, included heat transfer for two-phase cryogen, hot-gas (200 R) impingement on turbine blades, gas impingement on rotating disks and parallel plate fluid flow. The ?thermal analyzer' program employed to develop this model was the TRW Systems Improved Numerical Differencing Analyzer (SINDA). This program uses finite differencing with lumped parameter representation for each node. Also discussed are model development, simulations of turbopump startup/shutdown operations, and the effects of varying turbopump parameters on the thermal performance.
The Zero-Boil-Off Tank (ZBOT) experiments are a series of small scale two-phase fluids physics experiments conducted aboard the International Space Station (ISS) in 2017. A goal of these experiments is to study various fundamental fluid flow, heat and mass transport, and phase change phenomena that control storage tank pressurization and pressure control in microgravity. In this paper, results of a state-of-the-art two-phase CFD storage tank model are presented and validated against jet mixing data provided by microgravity ZBOT experiments. Jet-ullage interaction was studied in the experiment for different jet speeds representing different flow regimes, and tank fill levels. Effects of various physical parameters, such as, jet angle, gravitational acceleration and surface tension value were investigated. Computational parameters, such as turbulence and surface tension models and interface capturing scheme, affecting jet-ullage interaction were also studied. Computational results were compared to experimental results for model validation.
The fluid mechanics associated with crystal growth processes on centrifuges is modeled using 2D and 3D models. Two-dimensional calculations show that flow bifurcations exist in such crystal growth configurations where the ampoule is oriented in the same direction as the resultant gravity vector and a temperature gradient is imposed on the melt. A scaling analysis is formulated to predict the flow transition point from the natural convection dominated regime to the Coriolis force dominated regime. Results of 3D calculations are presented for two thermal configurations of the crystal growth cell: top heated and bottom heated with respect to the centrifugal acceleration. In the top heated configuration, a substantial reduction in the convection intensity within the melt can be attained by centrifuge operations, and close to steady diffusion-limited thermal conditions can be achieved over a narrow range of the imposed microgravity level. In the bottom heated configuration the Coriolis force has a stabilizing effect on fluid motion by delaying the onset of unsteady convection.
A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.
In view of experimental considerations, we give a model-independent argument that the novel tricritical points in multicomponent fluid mixtures, where three phases simultaneously become critical, are points on the boundary of a single two-dimensional surface of critical points. This result is corroborated by the Landau model suggested by Griffiths. The relationship between these tricritical points and the complex 'higher-order' critical points proposed to exist in certain magnetic systems is elucidated.
The acoustic slot elements, CSLOTi, are applied to analyze two-dimensional enclosures with fixed or moving boundaries. The capability utilized to compute (a) the acoustic natural modes and frequencies of a rigid walled enclosure and (b) the sound pressure at any point inside an enclosure when the surrounding walls are forced to vibrate. Applications to an automobile passenger compartment illustrate the technique. The axisymmetric fluid elements, CFLUIDi, are used in conjunction with a suitable choice of symmetry planes and a model of the surrounding structure to approximate a two-dimensional enclosure with flexible walls. The enclosure walls are modeled using finite elements or structural modes. Illustrative examples include a comparison of rectangular cavity modes with those calculated using the acoustic slot element and the free vibration modes of two enclosures coupled through a flexible rectangular panel.
Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.
An inviscid model of the surface dynamics at the interface between two parallel-flowing insulating dielectric fluids, subject to a tangential electric field with a weak gradient, is developed. The two flows are oppositely directed, resulting in an electrohydrodynamically coupled Kelvin-Hemholtz instability. A simple stability criterion is used successfully to correlate experimental data obtained with a test configuration.