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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Experiments planned specifically for developing turbulence models in computations of flow fields around aerodynamic shapes

Building block experiments and companion numerical simulations intended to verify and guide turbulence modeling are described. A series of experiments and computations being used to enhance modeling development for the shock wave turbulent boundary layer interaction problem is emphasized. Results are given for transonic flow over a circular arc airfoil undergoing shock wave induced, boundary layer separation for supersonic flow along a tube wall undergoing normal shock wave induced, boundary layer separation. Experimental data which use the complete Navier-Stokes equations are discussed.

Marvin, J. G.↗

NASA NTRS1976

Atmospheric flow over two-dimensional bluff surface obstructions

The phenomenon of atmospheric flow over a two-dimensional surface obstruction, such as a building (modeled as a rectangular block, a fence or a forward-facing step), is analyzed by three methods: (1) an inviscid free streamline approach, (2) a turbulent boundary layer approach using an eddy viscosity turbulence model and a horizontal pressure gradient determined by the inviscid model, and (3) an approach using the full Navier-Stokes equations with three turbulence models; i.e., an eddy viscosity model, a turbulence kinetic-energy model and a two-equation model with an additional transport equation for the turbulence length scale. A comparison of the performance of the different turbulence models is given, indicating that only the two-equation model adequately accounts for the convective character of turbulence. Turbulence flow property predictions obtained from the turbulence kinetic-energy model with prescribed length scale are only insignificantly better than those obtained from the eddy viscosity model. A parametric study includes the effects of the variation of the characteristics parameters of the assumed logarithmic approach velocity profile. For the case of the forward-facing step, it is shown that in the downstream flow region an increase of the surface roughness gives rise to higher turbulence levels in the shear layer originating from the step corner.

Bitte, J.↗

NASA NTRS1976

Development of relaxation turbulence models

Relaxation turbulence models have been intensively studied. The complete time dependent mass averaged Navier-Stokes equations have been solved for flow into a two dimensional compression corner. A new numerical scheme has been incorporated into the developed computed code with an attendant order of magnitude reduction in computation time. Computed solutions are compared with experimental measurements of Law for supersonic flow. Details of the relaxation process have been studied; several different relaxation models, including different relaxation processes and varying relaxation length, are tested and compared. Then a parametric study has been conducted in which both Reynolds number and wedge angle are varied. To assess effects of Reynolds number and wedge angle, the parametric study includes the comparison of computed separation location and upstream extent of pressure rise; numerical results are also compared with the measurements of surface pressure, skin friction and mean velocity field.

Hung, C. M.↗

NASA NTRS1977

A numerical investigation of the finite element method in compressible primitive variable Navier-Stokes flow

The results of a comprehensive numerical investigation of the basic capabilities of the finite element method (FEM) for numerical solution of compressible flow problems governed by the two-dimensional and axis-symmetric Navier-Stokes equations in primitive variables are presented. The strong and weak points of the method as a tool for computational fluid dynamics are considered. The relation of the linear element finite element method to finite difference methods (FDM) is explored. The calculation of free shear layer and separated flows over aircraft boattail afterbodies with plume simulators indicate the strongest assets of the method are its capabilities for reliable and accurate calculation employing variable grids which readily approximate complex geometry and capably adapt to the presence of diverse regions of large solution gradients without the necessity of domain transformation.

Cook, C. H.↗

NASA NTRS1977

Reynolds number effects on shock-wave turbulent boundary-layer interactions - A comparison of numerical and experimental results

An experiment is described that tests and guides computations of a shock-wave turbulent boundary-layer interaction flow over a 20-deg compression corner at Mach 2.85. Numerical solutions of the time-averaged Navier-Stokes equations for the entire flow field, employing various turbulence models, are compared with the data. Each model is critically evaluated by comparisons with the details of the experimental data. Experimental results for the extent of upstream pressure influence and separation location are compared with numerical predictions for a wide range of Reynolds numbers and shock-wave strengths.

Horstman, C. C.↗

NASA NTRS1977

Numerical solution of axisymmetric boattail flow fields with plume simulators

Turbulent separating flows over axisymmetric afterbody-boattail configurations with solid sting plume simulators are computed with a time-dependent finite-difference method to solve the compressible Navier-Stokes equations. The Reynolds stress terms are replaced with a two-layer eddy viscosity model including a relaxation formula to model the nonequilibrium effects of the separated flow. The mesh is alined with the boattail body through an analytic transformation which accommodates a wide variety of boattail geometries. Numerical results for a series of boattail geometries over a wide range of Reynolds number (140,000 to 140 million) are presented and, when possible, compared with experimental data or independent numerical results.

Holst, T. L.↗

NASA NTRS1977

Numerical investigation of turbulence models for shock separated boundary-layer flows

Numerical solutions of the Navier-Stokes equations for shock separated turbulent boundary-layer flows are presented. Several turbulence models are investigated and assessed by their ability to predict the physical phenomena associated with two extensively documented experiments. The experimental flows consist of shock-wave boundary-layer interactions in axisymmetric internal and external geometries at Mach numbers of 1.5 and 7, respectively. Algebraic and one-equation eddy viscosity models are used to describe the Reynolds shear stress. Calculated values of skin friction, wall pressure distribution, kinetic energy of turbulence, and heat transfer are compared with measurements.

Viegas, J. R.↗

NASA NTRS1976

A numerical study of flow phenomena within a vortex sink rate sensor

A numerical investigation of a viscous incompressible fluid with constant properties within a vortex sink rate sensor is undertaken. A computer program was developed, for obtaining the numerical solution of the Navier-Stokes equations, assuming laminar flow, having general prescribed inlet conditions and axisymmetric boundary conditions. Computational results for various viscous flows and assorted boundary conditions have been obtained. In the vortex chamber, the radial momentum is conserved, whereas the axial momentum is negligible. In the sink region, the radial and tangential velocities, as well as the tangential vorticities, have significantly high values. In the sink tube, the radial momentum is negligible. The radial momentum within the vortex chamber is converted to the axial momentum in the sink tube. A change in the flow rate had an appreciable effect within the sensor and particularly near the sink tube entrance. As the flow rate was increased, both the tangential velocity and tangential vorticity increased rapidly. At the higher flow rates, vortices were produced at the corner of the entrance section of the sink tube.

Patel, D. K.↗

NASA NTRS1976

Recent computation of viscous effects in transonic flow

The paper describes some improvements to a method of solving the time-dependent Reynolds-averaged Navier-Stokes equations for two-dimensional compressible flows for describing the coupling between viscous and inviscid regions in transonic flows. Mesh patching is implemented at the sublayer-boundary layer interfaces and the boundary layer-inviscid flow interface. Mesh stretching is geometric and is small in regions of large gradients and large in regions of small gradients, thus preserving second-order accuracy. Mesh lines are arranged so as to accommodate MacCormack's rapid solver for hyperbolic systems and a wider class of turbulent transport models. Results of computations of transonic flows over lifting airfoils are compared with experimental data and other solutions.

Deiwert, G. S.↗

NASA NTRS1977

Numerical solution of three-dimensional free turbulent shear flows

Numerical results are presented for three-dimensional compressible turbulent jet and wake flows. An alternating direction implicit numerical procedure is used to solve the finite-difference form of the parabolic elliptic Navier-Stokes equations. A coordinate transformation maps the boundaries at infinity into a finite computational domain in order to properly specify infinity boundary conditions as well as contain the downstream growth of the viscous flow field in a fixed computational grid. Turbulence closure is achieved through an algebraic mixing length eddy viscosity model. Numerical results for supersonic flow are presented for an axisymmetric jet, an elliptical jet, an elliptical wake, and two interacting rectangular jets. Experimental data were not available for comparison with the numerical results. However, the results compare well with empirical results for free shear flows.

Oh, Y. H.↗

NASA NTRS1977

An experimental and computational investigation of the steady and unsteady transonic flow fields about an airfoil in a solid-wall test channel

An experimental and computational investigation of the steady and unsteady transonic flow field about a thick airfoil is described. An operational computer code for solving the two-dimensional, compressible Navier-Stokes equations for flow over airfoils was modified to include solid-wall, slipflow boundary conditions to properly assess the code and help guide the development of improved turbulence models. Steady and unsteady flow fields about an 18% thick circular arc airfoil at Mach numbers of 0.720, 0.754, and 0.783 and a chord Reynolds number of 11 x 10 to the 6th are predicted and compared with experiment. For the first time, computed results for unsteady turbulent flows with separation caused by a shock wave were obtained which qualitatively reproduce the time-dependent aspects of experiments. Features such as the intensity and reduced frequency of airfoil surface-pressure fluctuations, oscillatory regions of trailing-edge and shock-induced separation, and the Mach number range for unsteady flows were all qualitatively reproduced.

Levy, L. L., Jr.↗

NASA NTRS1977

Numerical solution of supersonic laminar flow over a three-dimensional compression corner

A newly developed, rapid numerical scheme is extended to three dimensions to solve the complete Navier-Stokes equations for a supersonic, laminar flow over a compression corner with sidewall effects. The program is coded so that it can solve for a general curved ramp surface geometry such as found in inlets and fuselage-wing-flap junctions. A test case of Mach 3.0 flow is calculated. In regions where three-dimensional effects are small, good agreement is obtained between the present calculation and previous two-dimensional solutions. In other regions, the results show complex three-dimensional flow-field interactions including shock-shock and shock/boundary-layer interactions

Hung, C. M.↗

NASA NTRS1977

Polynomial interpolation methods for viscous flow calculations

Higher-order collocation procedures which result in block-tridiagonal matrix systems are derived from (1) Taylor series expansions and from (2) polynomial interpolation, and the relationships between the two formulations, called respectively Hermite and spline collocation, are investigated. A Hermite block-tridiagonal system for a nonuniform mesh is derived, and the Hermite approach is extended in order to develop a variable-mesh sixth-order block-tridiagonal procedure. It is shown that all results obtained by Hermite development can be recovered by appropriate spline polynomial interpolation. The additional boundary conditions required for these higher-order procedures are also given. Comparative solutions using second-order accurate finite difference and spline and Hermite formulations are presented for the boundary layer on a flat plate, boundary layers with uniform and variable mass transfer, and the viscous incompressible Navier-Stokes equations describing flow in a driven cavity.

Rubin, S. G.↗

NASA NTRS1977

Numerical solutions for viscous and potential flow about arbitrary two-dimensional bodies using body-fitted coordinate systems

A procedure for numerical solution of the time-dependent, incompressible Navier-Stokes equations for the flow about arbitrarily shaped two-dimensional bodies is given. This solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having a coordinate line coincident with the body contour regardless of its shape. The implicit solution utilizes the vorticity-stream function formulation with a false-position iterative adjustment of the surface vorticity in satisfaction of the no-slip boundary condition. Excellent agreement with the Blasius boundary layer solution is obtained for a semi-infinite flat plate. Results are presented for Reynolds numbers up to 2000 for several airfoils and a cambered rock.

Thames, F. C.↗

NASA NTRS1977

Computation of viscous transonic flow about a lifting airfoil

The Reynolds averaged Navier-Stokes equations are solved numerically for the viscous transonic flow about a stationary NACA 64A010 airfoil in free air. This paper presents descriptions of the numerical method, turbulence models employed, and boundary conditions appropriate to simulation of free-air flight. Computed results are presented for the airfoil at a free-stream Mach number of 0.8, angles of attack of 0 and 2 deg, and a Reynolds number based on a chord of 4 x 10 to the 6th. For the lifting case, unsteady periodic motion was calculated along the aft portion of the airfoil and in its wake. Recent experimental results obtained by Johnson indicate periodicity aft of the shock closely approximates the computed frequency, but the amplitude of the disturbances was significantly less than the calculated amplitude.

Walitt, L.↗

NASA NTRS1977

Neutrally stable atmospheric flow over a two-dimensional rectangular block

The phenomena of atmospheric flow over a two dimensional surface obstruction such as a building modeled as a rectangular block are analyzed by an approach using the Navier-Stokes equations with a two equation model of turbulence. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale are solved by a finite difference technique. The predicted results are in agreement with the limited experimental data available. Current computed results show that the separation bubble originates from the upper front corner of the block and extends approximately 11.5 block heights behind the block. The decay of the mean velocity along the wake center line coincides almost perfectly with the experimental data. The vertical profiles of the mean velocity defect are also in reasonable agreement with wind tunnel results. Velocity profiles in the mixing region are shown to agree with the error function profile typically found in the shear layer. Details of the behavior of the turbulence kinetic energy and the turbulence length scale are also discussed.

Shieh, C. F.↗

NASA NTRS1977

Computations of supersonic viscous flow over a finite-width plate

Finite difference methods were applied to solve the parabolic Navier-Stokes equations for the flow over a finite width plate at 0 deg and 10 deg angles of attack. The methods were developed on the basis of the operator factorization concept resulting in the split of a three dimensional equation into successive two dimensional equations. Backward and centered implicit factorization schemes, were used and their results were compared. Available numerical solutions and experimental data obtained at low Reynolds number conditions were also used for comparison. The backward implicit method provides a more successful solution, which ranges from the merged layer to the strong interaction regimes. Detailed structures were revealed of the shear layer around and behind the side edge.

Li, C. P.↗

NASA NTRS1978

Computational aerodynamics and the numerical aerodynamic simulation facility

Technical and economic reasons for accelerating the maturation of the discipline of computational aerodynamics include the cost of conducting the experiments required to provide the empirical data base for new aeronautical vehicles and the limitations in test facilities (Reynolds number, wall and support interferences, aeroelastic distortions, real-gas effects, etc.) for simulating the full-scale vehicle environment. General purpose computers do not have the necessary capability for the next stage of development. Solution of the three dimensional Reynolds averaged Naiver-Stokes equations in a short time to be practical for design purposes will require 40 times the power of current supercomputers. However, it is feasible to construct a special purpose processor that will meet these requirements to enhance the nation's aerodynamic design capability in the 1980's.

Peterson, V. L.↗