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 1,081 records · Page 60

Comparison of Coupled Radiative Navier-Stokes Flow Solutions with the Project Fire II Flight Data

A nonequilibrium, axisymmetric, Navier-Stokes flow solver with coupled radiation has been developed to use in the design of thermal protection systems for vehicles where radiation effects are important. The present method has been compared with an existing flow and radiation solver and with the Project Fire II experimental data. Very good agreement has been obtained over the entire Fire II trajectory with the experimentally determined values of the stagnation radiation intensity in the .2 to 6.2 eV range and with the total stagnation heating. The agreement was significantly better than previous numerical predictions. The effects of a number of flow models are examined to determine which combination of physical models produces the best agreement with the experimental data. These models include radiation coupling, multi-temperature thermal models, finite-rate chemistry, and a quasi-steady-state or Boltzmann assumption for the calculation of the excited electronic states. Finally, the computational efficiency of the present model is evaluated. The radiation properties model developed for this study is shown to offer significant computational savings compared to existing codes.

Olynick, David R.↗

Navier-Stokes Flowfield Simulation of Boeing 747-200 as Platform for SOFIA

Steady and unsteady viscous, three-dimensional flowfields are calculated using a thin layer approximation of Navier-Stokes equations in conjunction with Chimera overset grids. The finite-difference numerical scheme uses structured grids and a pentadiagonal flow solver called "OVERFLOW". The configuration of Boeing 747-200 has been chosen as one of configurations to be used as a platform for the SOFIA (Stratospheric Observatory For Infrared Astronomy). Initially, the steady flowfield of the full aircraft is calculated for the clean configuration (without a cavity to house telescope). This solution is then used to start the unsteady flowfield of a configuration containing cavity housing the observation telescope and its peripheral units. Analysis of unsteady flowfield in the cavity and its influence on the tail empennage, as well as the noise due to turbulence and optical quality of the flow are the main focus of this study. For the configuration considered here, the telescope housing cavity is located slightly downstream of the portwing. The entire flow-field is carefully constructed using 45 overset grids and consists of nearly 4 million grid points. All the computations axe done at one freestream flow condition of M(sub infinity) = 0.85, alpha = 2.5deg, and a Reynolds of Re = 1.85x10deg

Srinivasan, G.R.↗

Navier-Stokes Calculations For a Venus Composition Probe With Ablation And Wake Radiation

A numerical investigation is carried out to determine the magnitude of wake radiation for a proposed Venus composition probe. One of the scientific goals of the mission is to determine the atmospheric composition of Venus by examining the intensity of scattered sunlight through the wake of the vehicle during planetary entry. In the wake of the vehicle, excited particles generated in the bow shock and boundary layers absorb and emit radiation. Thus, the purpose of this study is to determine if the radiation sensor will be able to sense the incoming solar radiative flux relative to the radiative flux generated in the wake. During portions of the entry trajectory the incident surface heat flux will be high enough to produce significant ablation. Ablation products such as CN are known to be strong radiators. Also, the ablation will be driven by strong radiation emanating from the bow shock. Thus, radiation and ablation will be coupled into the Navier-Stokes flow solutions.

Olynick, David R.↗

A Data Parallel Multizone Navier-Stokes Code

We have developed a data parallel multizone compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the "chimera" approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. The design choices can be summarized as: 1. finite differences on structured grids; 2. implicit time-stepping with either distributed solves or data motion and local solves; 3. sequential stepping through multiple zones with interzone data transfer via a distributed data structure. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran (HPF). One interesting feature is the issue of turbulence modeling, where the architecture of a parallel machine makes the use of an algebraic turbulence model awkward, whereas models based on transport equations are more natural. We will present some performance figures for the code on the CM-5, and consider the issues involved in transitioning the code to HPF for portability to other parallel platforms.

Jespersen, Dennis C.↗

Fractional Step and Pseudocompressibility Methods for the Solution of Incompressible Navier-Stokes Equations

The fractional step and the pseudocompressibility methods for the solution of the incompressible Navier-Stokes equations are outlined. The fractional step method is based on finite-volume formulation and uses the pressure and the volume fluxes across the faces of each cell as dependent variables. The momentum equations are solved implicitly and the Poisson equation for the pressure is solved by using the multigrid method. The pseudocompressibility approach uses an implicit-higher-order-upwind differencing scheme for the convective terms together with the Gauss-Seidel line relaxation method. The dependent variables in the pseudocompressibility approach are the pressure and the cartesian velocity components in unstaggered mesh orientation. The 90-degree square duct flow, the wing-tip vortex wake flow and unsteady turbulent flows over an oscillating NACA 0015 airfoil are computed using both the fractional step and the pseudocompressibility methods. The results obtained from two different schemes are compared against experimental measurements.

Kiris, Cetin↗

Progress in Navier-Stokes Simulation of Delta Wing Maneuver Aerodynamics

This paper will summarize recent progress in the numerical simulation of high incidence vortical flow about a generic 65 degree sweep delta wing using the three dimensional, time-dependent, Reynolds averaged, Navier-Stokes (RANS) equations. Computations have been carried out at 15 and 30 degrees angle of attack under subsonic turbulent flow conditions, and compared with experimental data provided by Hanff, Jenkins, and their colleagues. This work has already been published elsewhere and widely disseminated. Computations carried out at 15 degrees angle of attack included static roll angles ranging up to 65 degrees, and a large-amplitude (40 degree), high rate (7 Hz), forced roll motion. There was very good agreement between computed and experimental forces and moments, and static surface pressures. There was a significant hysteresis of the dynamic rolling moment due to the high rate of roll motion. At this angle of attack, no vortex breakdown was observed in the computations or experiment. Computations were also carried out at 30 degrees angle of attack, where vortex breakdown was present in both the computations and experiment. There was overall good agreement in the computed and experimental forces and moments. The static rolling moment varied with roll angle in a highly nonlinear manner, and exhibited three stable trim points and two unstable trim points. This behavior was attributed to the presence of vortex breakdown. Two large-amplitude (30 degrees), high-rate (10 Hz) forced roll motions were computed. The dynamics of the vortex breakdown motion was dramatically visualized by tracking the time-dependent motion of particles released near the delta wing apex. This numerical visualization is analogous to experimental smoke flow techniques. In one of the dynamic cases the breakdown was found to move off the wing, convected downwind of the trailing edge, and later reformed near the trailing edge through an instability of the vortex core. A damped free-to-roll motion was also computed by releasing the wing from rest at 40 degrees of roll. The wing went to the same trim point as in the experiment.

Chaderjian, Neal M.↗

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Chapter 1 briefly reviews several related topics associated with the symmetrization of systems of conservation laws and quasi-conservation laws: (1) Basic Entropy Symmetrization Theory; (2) Symmetrization and eigenvector scaling; (3) Symmetrization of the compressible Navier-Stokes equations; and (4) Symmetrization of the quasi-conservative form of the magnetohydrodynamic (MHD) equations. Chapter 2 describes one of the best known tools employed in the study of differential equations, the maximum principle: any function f(x) which satisfies the inequality f(double prime)>0 on the interval [a,b] attains its maximum value at one of the endpoints on the interval. Chapter three examines the upwind finite volume schemes for scalar and system conservation laws. The basic tasks in the upwind finite volume approach have already been presented: reconstruction, flux evaluation, and evolution. By far, the most difficult task in this process is the reconstruction step.

Bart, Timothy J.↗

Incompressible Navier-Stokes Simulation Procedure for a Wingtip Vortex Flow Analysis

The pacing items to reach the highly desirable goal of obtaining computationally accurate flow simulation of a wingtip vortex include; super-computer development, solver accuracy, grid generation and turbulence modeling. In these four areas, many advances have been made but the fact remains that most wing computations are, at best, five percent accurate (in drag coefficient, for example). This level of accuracy has been sufficient for many purposes such as airfoil design, rudimentary wing design, and some forms of optimization. However, this accuracy level will not allow commercial aircraft designers to extract the remaining few percent of efficiency theoretically possible for conventional aircraft configurations. Thus further research is needed, particularly in the areas of solver development and turbulence modeling, to advance the state of the art of viscous computational techniques as applied to problems in aerodynamics. During the course of this study, a substantial amount of measured and computed results have been acquired. In this paper, only a small selection of experimental and computational results will be presented. This paper will outline and discuss a simulation procedure for a wingtip vortex flow analysis using the method of artificial compressibility to solve the three-dimensional, incompressible, Navier-Stokes equations (INS3D-UP).

Jennifer Dacles-Mariani↗

Incompressible Navier-Stokes, GMRES, and Multi-Element Airfoils

The talk will discuss current work in developing and applying the INS2D and INS3D flow solvers to problems in high-lift aerodynamics, primarily multi-element airfoils. These flow solvers solve the incompressible Navier-Stokes equations using the method of artificial compressibility. High-lift system flowfields are perhaps the most challenging aeronautical configurations for CFD. Difficulties with convergence for multi-element configurations using fine-grid overset meshes has led to research in using Krylov-space iterative matrix solvers. In particular, a Generalized Minimum Residual (GMRES) solver has found to be a dramatic improvement over a previously used Gauss-Seidal line-relaxation solver.

Rogers, Stuart↗

RANS-MP: A Portable Parallel Navier-Stokes Solver

RANS-MP, a new implementation of a single-grid Navier-Stokes solver using the diagonalized Beam-Warming approximate-factorization scheme, is presented. This first release of the completely rewritten solver employs the following optimizations: (1) Bi-directional multi-partition method for the ADI solver part; this improves granularity and load balance; (2) Improved cache usage through elimination of non-unit-stride array access (possible in part due to multi-partitioning); (3) Preprocessing of communicating boundary conditions to streamline logic during time stepping; (4) Truly parallel, high-performance I/O using the newly-developed MPI-IO library; (5) Elimination of large amounts of redundant operations through efficient use of workspace. Results of some realistic wing computations on the IBM SP2 computer will be presented. We will demonstrate that excellent absolute performance and scalability are obtained with RANS-MP, even for relatively small grid sizes. Besides high performance, an outstanding feature of RANS-MP is its true portability, due to the use of the portable message passing and I/O libraries MPI and MPI-IO.

VanderWijngaart, Rob F.↗

Navier-Stokes Simulation of the Canard-Wing-Body Longitudinal Dynamic Stability Characteristics

Many modern aircraft are canard-configured for aircraft control and improved aerodynamic performance. Canards can often enhance aircraft cruise performance, maneuverability and agility. For close-coupled canard configurations, the aerodynamic interaction between the canard and wing significantly changes the flow characteristics of the wing. In unsteady flow, such changes in the flow structure and performance of wings can be quite pronounced. Accurate modeling of the unsteady aerodynamics is essential for potential CFD design and analysis of such configurations. A time-accurate numerical simulation is performed to study the unsteady aerodynamic interaction between a canard and wing with emphasis on the effects of the canard on the configuration's dynamic response characteristics. The thin-layer Reynolds-averaged Navier-Stokes Equations with various turbulence models are used in this study. Computations are made on a generic, analytically-defined, close-coupled canard-wing-body configuration which has been the subject of numerous previously published experimental studies during the 1970's to mid-80's. More recently, a series of steady-flow simulations has been performed and published by the author. In the current study, the configuration is given prescribed ramp and oscillatory motions in order to predict characteristics such as the damping-in-pitch and oscillatory longitudinal stability parameters. The current computations are made at high-subsonic and transonic Mach numbers, moderate angles-of- attack from -4 to 20 degrees, and at various pitch rates and reduced frequencies. Comparisons of pressures and integrated force quantities (e.g. lift, drag, pitching moment and selected dynamic parameters) are made with other published computational results and available experimental data. Results showing the unsteady effects of the canard on surface pressures, integrated forces, canard-wing vortex interaction and vortex breakdown will be presented.

Tu, Eugene L.↗

Incompressible Navier-Stokes Computations in Aerospace and Beyond

Challenges in applying CFD to real world problems will be discussed using incompressible flow examples. Despite the remarkable progress made in CFD technology during the past two decades, new challenges are ahead of us. Some of these will be discussed in conjunction with the incompressible Navier-Stokes solver development and their applications. Computed examples include wing tip vortex formation and propagation, flow simulation of an advanced rocket pump, and an extension of the similar technology to biofluid analysis and design. Numerical issues dealing with these applications will be discussed both from algorithm development and from application point of view. The material will be printed in a bound volume after the symposium.

Kwak, Dochan↗

Simulation of Rotary-Wing Near-Wake Vortex Structures Using Navier-Stokes CFD Methods

This paper will use high-resolution Navier-Stokes computational fluid dynamics (CFD) simulations to model the near-wake vortex roll-up behind rotor blades. The locations and strengths of the trailing vortices will be determined from newly-developed visualization and analysis software tools applied to the CFD solutions. Computational results for rotor nearwake vortices will be used to study the near-wake vortex roll up for highly-twisted tiltrotor blades. These rotor blades typically have combinations of positive and negative spanwise loading and complex vortex wake interactions. Results of the computational studies will be compared to vortex-lattice wake models that are frequently used in rotorcraft comprehensive codes. Information from these comparisons will be used to improve the rotor wake models in the Tilt-Rotor Acoustic Code (TRAC) portion of NASA's Short Haul Civil Transport program (SHCT). Accurate modeling of the rotor wake is an important part of this program and crucial to the successful design of future civil tiltrotor aircraft. The rotor wake system plays an important role in blade-vortex interaction noise, a major problem for all rotorcraft including tiltrotors.

Kenwright, David↗

Time-Accurate Solutions of Incompressible Navier-Stokes Equations for Potential Turbopump Applications

Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two method are compared by obtaining unsteady solutions for the evolution of twin vortices behind a at plate. Calculated results are compared with experimental and other numerical results. For an un- steady ow which requires small physical time step, pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kiris, Cetin↗

Parametric Study of a YAV-8B Harrier in Ground Effect Using Time-Dependent Navier-Stokes Computations

Flow simulations using the time-dependent Navier-Stokes equations remain a challenge for several reasons. Principal among them are the difficulty to accurately model complex flows, and the time needed to perform the computations. A parametric study of such complex problems is not considered practical due to the large cost associated with computing many time-dependent solutions. The computation time for each solution must be reduced in order to make a parametric study possible. With successful reduction of computation time, the issue of accuracy, and appropriateness of turbulence models will become more tractable.

Shishir, Pandya↗

Numerical Experiment with Time and Spatial Accuracy of Navier-Stokes Computation For Helicopter Problems

Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other blades. This interaction, known as blade-vortex interactions (BVI), results in unsteady loading of the blades and can cause a distinctive acoustic signature. Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interactions can help rotor designers and engineers to predict rotor performance and to develop designs for low acoustic signature. Such a predictive method must preserve a blade's shed vortex for several blade revolutions before being dissipated. A number of researchers have explored the requirements for this task. This paper will outline some new capabilities that have been added to the NASA Ames' OVERFLOW code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these improvements, a number of case studies will be presented. These case studies consist of free convection of a 2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a full 3-D unsteady viscous solution of a helicopter rotor in forward flight In this study both central and upwind difference schemes are modified to be more accurate. Central difference scheme is chosen for this simulation because the flowfield is not dominated by strong shocks. The feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interaction. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in fully-implicit manner at each time-step. The spatial accuracy is either second and fourth-order central difference or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and demonstrate the methods for several sample cases and for a helicopter rotor. Preliminary computations on a rotor were performed by using this method and are in the process of documentation.

Ahmad, Jasim↗

Turbulent Navier-Stokes Flow Analysis of an Advanced Semispan Diamond-Wing Model in Tunnel and Free Air at High-Lift Conditions

Turbulent Navier-Stokes computational results are presented for an advanced diamond wing semispan model at low-speed, high-lift conditions. The numerical results are obtained in support of a wind-tunnel test that was conducted in the National Transonic Facility at the NASA Langley Research Center. The model incorporated a generic fuselage and was mounted on the tunnel sidewall using a constant-width non-metric standoff. The computations were performed at to a nominal approach and landing flow conditions.The computed high-lift flow characteristics for the model in both the tunnel and in free-air environment are presented. The computed wing pressure distributions agreed well with the measured data and they both indicated a small effect due to the tunnel wall interference effects. However, the wall interference effects were found to be relatively more pronounced in the measured and the computed lift, drag and pitching moment. Although the magnitudes of the computed forces and moment were slightly off compared to the measured data, the increments due the wall interference effects were predicted reasonably well. The numerical results are also presented on the combined effects of the tunnel sidewall boundary layer and the standoff geometry on the fuselage forebody pressure distributions and the resulting impact on the configuration longitudinal aerodynamic characteristics.

Ghaffari, Farhad↗

Incompressible Navier-Stokes Solvers in Primative Variables and their Applications to Steady and Unsteady Flow Simulations

This paper reviews recent progress made in incompressible Navier-Stokes simulation procedures and their application to problems of engineering interest. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. A summary of efforts in flow solver development is given followed by numerical studies of a few example problems of current interest. Both steady and unsteady solution algorithms and their salient features are discussed. Solvers discussed here are based on a structured-grid approach using either a finite -difference or a finite-volume frame work. As a grand-challenge application of these solvers, an unsteady turbopump flow simulation procedure has been developed which utilizes high performance computing platforms. In the paper, the progress toward the complete simulation capability of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of two parallel computing algorithms that have been implemented in the INS3D code. The relative motion of the grid systems for the rotorstator interaction was obtained using overact grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on SCSI Origin 3000 systems at NASA Ames Research Center. The same procedure has been extended to the development of NASA-DeBakey Ventricular Assist Device (VAD) that is based on an axial blood pump. Computational, and clinical analysis of this device are presented.

Kiris, Cetin C.↗