Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes solver”

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.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Creating a Simple Single Computational Approach to Modeling Rarefied and Continuum Flow About Aerospace Vehicles

We proposed to create a single computational code incorporating methods that can model both rarefied and continuum flow to enable the efficient simulation of flow about space craft and high altitude hypersonic aerospace vehicles. The code was to use a single grid structure that permits a smooth transition between the continuum and rarefied portions of the flow. Developing an appropriate computational boundary between the two regions represented a major challenge. The primary approach chosen involves coupling a four-speed Lattice Boltzmann model for the continuum flow with the DSMC method in the rarefied regime. We also explored the possibility of using a standard finite difference Navier Stokes solver for the continuum flow. With the resulting code we will ultimately investigate three-dimensional plume impingement effects, a subject of critical importance to NASA and related to the work of Drs. Forrest Lumpkin, Steve Fitzgerald and Jay Le Beau at Johnson Space Center. Below is a brief background on the project and a summary of the results as of the end of the grant.

Goldstein, David B.↗
NASA NTRS1998

Parallel PAB3D: Experiences with a Prototype in MPI

PAB3D is a three-dimensional Navier Stokes solver that has gained acceptance in the research and industrial communities. It takes as computational domain, a set disjoint blocks covering the physical domain. This is the first report on the implementation of PAB3D using the Message Passing Interface (MPI), a standard for parallel processing. We discuss briefly the characteristics of tile code and define a prototype for testing. The principal data structure used for communication is derived from preprocessing "patching". We describe a simple interface (COMMSYS) for MPI communication, and some general techniques likely to be encountered when working on problems of this nature. Last, we identify levels of improvement from the current version and outline future work.

Guerinoni, Fabio↗
NASA NTRS1998

Computational Study of Surface Tension and Wall Adhesion Effects on an Oil Film Flow Underneath an Air Boundary Layer

The fringe-imaging skin friction (FISF) technique, which was originally developed by D. J. Monson and G. G. Mateer at Ames Research Center and recently extended to 3-D flows, is the most accurate skin friction measurement technique currently available. The principle of this technique is that the skin friction at a point on an aerodynamic surface can be determined by measuring the time-rate-of-change of the thickness of an oil drop placed on the surface under the influence of the external air boundary layer. Lubrication theory is used to relate the oil-patch thickness variation to shear stress. The uncertainty of FISF measurements is estimated to be as low as 4 percent, yet little is known about the effects of surface tension and wall adhesion forces on the measured results. A modified version of the free-surface Navier-Stokes solver RIPPLE, developed at Los Alamos National Laboratories, was used to compute the time development of an oil drop on a surface under a simulated air boundary layer. RIPPLE uses the volume of fluid method to track the surface and the continuum surface force approach to model surface tension and wall adhesion effects. The development of an oil drop, over a time period of approximately 4 seconds, was studied. Under the influence of shear imposed by an air boundary layer, the computed profile of the drop rapidly changes from its initial circular-arc shape to a wedge-like shape. Comparison of the time-varying oil-thickness distributions computed using RIPPLE and also computed using a greatly simplified numerical model of an oil drop equation which does not include surface tension and wall adhesion effects) was used to evaluate the effects of surface tension on FISF measurement results. The effects of surface tension were found to be small but not necessarily negligible in some cases.

Celic, Alan↗
NASA NTRS1998

Neural Net-Based Redesign of Transonic Turbines for Improved Unsteady Aerodynamic Performance

A recently developed neural net-based aerodynamic design procedure is used in the redesign of a transonic turbine stage to improve its unsteady aerodynamic performance. The redesign procedure used incorporates the advantages of both traditional response surface methodology (RSM) and neural networks by employing a strategy called parameter-based partitioning of the design space. Starting from the reference design, a sequence of response surfaces based on both neural networks and polynomial fits are constructed to traverse the design space in search of an optimal solution that exhibits improved unsteady performance. The procedure combines the power of neural networks and the economy of low-order polynomials (in terms of number of simulations required and network training requirements). A time-accurate, two-dimensional, Navier-Stokes solver is used to evaluate the various intermediate designs and provide inputs to the optimization procedure. The optimization procedure yields a modified design that improves the aerodynamic performance through small changes to the reference design geometry. The computed results demonstrate the capabilities of the neural net-based design procedure, and also show the tremendous advantages that can be gained by including high-fidelity unsteady simulations that capture the relevant flow physics in the design optimization process.

Madavan, Nateri K.↗
NASA NTRS1998

Parallel 3D Multi-Stage Simulation of a Turbofan Engine

A 3D multistage simulation of each component of a modern GE Turbofan engine has been made. An axisymmetric view of this engine is presented in the document. This includes a fan, booster rig, high pressure compressor rig, high pressure turbine rig and a low pressure turbine rig. In the near future, all components will be run in a single calculation for a solution of 49 blade rows. The simulation exploits the use of parallel computations by using two levels of parallelism. Each blade row is run in parallel and each blade row grid is decomposed into several domains and run in parallel. 20 processors are used for the 4 blade row analysis. The average passage approach developed by John Adamczyk at NASA Lewis Research Center has been further developed and parallelized. This is APNASA Version A. It is a Navier-Stokes solver using a 4-stage explicit Runge-Kutta time marching scheme with variable time steps and residual smoothing for convergence acceleration. It has an implicit K-E turbulence model which uses an ADI solver to factor the matrix. Between 50 and 100 explicit time steps are solved before a blade row body force is calculated and exchanged with the other blade rows. This outer iteration has been coined a "flip." Efforts have been made to make the solver linearly scaleable with the number of blade rows. Enough flips are run (between 50 and 200) so the solution in the entire machine is not changing. The K-E equations are generally solved every other explicit time step. One of the key requirements in the development of the parallel code was to make the parallel solution exactly (bit for bit) match the serial solution. This has helped isolate many small parallel bugs and guarantee the parallelization was done correctly. The domain decomposition is done only in the axial direction since the number of points axially is much larger than the other two directions. This code uses MPI for message passing. The parallel speed up of the solver portion (no 1/0 or body force calculation) for a grid which has 227 points axially.

Turner, Mark G.↗
NASA NTRS1999

A Numerical Study of Low-Reynolds-Number Separation Bubbles

The present study uses two dimensional numerical simulations to study unsteady low-Reynolds-number separation bubbles. The numerical study is in two parts: (1) a two dimensional time-accurate Navier-Stokes solver is used to simulate flows over the APEX airfoil, and (2) a numerical procedure is developed for localized simulations of transitional separation bubbles. The 2-D computations of flow over the APEX airfoil show that the flow is unsteady with periodic vortex shedding. A linear stability analysis of the separated flow shows that the vortex shedding is caused due to the instability of the separated flow. For transonic flows over the APEX airfoil the vortex shedding is additionally influenced by the presence of shocks. The flowfield has two characteristic time scales, one corresponding to the vortex shedding and another corresponding to the movement of the shocks. The two dimensional (2-D) airfoil simulations also showed the presence of nonlinear effects in the separated region. To better understand the characteristics of separation bubbles a numerical procedure has been developed for localized separation bubble calculations. This procedure is used to perform computations for a flat plate separation bubble test case. The separation bubble is induced by specifying a velocity gradient in the freestream. The growth of disturbances in the separation bubble is analyzed by introducing disturbances upstream of the separation bubble.

Tatineni, Mahidhar↗
NASA NTRS1999

Numerical Simulation of Complex Turbomachinery Flows

An unsteady, multiblock, Reynolds Averaged Navier Stokes solver based on Runge-Kutta scheme and Pseudo-time step for turbo-machinery applications was developed. The code was validated and assessed against analytical and experimental data. It was used to study a variety of physical mechanisms of unsteady, three-dimensional, turbulent, transitional, and cooling flows in compressors and turbines. Flow over a cylinder has been used to study effects of numerical aspects on accuracy of prediction of wake decay and transition, and to modify K-epsilon models. The following simulations have been performed: (a) Unsteady flow in a compressor cascade: Three low Reynolds number turbulence models have been assessed and data compared with Euler/boundary layer predictions. Major flow features associated with wake induced transition were predicted and studied; (b) Nozzle wake-rotor interaction in a turbine: Results compared to LDV data in design and off-design conditions, and cause and effect of unsteady flow in turbine rotors were analyzed; (c) Flow in the low-pressure turbine: Assessed capability of the code to predict transitional, attached and separated flows at a wide range of low Reynolds numbers and inlet freestream turbulence intensity. Several turbulence and transition models have been employed and comparisons made to experiments; (d) leading edge film cooling at compound angle: Comparisons were made with experiments, and the flow physics of the associated vortical structures were studied; and (e) Tip leakage flow in a turbine. The physics of the secondary flow in a rotor was studied and sources of loss identified.

Chernobrovkin, A. A.↗
NASA NTRS1996

Improved Models for the Porous Surface with Passive Control

A computational investigation of afterbody flow using a passive control method is conducted. The passive control method consists of a porous surface placed over a plenum. The purpose of the passive control method is to exploit the adverse pressure gradient present in afterbody flow in an attempt to reduce boundary layer separation and afterbody drag. Four different porous wall models are used to model the transpiration velocity in the region of passive control. A three-dimensional, time-dependent, Reynolds-averaged, simplified Navier-Stokes solver, PAB3D, is used to simulate afterbody flow with and without passive control. Three afterbody configurations with boat-tail angles of 10, 20, and 30 deg. are used to obtain two-dimensional solutions with a freestream Mach number of 0.6 and nozzle pressure ratio of 6. The region of passive control was initially placed from 20-60% of the nozzle length. The effect of the porous placement and porous extent is also studied. Baseline (no porosity) two-dimensional solutions are qualitatively similar to experimental data but under-predict the magnitude of the pressure recovery. Results for the subsonic solutions show losses in the pressure recovery for some cases with passive control. Three-dimensional effects are also investigated and seen to be very significant. Three-dimensional baseline solutions, for both sub- and super-critical freestream Mach numbers, compare very favorably with the experimental data in comparison to the two-dimensional solution. Future work is required to examine three-dimensional afterbody flows with passive porosity.

Kraushaar, Sandra L.↗
NASA NTRS1998

Numerical Simulations of Unsteady Low-Reynolds-Number Flows Over the APEX Airfoil

Laminar and transitional separation bubbles are an important feature of low-Reynolds-number flows over airfoils. The separation bubbles are unsteady and have a significant impact on the aerodynamic properties of the airfoils. In this paper unsteady low-Reynolds-number separated flows over the APEX airfoil are calculated using a Navier-Stokes solver. The numerical results show the presence of unsteady separation bubbles in the flowfield. An analysis of the numerical results shows that flowfield disturbances are amplified significantly in the separation bubble, leading to periodic vortex shedding. A linear stability analysis of the separated boundary layer is performed and the results show that the dominant wavenumber and frequency in the numerical simulations agree with the most unstable wavenumber and frequency from the linear stability analysis. The numerical results also show the growth and interaction of disturbance waves in the separation bubble. For transonic flows over the APEX airfoil, the calculations show that the presence of shocks causes significant changes in the separation location and consequently, the overall flowfield.

Tatineni, Mahidhar↗
NASA NTRS1999

A Numerical Analysis of Heat Transfer and Effectiveness on Film Cooled Turbine Blade Tip Models

A computational study has been performed to predict the distribution of convective heat transfer coefficient on a simulated blade tip with cooling holes. The purpose of the examination was to assess the ability of a three-dimensional Reynolds-averaged Navier-Stokes solver to predict the rate of tip heat transfer and the distribution of cooling effectiveness. To this end, the simulation of tip clearance flow with blowing of Kim and Metzger was used. The agreement of the computed effectiveness with the data was quite good. The agreement with the heat transfer coefficient was not as good but improved away from the cooling holes. Numerical flow visualization showed that the uniformity of wetting of the surface by the film cooling jet is helped by the reverse flow due to edge separation of the main flow.

Ameri, A. A.↗
NASA NTRS1999

Turbulence Model Comparisons for a High-Speed Aircraft

Four turbulence models are described and evaluated for transonic flows over the High-Speed Research/industry baseline configuration known as Reference H by using the thin-layer, upwind, Navier-Stokes solver known as CFL3D. The turbulence models studied are the equilibrium model of Baldwin-Lomax (B-L) with the Degani-Schiff (D-S) modifications, the one-equation Baldwin-Barth (B-B) model, the one-equation Spalart-Allmaras (S-A) model, and Menter's two-equation Shear Stress Transport (SST) model. The flow conditions, which correspond to tests performed in the National Transonic Facility (NTF) at Langley Research Center, are a Mach number of 0.90 and a Reynolds number of 30 x 10 (exp. 6) based on mean aerodynamic chord for angles of attack of 1 deg., 5 deg., and 10 deg. The effects of grid topology and the representation of the actual wind tunnel model geometry are also investigated. Computed forces and surface pressures compare reasonably well with the experimental data for all four turbulence models.

Rivers, Melissa B.↗
NASA NTRS1999

Development of a Multiobjective Optimization Procedure for Sonic Boom Minimization

A design optimization procedure for improved sonic boom and aerodynamic performance of high speed aircraft is presented. The multiobjective optimization procedure simultaneously minimizes the sonic boom at a given distance from the aircraft and the drag-to-lift ratio (C(sub D)/C(sub L)) Of the aircraft. Upper and lower bounds are also imposed on the lift coefficient. The Kreisselmeier - Steinhauser function is used for the multiobjective optimization formulation. A discrete semi-analytical aerodynamic sensitivity analysis procedure coupled with an analytical grid sensitivity analysis technique is used for evaluating design sensitivities. The use of the semi-analytical sensitivity analysis techniques results in significant computational savings. The flow equations are solved using a three-dimensional parabolized Navier-Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. A nonlinear programming technique and an approximate analysis procedure are used for the optimization. The optimization procedure developed is applied to the design of two high speed configurations, namely, a doubly swept wing-body configuration and a delta wing-body configuration. For the two sweep case only, minimization of the first peak in the pressure signature is performed first by optimizing only the nose radius and length of the aircraft. Minimization of the second peak in the pressure signature is performed next by optimizing only the wing geometric parameters. Significant improvements are obtained in the sonic boom characteristics and the aerodynamic performance of the wing-body configurations.

Narayan, J. R.↗
NASA NTRS1999

Viscous Analysis of Ref. H Based Wing/Bodies

Computations have been performed on the baseline Reference H wing/body configuration, as well as the Wing 704 configuration, an optimized wing and fuselage combination derived from Ref. H through automated optimization. The parabolized Navier-Stokes solver UPS was employed with viscous terms in two directions in an effort to understand the source and level of potential viscous/inviscid interactions. The paper briefly describes the UPS code and the grids used to obtain the solutions before the discussion of results. Results of these computations indicate that viscous/inviscid interaction can contribute increments to both the pressure- and friction-related drag. Computations were performed for wind tunnel conditions-1.675% scale models at a Reynolds number of 4 million per foot. Turbulent flow results were obtained using the Baldwin-Lomax algebraic turbulence model and were compared with laminar flow results. The laminar flow fields were used to obtain upper bounds on potential interaction effects.

Lawrence, Scott L.↗
NASA NTRS2000

Towards Flange-to-Flange Turbopump Simulations for Liquid Rocket Engines

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbopump geometry through numerical simulation will be of significant value toward design. This will help to improve safety of future space missions. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the MPI and MLP versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology. INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbopump simulations, moving boundary capability and an efficient time-accurate integration methods are build in the flow solver. To handle the geometric complexity and moving boundary problems, overset grid scheme is incorporated with the solver that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. The performance of the two time integration schemes for time-accurate computations is investigated. For an unsteady flow which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive. The current geometry for the LOX boost turbopump has various rotating and stationary components, such as inducer, stators, kicker, hydrolic turbine, where the flow is extremely unsteady. Figure 1 shows the geometry and computed surface pressure of the inducer. The inducer and the hydrolic turbine rotate in different rotational speed.

Kiris, Cetin↗
NASA NTRS1999

Aeroperformance and Acoustics of the Nozzle with Permeable Shell

Several simple experimental acoustic tests of a spraying system were conducted at the NASA Langley Research Center. These tests have shown appreciable jet noise reduction when an additional cylindrical permeable shell was employed at the nozzle exit. Based on these results, additional acoustic tests were conducted in the anechoic chamber AK-2 at the Central Aerohydrodynamics Institute (TsAGI, Moscow) in Russia. These tests examined the influence of permeable shells on the noise from a supersonic jet exhausting from a round nozzle designed for exit Mach number, M (sub e)=2.0, with conical and Screwdriver-shaped centerbodies. The results show significant acoustic benefits of permeable shell application especially for overexpanded jets by comparison with impermeable shell application. The noise reduction in the overall pressure level was obtained up to approximately 5-8%. Numerical simulations of a jet flow exhausting from a convergent-divergent nozzle designed for exit Mach number, M (sub e)=2.0, with permeable and impermeable shells were conducted at the NASA LaRC and Hampton University. Two numerical codes were used. The first is the NASA LaRC CFL3D code for accurate calculation of jet mean flow parameters on the basis of a full Navier-Stokes solver (NSE). The second is the numerical code based on Tam's method for turbulent mixing noise (TMN) calculation. Numerical and experimental results are in good qualitative agreement.

Gilinsky, M.↗
NASA NTRS2000

Force, Surface Pressure and Flowfield Measurements on Slender Missile Configurations at Supersonic Speeds

This paper describes a series of wind tunnel experiments carried out with the aim of providing data suitable for evaluating the performance of Computational Fluid Dynamics (CFD) codes. The configurations and flow conditions studied are most relevant to slender supersonic missiles. However, the data obtained, which includes forces and moments, surface pressures, flowfield surveys and a selection of flow visualization images, should he of interest to other CFD practitioners. Results for three test cases are presented and discussed in this paper. These cases have been the subject of a collaborative study concerned with the evaluation of Navier-Stokes solvers for missiles, carried out under the auspices of The Technical Cooperation Programme (TTCP).

Birch, T. J.↗
NASA NTRS2000

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

In this project, we continue to develop the previous joint research between the Fluid Mechanics and Acoustics Laboratory (FM&AL) at Hampton University (HU) and the Jet Noise Team (JNT) at the NASA Langley Research Center (NASA LaRC). The FM&AL was established at Hampton University in June of 1996 and has conducted research under two NASA grants: NAG-1-1835 (1996-99), and NAG-1-1936 (1997-00). In addition, the FM&AL has jointly conducted research with the Central AeroHydrodynamics Institute (TsAGI, Moscow) in Russia under a Civilian Research and Development Foundation (CRDF) grant #RE2-136 (1996-99). The goals of the FM&AL programs are twofold: (1) to improve the working efficiency of the FM&AUs team in generating new innovative ideas and in conducting research in the field of fluid dynamics and acoustics, basically for improvement of supersonic and subsonic aircraft engines, and (2) to attract promising minority students to this research and training and, in cooperation with other HU departments, to teach them basic knowledge in Aerodynamics, Gas Dynamics, and Theoretical and Experimental Methods in Aeroacoustics and Computational Fluid Dynamics (CFD). The research at the HU FM&AL supports reduction schemes associated with the emission of engine pollutants for commercial aircraft and concepts for reduction of observables for military aircraft. These research endeavors relate to the goals of the NASA Strategic Enterprise in Aeronautics concerning the development of environmentally acceptable aircraft. It is in this precise area, where the US aircraft industry, academia, and Government are in great need of trained professionals and which is a high priority goal of the Minority University Research and Education (MUREP) Program, that the HU FM&AL can make its most important contribution. The main achievements for the reporting period in the development of concepts for noise reduction and improvement in efficiency for jet exhaust nozzles and inlets for aircraft engines are as follows: (1) Publications- The AIAA Paper #99-1924 has been presented at the 5th AIAA/CEAS Aeroacoustics Conference, May 10-12, 1999, Seattle, WA; the AIAA Paper #00-3315 has been accepted for the 36th AIAA/ASME/ SAE/ASEE Joint Propulsion Conference, 17-19 July, 2000, Huntsville, AL; and another paper has been accepted for the International Environmental Congress, 14-16 June, 2000, St.-Petersburg, Russia. (2) Two patents were granted on July 20, 1999, and January 12, 2000. (3) Three reports/presentations at the NASA LaRC and GRC (06/22199, 09/26/ 99, and 06/25/00). (4) Grants and Proposals: Four proposals were submitted to the NASA and CRDF; a NASA Faculty Award was granted on January, 2000. A CRDF Young Investigator Program Award was granted for a 3 months visit of the Russian scientist to the HU FM&AL (03/99-05/99). (5) 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 codes based on full Euler and Navier-Stokes solvers: CFL3D, CRAFT, GODUNOV, and others. New approach for environmental monitoring via infrasound. (6) Experimental Tests: Experimental acoustic tests at the TsAGI, Moscow, with nozzles having Screwdriver or Axisymmetric Plug and Permeable Shells. A small scale working model of the NASA Low Speed Wind Tunnel (LSWT) has been installed in the Experimental Hall of the HU FM&AL (June, 1999). Preliminary preparations for experimental tests were made. (7) Students Research Activity: Involvement of the two graduate students as research assistants in the current research project.

Gilinsky, Mikhail↗
NASA NTRS1998

Evaluation of Turbulence-Model Performance as Applied to Jet-Noise Prediction

The accurate prediction of jet noise is possible only if the jet flow field can be predicted accurately. Predictions for the mean velocity and turbulence quantities in the jet flowfield are typically the product of a Reynolds-averaged Navier-Stokes solver coupled with a turbulence model. To evaluate the effectiveness of solvers and turbulence models in predicting those quantities most important to jet noise prediction, two CFD codes and several turbulence models were applied to a jet configuration over a range of jet temperatures for which experimental data is available.

Woodruff, S. L.↗