On the stability of Poiseuille flow in a pipe
Poiseuille pipe flow stability from finite difference equations approximation to nonlinear axisymmetric Navier-Stokes equations under stream function perturbation
SEARCH · Search NASA
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.
Poiseuille pipe flow stability from finite difference equations approximation to nonlinear axisymmetric Navier-Stokes equations under stream function perturbation
For computational simulation of the convection and mixing of stratified fluid in a rotating tank (such as used in Apollo flights) with time-dependent rotation, the Navier-Stokes convection problem was formulated for a circular tank configuration. The final equations results from a general approximate theory for combined forced and contained natural convection in a time-dependent rotating system. The equations are cast in terms of vorticity and stream function in a form convenient for computation, with a transformed coordinate system, and appropriate boundary conditions are derived. Accurate representations for the cryogenic supercritical oxygen thermodynamic properties are used in the computations, and an efficient numerical finite difference scheme and computational method are employed.
Rapid iterative (or semidirect) computation methods are developed for the finite-difference solution of the nonlinear equations of subsonic and transonic aerodynamics. At each iteration, a fast, direct elliptic algorithm solves the entire computation field. In an application to subsonic flow over a lifting airfoil, the full nonlinear stream-function equation is solved. Finally, a direct Cauchy-Riemann solver is used for the nonlinear transonic small-disturbance equations for a biconvex airfoil. At M = 0.7, t/c = 0.1 (subcritical), three iterations on a 39 x 32 mesh (totaling 2.45 sec on an IBM 360/67 computer) obtain convergence within 0.1%. A slightly supercritical case requires seven iterations (6.75 sec) for convergence within 1%.
Two computers programs, known as Matrix Through-Flow and Matrix Blade-To-Blade, for analyzing the meridional and blade-to-blade flow patterns are described. The numerical solutions are obtained by finite difference approximations to the governing Poisson-type differential equations for the stream function. Solutions for several turbomachines, giving flow patterns and velocity distributions, are included.
The application of the elements of quantum (wave) mechanics to some special problems in the field of macroscopic fluid dynamics is discussed. Emphasis is placed on the flow of a viscous, incompressible fluid around a circular cylinder. The following subjects are considered: (1) the flow of a nonviscous fluid around a circular cylinder, (2) the restrictions imposed the stream function by the number of dimensions of space, and (3) the flow past three dimensional bodies in a viscous fluid, particularly past a circular cylinder in the symmetrical case.
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.
A hydrodynamic model of aircraft vortex wakes in an irregular wind shear field near the ground is developed and used as a basis for modeling the characteristics of a laser Doppler detection and vortex location system. The trailing vortex sheet and the wind shear are represented by discrete free vortices distributed over a two-dimensional grid. The time dependent hydrodynamic equations are solved by direct numerical integration in the Boussinesq approximation. The ground boundary is simulated by images, and fast Fourier Transform techniques are used to evaluate the vorticity stream function. The atmospheric turbulence was simulated by constructing specific realizations at time equal to zero, assuming that Kolmogoroff's law applies, and that the dissipation rate is constant throughout the flow field. The response of a simulated laser Doppler velocimeter is analyzed by simulating the signal return from the flow field as sensed by a simulation of the optical/electronic system.
Abstracts are presented on a method of numerical solution of the Navier-Stokes equation for the flow about arbitrary airfoils, using a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Results of continuing research are reported and include: application of the Navier-Stokes solution in the vorticity-stream function formulation to a number of single airfoils at Reynolds numbers up to 2000; programming of the Navier-Stokes solution for multiple airfoils in the primitive variable formulation; testing of the potential flow solution of multiple bodies; and development of a generalized coordinate system program.
Rotor and propeller performance and induced potential flowfields were studied on the basis of a rotating actuator disk concept, with special emphasis on rotors hovering out of ground effect. A new theory for the optimum performance of rotors hovering OGE is developed and presented. An extended theory for the optimum performance of rotors and propellers in axial motion is also presented. Numerical results are presented for the optimum distributions of blade-bound circulation together with axial inflow and ultimate wake velocities for the hovering rotor over the range of thrust coefficient of interest in rotorcraft applications. Shapes of the stream tubes and of the velocities in the slipstream are obtained, using available methods, for optimum and off-optimum circulation distributions for rotors hovering in and out of ground effect. A number of explicit formulae useful in computing rotor and propeller induced flows are presented for stream functions and velocities due to distributions of circular vortices over axi-symmetric surfaces.
A cubic spline approximation is presented which is suited for many fluid-mechanics problems. This procedure provides a high degree of accuracy, even with a nonuniform mesh, and leads to an accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several implicit and explicit integration schemes are presented. For two-dimensional flows, a spline-alternating-direction-implicit method is evaluated. The spline procedure is assessed, and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.
The stalling characteristics of an airfoil in a laminar viscous incompressible fluid are investigated. The governing equations in terms of the vorticity and stream function are solved using an implicit finite-difference scheme and point successive relaxation procedure. The development of the impulsively started flow, the initial generation of circulation, and the behavior of the forces at large times are studied. Following the impulsive start, the lift is at first very large and then it rapidly drops. The subsequent growth of circulation and lift is associated with the starting vortex. After incipient separation, the lift increases owing to enlargement of the separation bubble and intensification of the flow rotation in it. The extension of this bubble beyond the trailing edge causes it to rupture and brings about the stalling characteristics of the airfoil. Subsequently, new bubbles are formed near the upper surface of the airfoil and are swept away. The behavior of the lift acting on the airfoil is explained in terms of the strength and sense of these bubbles.
Results are presented which were obtained with the Winge-Coleman model for theoretical predictions of latitudinal dependencies in the solar wind. A first-order expansion is described which allows analysis of first-order latitudinal variations in the coronal boundary conditions and results in a second-order partial differential equation for the perturbation stream function. Latitudinal dependencies are analytically separated out in the form of Legendre polynomials and their derivative, and are reduced to the solution of radial differential equations. This analysis is shown to supply an estimate of how large the coronal variation in latitude must be to produce an 11 km/sec/deg gradient in the radial velocity of the solar wind, assuming steady-state processes.
A cubic spline approximation is used for the solution of several problems in fluid mechanics. This procedure provides a high degree of accuracy even with a nonuniform mesh, and leads to a more accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several typical integration schemes are presented. For two-dimensional flows a spline-alternating-direction-implicit (SADI) method is evaluated. The spline procedure is assessed and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.
A method based on the Navier-Stokes equations was developed for analyzing the unsteady incompressible viscous flow around oscillating airfoils at high Reynolds numbers. The Navier-Stokes equations have been integrated in their classical Helmholtz vorticity transport equation form, and the instantaneous velocity field at each time step was determined by the solution of Poisson's equation. A refined finite element was utilized to allow for a conformable solution of the stream function and its first space derivatives at the element interfaces. A corresponding set of accurate boundary conditions was applied; thus obtaining a rigorous solution for the velocity field. The details of the computational procedure and examples of computed results describing the unsteady flow characteristics around the airfoil are presented.
A comparative study of the efficiency of some finite difference methods for the solution of the Navier-Stokes equations was conducted. The study was restricted to the two-dimensional steady, uniform property vorticity-stream function equations. The comparisons were drawn by recording the CPU time required to obtain a solution as well as the accuracy of this solution using five numerical methods: central differences, first order upwind differences, second order upwind differences, exponential differences, and an ADI solution of the central difference equations. Solutions were obtained for two test cases: a recirculating eddy inside a square cavity with a moving top, and an impinging jet flow. The results show that whenever the central difference method is stable it generates results with a given accuracy for less CPU time than any other method.
Numerical solutions of the laminar, incompressible boundary layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified; the pressure is deduced from the solution. A vorticity-stream-function formulation of the boundary layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite difference equations, in which the flow direction is taken into account, are obtained by several global iteration schemes which are stable and have unconditional diagonal dominance. Results are compared with Klineberg and Steger's separated boundary layer calculations, and with Briley's solution of Navier-Stokes equations for a separated region. In addition, an approximate technique is presented in which the streamwise convection of vorticity is set equal to zero in the reversed flow region; such a technique results in a quick forward marching procedure for separated flows.
An integrodifferential method, previously formulated in terms of velocity and vorticity vectors, is reformulated in terms of stream function and vorticity for two-dimensional incompressible viscous flows. The reformulated integrodifferential method is shown to retain the distinguishing feature of the previous formulation in permitting the confinement of the solution field to the viscous region of the flow and consequently offers great computational advantages. The application of this procedure in a study of an incompressible flow around an impulsively started 9% thick symmetric Joukowski airfoil at an angle of attack of 15 deg and a Reynolds number of 1000 is discussed. Numerical results are presented and compared with available finite-difference results.
A computer code has been developed for the direct calculation of shockless transonic airfoils whose pressure distributions can be assigned within reasonable limits. The partial differential equations of two-dimensional inviscid gas dynamics are solved by analytic continuation into the domain of two independent complex characteristic coordinates. The domain of integration is mapped conformally onto the unit circle in the hodograph plane of one of these coordinates. It is possible to formulate a boundary value problem on this circle for the stream function that is well posed in the case of transonic flow. This enables the formulation of a procedure for the calculation of an airfoil on which the speed is prescribed as a function of the arc length