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 1,657 records · Page 92

Navier-Stokes solutions using Stetter's method

Compressible Navier-Stokes equations for quasi-one-dimensional flow in a converging-diverging nozzle have been solved using Stetter's three-step predictor-corrector technique. Particular emphasis is given to the minimum iterative step feature for steady-state solutions. It is found that for the nonviscous-dominated case, Stetter's method attains the steady-state solution in the fewest steps when compared to four other currently used techniques.

Graves, R. A., Jr.↗

Heating analysis of bent-nose biconics at high angles of attack using the parabolized Navier-Stokes equations

A method based on the Parabolized Navier-Stokes equations is used to calculate the flow field and heat transfer of lifting entry vehicles. The method is based on the Bean and Warming implicit algorithm and uses a new procedure for preventing departure solutions. Calculations are carried out for blunt on-axis and bent biconics, assuming a perfect gas and laminar flow, and compared with available heat transfer, surface pressure and shock shape measurements for a range of Mach numbers and angles of attack. In all calculations presented here, the starting solution is obtained from available inviscid and boundary layer codes. Good agreement with experiment is indicated. Thus, the method provides an accurate and rather inexpensive procedure for calculating three-dimensional flows at supersonic Mach numbers.

Stephenson, B. L.↗

Numerical simulation of the leading-edge separation vortex for a wing and strake-wing configuration

In the present investigation, the 'thin layer' Navier-Stokes equations are used to compute the flow about a delta wing and a strake-delta wing configuration. Both configurations possess blunt noses and rounded leading edges. The computational grid about these configurations is generated using a newly developed three-dimensional grid-generation code which solves a set of Poisson equations with spherical coordinate variables using an alternating direction implicit (ADI) scheme. Computational results are obtained for an isolated wing with 60 deg sweep and a strake-wing configuration with 80-60 deg sweep. Computations for angles of attack in the range from 6 to 30 deg are presented. Attention is given to the effect of the angle of attack, the Mach number, and the Reynolds number. Theoretical results are compared with experimental data.

Fujii, K.↗

Computations of two-dimensional airfoil-vortex interactions

A procedure for calculating the interaction of a vortex with a two dimensional airfoil in a uniform free stream is presented along with results for several test cases. A Lamb like anaytical vortex having a finite core and convect in a uniform free stream interacts with the flowfield of NACA 0012 or NACA 64A006 airfoil in transonic and subsonic flow. Euler and thin layer Navier-Stokes solutions are computed and the results are compared with the results from transonic small disturbance code and available experimental results. For interactions within the limits of transonic small disturbance assumptions, the three methods gave qualitatively similar results of a two bladed helicopter rotor and suggest that the time lag effects of the free stream velocity approaching the blade may be important and should be considered in the analysis. In general, the results show a tremendous influence of the interacting vortex on the flowfield around the airfoil. This is particularly true when the vortex is stationary. For a convecting vortex, the most dramatic changes in the flowfield seem to occur when the vortex is within one chord of the airfoil.

Srinivasan, G. R.↗

Numerical simulation and comparison with experiment for self-excited oscillations in a diffuser flow

This paper describes numerical simulations of self-excited oscillations in a two-dimensional transonic diffuser flow obtained by solving the Navier-Stokes equations with a two-equation turbulence model. Comparisons were made between the computational results and experimental data. For the mean flowfields, the agreement between computation and experiment is good for the wall pressures, shock location, and the separation and reattachment points. However, the thickness of the computed recirculation zone is about 50 percent of the measured thickness. For the fluctuating flowfields, a great deal of qualitative similarity exists between the computation and experiment; however, the predicted oscillation frequency is about 50 percent higher than the measured value. The formation of a succession of downstream-traveling counter-rotating vortices, as seen experimentally, is also vividly displayed in the numerical results.

Hsieh, T.↗

Developments in the simulation of separated flows using finite difference methods

Compressible viscous flow simulation using finite difference Navier-Stokes and viscous-inviscid interaction methods is described. Recent developments are reviewed that significantly improve the computational efficiency of approximately factored implicit Navier-Stokes algorithms. Compared to Navier-Stokes codes, modern viscous-inviscid interaction codes are more computationally efficient, but have restricted application and are more complicated to program. Therefore, less efficient but more general viscous-inviscid interaction methods are investigated that use forcing functions instead of boundary condition matching, and a simple, direct/inverse, three-dimensional, finite-difference, boundary layer code is presented.

Steger, J. L.↗

Steady supersonic Navier-Stokes solutions of a 75 deg delta wing

Steady solutions about a slender sharp edged delta wing in a supersonic freestream for moderate and high angles of attack are obtained numerically by time integration of the unsteady compressible three dimensional laminar Navier-Stokes equations. The main features of the flow, including primary and secondary separation, and vortex position and strength, are adequately simulated in the numerical solutions. Improved resolution of the computational grid in the leading edge region from a previous solution had considerable effect on the accuracy of the solutions. Good agreement between numerical solutions and experimental data was obtained for two cases. A local timestepping procedure is used to speed convergence by approximately a factor of two.

Buter, T. A.↗

Direct numerical simulation of equilibrium turbulent boundary layers

This paper describes the simulation of turbulent boundary layers by direct numerical solution of the three-dimensional, time-dependent Navier-Stokes equations, using a spectral method. The flow is incompressible, with Re sub delta = 1000. The equations are written in the self-similar coordinate system and periodic streamwise and spanwise boundary conditions are imposed. A family of nine 'equilibrium' boundary layers, from the strongly accelerated 'sink' flow to Stratford's separating boundary layer is treated. Good general agreement with experiments is observed. The effects of pressure gradients on the structures and statistics, both in the wall and wake regions are discussed.

Spalart, P. R.↗

On the instability of a three-dimensional attachment-line boundary layer - Weakly nonlinear theory and a numerical approach

The instability of a three-dimensional attachment-line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite-amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time-dependent Navier-Stokes equations for the attachment-line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite-amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment-line boundary layer is also investigated.

Hall, P.↗

Implicit hybrid schemes for the flux-difference split, three-dimensional Navier-Stokes equations

Implicit hybrid algorithms employing symmetric planar Gauss-Seidel (SPGS) relaxation and either block-tridiagonally structured coefficient matrices (AF-SPGS) or block-triangular coefficient matrices (LU-SPGS) are derived to solve the flux-difference-split Navier-Stokes equations for three-dimensional incompressible flow in an upwind scheme. The physical basis of the approach is discussed, and results for problems involving vortex flow around a thin delta wing at Reynolds numbers 900,000 and 10,000 are presented graphically. It is found that AF-SPGS converges faster on vector computers which depend on long vector lengths to achieve optimum performance, whereas LU-SPGS is preferable on sequentially operating machines and vector computers using shorter vector lengths.

Hartwich, P.-M.↗

Simulation of transonic viscous wing and wing-fuselage flows using zonal methods

The thin-layer Navier-Stokes equations are coupled with a zonal scheme (or domain-decomposition method) to develop the Transonic Navier-Stokes (TNS) wing-alone code. The TNS has a total of 4 zones and is extended to a total of 16 zones for the wing-fuselage version of the code. Results are compared on the Cray X-MP-48 and compared with experimental data.

Flores, Jolen↗

Computational fluid dynamics in a marine environment

The introduction of the supercomputer and recent advances in both Reynolds averaged, and large eddy simulation fluid flow approximation techniques to the Navier-Stokes equations, have created a robust environment for the exploration of problems of interest to the Navy in general, and the Naval Underwater Systems Center in particular. The nature of problems that are of interest, and the type of resources needed for their solution are addressed. The goal is to achieve a good engineering solution to the fluid-structure interaction problem. It is appropriate to indicate that a paper by D. Champman played a major role in developing the interest in the approach discussed.

Carlson, Arthur D.↗

Migration of the separation point on a deforming cylinder

An iterative scheme of solving the Navier-Stokes equations for unsteady two-dimensional flow about a deforming and translating cylinder is given. In the kth iteration, the convected vorticity appears as the source term in an equation of transient diffusion of vorticity. A novel integral transform is used to reduce this transient vorticity equation to a k-dimensional heat equation. The bounded solution of this equation is obtained with a general method of superposition for problems involving a moving boundary. An equation describing the migration of the separation point on a deforming cylinder in unsteady cross flows is derived from the analytically obtained velocity field. The radial expansion of the cylinder surface is shown to hasten the separation time and to increase the separation angle. The results imply that in the steady flow past a body that is moving forward and sinking at constant rates, the locus of points at which the azimuthal component of skin friction changes sign originates on the leeward ray at a point downstream of the front tip.

Lin, S. P.↗

Nearly two-dimensional solutions of Euler's equations

Motivated by a stretched version of the Taylor-Green initial-value problem for the Euler and Navier-Stokes equations, nearly two-dimensional incompressible flows are considered involving a single direction of slow spatial variation. A multiple-scale formulation of the problem leads by contour averaging to a system that determines the slow development in space and time. It is shown that the latter system is equivalent to an axisymmetric problem with nonstandard connection between circulation, cylindrical radius, and angular velocity component. Special solutions of the system in the exact axisymmetric case suggest that near-two-dimensionality is lost in finite time.

Childress, Stephen↗

Turbulence statistics in fully developed channel flow at low Reynolds number

A direct numerical simulation of a turbulent channel flow is performed. The unsteady Navier-Stokes equations are solved numerically at a Reynolds number of 3300, based on the mean centerline velocity and channel half-width, with about 4 million grid points. All essential turbulence scales are resolved on the computational grid and no subgrid model is used. A large number of turbulence statistics are computed and compared with the existing experimental data at comparable Reynolds numbers. Agreements as well as discrepancies are discussed in detail. Particular attention is given to the behavior of turbulence correlations near the wall. A number of statistical correlations which are complementary to the existing experimental data are reported for the first time.

Kim, John↗

Explicit upwind algorithm for the parabolized Navier-Stokes equations

A new explicit upwind algorithm based on Roe's flux-difference splitting (FDS) method has been developed for the three-dimensional Parabolized Navier-Stokes (PNS) equations. For three-dimensional flows, FDS's are determined separately for the two nonmarching directions and modified to account for the calculated shock angle in the crossflow plane. Second-order FDS is applied to the pressure and convection terms with the streamwise pressure gradient limited in the subsonic region to maintain a hyperbolic inviscid equation set. Second-order central differencing is obtained in the two-step algorithm for the shear and heat flux terms. The new algorithm is demonstrated for three laminar flow test cases: supersonic flow over a flat plate, hypersonic flow over a 15 deg ramp, and hypersonic flow past a 10 deg cone at a 24 deg angle of attack. The computed results agree well with experimental measurements.

Korte, John J.↗

Numerical simulation of a subsonic jet in a crossflow

The aerodynamic/propulsive interaction between a subsonic jet exhausting perpendicularly through a flat plate into a crossflow is investigated numerically using an approximately factored, partially flux-split, implicit solver for the three-dimensional, thin-layer Navier-Stokes equations. This algorithm is applied to flows with a range of jet-to-crossflow velocity ratios between 4 and 8. The computations model the jet trajectory, the contrarotating vortex pair and the wake region near the plate downstream of the jet orifice. Both qualitative and quantitative agreement with the existing experimental database are demonstrated. Flow visualization is instructive for understanding the physics of this flowfield.

Roth, Karlin R.↗

Applied transonics at Grumman

A review of several applications of Computational Fluid Dynamics (CFD) to various aspects of aerodynamic design recently carried out at Grumman is presented. The emphasis is placed on project-oriented applications where the ease of use of the methods and short start-to-completion times are required. Applications cover transonic wing design/optimization, wing mounted stores load prediction, transonic buffet alleviation, fuselage loads estimation, and compact offset diffuser design for advanced aircraft configurations. Computational methods employed include extended transonic small disturbance (automatic grid embedding) formulation for analysis/design/optimization and a thin layer Navier-Stokes formulation for both external and internal flow analyses.

Davis, W. H.↗