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,009 records · Page 56

Calculations Of Transonic Flow About A Wing

Report describes calculations of transonic airflows about wing in wind tunnel. Basic equations of flow used in study are Reynolds-averaged Navier-Stokes equations in strong conservation-law form. Equations of flow incorporated into finite-difference computer code called TNS (Transonic Navier-Stokes). Computational grid generated by solution of partial differential equations yielding smooth meshes conforming to surfaces of wing and wind tunnel.

Holst, Terry L.↗

Development and Validation of a Fast, Accurate and Cost-Effective Aeroservoelastic Method on Advanced Parallel Computing Systems

Progress to date towards the development and validation of a fast, accurate and cost-effective aeroelastic method for advanced parallel computing platforms such as the IBM SP2 and the SGI Origin 2000 is presented in this paper. The ENSAERO code, developed at the NASA-Ames Research Center has been selected for this effort. The code allows for the computation of aeroelastic responses by simultaneously integrating the Euler or Navier-Stokes equations and the modal structural equations of motion. To assess the computational performance and accuracy of the ENSAERO code, this paper reports the results of the Navier-Stokes simulations of the transonic flow over a flexible aeroelastic wing body configuration. In addition, a forced harmonic oscillation analysis in the frequency domain and an analysis in the time domain are done on a wing undergoing a rigid pitch and plunge motion. Finally, to demonstrate the ENSAERO flutter-analysis capability, aeroelastic Euler and Navier-Stokes computations on an L-1011 wind tunnel model including pylon, nacelle and empennage are underway. All computational solutions are compared with experimental data to assess the level of accuracy of ENSAERO. As the computations described above are performed, a meticulous log of computational performance in terms of wall clock time, execution speed, memory and disk storage is kept. Code scalability is also demonstrated by studying the impact of varying the number of processors on computational performance on the IBM SP2 and the Origin 2000 systems.

Goodwin, Sabine A.↗

A Parallel Incompressible Navier-Stokes Solver with Multigrid Iterations

We developed a parallel, numerically accurate and stable, and computationally efficient finate-difference incompressible Navier-Stokes (N-S) fluid flow solver. The solver runs on both sequential and massively parallel computers. The numerical method used here is a second-order projection method on a staggered grid. The code is highly modular and it can be used either as a stand-alone flow solver and or a template code which can be adapted or expanded to a specific application. Numerical results and parallel performances of our code on Intel Delta and Paragon are reported.

Navier-Stokes solver↗

Applications of the contravariant form of the Navier-Stokes equations

The contravariant Navier-Stokes equations in weak conservation form are well suited to certain fluid flow analysis problems. Three dimensional contravariant momentum equations may be used to obtain Navier-Stokes equations in weak conservation form on a nonplanar two dimensional surface with varying streamsheet thickness. Thus a three dimensional flow can be simulated with two dimensional equations to obtain a quasi-three dimensional solution for viscous flow. When the Navier-Stokes equations on the two dimensional nonplanar surface are transformed to a generalized body fitted mesh coordinate system, the resulting equations are similar to the equations for a body fitted mesh coordinate system on the Euclidean plane. Contravariant momentum components are also useful for analyzing compressible, three dimensional viscous flow through an internal duct by parabolic marching. This type of flow is efficiently analyzed by parabolic marching methods, where the streamwise momentum equation is uncoupled from the two crossflow momentum equations. This can be done, even for ducts with a large amount of turning, if the Navier-Stokes equations are written with contravariant components.

Katsanis, T.↗

Investigation of nose bluntness and angle of attack effects on slender bodies in viscous hypersonic flows

Hypersonic flows over cones and straight biconic configurations are calculated for a wide range of free stream conditions in which the gas behind the shock is treated as perfect. Effect of angle of attack and nose bluntness on these slender cones in air is studied extensively. The numerical procedures are based on the solution of complete Navier-Stokes equations at the nose section and parabolized Navier-Stokes equations further downstream. The flow field variables and surface quantities show significant differences when the angle of attack and nose bluntness are varied. The complete flow field is thoroughly analyzed with respect to velocity, temperature, pressure, and entropy profiles. The post shock flow field is studied in detail from the contour plots of Mach number, density, pressure, and temperature. The effect of nose bluntness for slender cones persists as far as 200 nose radii downstream.

Sehgal, A. K.↗

Combined effect of nose bluntness and angle of attack on slender bodies in viscous hypersonic flows

Hypersonic flows over cones and straight biconic configurations are calculated for a wide range of free stream conditions in which the gas behind the shock is treated as perfect. Effect of angle of attack and nose bluntness on these slender cones in air is studied extensively. The numerical procedures are based on the solution of complete Navier-Stokes equations at the nose section and parabolized Navier-Stokes equations further downstream. The flow field variables and surface quantities show significant differences when the angle of attack and nose bluntness are varied. The complete flow field is thoroughly analyzed with respect to velocity, temperature, pressure, and entropy profiles. The post shock flow field is studied in detail from the contour plots of Mach number, density, pressure, and temperature. The effect of nose bluntness for slender cones persists as far as 200 nose radii downstream.

Tiwari, S. N.↗

Three-dimensional Navier-Stokes heat transfer predictions for turbine blade rows

Results are shown for a three-dimensional Navier-Stokes analysis of both the flow and the surface heat transfer for turbine applications. Heat transfer comparisons are made with the experimental shock-tunnel data of Dunn and Kim, and with the data of Blair for the rotor of the large scale rotating turbine. The analysis was done using the steady-state, three-dimensional, thin-layer Navier-Stokes code developed by Chima, which uses a multistage Runge-Kutta scheme with implicit residual smoothing. An algebraic mixing length turbulence model is used to calculate turbulent eddy viscosity. The variation in heat transfer due to variations in grid parameters is examined. The effects of rotation, tip clearance, and inlet boundary layer thickness variation on the predicted blade and endwall heat transfer are examined.

Boyle, R. J.↗

Three-dimensional Navier-Stokes heat transfer predictions for turbine blade rows

Results are shown for a three-dimensional Navier-Stokes analysis of both the flow and the surface heat transfer for turbine applications. Heat transfer comparisons are made with the experimental shock-tunnel data of Dunn and Kim, and with the data of Blair for the rotor of the large scale rotating turbine. The analysis was done using the steady-state, three-dimensional, thin-layer Navier-Stokes code developed by Chima, which uses a multistage Runge-Kutta scheme with implicit residual smoothing. An algebraic mixing length turbulence model is used to calculate turbulent eddy viscosity. The variation in heat transfer due to variations in grid parameters is examined. The effects of rotation, tip clearance, and inlet boundary layer thickness variation on the predicted blade and endwall heat transfer are examined.

Boyle, Robert J.↗

Numerical Simulations of Plasma Based Flow Control Applications

A mathematical model was developed to simulate flow control applications using plasma actuators. The effects of the plasma actuators on the external flow are incorporated into Navier Stokes computations as a body force vector. In order to compute this body force vector, the model solves two additional equations: one for the electric field due to the applied AC voltage at the electrodes and the other for the charge density representing the ionized air. The model is calibrated against an experiment having plasma-driven flow in a quiescent environment and is then applied to simulate a low pressure turbine flow with large flow separation. The effects of the plasma actuator on control of flow separation are demonstrated numerically.

Suzen, Y. B.↗

A 3D-PNS computer code for the calculation of supersonic combusting flows

A computer code has been developed based on the three-dimensional parabolized Navier-Stokes (PNS) equations which govern the supersonic combusting flow of the hydrogen-air system. The finite difference algorithm employed was a hybrid of the Schiff-Steger algorithm and the Vigneron, et al., algorithm which is fully implicit and fully coupled. The combustion of hydrogen and air was modeled by the finite-rate two-step combustion model of Rogers-Chinitz. A new dependent variable vector was introduced to simplify the numerical algorithm. Robustness of the algorithm was considerably enhanced by introducing an adjustable parameter. The computer code was used to solve a premixed shock-induced combustion problem and the results were compared with those of a full Navier-Stokes code. Reasonably good agreement was obtained at a fraction of the cost of the full Navier-Stokes procedure.

Chitsomboon, Tawit↗

Source Term Model for Steady Micro Jets in a Navier-Stokes Computer Code

A source term model for steady micro jets was implemented into a non-proprietary Navier-Stokes computer code, OVERFLOW. The source term models the mass flow and momentum created by a steady blowing micro jet. The model is obtained by adding the momentum and mass flow created by the jet to the Navier-Stokes equations. The model was tested by comparing with data from numerical simulations of a single, steady micro jet on a flat plate in two and three dimensions. The source term model predicted the velocity distribution well compared to the two-dimensional plate using a steady mass flow boundary condition, which was used to simulate a steady micro jet. The model was also compared to two three-dimensional flat plate cases using a steady mass flow boundary condition to simulate a steady micro jet. The three-dimensional comparison included a case with a grid generated to capture the circular shape of the jet and a case without a grid generated for the micro jet. The case without the jet grid mimics the application of the source term. The source term model compared well with both of the three-dimensional cases. Comparisons of velocity distribution were made before and after the jet and Mach and vorticity contours were examined. The source term model allows a researcher to quickly investigate different locations of individual or several steady micro jets. The researcher is able to conduct a preliminary investigation with minimal grid generation and computational time.

Waithe, Kenrick A.↗

Numerical solution of supersonic laminar flow over an inclined body of revolution

A mixed explicit-implicit scheme is used to solve the time-dependent thin-layer approximation of the Navier-Stokes equations for a supersonic laminar flow over an inclined body of revolution. Test cases for Mach 2.8 flow over a cylinder with 15 deg flare angle at angles of attack of 0, 1, and 4 deg are calculated. Good agreement is obtained between the present computed results and experimental measurements of surface pressure. A pair of vortices on the leeward and a peak in the normal force distribution near the flared juncture are predicted; the role of circumferential communication is discussed.

Hung, C. M.↗

Numerical Solution of Supersonic Laminar Flow Over an Inclined Body of Revolution

A mixed explicit-implicit scheme is used to solve the time-dependent thin-layer approximation of the Navier-Stokes equations for a supersonic laminar flow over an inclined body of revolution. Test cases for Mach 2.8 flow over a cylinder with 15-deg flare angle at angles of attack of 0,1, and 4 deg are calculated. Good agreement is obtained between the present computed results and experimental measurements of surface pressure. A pair of vortices on the leeward and a peak in the normal force distribution near the flared juncture are predicted; the role of circumferential communication is discussed.

Hung, C. M.↗

Hydrodynamic stability problem formulated by numerical solutions of the Navier-Stokes equations

The problem of hydrodynamic stability and the transition from laminar to turbulent flows are reformulated by seeking numerical solutions of the full, unapproximated Navier-Stokes equations. This method differs significantly from the well known Orr-Sommerfeld equation approach. The oncoming laminar flow is disturbed by forced, time dependent perturbations. The magnitudes of these perturbations are arbitrary. Then, the ensuing spatial and temporal development of the imposed perturbations on the basic flow is calculated by direct numerical solutions of the time dependent Navier-Stokes equations. One of the main advantages of this method is its ability to simulate nonlinear processes. As a specific application of this technique to SSME (Space Shuttle Main Engine) flow configurations, computer programs have been written for the two dimensional flow over a backward facing step. This numerical code will be tested for operational use as part of a continued research collaboration effort with the NASA/MSFC counterparts.

Hyun, J. M.↗

Effect of Plasma Sheaths on Earth Re-entry MHD Processes

The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. Results obtained indicate that the cathode sheath can exert significant resistance to current flow and thus reduce the MHD-generated force by 40% or more. The effect of the sheath on the MHD force is found to be more pronounced at lower applied magnetic field.

Earth Re-entry↗

Effect of Plasma Sheaths on Earth Re-entry MHD Processes

The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. The effect of the electrode material on the MHD process is studied. Using thoriated tungsten instead of graphite leads to a four-fold increase in the Lorentz forces and also leads to significantly reduced heat fluxes on the cathode. This is attributed to the much higher electrical conductivity of the thoriated tungsten sheath reducing by orders of magnitude the plasma electrical resistance near the surfaces.

Bernard Parent↗

Recent computation of viscous effects in transonic flow

The paper describes some improvements to a method of solving the time-dependent Reynolds-averaged Navier-Stokes equations for two-dimensional compressible flows for describing the coupling between viscous and inviscid regions in transonic flows. Mesh patching is implemented at the sublayer-boundary layer interfaces and the boundary layer-inviscid flow interface. Mesh stretching is geometric and is small in regions of large gradients and large in regions of small gradients, thus preserving second-order accuracy. Mesh lines are arranged so as to accommodate MacCormack's rapid solver for hyperbolic systems and a wider class of turbulent transport models. Results of computations of transonic flows over lifting airfoils are compared with experimental data and other solutions.

Deiwert, G. S.↗

Flow-thermal-structural study of aerodynamically heated leading edges

A finite element approach for integrated fluid-thermal-structural analysis of aerodynamically heated leading edges is presented. The Navier-Stokes equations for high speed compressible flow, the energy equation, and the quasi-static equilibrium equations for the leading edge are solved using a single finite element approach in one integrated, vectorized computer program called LIFTS. The fluid-thermal-structural coupling is studied for Mach 6.47 flow over a 3-inch diameter cylinder for which the flow behavior and the aerothermal loads are calibrated by experimental data. Issues of the thermal-structural response are studied for hydrogen cooled, super thermal conducting leading edges subjected to intense aerodynamic heating.

Dechaumphai, Pramote↗