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 649 records · Page 36

Numerical solution of unsteady rotational flow past fixed and rotary wing configurations

The application of unsteady 3-D Euler and Navier-Stokes equations to transonic flow past rotor blades, and wing-alone configurations is described. A promising approach for the numerical solution of these equations is examined. Additional work is needed for improving the efficiency of the present procedure. It is hoped that the techniques presented will find use in fixed and rotary wing aircraft analysis.

Sankar, N. L.↗

Navier-Stokes and viscous-inviscid interaction

Some considerations toward developing numerical procedures for simulating viscous compressible flows are discussed. Both Navier-Stokes and boundary layer field methods are considered. Because efficient viscous-inviscid interaction methods have been difficult to extend to complex 3-D flow simulations, Navier-Stokes procedures are more frequently being utilized even though they require considerably more work per grid point. It would seem a mistake, however, not to make use of the more efficient approximate methods in those regions in which they are clearly valid. Ideally, a general purpose compressible flow solver that can optionally take advantage of approximate solution methods would suffice, both to improve accuracy and efficiency. Some potentially useful steps toward this goal are described: a generalized 3-D boundary layer formulation and the fortified Navier-Stokes procedure.

Steger, Joseph L.↗

Incompressible Navier-Stokes Solvers in Primative Variables and their Applications to Steady and Unsteady Flow Simulations

This paper reviews recent progress made in incompressible Navier-Stokes simulation procedures and their application to problems of engineering interest. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. A summary of efforts in flow solver development is given followed by numerical studies of a few example problems of current interest. Both steady and unsteady solution algorithms and their salient features are discussed. Solvers discussed here are based on a structured-grid approach using either a finite -difference or a finite-volume frame work. As a grand-challenge application of these solvers, an unsteady turbopump flow simulation procedure has been developed which utilizes high performance computing platforms. In the paper, the progress toward the complete simulation capability of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of two parallel computing algorithms that have been implemented in the INS3D code. The relative motion of the grid systems for the rotorstator interaction was obtained using overact grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on SCSI Origin 3000 systems at NASA Ames Research Center. The same procedure has been extended to the development of NASA-DeBakey Ventricular Assist Device (VAD) that is based on an axial blood pump. Computational, and clinical analysis of this device are presented.

Kiris, Cetin C.↗

Steady flow past sudden expansions at large Reynolds number. II - Navier-Stokes solutions for the cascade expansion

The equations of motion in the steady laminar flow past a sudden expansion at large Reynolds number, R, reduce to the boundary layer equations as R tends to the freestream valve, when the longitudinal length scale of the separated eddy increases linearly and indefinitely with R. A global Newton method is presently used to obtain finite difference solutions to the steady Navier-Stokes equations up to R of 1000 for a uniform inflow past a cascade of sudden expansions. For large expansion ratio values, the eddy length increases linearly with R; for smaller values, however, where boundary layer equation solutions could not be found, the steady solutions to the Navier-Stokes equations approach the limit of an inviscid eddy in length with increasing R.

Milos, Frank S.↗

A simple numerical orthogonal coordinate generator for fluid dynamic applications

An application of a simple numerical technique which allows for the rapid construction of orthogonal coordinate systems about two dimensional and axisymmetric bodies is presented. The technique is based on a predictor corrector numerical method. It can be used to generate orthogonal meshes which have unequally spaced points in two directions. These orthogonal meshes in their transformed computational plane are, however, equally spaced so that the differencing for the metric coefficients and the fluid dynamic equation terms can be easily determined using equally spaced central finite differences. Solutions to the Navier-Stokes equations for flow over blunt bodies with reverse curvature are presented. The coupling of the time dependent fluid dynamic equations and the coordinate generator worked well with no undesirable effects noted.

Graves, R. A., Jr.↗

Code Solves Three-Dimensional Navier-Stokes Equations

Set of computer codes solves three-dimensional Navier-Stokes equations for flow over nonaxisymmetric nozzles. Codes compute internal and external viscous flowfield about isolated nozzle, so flow characteristics and performance of three-dimensional jet engine exhaust nozzles can be predicted. Programs written in FORTRAN IV and ASSEMBLER.

Thomas, P.↗

Airfoil computation at high angles of attack, inviscid and viscous phenomena

An implicit central difference code is used to calculate two dimensional inviscid and thin-layer Navier-Stokes solutions for flow about an NACA0012 airfoil at high angles of attack. Among the issues addressed are whether separation can occur in an inviscid calculation and what the causes would be of such separation. Examples are shown of inviscid shocked flow with and without separation and shock-free flow with separation. An Euler solution with self-induced oscillation and separation driven by a strong shock is contrasted with a shock-free solution whose separation is caused by numerical error. Computed solutions to the Euler equations are compared to those of the potential equations. Comparisons are also made between experimental data from wind tunnel tests and viscous calculations at similar conditions.

Barton, J. T.↗

Numerical study of large-eddy breakup and its effect on the drag characteristics of boundary layers

The break-up of a field of eddies by a flat-plate obstacle embedded in a boundary layer is studied using numerical solutions to the two-dimensional Navier-Stokes equations. The flow is taken to be incompressible and unsteady. The flow field is initiated from rest. A train of eddies of predetermined size and strength are swept into the computational domain upstream of the plate. The undisturbed velocity profile is given by the Blasius solution. The disturbance vorticity generated at the plate and wall, plus that introduced with the eddies, mix with the background vorticity and is transported throughout the entire flow. All quantities are scaled by the plate length, the unidsturbed free-stream velocity, and the fluid kinematic viscosity. The Reynolds number is 1000, the Blasius boundary layer thickness is 2.0, and the plate is positioned a distance of 1.0 above the wall. The computational domain is four units high and sixteen units long.

Kinney, R. B.↗

Hypersonic turbulent wall boundary layer computations

The Baldwin-Lomax algebraic turbulence model was modified for hypersonic flow conditions. Two coefficients in the outer layer eddy viscosity model were determined as functions of Mach number and temperature ratio. By matching the solutions from the Baldwin-Lomax model to those from the Cebeci-Smith model for a flat plate at hypersonic speed, the new values of the coefficient were obtained. The results show that the values of C sub cp and C sub kleb are functions of both Mach number and wall temperature ratio. The C sub cp and C sub kleb variations with Mach number and wall temperature were used for the calculations of both a 4 deg wedge flow at Mach 18 and an axisymmetric Mach 20 nozzle flow. The Navier-Stokes equations with thin layer approximation were solved for the above hypersonic flow conditions and the results were compared with existing experimental data. The agreement between the numerical solutions and the existing experimental data were good. The modified Baldwin-Lomax model thus is useful in the computations of hypersonic flows.

Kim, S. C.↗

Direct numerical study of leading-edge contamination

Instability, turbulence, and relaminarization in the attachment-line region of swept and unswept cylindrical bodies are studied by numerical solution of the full Navier-Stokes equations. The flow is simulated over a strip containing the attachment-line and treated as homogeneous in the spanwise direction; the disturbances decay exponentially upstream. Transpiration through the wall may be prescribed. The new method, which admits completely general disturbance, agrees with published linear-stability results, which were limited to an apparently restrictive form of disturbance. Fully developed turbulent solutions with sweep are generated and compare well with the experiment. The turbulence is subcritical (except for blowing), resulting in large hysteresis loops. By lowering the sweep Reynolds number, or increasing the suction, the turbulent flow is made to relaminarize. The relaminarization Reynolds number is much less sensitive to suction than the linear-stability Reynolds number. Extensive attempts to detect the postulated nonlinear instability of the unswept flow failed, suggesting that this flow is linearly and nonlinearly stable.

Spalart, P. R.↗

CFD code evaluation for internal flow modeling

Research on the computational fluid dynamics (CFD) code evaluation with emphasis on supercomputing in reacting flows is discussed. Advantages of unstructured grids, multigrids, adaptive methods, improved flow solvers, vector processing, parallel processing, and reduction of memory requirements are discussed. As examples, researchers include applications of supercomputing to reacting flow Navier-Stokes equations including shock waves and turbulence and combustion instability problems associated with solid and liquid propellants. Evaluation of codes developed by other organizations are not included. Instead, the basic criteria for accuracy and efficiency have been established, and some applications on rocket combustion have been made. Research toward an ultimate goal, the most accurate and efficient CFD code, is in progress and will continue for years to come.

Chung, T. J.↗

High alpha aerodynamic control by tangential fuselage blowing

This work explores the effect of tangential blowing on the vortical structures that develop around a tangent-ogive cylinder configuration at high angle of attack. A lateral force results if blowing is applied asymmetrically. The study is conducted numerically by solving the three-dimensional, compressible-flow Navier-Stokes equations. The computation was done for a Reynolds number of 52,000, Mach number of 0.2, blowing momentum coefficients of 0.0, 0.1, 0.2, and 0.4, and angle of attack of 10, 30, and 45 deg. Only asymmetrical blowing was considered. Surface streamlines, helicity contours and pressure distributions are obtained for each case.

Font, G. I.↗

Computer Code Generates Homotopic Grids

HOMAR is computer code using homotopic procedure to produce two-dimensional grids in cross-sectional planes, which grids then stacked to produce quasi-three-dimensional grid systems for aerospace configurations. Program produces grids for use in both Euler and Navier-Stokes computation of flows. Written in FORTRAN 77.

Moitra, Anutosh↗

Comparison Of Two Viscous-Flow Computer Codes

Two viscous-flow computer codes compared by applying them to five test cases of steady-state transonic viscous flows about transonic airfoils. Two codes were: FLOMG, which solves Navier-Stokes differential equations of flow by implementing explicit, Runge-Kutta, finite-volume, multigrid numerical-integration algorithm; and ARC2D, which implements implicit, finite-difference, approximate-factorization, eigenvector-diagonalization numerical-integration algorithm.

Maksymiuk, C. M.↗

Parallelization of a Parabolized Navier-Stokes Solver with a Design Optimizer

The design of future supersonic aircraft, such as the High Speed Civil Transport (HSCT), will rely heavily on computational methods for aircraft design and the prediction of the complex aerodynamics encountered in flight. Parabolized Navier-Stokes (PNS) equation flow solvers are recognized as efficient and accurate computational tools for the solution of supersonic and hypersonic flow-fields, while design optimizers have the potential to be valuable tools within the overall design process. Presently, however, the execution of the flow solver in conjunction with a design optimizer presents a computationally intensive and formidable problem. To meet the challenges and increasing demand for multidisciplinary numerical tools which are faster, more robust and provide greater functionality, alternative strategies are explored to increase computational throughput by coupling a design optimizer and flow solver in a parallel processing environment. To address this problem the parallel processing of a PNS flow solver with a nonlinear constraint design optimizer is investigated as an alternative computational method.

Pallis, J. M.↗

Direct Numerical Simulation of a Temporally Evolving Incompressible Plane Wake: Effect of Initial Conditions on Evolution and Topology

Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Sondergaard, R.↗