Search

The method of pseudocompressibility is tested for its accuracy in solving the incompressible Navier-Stokes equations. An implicit, finite difference computer code is used to solve the equations in a three-dimensional, curvilinear coordinate system. The code employs artificial compressibiity for solving the pressure field, coupled with an implicit, approximate factorization scheme. This coupling is known as the pseudocompressibility method and introduces pressure waves of finite speed into the fluid medium that would otherwise have an infinite sound speed. The waves die out as the solution converges, and the steady state solution approaches a divergence-free condition. These waves, however, limit the time accuracy of the computations. The effects of these waves are analysed and criteria are set for choosing the pseudocompressibility parameter that governs the pressure wave speed; test cases are presented that verify these criteria. The code is tested by computing laminar flow over a two-dimensional, backward-facing step and over a two-dimensional, impulsively started circular cylinder.

The laminar, steady juncture flow around single or multiple posts mounted between two flat plates is simulated using the three dimensional incompressible Navier-Stokes code, INS3D. The three dimensional separation of the boundary layer and subsequent formation and development of the horseshoe vortex is computed. The computed flow compares favorably with the experimental observation. The recent numerical study to understand and quantify the juncture flow relevant to the Space Shuttle main engine power head is summarized.

The incompressible laminar flow around single and double rows of an infinite number of posts between two flat plates has been calculated numerically. A C-grid with periodic boundary conditions has been used. The angle of attack, measured from the line normal to the row of posts, was varied between zero and less than 90 deg. The pressure loading was computed for all of the posts in each of these cases. Most of these calculations have been carried out on the Numerical Aerodynamic Simulation Cray 2 at NASA Ames Research Center.

The method of pseudocompressibility is being tested for its accuracy in solving the incompressible Navier-Stokes equations. An implicit, finite-difference computer code is used to solve the equations in a three-dimensional, curvilinear coordinate system. The code employs artificial compressibility for solving the pressure field, coupled with an implicit, approximate-factorization scheme. This coupling is known as the pseudocompressibility method. The pseudocompressibility method introduces pressure waves of finite speed into the fluid medium that would otherwise have an infinite sound speed. The waves die out as the solution converges, and the steady state solution approaches a divergence-free condition. However, these waves limit the time-accuracy of the computations. The effects of these waves are analyzed and criteria are set for choosing the pseudocompressibility parameters that govern the pressure wave speed. Test cases are presented that verify these criteria. The code is tested by computing laminar flow over a two-dimensional, backward-facing step and over a two-dimensional, impulsively started circular cylinder.

An implicit two-equation turbulence solver in generalized coordinates has been developed and is used in conjunction with the three-dimensional incompressible Navier-Stokes solver, INS3D, to calculate the internal flow in one channel and in an additional channel with a sudden 2:3 expansion. A new and consistent boundary procedure for a low Reynolds number form of the kappa-epsilon turbulence model is chosen to integrate the equations up to the wall. The high Reynolds number form of the equations are integrated using wall functions. The latter approach yields a faster convergence to the steady state solution than the former. For the case of channel flow, both the wall function and wall boundary condition approaches yield results in good agreement with the experimental data. the back-step (sudden expansion) flow is calculated using wall function approach. The predictions are in reasonable agreement with the experimental data.

An implicit finite difference code was developed which solves the transport equations for the turbulence kinetic energy and its dissipation rate in generalized coordinates in three dimensions. The finite difference equations are solved using the Beam-Warming algorithm. The kinetic energy-dissipation code, KEM, provides the closure; i.e., the turbulent viscosity for calculation of either compressible or incompressible flows. Turbulent internal flow over a backward-facing step has been calculated using the present code in conjunction with the Incompressible Navier-Stokes Code, INS3D. The results are in good agreement with experiments and two dimensional computations of other researchers.

The flow into an open return wind tunnel inlet was simulated using Euler equations. An explicit predictor-corrector method was employed to solve the system. The calculation is time-accurate and was performed to achieve a steady-state solution. The predictions are in reasonable agreement with the experimental data. Wall pressures are accurately predicted except in a region of recirculating flow. Flow-field surveys agree qualitatively with laser velocimeter measurements. The method can be used in the design process for open return wind tunnels.

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.

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching-step sizes with the implicit shock fitting procedure. A correlation is found between the marching-step size, Reynolds number and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic, and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching step sizes with the implicit shock fitting procedure. A correlation is found between the marching step size, Reynolds number, and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.