Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes”

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 721 records · Page 40

Spectral solution of the incompressible Navier-Stokes equations on the Connection Machine 2

The issue of solving the time-dependent incompressible Navier-Stokes equations on the Connection Machine 2 is addressed, for the problem of transition to turbulence of the steady flow in a channel. The spectral algorithm used serially requires O(N(4)) operations when solving the equations on an N x N x N grid; using the massive parallelism of the CM, it becomes an O(N(2)) problem. Preliminary timings of the code, written in LISP, are included and compared with a corresponding code optimized for the Cray-2 for a 128 x 128 x 101 grid.

Tomboulian, Sherryl↗

Least-squares finite elements for Stokes problem

A least-squares method based on the first-order velocity-pressure-vorticity formulation for the Stokes problem is proposed. This method leads to a minimization problem rather than to a saddle-point problem. The choice of the combinations of elements is thus not subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) condition. Numerical results are given for the optimal rate of convergence for equal-order interpolations.

Jiang, Bo-Nan↗

On the Nature of Navier-stokes Turbulence

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, largest Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The solutions, separate exponentially with time, having a positive Liapunov exponent. Thus the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number is reduced. Several qualitatively different flows are noted. These are, fully chaotic, complex period, weakly chaotic, simple periodic, and fixed-point. Of these, only the fully chaotic flows are classified as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, is both time-dependent and aperiodic. Turbulent solutions are obtained in which energy cascades from large to small-scale motions. In general, the spectral energy transfer takes place between wavenumber bands that are considerably separated. The special transfer can occur either as a result of nonlinear turbulence self-interaction or by interaction of turbulence with mean gradients. Turbulent systems are compared with those studied in kinetic theory. The two types of systems are fundamentally different (continuous and dissipative as opposed to discrete and conservative), but there are similarities. For instance, both are nonlinear and show sensitive dependence on initial conditions. Also, the turbulent and molecular stress tensors are identical if the macroscopic velocities for the turbulent stress are replaced by molecular velocities.

Deissler, Robert G.↗

Development of a Navier-Stokes algorithm for parallel-processing supercomputers

An explicit flow solver, applicable to the hierarchy of model equations ranging from Euler to full Navier-Stokes, is combined with several techniques designed to reduce computational expense. The computational domain consists of local grid refinements embedded in a global coarse mesh, where the locations of these refinements are defined by the physics of the flow. Flow characteristics are also used to determine which set of model equations is appropriate for solution in each region, thereby reducing not only the number of grid points at which the solution must be obtained, but also the computational effort required to get that solution. Acceleration to steady-state is achieved by applying multigrid on each of the subgrids, regardless of the particular model equations being solved. Since each of these components is explicit, advantage can readily be taken of the vector- and parallel-processing capabilities of machines such as the Cray X-MP and Cray-2.

Swisshelm, Julie M.↗

Aerodynamic applications of an efficient incompressible Navier-Stokes solver

A robust and accurate upwind scheme based on flux-difference splitting for three-dimensional, incompressible flows in general coordinates is described. Viscous effects are included through a thin-layer approximation to the Navier-Stokes equations. Turbulent flows are simulated with extensions to the Baldwin-Lomax (1978) turbulence model. The effects of accuracy and gridding on the computed flow field results are assessed. The results obtained from the code for vortical flows over a delta wing, a double delta wing, and a tangent-ogive forebody compare well with experimental data.

Hartwich, Peter M.↗

Navier-Stokes calculation of solid-propellant rocket motor internal flowfields

A comprehensive numerical analysis has been carried out to study the detailed physical and chemical processes involved in the combustion of homogeneous propellant in a rocket motor. The formulation is based on the time-dependent full Navier-Stokes equations, with special attention devoted to the chemical reactions in both gas and condensed phases. The turbulence closure is achieved using both the Baldwin-Lomax algebraic model and a modified k-epsilon two-equation scheme with a low Reynolds number and near-wall treatment. The effects of variable thermodynamic and transport properties are also included. The system of governing equations are solved using a multi-stage Runge-Kutta shceme with the source terms treated implicitly. Preliminary results clearly demonstrate the presence of various combustion regimes in the vicinity of propellant surface. The effects of propellant combustion on the motor internal flowfields are investigated in detail.

Hsieh, Kwang-Chung↗

Unsteady solution of incompressible Navier-Stokes equations

The numerical scheme used by the present time-accurate FEM numerical method for incompressible Navier-Stokes equations, using primitive variables as the unknowns, is a Crank-Nicholson implicit treatment of all equation terms with central differencing for space derivatives. The introduction of a continuous auxilliary system in pseudo-time, with artificial compressibility, yields the incompressible solution at the advanced time level; time-accurate solutions are thereby obtained for two-dimensional fluid flows in a square cavity, in the cases of both an impulsively starting lid and an oscillating lid.

Soh, W. Y.↗

Large-scale turbulence simulation on the Navier-Stokes Computer

Analytical techniques for turbulent air or water flow over vehicles and for the development of active control methods are described, with a focus on simulations being implemented on the Navier-Stokes Computer (NSC) at Princeton University. Consideration is given to the general NSC architecture (the central control, memory, and processing units; the switching network; and the prototype NSC), the algorithm for large-scale MHD turbulence (Orszag and Tang, 1979), and the pipeline configurations. Extensive diagrams and flow charts are provided.

Hayder, M. E.↗

Explicit and implicit solution of the Navier-Stokes equations on a massively parallel computer

The design, implementation, and performance of a two-dimensional time-accurate Navier-Stokes solver for the CM2 supercomputer are described. The program uses a single processor for each grid point. Two different time-stepping methods have so far been implemented: an explicit third-order Runge-Kutta method and an implicit approximation-factorization method. The CM2 results are checked against those of a mature well-vectorized Cray 2 program, both for correctness and performance. The code is found to be correct, and the performance in some cases is up to several times that of the Cray 2.

Levit, Creon↗

An implicit flux-split algorithm for the compressible Navier-Stokes equations

An implicit upwind scheme for the compressible Navier-Stokes equations is described and applied to the internal flow in a dual-throat nozzle. The method is second-order accurate spatially and naturally dissipative. A spatially-split approximate factorization method is used to obtain efficient steady-state solutions on the NASA Langley VPS-32 (CYBER 205) supercomputer.

Thomas, James L.↗

Navier-Stokes calculations of hovering rotor flowfields

Unsteady, thin-layer Navier-Stokes equations written in rotor coordinates are solved using a partially flux-split, implicit numerical algorithm to calculate the flowfields of a hovering rotor blade at subsonic and transonic conditions. Numerical results are in good agreement with experimental data for both nonlifting and lifting rectangular blades. For the lifting case, the rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blades. Alternate methods of calculating the hovering rotor flowfields in blade-fixed mode that have the same circulation distribution as a hovering blade are explored. All of the results presented in this paper were computed on a Cray-2 supercomputer.

Srinivasan, G. R.↗

Navier-Stokes calculations using Cartesian grids. I - Laminar flows

A finite-volume formulation for the Navier-Stokes equations using Cartesian grids is used to study flows past airfoils. In addition to the solution of the complete equations, solutions for two simplified versions of the governing equations were obtained and compared with those using body-fitted grids. Results are presented for two airfoil sections, NACA 0012 and RAE 2822, for a range of Mach numbers, angles of attack, and Reynolds numbers. It is shown that the results are highly dependent on the smoothness of the surface grid. Without such smoothness, the skin friction and pressure converge to nonuniform distributions. On the other hand, when surface cells with smoothly varying areas are used, the results compared favorably with calculations employing body-fitted grids.

Frymier, Paul D., Jr.↗

Splitting methods for low Mach number Euler and Navier-Stokes equations

Examined are some splitting techniques for low Mach number Euler flows. Shortcomings of some of the proposed methods are pointed out and an explanation for their inadequacy suggested. A symmetric splitting for both the Euler and Navier-Stokes equations is then presented which removes the stiffness of these equations when the Mach number is small. The splitting is shown to be stable.

Abarbanel, Saul↗

The solution of the Navier-Stokes equations using Gauss-Seidel line relaxation

The Navier-Stokes equations in an implicit flux-split difference formulation are solved numerically using a Gauss-Seidel line-relaxation procedure. Particular attention is given to the selection of flux-vector splitting method and flux splitting in boundary layers. Results for sample problems involving (1) turbulent supersonic flow over a cone and (2) the viscous hypersonic flow of a chemically reacting gas in thermal nonequilibrium past a blunted cone are presented in extensive graphs and briefly characterized. The present flux-split procedures are shown to provide accurate shear-layer calculations.

Maccormack, Robert W.↗

Unsteady Navier-Stokes computations past oscillating delta wing at high incidence

The unsteady, thin-layer, compressible Navier-Stokes equations, written in the moving frame of reference for the flow relative motion, is solved for the steady and unsteady supersonic flow around a round-edged delta wing. For supersonic flow, local conical flow solution has been obtained from the three-dimensional equations. Pseudotime stepping is used for the steady flow, while time-accurate stepping is used for the unsteady flow. The computational scheme is an implicit approximately-factored finite volume scheme which uses explicit and implicit dissipation terms. The scheme is verified for the steady flow solution. The scheme is then applied to a delta wing undergoing rolling oscillation at a reduced frequency of 1.137 with 15- deg maximum amplitude about a mean angle of attack of 10 deg for a Mach number of 2 and a Reynolds number of 500,000.

Kandil, Osama A.↗

Multigrid solution of the Navier-Stokes equations on triangular meshes

A new Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite-element Galerkin approximation, which can be shown to be equivalent to a finite-volume approximation for regular equilateral triangular meshes. Integration to steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. The convergence efficiency of the present method is shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, D. J.↗

Upwind Navier-Stokes solutions for leading-edge vortex flows

An incompressible Navier-Stokes solver using an upwind finite-difference algorithm is employed to investigate low-speed, three-dimensional, laminar, leading-edge vortex flows over three round-edged low-aspect-ratio wings. The effects of grid density, Reynolds number, and wing planform on the flowfield structures and integral values are studied. Computed results show good qualitative and quantitative agreement with the available experimental data.

Hsu, C.-H.↗

Comparison of LDV measurements and Navier-Stokes solutions in a two-dimensional 180-degree turn-around duct

Results from an experimental and numerical investigation of turbulent subsonic flow inside a two-dimensional, strongly curved, 180 deg turn-around duct are presented. Data measured with a two-component, two-color laser Doppler velocimeter include profiles of mean axial velocity and local flow angle. Static pressure distributions are also measured. Results are obtained at a Mach number of 0.1, and at Reynolds numbers of 100,000 and a million based on channel height. Numerical calculations are performed using the incompressible Navier-Stokes equations with a Prandtl mixing length zero-equation turbulence model modified for internal flows. Theory and experiment are compared to evaluate the ability of the turbulence model to predict this class of internal flows with strong curvature.

Monson, Daryl J.↗