Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes equations”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Nonlinear stability of oscillatory core-annular flow: A generalized Kuramoto-Sivashinsky equation with time periodic coefficients

In this paper the nonlinear stability of two-phase core-annular flow in a pipe is examined when the acting pressure gradient is modulated by time harmonic oscillations and viscosity stratification and interfacial tension is present. An exact solution of the Navier-Stokes equations is used as the background state to develop an asymptotic theory valid for thin annular layers, which leads to a novel nonlinear evolution describing the spatio-temporal evolution of the interface. The evolution equation is an extension of the equation found for constant pressure gradients and generalizes the Kuramoto-Sivashinsky equation with dispersive effects found by Papageorgiou, Maldarelli & Rumschitzki, Phys. Fluids A 2(3), 1990, pp. 340-352, to a similar system with time periodic coefficients. The distinct regimes of slow and moderate flow are considered and the corresponding evolution is derived. Certain solutions are described analytically in the neighborhood of the first bifurcation point by use of multiple scales asymptotics. Extensive numerical experiments, using dynamical systems ideas, are carried out in order to evaluate the effect of the oscillatory pressure gradient on the solutions in the presence of a constant pressure gradient.

Coward, Adrian V.↗
NASA NTRS1996

On buffer layers as non-reflecting computational boundaries

We examine an absorbing buffer layer technique for use as a non-reflecting boundary condition in the numerical simulation of flows. One such formulation was by Ta'asan and Nark for the linearized Euler equations. They modified the flow inside the buffer zone to artificially make it supersonic in the layer. We examine how this approach can be extended to the nonlinear Euler equations. We consider both a conservative and a non-conservative form modifying the governing equations in the buffer layer. We compare this with the case that the governing equations in the layer are the same as in the interior domain. We test the effectiveness of these buffer layers by a simulation of an excited axisymmetric jet based on a nonlinear compressible Navier-Stokes equations.

Hayder, M. Ehtesham↗
NASA NTRS1996

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗
NASA NTRS2008

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗
NASA NTRS1991

Investigation of vortex development on a pitching slender body of revolution

A computational study of the unsteady flow about a pitching 3.5 caliber tangent ogive forebody is presented. The flow is simulated using the full 3D unsteady Navier-Stokes equations and a time-accurate implicit algorithm. Comparison to available experimental data for a steady 20 deg case is presented as validation. Effects of grid resolution and a comparison of solutions using full Navier-Stokes and the thin-layer approximation are included. The forebody is simulated in a 'pitch-up to 20 deg and hold' maneuver, and two different pitch axis locations are used in the study. Examination of the unsteady vorticity field for the pitch-up cases reveals the formation of strong shear layers as the body decelerates, and their roll-up into vortical structures in a process similar to that observed in airfoil dynamic stall. Pronounced vortex/surface interactions are seen which produce multiple secondary separation regions, ejection of vorticity from the surface, and embedded regions of high suction.

Stanek, M. J.↗
NASA NTRS1978

Stochastic partial differential equations in turbulence related problems

The theory of stochastic partial differential equations (PDEs) and problems relating to turbulence are discussed by employing the theories of Brownian motion and diffusion in infinite dimensions, functional differential equations, and functional integration. Relevant results in probablistic analysis, especially Gaussian measures in function spaces and the theory of stochastic PDEs of Ito type, are taken into account. Linear stochastic PDEs are analyzed through linearized Navier-Stokes equations with a random forcing. Stochastic equations for waves in random media as well as model equations in turbulent transport theory are considered. Markovian models in fully developed turbulence are discussed from a stochastic equation viewpoint.

Chow, P.-L.↗
NASA NTRS1982

Solution of 3-dimensional time-dependent viscous flows. Part 3: Application to turbulent and unsteady flows

A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.

Weinberg, B. C.↗
NASA NTRS1986

Numerical simulation of electrophoresis separation processes

A new Petrov-Galerkin finite element formulation has been proposed for transient convection-diffusion problems. Most Petrov-Galerkin formulations take into account the spatial discretization, and the weighting functions so developed give satisfactory solutions for steady state problems. Though these schemes can be used for transient problems, there is scope for improvement. The schemes proposed here, which consider temporal as well as spatial discretization, provide improved solutions. Electrophoresis, which involves the motion of charged entities under the influence of an applied electric field, is governed by equations similiar to those encountered in fluid flow problems, i.e., transient convection-diffusion equations. Test problems are solved in electrophoresis and fluid flow. The results obtained are satisfactory. It is also expected that these schemes, suitably adapted, will improve the numerical solutions of the compressible Euler and the Navier-Stokes equations.

Ganjoo, D. K.↗
NASA NTRS1989

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗
NASA NTRS1993

Grid and aerodynamic sensitivity analyses of airplane components

An algorithm is developed to obtain the grid sensitivity with respect to design parameters for aerodynamic optimization. The procedure is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B-Splines) for defining the wing-section geometry. An interactive algebraic grid generation technique, known as Two-Boundary Grid Generation (TBGG) is employed to generate C-type grids around wing-sections. The grid sensitivity of the domain with respect to geometric design parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations such as a wing or fuselage. The aerodynamic sensitivity coefficients are obtained by direct differentiation of the compressible two-dimensional thin-layer Navier-Stokes equations. An optimization package has been introduced into the algorithm in order to optimize the wing-section surface. Results demonstrate a substantially improved design due to maximized lift/drag ratio of the wing-section.

Sadrehaghighi, Ideen↗
NASA NTRS1995

Grid Sensitivity and Aerodynamic Optimization of Generic Airfoils

An algorithm is developed to obtain the grid sensitivity with respect to design parameters for aerodynamic optimization. The procedure is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B- Splines) for defining the airfoil geometry. An interactive algebraic grid generation technique is employed to generate C-type grids around airfoils. The grid sensitivity of the domain with respect to geometric design parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations such as a wing or fuselage. The aerodynamic sensitivity coefficients are obtained by direct differentiation of the compressible two-dimensional thin-layer Navier-Stokes equations. An optimization package has been introduced into the algorithm in order to optimize the airfoil surface. Results demonstrate a substantially improved design due to maximized lift/drag ratio of the airfoil.

Sadrehaghighi, Ideen↗
NASA NTRS1972

Compressible boundary layer with normal pressure gradients: Quasi-similarity equations - Their properties at the wall and at sharp and blunt leading edges.

The failure of most viscous-inviscid interaction methods at strong interactions is attributed to the presence of a normal pressure gradient. A new theory is proposed for supersonic laminar boundary layers that can generate normal pressure gradients. The Navier-Stokes equations are reexamined by an order of magnitude analysis and all first and second order terms are retained. The approximation is found to be dependent not only on the boundary layer thickness but also on the ratio of the dimensionless viscosity and density. The equations are transformed into two quasi-similar, nonlinear, third order, ordinary integro-differential equations for the velocity and pressure as functions of a single transverse variable. The properties of the equations at the boundaries are discussed.

Seginer, A.↗
NASA NTRS1981

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 1: Theory document (version 1.1)

Panel aerodynamics (PAN AIR) is a system of computer programs designed to analyze subsonic and supersonic inviscid flows about arbitrary configurations. A panel method is a program which solves a linear partial differential equation by approximating the configuration surface by a set of panels. An overview of the theory of potential flow in general and PAN AIR in particular is given along with detailed mathematical formulations. Fluid dynamics, the Navier-Stokes equation, and the theory of panel methods were also discussed.

Magnus, A. E.↗
NASA NTRS1992

Three-dimensional boundary layer stability and transition

Nonparallel and nonlinear stability of a three-dimensional boundary layer, subject to crossflow instability, is investigated using parabolized stability equations (PSEs). Both traveling and stationary disturbances are considered and nonparallel effect on crossflow instability is found to be destabilizing. Our linear PSE results for stationary disturbances agree well with the results from direct solution of Navier-Stokes equations obtained by Spalart (1989). Nonlinear calculations have been carried out for stationary vortices and the computed wall vorticity pattern results in streamwise streaks which resemble remarkably well with the surface oil-flow visualizations in swept-wing experiments. Other features of the stationary vortex development (half-mushroom structure, inflected velocity profiles, vortex doubling, etc.) are also captured in our nonlinear calculations. Nonlinear interaction of the stationary amplitude of the stationary vortex is large as compared to the traveling mode, and the stationary vortex dominates most of the downstream development. When the two modes have the same initial amplitude, the traveling mode dominates the downstream development owing to its higher growth rate, and there is a tendency for the stationary mode to be suppressed. The effect of nonlinear wave development on the skin-friction coefficient is also computed.

Malik, M. R.↗
NASA NTRS1993

Numerical simulation of crossing/turbulent boundary layer interaction at Mach 8.3 comparison of zero and two-equation turbulence models

A 3D hypersonic crossing shock wave/turbulent boundary layer interaction is examined numerically. The test geometry consists of a pair of opposing sharp fins of angle alpha = 15 deg mounted on a flat plate. The freestream Mach number is 8.28. Two theoretical models are evaluated. The full 3D Reynolds-averaged Navier-Stokes equations are solved using the Baldwin-Lomax algebraic turbulent eddy viscosity model and the Rodi turbulence model. Computed results for both cases show good agreement with experiment for flat plate surface pressure and for pitot pressure and yaw angle profiles in the flowfield. General agreement is obtained for surface flow direction. Fair to poor agreement is obtained for surface heat transfer, indicating a need for more accurate turbulence models. The overall flowfield structure is similar to that observed in previous crossing shock interaction studies.

Narayanswami, N.↗
NASA NTRS1977

Over-relaxation applied to the MacCormack finite-difference scheme

An over-relaxation procedure is applied to the MacCormack finite-difference scheme in order to reduce the computation time required to obtain a steady-state solution. The implementation of this acceleration procedure to an existing computer program using the regular MacCormack method is extremely simple and does not require additional storage. The over-relaxation procedure does not alter the steady-state solution, which is second-order accurate. The method is first applied to Burgers' equation. A stability condition and an expression for the increase in the rate of convergence are derived. The method is then applied to the calculation of the hypersonic viscous flow over a flat plate, using the complete Navier-Stokes equations, and the inviscid flow over a wedge. Reductions in computing time by factors of 3 and 1.5, respectively, are obtained by over-relaxation.

Desideri, J.-A.↗
NASA NTRS1985

Application of a viscous-inviscid interaction method to predict transonic separated flows

A viscous-inviscid interaction numerical procedure is used to predict turbulent transonic separated flows over axisymmetric bodies. The viscous flow solutions are obtained by solving the boundary-layer equations inversely in a coupled manner by a finite-difference scheme; the inviscid flow is determined by a finite-difference solution to the full-potential equation in conservative form. Parameters affecting the convergence rate of the interaction procedure are discussed. The predictions of a new turbulence model proposed by Johnson and King are compared with experimental data and numerical solutions to the Navier-Stokes equations for flows with large regions of separation.

Lee, D. S.↗
NASA NTRS1985

The collisional dynamics of particulate disks

It is shown that the use of the Krook equation greatly simplifies the discussion of the collisional dynamics of particulate disks. Collisions are assumed to be inelastic from the outset. A simple heuristic argument is used to compute the effective collision rate in a disk of spherical particles with a power-law distribution of sizes. For Saturn's rings, the effective collision rate for momentum transport is substantially lower than that conventionally estimated on the basis of an observed optical depth at visual wavelengths. The closure of the vertically integrated set of moment equations without discarding the third-order moments at the outset is discussed, allowing for the possibility of a bent disk. In the limit that the collision frequency is much larger than the orbit frequency, the usual Navier-Stokes equations of viscous hydrodynamics for a thin disk are recovered with an explicit expression for the shear viscosity. The method is generalized to include the effects of gravitational scattering.

Shu, F. H.↗