Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes 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 1,315 records · Page 73

Computational flow development for unsteady viscous flows: Foundation of the numerical method

A procedure is presented for effective consideration of viscous effects in computational development of high Reynolds number flows. The procedure is based on the interpretation of the Navier-Stokes equations as vorticity transport equations. The physics of the flow was represented in a form suitable for numerical analysis. Lighthill's concept for flow development for computational purposes was adapted. The vorticity transport equations were cast in a form convenient for computation. A statement for these equations was written using the method of weighted residuals and applying the Galerkin criterion. An integral representation of the induced velocity was applied on the basis of the Biot-Savart law. Distribution of new vorticity, produced at wing surfaces over small computational time intervals, was assumed to be confined to a thin region around the wing surfaces.

Bratanow, T.↗

Calculation of hypersonic forebody/inlet flow fields

The National Aero-Space Plane (NASP) is unique in its integration of the airframe and propulsion system. The performance of the NASP propulsion system (scramjet) is dependent upon the quality of the flow into the scramjet inlet. Hence, the calculation of the inlet flow field must be coupled to the forebody flow field for accurate performance predictions. The Compressible Navier-Stokes code has been extended to perform both external and internal hypersonic flow fields. Computations are presented that validate these changes and demonstrate the capability to perform coupled forebody/inlet computations for hypersonic flows.

Bennett, Bradford C.↗

Mathematical analysis of the Navier-Stokes equations with non standard boundary conditions

One of the major applications of the domain decomposition time marching algorithm is the coupling of the Navier-Stokes systems with Boltzmann equations in order to compute transitional flows. Another important application is the coupling of a global Navier-Stokes problem with a local one in order to use different modelizations and/or discretizations. Both of these applications involve a global Navier-Stokes system with nonstandard boundary conditions. The purpose of this work is to prove, using the classical Leray-Schauder theory, that these boundary conditions are admissible and lead to a well posed problem.

Tidriri, M. D.↗

Navier-Stokes Calculations With Deforming Grid

Scheme for numerical simulation of transonic, compressible, viscous flow about two-dimensional cascade of oscillating rigid airfoils involves finite-difference solution of Navier-Stokes equations on coordinate grid deforming with motions of airfoils. Contributes to understanding of flows around advanced turboprop airfoil sections. Cascade represented by periodic grid conforming to airfoils and to outer boundaries. Oscillation of airfoil represented by deformation of grid. Applications of deforming-grid technique includes turbomachinery and aeroelastic modeling of nonrigid airfoils.

Huff, Dennis L.↗

Unsteady stator/rotor interaction

The major thrust of the computational analysis of turbomachinery to date has been the steady-state solution of isolated blades using mass-averaged inlet and exit conditions. Unsteady flows differ from the steady solution due to interaction of pressure waves and wakes between blade rows. To predict the actual complex flow conditions one must look at the time accurate solution of the entire turbomachine. Three quasi-three-dimensional Euler and thin layer Navier-Stokes equations are solved for unsteady turbomachinery flows.

Jorgenson, Philip C. E.↗

Numerical simulation of flows in curved diffusers with cross-sectional transitioning using a three-dimensional viscous analysis

A three dimensional analysis for fully viscous, subsonic, compressible flow is evaluated. An approximate form of the Navier Stokes equations is solved by an implicit spatial marching technique. Calculations were made for flow in a circular S duct and in the F 16 inlet duct. The computed total pressure contours and secondary flow velocity vectors are presented. Qualitative comparisons with experiment are shown for both ducts. The analysis is used to show how the cross section transitioning in the F 16 inlet suppresses the development of a secondary flow vortex.

Towne, C. E.↗

Investigation of turbulent shear flows by numerical simulation

A data base (Kim et al., 1986) generated by direct numerical integration of the Navier-Stokes equations for unsteady three-dimensional turbulent channel flow at low Reynolds numbers is applied to investigate the physical processes acting in wall-bounded shear flows. The techniques used to construct the data base are outlined; the modeling of turbulence, streaks and streamwise vortices, horseshoe and inverted-horseshoe vortices, and vortex rings is described; the conditional sampling techniques employed are discussed; and typical results are presented graphically. An organized vorticity-field structure is revealed and found to be associated with the bursting process and a local peak in the wall pressure, indicating a possible approach to reducing viscous drag.

Kim, J.↗

Numerical studies of unsteady transonic flow over an oscillating airfoil

A finite-difference solution to the Navier-Stokes equations combined with a time-varying grid-generation technique was used to compute unsteady transonic flow over an oscillating airfoil. These computations were compared with experimental data (obtained at Ames Research Center) which form part of the AGARD standard configuration for aeroelastic analysis. A variety of approximations to the full Navier-Stokes equations was used to determine the effect of frequency, shock-wave motion, flow separation, and airfoil geometry on unsteady pressures and overall air loads. Good agreement is shown between experiment and theory with the limiting factor being the lack of a reliable turbulence model for high-Reynolds-number, unsteady transonic flows.

Chyu, W. J.↗

Triangular spectral elements for incompressible fluid flow

We discuss the use of triangular elements in the spectral element method for direct simulation of incompressible flow. Triangles provide much greater geometric flexibility than quadrilateral elements and are better conditioned and more accurate when small angles arise. We employ a family of tensor product algorithms for triangles, allowing triangular elements to be handled with comparable arithmetic complexity to quadrilateral elements. The triangular discretizations are applied and validated on the Poisson equation. These discretizations are then applied to the incompressible Navier-Stokes equations and a laminar channel flow solution is given. These new triangular spectral elements can be combined with standard quadrilateral elements, yielding a general and flexible high order method for complex geometries in two dimensions.

Mavriplis, C.↗

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri=-2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri = -2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗

A finite difference scheme for three-dimensional steady laminar incompressible flow

A finite difference scheme for three-dimensional steady laminar incompressible flows is presented. The Navier-Stokes equations are expressed conservatively in terms of velocity and pressure increments (delta form). First order upwind differences are used for first order partial derivatives of velocity increments resulting in a diagonally dominant matrix system. Central differences are applied to all other terms for second order accuracy. The SIMPLE pressure correction algorithm is used to satisfy the continuity equation. Numerical results are presented for cubic cavity flow problems for Reynolds numbers up to 2000 and are in good agreement with other numerical results.

Hwang, Danny P.↗

A finite difference scheme for three-dimensional steady laminar incompressible flow

A finite difference scheme for three-dimensional steady laminar incompressible flows is presented. The Navier-Stokes equations are expressed conservatively in terms of velocity and pressure increments (delta form). First order upwind differences are used for first order partial derivatives of velocity increments resulting in a diagonally dominant matrix system. Central differences are applied to all other terms for second order accuracy. The SIMPLE pressure correction algorithm is used to satisfy the continuity equation. Numerical results are presented for cubic cavity flow problems for Reynolds numbers up to 2000 and are in good agreement with other numerical results.

Hwang, Danny P.↗

Triangular spectral elements for incompressible fluid flow

We discuss the use of triangular elements in the spectral element method for direct simulation of incompressible flow. Triangles provide much greater geometric flexibility and are better conditioned and more accurate than quadrilateral elements when small angles arise. We employ a family of tensor product algorithms for triangles, allowing triangular elements to be handled with comparable arithmetic complexity to quadrilateral elements. The triangular discretizations are applied to the Poisson equation and are validated. The triangular discretizations are then applied to the incompressible Navier-Stokes equations, and a laminar channel flow solution is given. The new triangular spectral elements can be combined with standard quadrilateral elements, yielding a general and flexible high order method for complex geometries in two dimensions.

Mavriplis, Catherine↗

A well-posed optimal spectral element approximation for the Stokes problem

A method is proposed for the spectral element simulation of incompressible flow. This method constitutes in a well-posed optimal approximation of the steady Stokes problem with no spurious modes in the pressure. The resulting method is analyzed, and numerical results are presented for a model problem.

Maday, Y.↗

A semi-elliptic analysis for 2-D viscous flows through cascade configurations

A semielliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is negligible but which are significantly influenced by upstream interactions. A two-step alternating-direction-explicit numerical scheme is developed to solve the resulting governing equations efficiently. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Various simple channel and cascade flow configurations are examined for several values of Re, thickness ratio and Mach number, in order to establish the applicability of the IPNS model for these flows. With this model, solutions to compressible as well as nearly incompressible flows are obtained without any modification either in the analysis or in the solution procedure.

Ramamurti, R.↗

Radiative interactions in chemically reacting compressible nozzle flows using Monte Carlo simulations

The two-dimensional spatially elliptic Navier-Stokes equations have been used to investigate the radiative interactions in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The radiative heat transfer term in the energy equation is simulated using the Monte Carlo method (MCM). The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The spectral correlation has been considered in the Monte Carlo formulations. Results obtained demonstrate that the effect of radiation on the flow field is minimal but its effect on the wall heat transfer is significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and the nozzle size on the radiative and conductive wall fluxes.

Liu, J.↗

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗