Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvers”

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 775 records · Page 43

High speed inlet calculations with real gas effects

A 2-D steady-state Navier-Stokes solver has been upgraded to include the effects of frozen and equilibrium air chemistry for applications to high speed flight vehicles. To provide a computationally economical first order approximation to the high temperature physics, variable thermodynamic data is used for the chemically frozen mode to allow for a variation with temperature of the air specific heats and enthalpy. For calculations involving air in chemical equilibrium, a specially modified version of the NASA Lewis Chemical Equilibrium Code, CEC, is used to compute the chemical composition and resultant thermochemical properties. The upgraded solver is demonstrated by comparing results from calorically perfect (C sub p=constant), thermally perfect (frozen) and equilibrium air calculations for a variety of geometries, and flight Mach numbers.

Coirier, William J.↗

Numerical simulation methods of incompressible flows and an application to the Space Shuttle main engine

This paper discusses incompressible Navier-Stokes solution methods with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow solver code has been used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.

Chang, J. L. C.↗

Vortex-dominated conical-flow computations using unstructured adaptively-refined meshes

A conical Euler/Navier-Stokes algorithm is presented for the computation of vortex-dominated flows. The flow solver involves a multistage Runge-Kutta time stepping scheme which uses a finite-volume spatial discretization on an unstructured grid made up of triangles. The algorithm also employs an adaptive mesh refinement procedure which enriches the mesh locally to more accurately resolve the vortical flow features. Results are presented for several highly-swept delta wing and circular cone cases at high angles of attack and at supersonic freestream flow conditions. Accurate solutions were obtained more efficiently when adaptive mesh refinement was used in contrast with refining the grid globally. The paper presents descriptions of the conical Euler/Navier-Stokes flow solver and adaptive mesh refinement procedures along with results which demonstrate the capability.

Batina, John T.↗

A natural flow wing design employing 3-D nonlinear analysis applied at supersonic speeds

A wing-design study has been conducted on a 65-deg-swept leading-edge delta wing in which a near-conical geometry was employed to take advantage of the naturally occurring conical flow which arises over such a wing in a supersonic flow field. Three-dimensional nonlinear analysis methods were used in the study. In preliminary design, wing planform, design conditions, and near-conical concept were derived and a baseline standard wing (conventional airfoil distribution) and a baseline near-conical wing were chosen. During the initial analysis, a full-potential solver was employed to determine the aerodynamic characteristics of the baseline standard delta wing and the near-conical delta wing. Modifications due to airfoil thickness, leading-edge radius, and camber were then applied to the baseline near-conical wing. The final design employed a Euler solver to analyze the best wing configurations found in the initial design, and to extend this study to develop a more refined wing. Benefits due to each modification are discussed, and a final natural flow wing geometry is chosen and its aerodynamic characteristics are compared with the baseline wings.

Bauer, Steven X. S.↗

An Euler analysis of a High-Speed Civil Transport (HSCT) concept at Mach 3

A marching Euler solver, GEM3D, was used to predict the Mach 3 flow field for the wing and body of a High-Speed Civil Transport concept. The analysis focused on a typical cruise lift coefficient of 0.1 at alpha = 3 deg. The Euler solution indicated that embedded shocks formed on the upper surface of the inboard wing panel and at the leading-edge of the outboard wing panel, due to its supersonic leading edge condition. According to a simple static-pressure criterion, the embedded wing upper-surface shocks were sufficiently strong to separate a turbulent boundary layer. Comparison of aerodynamic coefficients from the Euler solver with those from linear theory shows that the linear theory estimates of lift and drag are optimistic, which would lead to optimistic estimates of cruise range.

Pittman, James L.↗

Computation of viscous incompressible flows

Incompressible Navier-Stokes solution methods and their applications to three-dimensional flows are discussed. A brief review of existing methods is given followed by a detailed description of recent progress on development of three-dimensional generalized flow solvers. Emphasis is placed on primitive variable formulations which are most promising and flexible for general three-dimensional computations of viscous incompressible flows. Both steady- and unsteady-solution algorithms and their salient features are discussed. Finally, examples of real world applications of these flow solvers are given.

Kwak, Dochan↗

Numerical simulation methods of incompressible flows and an application to the space shuttle main engine

Incompressible Navier-Stokes solution methods are discussed with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow-solver code was used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.

Chang, J. L. C.↗

Application of the TRANAIR rectangular grid approach to the aerodynamic analysis of complex configurations

A numerical method is described which uses a rectangular grid to solve the nonlinear full potential equation about complex configurations. The grid is locally refined to resolve high velocity gradients arising from leading edge expansions or shock waves. The grid penetrates the boundary (described by networks of quadrilateral panels) and is generated automatically. Discrete operators are constructed using the finite element method. The system of nonlinear discrete equations is solved iteratively using a Krylov subspace method preconditioned by an exterior Poisson solver and a direct sparse solver. The primary emphasis is to provide design engineers with an aerodynamic analysis tool (the TRANAIR code) which is accurate, reliable, economical, and flexible to use. Computational results for many interesting configurations are presented.

Johnson, Forrester T.↗

Towards a genuinely multi-dimensional upwind scheme

Methods of incorporating multi-dimensional ideas into algorithms for the solution of Euler equations are presented. Three schemes are developed and tested: a scheme based on a downwind distribution, a scheme based on a rotated Riemann solver and a scheme based on a generalized Riemann solver. The schemes show an improvement over first-order, grid-aligned upwind schemes, but the higher-order performance is less impressive. An outlook for the future of multi-dimensional upwind schemes is given.

Powell, Kenneth G.↗

The aerodynamic design of the oblique flying wing supersonic transport

The aerodynamic design of a supersonic oblique flying wing is strongly influenced by the requirement that passengers must be accommodated inside the wing. It was revealed that thick oblique wings of very high sweep angle can be efficient at supersonic speeds when transonic normal Mach numbers are allowed on the upper surface of the wing. The goals were motivated by the ability to design a maximum thickness, minimum size oblique flying wing. A 2-D Navier-Stokes solver was used to design airfoils up to 16 percent thickness with specified lift, drag and pitching moment. A new method was developed to calculate the required pressure distribution on the wing based on the airfoil loading, normal Mach number distribution and theoretical knowledge of the minimum drag of oblique configurations at supersonic speeds. The wing mean surface for this pressure distribution was calculated using an inverse potential flow solver. The lift to drag ratio of this wing was significantly higher than that of a comparable delta wing for cruise speeds up to Mach 2.

Vandervelden, Alexander J. M.↗

A package for 3-D unstructured grid generation, finite-element flow solution and flow field visualization

A set of computer programs for 3-D unstructured grid generation, fluid flow calculations, and flow field visualization was developed. The grid generation program, called VGRID3D, generates grids over complex configurations using the advancing front method. In this method, the point and element generation is accomplished simultaneously, VPLOT3D is an interactive, menudriven pre- and post-processor graphics program for interpolation and display of unstructured grid data. The flow solver, VFLOW3D, is an Euler equation solver based on an explicit, two-step, Taylor-Galerkin algorithm which uses the Flux Corrected Transport (FCT) concept for a wriggle-free solution. Using these programs, increasingly complex 3-D configurations of interest to aerospace community were gridded including a complete Space Transportation System comprised of the space-shuttle orbitor, the solid-rocket boosters, and the external tank. Flow solutions were obtained on various configurations in subsonic, transonic, and supersonic flow regimes.

Parikh, Paresh↗

A parallel finite-difference method for computational aerodynamics

A finite-difference scheme for solving complex three-dimensional aerodynamic flow on parallel-processing supercomputers is presented. The method consists of a basic flow solver with multigrid convergence acceleration, embedded grid refinements, and a zonal equation scheme. Multitasking and vectorization have been incorporated into the algorithm. Results obtained include multiprocessed flow simulations from the Cray X-MP and Cray-2. Speedups as high as 3.3 for the two-dimensional case and 3.5 for segments of the three-dimensional case have been achieved on the Cray-2. The entire solver attained a factor of 2.7 improvement over its unitasked version on the Cray-2. The performance of the parallel algorithm on each machine is analyzed.

Swisshelm, Julie M.↗

Staggered solution procedures for multibody dynamics simulation

A computational procedure is presented for the direct integration of the multibody dynamical (MBD) equations with constraints. Numerical experiments conducted with this MBD procedure indicate that it yields robust solutions in cases where the step size furnishes more than 20 samples for the period of highest apparent response frequency in a given multibody system. The MBD solution procedure is implemented in two separate modules: the generalized coordinated solver, 'CINT', and the constraint Lagrange multiplier solver, 'LINT'.

Park, K. C.↗

Experimental and predicted aerodynamic characteristics of a proposed Assured Crew Return Vehicle (ACRV) lifting-body configuration at Mach 6 and 10

The effects of Mach number, Reynolds number, and ratio of specific heats on the aerodynamic characteristics of a proposed Assured Crew Return Vehicle (ACRV) lifting-body configuration were examined for a range of angles of attack from -5 deg to 50 deg. Predictions made with a Langley-developed, three-dimensional Navier Stokes solver known as LAURA, which was exercised as an Euler solver for the present study, are compared with the experimental results. Unlike the Shuttle Orbiter, which experienced a significant nose-up increment in pitching moment with decreasing specific heat ratio (i.e., real gas effects), the aerodynamic characteristics of this lifting-body configuration are insensitive to changes in specific heat ratio. The maximum trimmed lift-to-drag ratio achieved was about 1.5. Predicted inviscid values of aerodynamic coefficients were generally in good agreement with measurement.

Micol, John R.↗

Modeling of high speed chemically reacting flow-fields

The SPARK3D and SPARK3D-PNS computer programs were developed to model 3-D supersonic, chemically reacting flow-fields. The SPARK3D code is a full Navier-Stokes solver, and is suitable for use in scramjet combustors and other regions where recirculation may be present. The SPARK3D-PNS is a parabolized Navier-Stokes solver and provides an efficient means of calculating steady-state combustor far-fields and nozzles. Each code has a generalized chemistry package, making modeling of any chemically reacting flow possible. Research activities by the Langley group range from addressing fundamental theoretical issues to simulating problems of practical importance. Algorithmic development includes work on higher order and upwind spatial difference schemes. Direct numerical simulations employ these algorithms to address the fundamental issues of flow stability and transition, and the chemical reaction of supersonic mixing layers and jets. It is believed that this work will lend greater insight into phenomenological model development for simulating supersonic chemically reacting flows in practical combustors. Currently, the SPARK3D and SPARK3D-PNS codes are used to study problems of engineering interest, including various injector designs and 3-D combustor-nozzle configurations. Examples, which demonstrate the capabilities of each code are presented.

Drummond, J. P.↗

CFD analysis for high speed inlets

The increased national interest in high speed flight has increased research for high speed propulsion components. The highly 3-D flows present in supersonic/hypersonic inlets are currently being studied at NASA-Lewis both experimentally and computationally using a family of steady Parabolized Navier-Stokes (PNS) and Navier-Stokes (NS) solvers and unsteady NS solvers. Some of the results of these efforts are presented with an emphasis on the comparison of the computational and experimental results. The flow in high speed inlets typically involves the interaction of compression shock waves and boundary layers on the internal surfaces. The fundamentals of these interactions have been studied experimentally for many years, while more recently, computations have been used to study these complex 3-D flow fields. Attempts to control the flow through boundary layer bleed are being investigated computationally prior to wind tunnel experiments. The ultimate goal is the higher performing inlets required for high speed flight.

Benson, Tom↗

Conical Euler simulation and active suppression of delta wing rocking motion

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling highly-swept delta wings, undergoing either forced or free-to-roll motions including active roll suppression. The flow solver of the code involves a multistage Runge-Kutta time-stepping scheme which uses a finite volume spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free-to-roll case, by including the rigid-body equation of motion for its simultaneous time integration with the governing flow equations. Results are presented for a 75 deg swept sharp leading edge delta wing at a freestream Mach number of 1.2 and at alpha equal to 10 and 30 deg angle of attack. A forced harmonic analysis indicates that the rolling moment coefficient provides: (1) a positive damping at the lower angle of attack equal to 10 deg, which is verified in a free-to-roll calculation; (2) a negative damping at the higher angle of attack equal to 30 deg at the small roll amplitudes. A free-to-roll calculation for the latter case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation. The wing rocking motion may be actively suppressed, however, through the use of a rate-feedback control law and antisymmetrically deflected leading edge flaps. The descriptions of the conical Euler flow solver and the free-to-roll analysis are presented. Results are also presented which give insight into the flow physics associated with unsteady vortical flows about forced and free-to-roll delta wings, including the active roll suppression of this wing-rock phenomenon.

Lee, Elizabeth M.↗

A three-dimensional application with the numerical grid generation code: EAGLE (utilizing an externally generated surface)

Program EAGLE (Eglin Arbitrary Geometry Implicit Euler) is a multiblock grid generation and steady-state flow solver system. This system combines a boundary conforming surface generation, a composite block structure grid generation scheme, and a multiblock implicit Euler flow solver algorithm. The three codes are intended to be used sequentially from the definition of the configuration under study to the flow solution about the configuration. EAGLE was specifically designed to aid in the analysis of both freestream and interference flow field configurations. These configurations can be comprised of single or multiple bodies ranging from simple axisymmetric airframes to complex aircraft shapes with external weapons. Each body can be arbitrarily shaped with or without multiple lifting surfaces. Program EAGLE is written to compile and execute efficiently on any CRAY machine with or without Solid State Disk (SSD) devices. Also, the code uses namelist inputs which are supported by all CRAY machines using the FORTRAN Compiler CF177. The use of namelist inputs makes it easier for the user to understand the inputs and to operate Program EAGLE. Recently, the Code was modified to operate on other computers, especially the Sun Spare4 Workstation. Several two-dimensional grid configurations were completely and successfully developed using EAGLE. Currently, EAGLE is being used for three-dimension grid applications.

Houston, Johnny L.↗