SEARCH · Engineering Papers
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.
The Navier-Stokes equation is transformed into a closed nonlinear equation for the convected curl operator. It is shown how Lie algebraic methods can be employed to solve the obtained double-commutator equation. The proposed operator equation provides a way to the further development of the Navier-Stokes fluid flow theory.
The unaveraged Navier-Stokes equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three dimensional cosine velocity fluctuations and periodic boundary conditions are used. No mean gradients are present. The three components of the mean square velocity fluctuations are equal for the initial conditions chosen. The resulting solution shows characteristics of turbulence, such as the nonlinear excitation of small scale fluctuations. For the higher Reynolds numbers the initially nonrandom flow develops into an apparently random turbulence.
Navier-Stokes equation in predicting solar wind flow pattern
Explicit, implicit, and characteristic finite-difference methods are applied to solve model equations representative of the compressible Navier-Stokes equations. An approach is then formulated for solving the Navier-Stokes equation at high Reynolds numbers. The approach has drastically reduced the computation time required to obtain viscous flow solutions. Computational results for shock wave separated flows are presented.
A grid generation technique called the two boundary technique is developed and applied for the solution of the three dimensional Navier-Stokes equations. The Navier-Stokes equations are transformed from a cartesian coordinate system to a computational coordinate system, and the grid generation technique provides the Jacobian matrix describing the transformation. The two boundary technique is based on algebraically defining two distinct boundaries of a flow domain and the distribution of the grid is achieved by applying functions to the uniform computational grid which redistribute the computational independent variables and consequently concentrate or disperse the grid points in the physical domain. The Navier-Stokes equations are solved using a MacCormack time-split technique. Grids and supersonic laminar flow solutions are obtained for a family of three dimensional corners and two spike-nosed bodies.
It is well known that the classical Navier-Stokes equation for three-dimensional incompressible viscous flow is invariant under Galilean transformations. In the present paper, it is shown that, although the general invariance group is no larger than Galilean-scale dilitation, the Navier-Stokes equation for incompressible viscous flow admits time-dependent rotational invariance for solutions of special types.
The adaptive solution of the Navier-Stokes equations depends upon the successful interaction of three key elements: (1) the ability to flexibly select grid length scales in composite grids, (2) the ability to efficiently control residual error in composite grids, and (3) the ability to define reliable, convenient error norms to guide the grid adjustment and optimize the residual levels relative to the local truncation errors. An initial investigation was conducted to explore how to approach developing these key elements. Conventional error assessment methods were defined and defect and deferred correction methods were surveyed. The one dimensional potential equation was used as a multigrid test bed to investigate how to achieve successful interaction of these three key elements.
The full Navier-Stokes equations for incompressible turbulent flow must be solved to accurately represent all flow phenomena which occur in a high Reynolds number incompressible flow. A two layer algebraic eddy viscosity turbulence model is used to represent the Reynolds stress in the primitive variable formulation. The development of the boundary-fitted coordinate systems makes the numerical solution of these equations feasible for arbitrarily shaped bodies. The nondimensional time averaged Navier-Stokes equations, including the turbulence mode, are represented by finite difference approximations in the transformed plane. The resulting coupled system of nonlinear algebraic equations is solved using a point successive over relaxation iteration. The test case considered was a NACA 64A010 airfoil section at an angle of attack of two degrees and a Reynolds number of 2,000,000.
The steady-state Navier-Stokes equations are solved for hypersonic flow about blunt axisymmetric bodies. The equations of motion are solved by successive approximations using an implicit finite-difference scheme. The results are compared with viscous shock-layer theory, experimental data, and time-dependent solutions of the Navier-Stokes equations. It is demonstrated that viscous shock-layer theory is sufficiently accurate for the range of flight conditions normally encountered by entry vehicles.
The steady-state Navier-Stokes equations are solved for hypersonic flow about blunt axisymmetric bodies. The equations of motion are solved by successive approximations using an implicit finite-difference scheme. The results are compared with viscous shock-layer theory, experimental data, and time-dependent solutions of the Navier-Stokes equations. It is demonstrated that viscous shock-layer theory is sufficiently accurate for the range of flight conditions normally encountered by entry vehicles.
Four explicit finite difference techniques designed to solve the time-dependent, compressible Navier Stokes equations are compared. These techniques are: (1) MacCormack, (2) modified Du Fort-Frankel, (3) modified hopscotch, and (4) Brailovskaya. The comparison was made numerically by solving the quasi-one dimensional Navier Stokes equations for the flow in a converging-diverging nozzle. Solutions with and without standing normal shock waves were computed for unit Reynolds numbers (based on total conditions) ranging from 45374 to 2269. The results indicate that all four techniques are comparable in accuracy; however, the modified hopscotch scheme is two to three times faster than the Brailovskaya and MacCormack schemes and three to six times faster than the modified Du Fort-Frankel scheme.
Steady solutions of the Navier-Stokes equations, in terms of velocity and pressure, for breakdown in an unconfined viscous vortex are obtained numerically using the artificial compressibility technique of Chorin combined with an ADI finite-difference scheme. Axisymmetry is assumed and boundary conditions are carefully applied at the boundaries of a large finite region in an axial plane while resolution near the axis is maintained by a coordinate transformation. The solutions, which are obtained for Reynolds numbers up to 200 based on the free-stream axial velocity and a characteristic core radius, show that breakdown results from the diffusion and convection of vorticity away from the vortex core which, because of the strong coupling between the circumferential and axial velocity fields in strongly swirling flows, can lead to stagnation and reversal of the axial flow near the axis.
The time dependent Navier-Stokes equations in mass averaged variables are solved for transonic flow over axisymmetric boattail plume simulator configurations. Numerical solution of these equations is accomplished with the unsplit explict finite difference algorithm of MacCormack. A grid subcycling procedure and computer code vectorization are used to improve computational efficiency. The two layer algebraic turbulence models of Cebeci-Smith and Baldwin-Lomax are employed for investigating turbulence closure. Two relaxation models based on these baseline models are also considered. Results in the form of surface pressure distribution for three different circular arc boattails at two free stream Mach numbers are compared with experimental data. The pressures in the recirculating flow region for all separated cases are poorly predicted with the baseline turbulence models. Significant improvements in the predictions are usually obtained by using the relaxation models.
An explicit-implicit technique for solving Navier-Stokes equations is described which, is much less complex than other implicit methods. It is used to solve a complex, two-dimensional, steady-state, supersonic-flow problem. The computational efficiency of the method and the quality of the solution obtained from it at high Courant-Friedrich-Lewy (CFL) numbers are discussed. Modifications are discussed and certain observations are made about the method which may be helpful in using it successfully.
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.
The proposed approach to the derivation of the Navier-Stokes equation is thought to be more plausible and easier to understand than other derivations that can be found in works on fluid mechanics. The tensor character of the stress is central to the derivation. In particular, a linear relation between stress and strain rate is assumed only for the shear, rather than for the full stress tensor as is done in most other derivations. An assumption for the shear is naturally simpler and easier to verify experimentally. The use of tensor analysis is shown to greatly simplify the derivation.
The numerical solution of the Navier-Stokes equations in general two dimensional domains is considered. A proper finite difference approximation to the governing equations, on nonstaggered grids, in a transformed plane is formulated. Several aspects of a multigrid method for the solution of the finite difference equations are described. The efficiency of some relaxation schemes and the transfer among the grids which are nonuniform in the physical plane are emphasized.
To study the nonlinear physics of uniform turbulent shear flow, the unaveraged Navier-Stokes equations are solved numerically. This extends our previous work in which mean gradients were absent. For initial conditions, modified three-dimensional-cosine velocity fluctuations are used. The boundary conditions are modified periodic conditions on a stationary three-dimensional numerical grid. A uniform mean shear is superimposed on the initial and boundary conditions. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. As in the case of no shear the initially nonrandom flow develops into an apparently random turbulence at higher Reynolds number. Thus, randomness or turbulence can apparently arise as a consequence of the structure of the Navier-Stokes equations. Except for an initial period of adjustment, all fluctuating components grow with time. The initial equality of the three intensity components is destroyed by the shear, the transverse components becoming smaller than the longitudinal one, in agreement with experiment. Also, the shear creates a small-scale structure in the turbulence. The nonlinear solutions are compared with linearized ones.