Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes flow”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Reduction of vortex-induced vibration in vane geometries

Computations using a time-accurate, compressible Navier-Stokes flow model were conducted to analyze both unsteady laminar and turbulent flows over two curved configurations of vanes which were shaped to treat a 4.0 kHz vibration problem which has occurred in several SSMEs. This problem involves vortex shedding from vanes which causes excessive vane vibration and cracking. The original vane configuration exhibited strong flow-induced vibrations at a Strouhal number near 0.19 for a the first bending mode excitation. Scalloping of the leading edge raised the frequency of the vane's first torsional mode, which in turn increased the onset flow velocity at lock-in. Beveling the vane's trailing edge eliminated vortex shedding at the trailing edge, which decoupled the flow oscillation from the vibrating wave. A modified vane configuration with a beveled trailing edge at a 30 deg angle was also studied.

Liu, B. L.↗
NASA NTRS1993

Computation of wake/exhaust mixing downstream of advanced transport aircraft

The mixing of engine exhaust with the vortical wake of high speed aircraft operating in the stratosphere can play an important role in the formation of chemical products that deplete atmospheric ozone. An accurate analysis of this type of interaction is therefore necessary as a part of the assessment of the impact of proposed High Speed Civil Transport (HSCT) designs on atmospheric chemistry. This paper describes modifications to the parabolic Navier-Stokes flow field analysis in the UNIWAKE unified aircraft wake model to accommodate the computation of wake/exhaust mixing and the simulation of reacting flow. The present implementation uses a passive chemistry model in which the reacting species are convected and diffused by the fluid dynamic solution but in which the evolution of the species does not affect the flow field. The resulting analysis, UNIWAKE/PCHEM (Passive CHEMistry) has been applied to the analysis of wake/exhaust flows downstream of representative HSCT configurations. The major elements of the flow field model are described, as are the results of sample calculations illustrating the behavior of the thermal exhaust plume and the production of species important to the modeling of condensation in the wake. Appropriate steps for further development of the UNIWAKE/PCHEM model are also outlined.

Quackenbush, Todd R.↗
NASA NTRS1993

Numerical study of 3-D inducer and impeller for pump model development

Current design of high-performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study is to develop a robust and effective CFD pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, two key components of the turbopump, the inducer and impeller, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne.

Cheng, G. C.↗
NASA NTRS1993

A hydrodynamic model of the oscillating screen viscometer

The viscometer consists of an oscillating screen immersed in a fluid and free to rotate about an axis in its plane. The viscosity can be determined from the measured ratio of the periodic driving force to the screen motion when an adequate hydrodynamical model of the immersed oscillator is available. The screen is formed by a square mesh of thin wire whose dimensions invite comparison with asymptotic results for narrow hollow bodies translating in Stokes flow. These indicate that the closed hole structure of the grid plays an important role in determining its motion. It is shown that this role diminishes as the frequency increases. The computed results, obtained from a system of linear equations, are consistent with experimental values over the appropriate range of frequency.

Davis, A. M. J.↗
NASA NTRS1994

Inducer analysis/pump model development

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study was to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, three key components of the turbopump, the inducer, impeller, and diffuser, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne. For the numerical calculation of inducer flows with tip clearance, the turbulence model and grid spacing are very important. Meanwhile, the development of the cross-stream secondary flow, generated by curved blade passage and the flow through tip leakage, has a strong effect on the inducer flow. Hence, the prediction of the inducer performance critically depends on whether the numerical scheme of the pump model can simulate the secondary flow pattern accurately or not. The impeller and diffuser, however, are dominated by pressure-driven flows such that the effects of turbulence model and grid spacing (except near leading and trailing edges of blades) are less sensitive. The present CFD pump model has been proved to be an efficient and robust analytical tool for pump design due to its very compact numerical structure (requiring small memory), fast turnaround computing time, and versatility for different geometries.

Cheng, Gary C.↗
NASA NTRS1994

Numerical design of advanced multi-element airfoils

The current study extends the application of computational fluid dynamics to three-dimensional high-lift systems. Structured, overset grids are used in conjunction with an incompressible Navier-Stokes flow solver to investigate flow over a two-element high-lift configuration. The computations were run in a fully turbulent mode using the one-equation Baldwin-Barth turbulence model. The geometry consisted of an unswept wing which spanned a wind tunnel test section. Flows over full and half-span Fowler flap configurations were computed. Grid resolution issues were investigated in two dimensional studies of the flapped airfoil. Results of the full-span flap wing agreed well with experimental data and verified the method. Flow over the wing with the half-span was computed to investigate the details of the flow at the free edge of the flap. The results illustrated changes in flow streamlines, separation locations, and surface pressures due to the vortex shed from the flap edge.

Mathias, Donovan L.↗
NASA NTRS1994

Applications of automatic differentiation in computational fluid dynamics

Automatic differentiation (AD) is a powerful computational method that provides for computing exact sensitivity derivatives (SD) from existing computer programs for multidisciplinary design optimization (MDO) or in sensitivity analysis. A pre-compiler AD tool for FORTRAN programs called ADIFOR has been developed. The ADIFOR tool has been easily and quickly applied by NASA Langley researchers to assess the feasibility and computational impact of AD in MDO with several different FORTRAN programs. These include a state-of-the-art three dimensional multigrid Navier-Stokes flow solver for wings or aircraft configurations in transonic turbulent flow. With ADIFOR the user specifies sets of independent and dependent variables with an existing computer code. ADIFOR then traces the dependency path throughout the code, applies the chain rule to formulate derivative expressions, and generates new code to compute the required SD matrix. The resulting codes have been verified to compute exact non-geometric and geometric SD for a variety of cases. in less time than is required to compute the SD matrix using centered divided differences.

Green, Lawrence L.↗
NASA NTRS1993

CFD Parametric Study of Consortium Impeller

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design impact in a productive manner. The main goal of this study is to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A Finite Difference Navier-Stokes flow solver, FDNS, which includes the extended k-epsilon turbulence model and appropriate moving interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. A second-order central difference scheme plus adaptive dissipation terms was employed in the FDNS code, along with a predictor plus multi-corrector pressure-based solution procedure. The multi-zone, multi-block capability allows the FDNS code to efficiently solve flow fields with complicated geometry. The FDNS code has been benchmarked by analyzing the pump consortium inducer, and it provided satisfactory results. In the present study, a CFD parametric study of the pump consortium impeller was conducted using the FDNS code. The pump consortium impeller, with partial blades, is a new design concept of the advanced rocket engines. The parametric study was to analyze the baseline design of the consortium impeller and its modification which utilizes TANDEM blades. In the present study, the TANDEM blade configuration of the consortium impeller considers cut full blades for about one quarter chord length from the leading edge and clocks the leading edge portion with an angle of 7.5 or 22.5 degrees. The purpose of the present study is to investigate the effect and trend of the TANDEM blade modification and provide the result as a design guideline. A 3-D flow analysis, with a 103 x 23 x 30 mesh grid system and with the inlet flow conditions measured by Rocketdyne, was performed for the baseline consortium impeller. The numerical result shows that the mass flow rate splits through various blade passages are relatively uniform. Due to the complexity of blade geometries, the TANDEM blade configurations were analyzed with the multi-zone grid structure. Both the 7.5 deg- and the 22.5 deg-clocking TANDEM blade cases utilized a 80K mesh system. The numerical result of two TANDEM blade modifications indicates the efficiency and the head are worse than those of the baseline case due to larger flow distortion. The gap between the TANDEM blade and the full blade allows the flow passes through and heavily loads the pressure side of the partial blade such that flow reversal occurs near the suction side of the splitter. The flow split at the exit of impeller blades is very non-uniform for TANDEM blade cases, and this will greatly induce the side load on the diffuser. Therefore, the TANDEM blade modification in the present CFD analysis does not improve the performance of the consortium impeller.

Cheng, Gary C.↗
NASA NTRS1995

Development of numerical methods for overset grids with applications for the integrated Space Shuttle vehicle

Algorithms and computer code developments were performed for the overset grid approach to solving computational fluid dynamics problems. The techniques developed are applicable to compressible Navier-Stokes flow for any general complex configurations. The computer codes developed were tested on different complex configurations with the Space Shuttle launch vehicle configuration as the primary test bed. General, efficient and user-friendly codes were produced for grid generation, flow solution and force and moment computation.

Chan, William M.↗
NASA NTRS1995

The effect of adding roughness and thickness to a transonic axial compressor rotor

The performance deterioration of a high speed axial compressor rotor due to surface roughness and airfoil thickness variations is reported. A 0.025 mm (0.001 in.) thick rough coating with a surface finish of 2.54-3.18 RMS microns (100-125 RMS microinches) is applied to the pressure and suction surface of the rotor blades. Coating both surfaces increases the leading edge thickness by 10% at the hub and 20% at the tip. Application of this coating results in a loss in efficiency of 6 points and a 9% reduction in the pressure ratio across the rotor at an operating condition near the design point. To separate the effect of thickness and roughness, a smooth coating of equal thickness is also applied to the blade. The smooth coating surface finish is 0.254-0.508 RMS microns (10-20 RMS microinches), compared to the bare metal blade surface finish of 0.508 RMS microns (20 RMS microinches). The smooth coating results in approximately half of the performance deterioration found from the rough coating. Both coatings are then applied to different portions of the blade surface to determine which portions of the airfoil are most sensitive to thickness/roughness variations. Aerodynamic performance measurements are presented for a number of coating configurations at 60%, 80%, and 100% of design speed. The results indicate that thickness/roughness over the first 10% of blade chord accounts for virtually all of the observed performance degradation for the smooth coating, compared to about 70% of the observed performance degradation for the rough coating. The performance deterioration is investigated in more detail at design speed using laser anemometer measurements as well as predictions generated by a quasi-3D Navier-Stokes flow solver which includes a surface roughness model. Measurements and analysis are performed on the baseline blade and the full-coverage smooth and rough coatings. The results indicate that coating the blade causes a thickening of the blade boundary layers. The interaction between the rotor passage shock and the thickened suction surface boundary layer then results in an increase in blockage which reduces the diffusion level in the rear half of the blade passage, thus reducing the aerodynamic performance of the rotor.

Suder, Kenneth L.↗
NASA NTRS1993

Mathematical description of complex chemical kinetics and application to CFD modeling codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗
NASA NTRS1995

Parallel CFD design on network-based computer

Combining multiple engineering workstations into a network-based heterogeneous parallel computer allows application of aerodynamic optimization with advanced computational fluid dynamics codes, which can be computationally expensive on mainframe supercomputers. This paper introduces a nonlinear quasi-Newton optimizer designed for this network-based heterogeneous parallel computing environment utilizing a software called Parallel Virtual Machine. This paper will introduce the methodology behind coupling a Parabolized Navier-Stokes flow solver to the nonlinear optimizer. This parallel optimization package is applied to reduce the wave drag of a body of revolution and a wing/body configuration with results of 5% to 6% drag reduction.

Cheung, Samson↗
NASA NTRS1996

Observations Regarding Use of Advanced CFD Analysis, Sensitivity Analysis, and Design Codes in MDO

Observations regarding the use of advanced computational fluid dynamics (CFD) analysis, sensitivity analysis (SA), and design codes in gradient-based multidisciplinary design optimization (MDO) reflect our perception of the interactions required of CFD and our experience in recent aerodynamic design optimization studies using CFD. Sample results from these latter studies are summarized for conventional optimization (analysis - SA codes) and simultaneous analysis and design optimization (design code) using both Euler and Navier-Stokes flow approximations. The amount of computational resources required for aerodynamic design using CFD via analysis - SA codes is greater than that required for design codes. Thus, an MDO formulation that utilizes the more efficient design codes where possible is desired. However, in the aerovehicle MDO problem, the various disciplines that are involved have different design points in the flight envelope; therefore, CFD analysis - SA codes are required at the aerodynamic 'off design' points. The suggested MDO formulation is a hybrid multilevel optimization procedure that consists of both multipoint CFD analysis - SA codes and multipoint CFD design codes that perform suboptimizations.

Newman, Perry A.↗
NASA NTRS1997

Parallelization of the Implicit RPLUS Algorithm

The multiblock reacting Navier-Stokes flow solver RPLUS2D was modified for parallel implementation. Results for non-reacting flow calculations of this code indicate parallelization efficiencies greater than 84% are possible for a typical test problem. Results tend to improve as the size of the problem increases. The convergence rate of the scheme is degraded slightly when additional artificial block boundaries are included for the purpose of parallelization. However, this degradation virtually disappears if the solution is converged near to machine zero. Recommendations are made for further code improvements to increase efficiency, correct bugs in the original version, and study decomposition effectiveness.

Orkwis, Paul D.↗
NASA NTRS1994

Off-Axis Drag of Dendrite Fragments at Low Reynolds Number

The aim of the present investigation is to characterize the motion of dendrite fragments falling under the influence of gravity in a uniform liquid medium at low Reynolds number. In an earlier study, Zakhem, Weidman and de Groh (1992) reported on the settling speed of model equiaxed dendrite grains released along their axis of symmetry. In this follow-up study uniaxial model dendrite grains were released off-axis to observe and document their motion at different orientations. It was hypothesized that the dendrite models might rotate when released off-axis in which case an attempt would be made to document the ensuing unsteady motion. This latter event turned out to be in fact true: at the small but finite Reynolds numbers that existed, each uniaxial dendrite slowly rotated towards its equilibrium orientation while failing under the influence of gravity. In addition to completing the original goal, we have made use of a beads-on-a shell Stokes flow code to numerically determine the drag coefficient for capsules, i.e.. uniaxial dendrites without arms. The drag on horizontally and vertically falling capsules are reported and compared with measurements.

Weidman, P. D.↗
NASA NTRS1998

CFL3D User's Manual (Version 5.0)

This document is the User's Manual for the CFL3D computer code, a thin-layer Reynolds-averaged Navier-Stokes flow solver for structured multiple-zone grids. Descriptions of the code's input parameters, non-dimensionalizations, file formats, boundary conditions, and equations are included. Sample 2-D and 3-D test cases are also described, and many helpful hints for using the code are provided.

Krist, Sherrie L.↗
NASA NTRS1999

Performance and Application of Parallel OVERFLOW Codes on Distributed and Shared Memory Platforms

The presentation discusses recent studies on the performance of the two parallel versions of the aerodynamics CFD code, OVERFLOW_MPI and _MLP. Developed at NASA Ames, the serial version, OVERFLOW, is a multidimensional Navier-Stokes flow solver based on overset (Chimera) grid technology. The code has recently been parallelized in two ways. One is based on the explicit message-passing interface (MPI) across processors and uses the _MPI communication package. This approach is primarily suited for distributed memory systems and workstation clusters. The second, termed the multi-level parallel (MLP) method, is simple and uses shared memory for all communications. The _MLP code is suitable on distributed-shared memory systems. For both methods, the message passing takes place across the processors or processes at the advancement of each time step. This procedure is, in effect, the Chimera boundary conditions update, which is done in an explicit "Jacobi" style. In contrast, the update in the serial code is done in more of the "Gauss-Sidel" fashion. The programming efforts for the _MPI code is more complicated than for the _MLP code; the former requires modification of the outer and some inner shells of the serial code, whereas the latter focuses only on the outer shell of the code. The _MPI version offers a great deal of flexibility in distributing grid zones across a specified number of processors in order to achieve load balancing. The approach is capable of partitioning zones across multiple processors or sending each zone and/or cluster of several zones into a single processor. The message passing across the processors consists of Chimera boundary and/or an overlap of "halo" boundary points for each partitioned zone. The MLP version is a new coarse-grain parallel concept at the zonal and intra-zonal levels. A grouping strategy is used to distribute zones into several groups forming sub-processes which will run in parallel. The total volume of grid points in each group are approximately balanced. A proper number of threads are initially allocated to each group, and in subsequent iterations during the run-time, the number of threads are adjusted to achieve load balancing across the processes. Each process exploits the multitasking directives already established in Overflow.

Djomehri, M. Jahed↗
NASA NTRS1999

OVERAERO-MPI: Parallel Overset Aeroelasticity Code

An overset modal structures analysis code was integrated with a parallel overset Navier-Stokes flow solver to obtain a code capable of static aeroelastic computations. The new code was used to compute the static aeroelastic deformation of an arrow-wing-body geometry and a complex, full aircraft configuration. For the simple geometry, the results were similar to the results obtained with the ENSAERO code and the PVM version of OVERAERO. The full potential of this code suite was illustrated in the complex, full aircraft computations.

Gee, Ken↗