Search NASA⌕ Search

SEARCH · Search NASA

Results for “Incompressible 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 379 records · Page 21

Numerical Simulations of Vortex Shedding in Hydraulic Turbines

Turbomachines for rocket propulsion applications operate with many different working fluids and flow conditions. Oxidizer boost turbines often operate in liquid oxygen, resulting in an incompressible flow field. Vortex shedding from airfoils in this flow environment can have adverse effects on both turbine performance and durability. In this study the effects of vortex shedding in a low-pressure oxidizer turbine are investigated. Benchmark results are also presented for vortex shedding behind a circular cylinder. The predicted results are compared with available experimental data.

Dorney, Daniel↗

Numerical simulation of shock/turbulent boundary layer interaction

Most flows of aerodynamic interest are compressible and turbulent. However, our present knowledge on the structures and mechanisms of turbulence is mostly based on incompressible flows. In the present work, compressibility effects in turbulent, high-speed, boundary layer flows are systematically investigated using the Direct Numerical Simulation (DNS) approach. Three-dimensional, time-dependent, fully nonlinear, compressible Navier-Stokes equations were numerically integrated by high-order finite-difference methods; no modeling for turbulence is used during the solution because the available resolution is sufficient to capture the relevant scales. The boundary layer problem deals with fully-turbulent compressible flows over flat geometries. Apart from its practical relevance to technological flows, turbulent compressible boundary layer flow is the simplest experimentally realizable turbulent compressible flow. Still, measuring difficulties prohibit a detailed experimental description of the flow, especially in the near-wall region. DNS studies provide a viable means to probe the physics of compressible turbulence in this region. The focus of this work is to explore the paths of energy transfer through which compressible turbulence is sustained. The structural similarities and differences between the incompressible and compressible turbulence are also investigated. The energy flow patterns or energy cascades are found to be directly related to the evolution of vortical structures which are generated in the near-wall region. Near-wall structures, and mechanisms which are not readily accessible through physical experiments are analyzed and their critical role on the evolution and the behavior of the flow is documented extensively.

Biringen, Sedat↗

Turbulent flow between a rotating disk and a parallel wall

A scalar effective viscosity method is used to calculate the incompressible flow in the gap between two infinite, parallel smooth disks, one of them rotating and the other stationary. The flow is calculated for two values of the gap Reynolds number omega(s squared)/nu, namely 296 and 2852, where omega is the disk angular velocity, s the gap width, and nu the kinematic viscosity of the fluid; and over a range of radial Reynolds number omega(r squared)/nu from 0 to 10 million, where r is the radius from the axis of rotation. A numerical shooting technique is used to match the two boundary-layer flows arising on the disk and wall. Gradual transition is assumed for radial Reynolds numbers from 160,000 to 250,000. The results are in good agreement with Daily and Nece's measurements of velocity profiles and disk friction drag for enclosed rotating disks, and at sufficiently large radius the local skin friction approaches that of turbulent Couette flow.

Cooper, P.↗

A computational study of saddle point separation and horseshoe vortex system

Incompressible flow around a cylinder-end wall junction has been simulated by solving the incompressible Navier-Stokes equations in three dimensions. The equations, cast in generalized curvilinear coordinates, are solved in time as a hyperbolic system by adding a pressure term in the continuity equation and are marched to a steady state. Various physical quantities associated with the saddle point of separation and the horseshoe vortex system are calculated. Computational and experimental results are generally consistent. The skin friction and the pressure distribution on the end wall are consistent with the physics of the problem. Secondary flows both in front of the cylinder and behind it are predicted that are in qualitative agreement with flow visualization results. The calculations also indicate a strongly nonuniform pressure loading along the length of the cylinder. A new mechanism for the existence of the recirculation bubbles behind the cylinder-end wall with relatively low ratio of cylinder height to the approaching boundary layer thickness is observed which is markedly different from its two-dimensional counterpart.

Kaul, U. K.↗

A method for analyzing dynamic stall of helicopter rotor blades

A model for each of the basic flow elements involved in the unsteady stall of a two-dimensional airfoil in incompressible flow is presented. The interaction of these elements is analyzed using a digital computer. Computations of the loading during transient and sinusoidal pitching motions are in good qualitative agreement with measured loads. The method was used to confirm that large torsional response of helicopter blades detected in flight tests can be attributed to dynamic stall.

Crimi, P.↗

Numerical simulation of flow over wavy walls

Spectral methods are used to solve the Navier-Stokes equations numerically for two-dimensional incompressible flow over wavy surfaces. Results obtained for high Reynolds number laminar flows and turbulent flows are compared to experiment and analytical theory. Overall, good agreement has been achieved.

Balasubramanian, R.↗

Computational analyses of crystal growth

Two important aspects of Hg/Cd/Te crystal growth processes are discussed. First, the thermal field and second, the fluid movement in the melt zone. The thermal analysis includes numerical calculation of axisymmetric heat conduction within the sample. It also includes a three-dimensional radiation model to calculate the radiative heat exchange between the furnace and the crystal as determined by the complex geometry of the furnace and the adiabatic shield. The thermal analysis also includes a crystal conductivity which is dependent on temperature and composition. To tackle the fluid flow aspect of the problem, an attempt was made to use a newly developed incompressible flow code based on the slight compressibility, and hence the finite sound speed, of all real fluids.

Dakhoul, Youssef M.↗

Computational flow simulation of liquid oxygen in a SSME preburner injector element LOX post

Liquid oxygen (LOX) is simulated as an incompressible flow through a Space Shuttle main engine fuel preburner injector element LOX post for the full range of operating conditions. Axial profiles of axial velocity and static pressure are presented. For each operating condition analyzed, the minimum pressure downstream of the orifice is compared to the vapor pressure to determine if cavitation could occur. Flow visualization is provided by velocity vectors and stream function contours. The results indicate that the minimum pressure is too high for cavitation to occur. To establish confidence in the CFD analysis, the simulation is repeated with water flow through a superscaled LOX post and compared with experimental results. The agreement between calculated and experimental results is very good.

Rocker, Marvin↗

Symmetry plane model for turbulent flows with vortex generators

An approximate procedure is proposed for predicting the performance of counterrotating vortex-generator installations in incompressible flow. An inviscid calculation that includes the motion of the vortices is used to obtain crossflow velocities at the boundary-layer edge as a function of initial position, spacing, and strength of the vortices, and local values of the spanwise gradient are then folded into an integral turbulent-boundary layer procedure applied in the plane of symmetry. Special attention is paid to the consistency of the approximations and equations used. The two-dimensional aerodynamics of vortex generator installations on a NACA 0016 airfoil at angle-of-attack are estimated in this manner, and the results compared with experiments carried out with a 30-cm chord wing mounted in a 2.4 x 3.6-m cross-section wind tunnel and tested at chord Reynolds numbers of 0.7 and 1.4 x 10 to the 6th. Agreement in the separation location is found for these complex flows for a range of conditions.

Arnaud, Gilles L.↗

Study Of Disturbances In Rotating Plane Poisseuille Flow

Report describes computer-simulation of instability and transition to turbulence in plane Poisseuille flow in rotating channel. Motions of fluid governed by Navier-Stokes equations of incompressible flow, using periodic boundary conditions in streamwise and spanwise directions. Principal parameters affecting flow are Reynolds and rotation numbers.

Kim, John↗

Investigation of Separation of the Turbulent Boundary Layer

An investigation was conducted on a turbulent boundary layer near a smooth surface with pressure gradients sufficient to cause flow separation. The reynolds number was high, but the speeds were entirely within the incompressible flow range. The investigation consisted of measurements of mean flow, three components of turbulence intensity, turbulent shearing stress, and correlations between two fluctuation components at a point and between the same component of different points. The results are given in the form of tables and graphs. The discussion deals first with separation and then with the more fundamental question of basic concepts of turbulent flow.

Schubauer, G B↗

FAS multigrid calculations of three dimensional flow using non-staggered grids

Grid staggering is a well known remedy for the problem of velocity/pressure coupling in incompressible flow calculations. Numerous inconveniences occur, however, when staggered grids are implemented, particularly when a general-purpose code, capable of handling irregular three-dimensional domains, is sought. In several non-staggered grid numerical procedures proposed in the literature, the velocity/pressure coupling is achieved by either pressure or velocity (momentum) averaging. This approach is not convenient for simultaneous (block) solvers that are preferred when using multigrid methods. A new method is introduced in this paper that is based upon non-staggered grid formulation with a set of virtual cell face velocities used for pressure/velocity coupling. Instead of pressure or velocity averaging, a momentum balance at the cell face is used as a link between the momentum and mass balance constraints. The numerical stencil is limited to 9 nodes (in 2D) or 27 nodes (in 3D), both during the smoothing and inter-grid transfer, which is a convenient feature when a block point solver is applied. The results for a lid-driven cavity and a cube in a lid-driven cavity are presented and compared to staggered grid calculations using the same multigrid algorithm. The method is shown to be stable and produce a smooth (wiggle-free) pressure field.

Matovic, D.↗

Thermal analysis of HGFQ using FIDAP(trademark): Solidification front motion

The High Gradient Furnace with Quench (HGFQ) is being designed by NASA/MSFC for flight on the International Space Station. The furnace is being designed specifically for solidification experiments in metal and metallic alloy systems. The HGFQ Product development Team (PDT) has been active since January 1994 and their effort is now in early Phase B. Thermal models have been developed both by NASA and Sverdrup (support contractor) to assist in the HGFQ design effort. Both these models use SINDA as a solution engine, but the NASA model was developed using PATRAN and includes more detail than the Sverdrup model. These models have been used to guide design decisions and have been validated through experimentation on a prototypical 'Breadboard' furnace at MSFC. One facet of the furnace operation of interest to the designers is the sensitivity of the solidification interface location to changes in the furnace setpoint. Specifically of interest is the motion (position and velocity) of the solidification front due to a small perturbation in the furnace temperature. FIDAP(TM) is a commercially available finite element program for analysis of heat transfer and fluid flow processes. Its strength is in solution of the Navier-Stokes equations for incompressible flow, but among its capabilities is the analysis of transient processes involving radiation and solidification. The models presently available from NASA and Sverdrup are steady-state models and are incapable of computing the motion of the solidification front. The objective of this investigation is to use FIDAP(TM) to compute the motion of the solidification interface due to a perturbation in the furnace setpoint.

Woodbury, Keith A.↗

On solving the compressible Navier-Stokes equations for unsteady flows at very low Mach numbers

The properties of a preconditioned, coupled, strongly implicit finite-difference scheme for solving the compressible Navier-Stokes equations in primitive variables are investigated for two unsteady flows at low speeds, namely the impulsively started driven cavity and the startup of pipe flow. For the shear-driven cavity flow, the computational effort was observed to be nearly independent of Mach number, especially at the low end of the range considered. This Mach number independence was also observed for steady pipe flow calculations; however, rather different conclusions were drawn for the unsteady calculations. In the pressure-driven pipe startup problem, the compressibility of the fluid began to significantly influence the physics of the flow development at quite low Mach numbers. The present scheme was observed to produce the expected characteristics of completely incompressible flow when the Mach number was set at very low values. Good agreement with incompressible results available in the literature was observed.

Pletcher, R. H.↗

On solving the compressible Navier-Stokes equations for unsteady flows at very low Mach numbers

The properties of a preconditioned, coupled, strongly implicit finite difference scheme for solving the compressible Navier-Stokes equations in primitive variables are investigated for two unsteady flows at low speeds, namely the impulsively started driven cavity and the startup of pipe flow. For the shear-driven cavity flow, the computational effort was observed to be nearly independent of Mach number, especially at the low end of the range considered. This Mach number independence was also observed for steady pipe flow calculations; however, rather different conclusions were drawn for the unsteady calculations. In the pressure-driven pipe startup problem, the compressibility of the fluid began to significantly influence the physics of the flow development at quite low Mach numbers. The present scheme was observed to produce the expected characteristics of completely incompressible flow when the Mach number was set at very low values. Good agreement with incompressible results available in the literature was observed.

Pletcher, R. H.↗

Sharply curved turn around duct flow predictions using spectral partitioning of the turbulent kinetic energy and a pressure modified wall law

Computational predictions of turbulent flow in sharply curved 180 degree turn around ducts are presented. The CNS2D computer code is used to solve the equations of motion for two-dimensional incompressible flows transformed to a nonorthogonal body-fitted coordinate system. This procedure incorporates the pressure velocity correction algorithm SIMPLE-C to iteratively solve a discretized form of the transformed equations. A multiple scale turbulence model based on simplified spectral partitioning is employed to obtain closure. Flow field predictions utilizing the multiple scale model are compared to features predicted by the traditional single scale k-epsilon model. Tuning parameter sensitivities of the multiple scale model applied to turn around duct flows are also determined. In addition, a wall function approach based on a wall law suitable for incompressible turbulent boundary layers under strong adverse pressure gradients is tested. Turn around duct flow characteristics utilizing this modified wall law are presented and compared to results based on a standard wall treatment.

Santi, L. Michael↗

A rapid numerical solution to subsonic flow over planar and axisymmetric profiles at an angle of attack of 0 deg

A numerical solution is presented for the incompressible flow over thin planar and axisymmetric profiles at an angle of attack of 0 degrees. The method uses a finite-difference field solution to the governing equation with a Gauss-Seidel successive overrelaxation scheme. However, the use of a simple Cartesian grid system restricts this method to slender profiles. Results are presented for a cambered airfoil, airfoil in wall effect (two-dimensional flowthrough inlet), body of revolution, and flowthrough nacelle. A computer program is presented which can be used for any of the previously mentioned cases with simple input changes. Results for compressible flow are available with the use of the appropriate two-dimensional or axisymmetric compressibility corrections. Computational time for a typical field calculation of 3000 grid points and 200 cycles through the field is less than 1 minute with less than 50,000 octal storage on the Control Data Corporation 6600 computing system.

Mascitti, V. R.↗

The obstacle block as a device to measure turbulent skin friction in compressible flow

The obstacle block, developed as an alternative to the Preston tube for indirectly measuring skin friction on smooth surfaces in incompressible flows, is examined as a device for compressible flows as well. The block, which is congruent with a surface static pressure orifice, has a geometry which is easily specified and thus has a universal calibration. Data from two independent studies are used to establish such a calibration using 'wall' variables, valid for Mach numbers up to about 3. Various aspects concerning practical application of the device are examined, such as sensitivity to yaw and the minimum permissible axial spacing between blocks. Several examples showing the utility of the device are given.

Elfstrom, G. M.↗