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,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1971

Effect of buoyancy on fuel containment in an open-cycle gas-core nuclear rocket engine.

Analysis aimed at determining the scaling laws for the buoyancy effect on fuel containment in an open-cycle gas-core nuclear rocket engine, so conducted that experimental conditions can be related to engine conditions. The fuel volume fraction in a short coaxial flow cavity is calculated with a programmed numerical solution of the steady Navier-Stokes equations for isothermal, variable density fluid mixing. A dimensionless parameter B, called the Buoyancy number, was found to correlate the fuel volume fraction for large accelerations and various density ratios. This parameter has the value B = 0 for zero acceleration, and B = 350 for typical engine conditions.

Putre, H. A.↗

NASA NTRS1971

The three-dimensional steady radial expansion of a viscous gas from a sonic source into a vacuum.

The three-dimensional steady radial expansion of a viscous, heat-conducting, compressible fluid from a spherical sonic source into a vacuum is analyzed using the Navier-Stokes equations as a basis. It is assumed that the model fluid is a perfect gas having constant specific heats, a constant Prandtl number of order unity, and viscosity coefficients varying as a power of the absolute temperature. Limiting forms for the flow variable solutions are studied for the Reynolds number based on the sonic source conditions, going to infinity and the Newtonian parameter both fixed and going to zero. For the case of the viscosity-temperature exponent between .5 and 1, it is shown that the velocity as well as the pressure approach zero as the radial distance goes to infinity. The formulations of the distinct regions that span the domain extending from the sonic source to the vacuum are presented.

Bush, W. B.↗

NASA NTRS1971

Prediction of flow profiles in arteries from local measurements.

This paper develops an approximate numerical method for calculating flow profiles in arteries. The theory takes into account the nonlinear terms of the Navier-Stokes equations as well as the large deformations of the arterial wall. The method, assuming axially symmetric flow, determines velocity distribution and wall shear at a given location from the locally measured values of the pressure, pressure gradient, and pressure-radius relation. The computed results agree well with the corresponding experimental data.

Ling, S. C.↗

NASA NTRS1972

Two-dimensional vortex motion and 'negative temperatures.'

Explanation of the novel phenomenon, tentatively identified as the 'ergodic boundary' in a space of initial conditions for turbulent flow, suggested by the recent numerical integration of the two-dimensional Navier-Stokes equations at high Reynolds numbers reported by Deem and Zabusky (1971). The proposed explanation is presented in terms of negative temperatures for a point vortex model.

Montgomery, D.↗

NASA NTRS1972

Numerical study of strong plasma shock waves produced in an electromagnetic shock tube.

A set of two-fluid Navier-Stokes equations with classical physical transport coefficients is used to compute the evolution and structure of collisional plasma shock waves in an electromagnetic shock tube. Shock speeds up to 200 cm/microsec and shocked-plasma temperatures of the order of kilovolts are studied. A strong transverse bias magnetic field is employed, which significantly alters the size and shape of the shock profiles, when compared with the zero bias field case. The wave structure is different from that of a two-fluid gasdynamic shock without transverse magnetic field. Especially significant is the effect of the small ion Larmor radius in reducing the gasdynamic shock thickness by at least an order of magnitude in the transverse shock case, which permits collisional shocks to have thicknesses much smaller than the post-shock mean free path. These collisional shock waves produce very hot ions, the ion temperature increasing with shock speed, although the formation time and distances also increase substantially with shock speed.

Schneider, S. H.↗

NASA NTRS1972

Simulation of the 'negative temperature' instability for line vortices.

In previous numerical solution to the continuum Navier-Stokes equations, a 'negative temperature' instability for the two-dimensional motions of interacting line vortices was observed. The experiment is repeated for a discrete vortex model, thus obtaining a numerical simulation of the 'negative temperature' instability for a large number of discrete line vortices. Typical results which are shown, are thought to lie above and below the energy threshold for negative temperature instability.

Joyce, G.↗

NASA NTRS1972

Numerical and Experimental Studies of the Natural Convection Flow Within a Horizontal Cylinder Subjected to a Uniformly Cold Wall Boundary Condition

Numberical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time dependent natural convection flow within a horizontal cylinder. The early time flow development and wall heat transfer is obtained after imposing a uniformly cold wall boundary condition on the cylinder. Solutions are also obtained for the case of a time varying cold wall boundary condition. Windware explicit differ-encing is 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. The results encompass a range of Grashof numbers from 8.34 times 10,000 to 7 times 10 to the 7th power which is within the laminar flow regime for gravitationally driven fluid flows. Experiments within a small scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer and also the decay of wall heat transfer with time.

Stewart, R. B.↗

NASA NTRS1973

Nonlinear evolution of a disturbance in an unbounded viscous fluid with uniform shear

The evolution of a disturbance in the presence of a uniform mean velocity gradient is calculated by a power-series solution of the incompressible Navier-Stokes equations. Terms through those in time cubed are retained in the solution. For the initial condition a three-dimensional cosine distribution with two harmonic terms is assumed. The nonlinear interaction of these harmonic terms produces new harmonics which in turn interact. For large velocity gradients the energy of the disturbance can grow with time. The results shed some light on the maintenance or growth of turbulence in a shear flow.

Deissler, R. G.↗

NASA NTRS1972

A nonlinear analysis of pulsatile flow in arteries.

An approximate numerical method for calculating flow profiles in arteries is developed. The theory takes into account the nonlinear terms of the Navier-Stokes equations as well as the nonlinear behaviour and large deformations of the arterial wall. Through the locally measured values of the pressure, pressure gradient, and pressure-radius function, the velocity distribution and wall shear at a given location along the artery can be determined. The computed results agree well with the corresponding experimental data.

Ling, S. C.↗

NASA NTRS1972

On the numerical solution of two-dimensional, laminar compressible flows with imbedded shock waves.

The complete, time-dependent Navier-Stokes equations are expressed in conservation form and solved by employing an explicit finite difference numerical technique which incorporates artificial viscosity terms of the form first suggested by Rusanov for numerical stability in the vicinity of shock waves. Surface boundary conditions are developed in a consistent and unique manner through the use of a physically oriented extrapolation procedure. From numerical experimentation an extended range for the explicit stability parameter is established. Also employed is an additional convergence parameter which relates incremental spatial steps. Convergence of the transient solution to a steady state flow was obtained after 400 to 500 time steps.

Goodrich, W. D.↗

NASA NTRS1973

Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.

An analytical procedure based on Navier-Stokes equations was developed for representing unsteady flow patterns around oscillating obstacles. A variational formulation of the Helmholtz vorticity equation was discretized in finite element form and integrated numerically. At each step of the numerical integration the velocity field around the obstacle was determined from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints at each time step of the numerical integration. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗

NASA NTRS1973

Laminar film boiling on inclined isothermal flat plates.

Laminar film boiling from an inclined flat plate has been investigated analytically. Using the singular perturbation scheme, the complete set of Navier-Stokes equations is solved. The zeroth-order perturbation coinciding with the boundary-layer equations for vertical flat plates governs the problem. The higher-order perturbations become important near the leading edge and for large values of the inclination angle from the vertical. The assumption of zero interfacial velocity shows that, except for fluids having large (rho x mu) ratios, the results can be predicted using the vertical flat plate results by defining a modified Grashof parameter containing a cos phi term. When the interfacial shear is considered, the solutions indicate that for fluids having large (rho x mu) ratios, the heat transfer rates will be larger (approximately 15% maximum) than those predicted by the simplified model using zero interfacial velocity. In general, the inclination decreases the rate of heat transfer as well as the rate of evaporation.

Nagendra, H. R.↗

NASA NTRS1973

Computational considerations in application of the finite element method for analysis of unsteady flow around airfoils.

Difficulties related to the application of the finite element method for numerical solution of the Navier-Stokes equation are discussed. The unsteady flow around a stationary and oscillating NACA 0012 airfoil was analyzed. In an effort to establish the range of applicability of the developed numerical method, the accuracy and stability problems in representing unsteady flow patterns and determining pressure distribution around the airfoil were investigated. Error involved in determining the velocity field and in the numerical integration of the discretized equations were analyzed in terms of geometry of the finite element gridwork, boundary conditions for stationary and oscillating airfoils, angles of attack and flow conditions.

Bratanow, T.↗

NASA NTRS1973

The global wind system in the thermosphere.

The upper atmosphere wind distribution is calculated by solving the Navier-Stokes equations with models of upper atmospheric densities and ion distribution. All non-linear terms and viscous terms are included. It is found that the velocities from the non-linear solutions are smaller than those obtained from a linear treatment and that in the equatorial region the magnitude of the average meridional velocity is strongly dependent upon the latitudinal coupling. No appreciable super-rotation is obtained by the inclusion of the nonlinear terms; the small value of the average zonal velocity is mainly due to the ion distribution.

Blum, P. W.↗

NASA NTRS1973

A new technique for calculating reentry base heating

The laminar base flow field of a two-dimensional reentry body has been studied by Telenin's method. The flow domain was divided into strips along the x-axis, and the flow variations were represented by Lagrange interpolation polynomials in the transformed vertical coordinate. The complete Navier-Stokes equations were used in the near wake region, and the boundary layer equations were applied elsewhere. The boundary conditions consisted of the flat plate thermal boundary layer in the forebody region and the near wake profile in the downstream region. The resulting two-point boundary value problem of 33 ordinary differential equations was then solved by the multiple shooting method. The detailed flow field and thermal environment in the base region are presented in the form of temperature contours, Mach number contours, velocity vectors, pressure distributions, and heat transfer coefficients on the base surface. The maximum heating rate was found on the centerline, and the two-dimensional stagnation point flow solution was adquate to estimate the maximum heating rate so long as the local Reynolds number could be obtained.

Meng, J. C. S.↗

NASA NTRS1974

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗

NASA NTRS1974

Numerical computation of two dimensional viscous blunt body flows with an impinging shock, part 2

Two-dimensional viscous blunt body flows with an impinging shock have been computed using a time-dependent finite-difference method which solves the complete set of Navier-Stokes equations for a compressible flow. For low Reynolds number flows, the entire flow field, including the bow shock and impinging shock, has been captured in the computation. For higher Reynolds number flows, the bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied, while the boundary layer and interaction regions are captured as before. Using this latter shock-fitting approach, a Type III shock interaction flow field has been computed with flow conditions corresponding to the space shuttle orbiter freestream conditions at 61 km (200,000 ft).

Holst, T. L.↗

NASA NTRS1974

Viscous, radiating hypersonic flow about a blunt body

The viscous, radiating hypersonic flow past an axisymmetric blunt body is analyzed based on the Navier-Stokes equations, plus a radiative equation of transfer derived from the Milne-Eddington differential approximation. The fluid is assumed to be a perfect gas with constant specific heats, a constant Prandtl number of order unity, a viscosity coefficient varying as a power of the temperature, and an absorption coefficient varying as the first power of the density and as a power of the temperature. The gray gas assumption is invoked, thereby making the absorption coefficient independent of the spectral frequency. Limiting forms of the solutions are studied as the freestream Mach number freestream Reynolds number and the temperature ratio across the shock wave, go to infinity, and as the Bouguer number and the density ratio across the shock wave go to zero. The method of matched asymptotic expansions is used in the analysis, and it is shown that there is a far-field precursor, composed of two regions, in which the fluid mechanics can be neglected for all practical purposes but included for completeness.

Passamaneck, R. S.↗