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,171 records · Page 65

Similitude requirements for hypersonic, rarefied, nonequilibrium flow

Similitude requirements for hypersonic, rarefied flow with nonequilibrium chemistry and vibration are presented. The full Navier-Stokes equations with catalytic or noncatalytic walls and with or without slip conditions are nondimensionalized. The heat transfer coefficient is written in terms of fourteen dimensionless parameters and reduced to four by making the binary scaling assumption. Duplication of blunt and sharp nose heat transfer requires the use of air over a geometrically similar model with the same free stream velocity, wall temperature and product of free stream density and characteristic length. Estimates of this heat transfer coefficient are also presented.

Hendricks, W. L.↗

Separated flow over bodies of revolution using an unsteady discrete-vorticity cross wake. Part 2: Computer program description

A method is developed to determine the flow field of a body of revolution in separated flow. The computer was used to integrate various solutions and solution properties of the sub-flow fields which made up the entire flow field without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the two-dimensional unsteady separated flow problem based upon an unsteady, discrete-vorticity wake. Data for the forces and moments on aerodynamic bodies at low speeds and high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data. In addition, three dimensional steady separated regions and wake vortex patterns are determined. The computer program developed to perform the numerical calculations is described.

Marshall, F. J.↗

Separated flow over bodies of revolution using an unsteady discrete-vorticity cross wake. Part 1: Theory and application

A method is developed to determine the flow field of a body of revolution in separated flow. The technique employed is the use of the computer to integrate various solutions and solution properties of the sub-flow fields which made up the entire flow field without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the required two-dimensional unsteady separated flow problem based upon an unsteady, discrete-vorticity wake. Data for the forces and moments on aerodynamic bodies at low speeds and high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data. In addition, three dimensional steady separation regions and wake vortex patterns are determined.

Marshall, F. J.↗

Convection and magnetic fields in an atmosphere with constant temperature gradient. I - Hydrodynamic flows

In order to deduce the fluid and field conditions deep in the supergranules from the fluid motions and magnetic fields observed in the solar photosphere, it is necessary to have some idea of the depth to which the supergranules extend and the form of their internal fluid velocities. Toward this end a number of analytic solutions of the Navier-Stokes equations are given for circumstances possessing various facets of the supergranule environment in the sun. These solutions, together with the earlier work by others, shows that the granules extend to depths in excess of their mean diameters, that the convective velocities probably do not diminish much, if any, with depth, and that across the tops of the supergranules at the photosphere the horizontal velocity is more or less independent of depth.

Parker, E. N.↗

A finite-difference solution for unsteady wave interactions with an application to the expansion tube

An analytical study of the unsteady wave interaction of a shock wave overtaking a previously initiated shock tube flow has been conducted for a real gas by solving a finite-difference representation of the one-dimensional Navier-Stokes equations. This analysis has been used to determine the effect on the available test time of opening the secondary diaphragm in the expansion-tube operating cycle prior to the arrival of the incident shock wave. Special emphasis is placed on determining expansion-tube test-time limitations resulting from a nonuniform entropy distribution occurring during the early test-gas flow due to the varying Mach number of the incident shock wave as it traverses the expansion fan caused by the preopening process. Time distributions of velocity, temperature, pressure, and entropy at various distances from the diaphragm are presented as a function of the time required to open the diaphragm.

Weilmuenster, K. J.↗

Analysis of three-dimensional unsteady flow around oscillating wings

A method based on the Navier-Stokes equations was developed for determining analytically the three-dimensional unsteady flow patterns around oscillating wings. The Helmholz vorticity transport equations were discretized in three-dimensional finite element form from a variational formulation and integrated numerically. At each time step of the numerical integration the velocity field was calculated from the representation of the three-dimensional wing by a system of optimized distribution of vortices in space. During the numerical integration of the vorticity transport equations the time-dependent boundary conditions on the wing were specified as external constraint conditions. Examples of obtained results describing the three-dimensional unsteady flow around a wing were presented.

Bratanow, T.↗

Finite element solution algorithm for incompressible fluid dynamics

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the transient motion of a viscous incompressible fluid, i.e., hydrodynamics. Dependent variable transformation renders the differential equation description uniformly elliptic. The finite element algorithm is established using the Galerkin criterion on a local basis within the Method of Weighted Residuals. It is unconstrained with respect to system linearity, computational mesh uniformity or solution domain closure regularity. The finite element matrices are established using a linear 'natural coordinate function' description. Computational solutions using the COMOC computer program illustrate the various features of the algorithm including recirculating flows.

Baker, A. J.↗

An experimental and analytical study of slip and catalytic boundary conditions applied to spheres in low Reynolds number arc jet flows

Measurements of low-density dissociated arc-heated flows around low and high surface catalycity spheres are compared with numerical solutions to the full Navier-Stokes equations for reacting gas flow fields, taking into account slip and finite surface catalytic boundary conditions. A time-dependent finite difference solution technique developed by Li (1973, 1974) is used in the investigation and a simple transport model is employed for the flows considered. The correctness of the solution approach is demonstrated by the very good agreement obtained between experimental and computational data.

Scott, C. D.↗

The Boltzmann- and the Klimontovich formalisms with reference to turbulence resolution

The Navier-Stokes equation expressed in terms of Klimontovich microscopic density is compared with the conventional one based on the Boltzmann formalism. Their difference is described by the two-point BBGKY equation. An implication is given that the averaging procedure from the Klimontovich to the Boltzmann formalism wipes out not only thermal agitation, but also a certain hydrodynamic turbulence. This is substantiated by agreement with measurements due to Schubauer and Klebanoff.

Tsuge, S.↗

Non-Newtonian effects on flow-generated cavitation and on cavitation in a pressure field

Observations are presented which show that the stresses in a flow field of very dilute polymer are not well enough described by the Navier-Stokes equations to accurately predict cavitation. The contitutive equation for the particular polymer and concentration used is needed. The second-order fluid form in which accelerations are relatively important appears capable of explaining observed cavitation suppression by changing the pressure field due to flow. Bubble dynamics in stationary dilute polymer solutions are also examined and found to be little different from those in water.

Ellis, A. T.↗

Numerical and experimental studies of the natural convection within a horizontal cylinder

Numerical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time-dependent natural convection within a horizontal cylinder. The early flow development and wall heat transfer are obtained after a uniformly cold wall is imposed as a boundary condition on the cylinder. Results are also obtained for a time-varying cold wall as a boundary condition with windward explicit differencing used for the numerical solutions. The viscous truncation error associated with this scheme is controlled so that first-order accuracy is maintained in time and space. Experiments within a small-scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer as well as the decay of wall heat transfer with time. Agreement between temperature distributions obtained experimentally and numerically was generally good. The time decay of the dimensionless ratio of the Nusselt number to the one-fourth power of the Grashof number is found both numerically and experimentally, and good agreement is obtained between these two results over most of the cylinder wall.

Stewart, R. B.↗

Two dimensional turbulence in inviscid fluids or guiding center plasmas

Analytic theory for two-dimensional turbulent equilibria for the inviscid Navier-Stokes equations is examined mathematically. Application of the technique to electrostatic guiding center plasma is discussed. A good fit is demonstrated for the approach to a predicted energy per Fourier mode obtained from a two-temperature canonical ensemble. Negative as well as positive temperature regimes are explored. Fluctuations about the mean energy per mode also compare well with theory. In the regime of alpha less than zero, beta greater than zero, with the minimum value of alpha plus beta times k squared near zero, contour plots of the stream function reveal macroscopic vortex structures similar to those seen previously in discrete vortex simulations. Eulerian direct interaction equations, which can be used to follow the approach to inviscid equilibrium, are derived.

Seyler, C. E., Jr.↗

Turbulence coefficients and stability studies for the coaxial flow or dissimiliar fluids

Analytical investigations of fluid dynamics problems of relevance to the gaseous core nuclear reactor program are presented. The vortex type flow which appears in the nuclear light bulb concept is analyzed along with the fluid flow in the fuel inlet region for the coaxial flow gaseous core nuclear reactor concept. The development of numerical methods for the solution of the Navier-Stokes equations for appropriate geometries is extended to the case of rotating flows and almost completes the gas core program requirements in this area. The investigations demonstrate that the conceptual design of the coaxial flow reactor needs further development.

Weinstein, H.↗

Investigation A: Numerical solutions of driven vortices of binary fluid in a cylindrical geometry

The steady state laminar motion of a viscous, incompressible and binary fluid is studied for a rotating flow in a cylindrical geometry. The mathematical model employed is a cold flow simulation of the fluid mechanics of the light-bulb concept of the gaseous core nuclear engine. A numerical treatment is developed for the rotating flow which includes a description of the nuclear fuel addition. The problem is formulated with the complete Navier-Stokes equations in order to show the interaction between the fuel addition, the main flow, and the boundary layer flow in an accurate manner. The results presented show holdup of the nuclear fuel for the case of steady fuel addition.

Shin, Y. W.↗

Numerical computation of viscous blunt body flows with a planar impinging shock

Two- and three-dimensional, viscous blunt body flows with planar impinging shocks are computed using an explicit, time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which planar shock waves of different strengths and orientations are allowed to impinge on the flow field surroundings an infinite cylinder resulting in two- and three-dimensional shock interference patterns. The numerical results are compared with experiment.

Holst, T. L.↗

Computational aerodynamics

The application of computer techniques for solving Navier-Stokes equations in support of wind tunnel tests is discussed. The ILLIAC IV computer is considered for this purpose and its limitations are analyzed. The author states that improved computers will make it possible to solve many aerodynamic problems and reduce the amount of wind tunnel testing required for adequate data processing.

Chapman, D. R.↗

Numerical computation of three-dimensional blunt body flow fields with an impinging shock

A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.

Holst, T. L.↗

Numerical simulation of turbulence in the presence of shear

The numerical calculations are presented of the large eddy structure of turbulent flows, by use of the averaged Navier-Stokes equations, where averages are taken over spatial regions small compared to the size of the computational grid. The subgrid components of motion are modeled by a local eddy-viscosity model. A new finite-difference scheme is proposed to represent the nonlinear average advective term which has fourth-order accuracy. This scheme exhibits several advantages over existing schemes with regard to the following: (1) the scheme is compact as it extends only one point away in each direction from the point to which it is applied; (2) it gives better resolution for high wave-number waves in the solution of Poisson equation, and (3) it reduces programming complexity and computation time. Examples worked out in detail are the decay of isotropic turbulence, homogeneous turbulent shear flow, and homogeneous turbulent shear flow with system rotation.

Shaanan, S.↗