Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes equations”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗
NASA NTRS1975

Numerical studies of incompressible viscous flow in a driven cavity

A series of project papers is presented in computational fluid dynamics. The work was performed during the 1973-74 academic year at Old Dominion University. Each paper briefly examines a numerical method(s) that can be applied to the Navier-Stokes equations governing incompressible flow in a driven cavity. Solutions obtained with a cubic spline procedure are also included.

Source record↗
NASA NTRS1975

Numerical calculation of the internal flow field in a centrifugal compressor impeller

An iterative numerical method has been developed for the calculation of steady, three-dimensional, viscous, compressible flow fields in centrifugal compressor impellers. The computer code, which embodies the method, solves the steady three dimensional, compressible Navier-Stokes equations in rotating, curvilinear coordinates. The solution takes place on blade-to-blade surfaces of revolution which move from the hub to the shroud during each iteration.

Walitt, L.↗
NASA NTRS1976

Comments on the role of diagonal dominance in implicit difference methods

Numerical tests were made for a model of the Navier-Stokes equations using a second-order accurate implicit scheme which guarantees diagonal dominance. The results suggest that the failure of implicit methods using large marching steps may not always be attributed to the lack of diagonal dominance in the coefficient matrix. In some cases the failure may be caused by a nonlinear instability associated with the solution method.

Rudy, D. H.↗
NASA NTRS1976

An experimental and computational investigation of the flow field about a transonic airfoil in supercritical flow with turbulent boundary-layer separation

A combined experimental and computational research program is described for testing and guiding turbulence modeling within regions of separation induced by shock waves incident in turbulent boundary layers. Specifically, studies are made of the separated flow the rear portion of an 18%-thick circular-arc airfoil at zero angle of attack in high Reynolds number supercritical flow. The measurements include distributions of surface static pressure and local skin friction. The instruments employed include highfrequency response pressure cells and a large array of surface hot-wire skin-friction gages. Computations at the experimental flow conditions are made using time-dependent solutions of ensemble-averaged Navier-Stokes equations, plus additional equations for the turbulence modeling.

Rubesin, M. W.↗
NASA NTRS1976

Numerical simulation of turbulent jet noise, part 2

Results on the numerical simulation of jet flow fields were used to study the radiated sound field, and in addition, to extend and test the capabilities of the turbulent jet simulation codes. The principal result of the investigation was the computation of the radiated sound field from a turbulent jet. In addition, the computer codes were extended to account for the effects of compressibility and eddy viscosity, and the treatment of the nonlinear terms of the Navier-Stokes equations was modified so that they can be computed in a semi-implicit way. A summary of the flow model and a description of the numerical methods used for its solution are presented. Calculations of the radiated sound field are reported. In addition, the extensions that were made to the fundamental dynamical codes are described. Finally, the current state-of-the-art for computer simulation of turbulent jet noise is summarized.

Metcalfe, R. W.↗
NASA NTRS1976

A normal shock-wave turbulent boundary-layer interaction at transonic speeds

Experimental results, including surveys of the mean and fluctuating flow, and measurements of surface pressure, skin friction, and separation length, are compared with solutions to the Navier-Stokes equations utilizing various algebraic eddy viscosity models to describe the Reynolds shear stresses. The experimental data, obtained at a free-stream Mach number of 1.5 and Reynolds numbers between 10 million and 80 million, show that a separated zone forms near the foot of the shock and that its length is proportional to the initial boundary-layer thickness; that a supersonic region forms downstream of the shock; and that the shear stress increases significantly through the interaction and subsequently decays downstream. The computations adequately represent the qualitative features of the flow field throughout the interaction but quantitatively underpredict the extent of separation and the downstream level of skin friction.

Mateer, G. G.↗
NASA NTRS1976

Comparison of a two-dimensional shock impingement computation with experiment

Results of computations of two-dimensional viscous blunt-body flowfields with an impinging shock wave, with a time-dependent finite-difference method employed to solve the complete set of Navier-Stokes equations, are compared with experimental results. The experimental results were obtained in a 20-inch hypersonic tunnel with a planar shock impinging on the cylindrical leading edge of a fin, hence with the shock parallel to the centerline of the leading edge, so that type III and type IV interference patterns were generated. Close agreement is found. The overall effects of smoothing and grid size on the calculations are determined. A 31 x 51 mesh is adequate for wall pressure values (except in peaked regions).

Tannehill, J. C.↗
NASA NTRS1976

On the flow in an annulus surrounding a whirling cylinder

When fluid in an annulus between two cylinders is set in motion by whirling movements of one or both of the cylinders, dynamic forces are imposed by the fluid on the cylinders. Knowledge of these forces is frequently important, indeed often critical, to the engineer designing rotor systems or journal bearings. Quite general solutions of the Navier-Stokes equations are presented for this problem and are limited only by restrictions on the amplitude of the whirl motion. From these solutions, the forces are derived under a wide variety of circumstances, including large and small annular widths, high and low Reynolds numbers, and the presence and absence of a mean flow created by additional net rotation of one or both of the cylinders.

Brennen, C.↗
NASA NTRS1976

Modifications of the law of the wall and algebraic turbulence modelling for separated boundary layers

Various modifications of the conventional algebraic eddy viscosity turbulence model are investigated for application to separated flows. Friction velocity is defined in a way that avoids singular behavior at separation and reattachment but reverts to the conventional definition for flows with small pressure gradients. This leads to a modified law of the wall for separated flows. The effect on the calculated flow field of changes in the model that affect the eddy viscosity at various distances from the wall are determined by (1) switching from Prandtl's form to an inner layer formula due to Clauser at various distances from the wall, (2) varying the constant in the Van Driest damping factor, (3) using Clauser's inner layer formula all the way to the wall, and (4) applying a relaxation procedure in the evaluation of the constant in Clauser's inner layer formula. Numerical solutions of the compressible Navier-Stokes equations are used to determine the effects of the modifications. Experimental results from shock-induced separated flows at Mach numbers 2.93 and 8.45 are used for comparison. For these cases improved predictions of wall pressure distribution and positions of separation and reattachment are obtained from the relaxation version of the Clauser inner layer eddy viscosity formula.

Baldwin, B. S.↗
NASA NTRS1976

Analysis of neutrally stable atmospheric flow over a two-dimensional forward facing step

A mathematical model describing ground-wind induced flow fields around surface obstructions such as buildings, bridges or other man-made structures is presented. The Navier-Stokes equations with a two-equation turbulence model are used to analyze the flow over a two-dimensional forward facing step. A study of the assumed logarithmic approach velocity profile shows that an increase of surface roughness produces a rise to higher turbulence levels in the shear layer originating from the step corner. This rise leads to higher Reynolds stress and to faster reattachment of the separated flow.

Bitte, J.↗
NASA NTRS1976

Numerical simulation of supersonic and hypersonic turbulent compression corner flows using relaxation models

Relaxation turbulence eddy viscosity models are incorporated to solve the complete Navier-Stokes equations for supersonic and hypersonic flows over a two-dimensional compression corner. The system of equations is solved by a time-split, second-order accurate numerical scheme. Details of relaxation process are studied, and several relaxation models are tested and compared. Good improvement in the prediction of upstream pressure propagation is obtained for a Mach number of 2.96, and Reynolds number of 10,000,000, with wedge angle of 25 deg. However, the application of the relaxation models to hypersonic flow at Mach number 8.66 and Reynolds number of 22,000,000, with highly cooled wall, shows unfavorable effects on heat transfer and skin friction in the reattachment and recompression regions.

Hung, C. M.↗
NASA NTRS1976

Circularly polarized inertial wave vectors in rotating fluids

The Navier-Stokes equations for a rotating fluid are harmonically analyzed for planar motion in an infinite half-space. All solutions are shown to be a sum of two inertial wave vectors, one circularly polarized to the left (CPL) and the other circularly polarized to the right (CPR). These basic solutions are therefore presented in the same nomenclature and form as that found useful by experimentalists in analyzing flow data (called 'rotary spectra'). The CPL wave acts counter to the Coriolis force and consequently has a slower phase speed and larger damping than the CPR wave. At resonance (forcing frequency = Coriolis frequency) the CPR wave has an infinite phase speed and no damping and is the important component leading to the singular nature of the solutions for certain boundary conditions. All possible resonant singularities are explicitly shown. The unsteady development of these unbounded (limited space structure), cyclic (no time origin or structure) flows is presented to show that with time structure the resonant singularities evolve in a self-similar manner.

Cogley, A. C.↗
NASA NTRS1976

Repeated-cascade theory of strong turbulence in a magnetized plasma

A two-dimensional Navier-Stokes equation of vorticity in fluid turbulence is used to model drift turbulence in a plasma with a strong constant magnetic field and a constant mean density gradient. The nonlinear eddy diffusivity is described by a time-integrated Lagrangian correlation of velocities, and the repeated-cascade method is employed to choose the rank accounting for nearest-neighbor interactions, to calculate the Lagrangian correlation, and to close the correlation hierarchy. As a result, the diffusivity becomes dependent on the plasma's induced diffusion and is represented by a memory chain that is cut off by similarity and inertial randomization. Spectral laws relating the kinetic-energy spectrum to the -5, -5/2, -3, and -11 powers of wavenumber are derived for the velocity subranges of production, approach to inertia, inertia, and dissipation, respectively. It is found that the diffusivity is proportional to some inverse power of the magnetic field, that power being 1, 2/3, 5/6, and 2, respectively, for the four velocity subranges.

Tchen, C. M.↗
NASA NTRS1976

An experimental and computational investigation of the flow field about a transonic airfoil in supercritical flow with turbulent boundary-layer separation

A combined experimental and computational research program for testing and guiding turbulence modeling within regions of separation induced by shock waves incident on turbulent boundary layers is described. Specifically, studies are made of the separated flow over the rear portion of an 18%-thick circular-arc airfoil at zero angle of attack in high Reynolds number supercritical flow. The measurements include distributions of surface static pressure and local skin friction. The instruments employed include high-frequency response pressure cells and a large array of surface hot-wire skin-friction gages. Computations at the experimental flow conditions are made using time-dependent solutions of ensemble-averaged Navier-Stokes equations, plus additional equations for the turbulence modeling.

Rubesin, M. W.↗
NASA NTRS1976

Numerical solutions for laminar and turbulent viscous flow over single and multi-element airfoils using body-fitted coordinate systems

The technique of body-fitted coordinate systems is applied in numerical solutions of the complete time-dependent compressible and incompressible Navier-Stokes equations for laminar flow and to the time-dependent mean turbulent equations closed by modified Kolmogorov hypotheses for turbulent flow. Coordinate lines are automatically concentrated near to the bodies at higher Reynolds number so that accurate resolution of the large gradients near the solid boundaries is achieved. Two-dimensional bodies of arbitrary shapes are treated, the body contour(s) being simply input to the program. The complication of the body shape is thus removed from the problem.

Thompson, J. F.↗
NASA NTRS1976

The lift force on a drop in unbounded plane Poiseuille flow

The lift force on a deformable liquid sphere moving in steady, plane Poiseuille-Stokes flow and subjected to an external body force is calculated. The results are obtained by seeking a solution to Stokes' equations for the motion of the liquids inside and outside the slightly perturbed sphere surface, as expansions valid for small values of the ratio of the Weber number to the Reynolds number. When the ratio of the drop and external fluid viscosities is small, the lift exerted on a neutrally buoyant drop is found to be approximately one-tenth of the magnitude of the force reported by Wohl and Rubinow acting on the same drop in unbounded Poiseuille flow in a tube. The resultant trajectory of the drop is calculated and displayed as a function of the external body force.

Wohl, P. R.↗
NASA NTRS1976

The inviscid pressure field on the tip of a semi-infinite wing and its application to the formation of a tip vortex

A method was developed for determining the aerodynamic loads on the tip of an infinitely thin, swept, cambered semi-infinite wing at an angle of attack which is operating subsonically in an inviscid medium and is subjected to a sinusoidal gust. Under the assumption of linearized aerodynamics, the loads on the tip are obtained by superposition of the steady aerodynamic results for angle of attack and camber, and the unsteady results for the response to the sinusoidal gust. The near field disturbance pressures in the fluid surrounding the tip are obtained by assuming a dipole representation for the loading on the tip and calculating the pressures accordingly. The near field pressures are used to drive a reduced form of the Navier-Stokes equations which yield the tip vortex formation. The combined viscid-inviscid analysis is applied to determining the pressures and examining the vortex rollup in the vicinity of an unswept, uncambered wing moving steadily at a Mach number of 0.2 at an angle of attack of 0.1 rad. The viscous tip flow calculation shows features expected in the tip flow such as the qualitatively proper development of boundary layers on both the upper and lower airfoil surfaces. In addition, application of the viscous solution leads to the generation of a circular type flow pattern above the airfoil suction surface.

Hall, G. F.↗