Search NASA⌕ Search

SEARCH · Search NASA

Results for “Euler characteristic”

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 109 records · Page 6

Numerical solution of the two-dimensional time-dependent incompressible Euler equations

A numerical method is presented for solving the artificial compressibility form of the 2D time-dependent incompressible Euler equations. The approach is based on using an approximate Riemann solver for the cell face numerical flux of a finite volume discretization. Characteristic variable boundary conditions are developed and presented for all boundaries and in-flow out-flow situations. The system of algebraic equations is solved using the discretized Newton-relaxation (DNR) implicit method. Numerical results are presented for both steady and unsteady flow.

Whitfield, David L.↗

A Three-Dimensional Linearized Unsteady Euler Analysis for Turbomachinery Blade Rows

A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic response characteristics of axial-flow turbomachinery blading. The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. In addition, a numerical model for linearized inviscid unsteady flow, which is based upon an existing nonlinear, implicit, wave-split, finite volume analysis, is described. These aerodynamic and numerical models have been implemented into an unsteady flow code, called LINFLUX. A preliminary version of the LINFLUX code is applied herein to selected, benchmark three-dimensional, subsonic, unsteady flows, to illustrate its current capabilities and to uncover existing problems and deficiencies. The numerical results indicate that good progress has been made toward developing a reliable and useful three-dimensional prediction capability. However, some problems, associated with the implementation of an unsteady displacement field and numerical errors near solid boundaries, still exist. Also, accurate far-field conditions must be incorporated into the FINFLUX analysis, so that this analysis can be applied to unsteady flows driven be external aerodynamic excitations.

Montgomery, Matthew D.↗

A One Dimensional, Time Dependent Inlet/Engine Numerical Simulation for Aircraft Propulsion Systems

The NASA Lewis Research Center (LeRC) and the Arnold Engineering Development Center (AEDC) have developed a closely coupled computer simulation system that provides a one dimensional, high frequency inlet/engine numerical simulation for aircraft propulsion systems. The simulation system, operating under the LeRC-developed Application Portable Parallel Library (APPL), closely coupled a supersonic inlet with a gas turbine engine. The supersonic inlet was modeled using the Large Perturbation Inlet (LAPIN) computer code, and the gas turbine engine was modeled using the Aerodynamic Turbine Engine Code (ATEC). Both LAPIN and ATEC provide a one dimensional, compressible, time dependent flow solution by solving the one dimensional Euler equations for the conservation of mass, momentum, and energy. Source terms are used to model features such as bleed flows, turbomachinery component characteristics, and inlet subsonic spillage while unstarted. High frequency events, such as compressor surge and inlet unstart, can be simulated with a high degree of fidelity. The simulation system was exercised using a supersonic inlet with sixty percent of the supersonic area contraction occurring internally, and a GE J85-13 turbojet engine.

Doug Garrard↗

Three-dimensional unsteady Euler equations solution using flux vector splitting

A method for numerically solving the three-dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of 2 and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far-field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three-dimensional Euler code and experimental data.

Whitfield, D. L.↗

Euler code evaluation of a conical vehicle at supersonic speeds

A study has been conducted on a conical wing-body transatmospheric vehicle at Mach numbers from 1.5 to 4.5. The objective of the study was to evaluate the application of a Euler code to a wing-body configuration across the supersonic Mach number range. The baseline wing-body fuselage consisted of a circular 5-deg half-angle conical forebody, a cylindrical midbody, and a 9-deg truncated conical boattail. The wing was a 4-percent-thick diamond airfoil wing of aspect ratio 1. Results show that the Euler code accurately predicted the lift, drag, and pitching moment coefficients as a function of Mach number and angle of attack. The longitudinal stability characteristics were predicted as a function of Mach number. The code also predicted the changes in lift, drag, and pitching moment for forward and rearward shift of the wing from the nominal position.

Mcgrath, Brian E.↗

An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case

Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.

Goodrich, John W.↗

Transformation of two and three-dimensional regions by elliptic systems

Efforts in transferring computational work from the LRC computer to the IRIS Graphics Workstation at MSU are reported and the computation of a conservative solution of a simple hyperbolic equation on an overlapping grid is discussed. Several conclusions concerning computations on overlapping grids are apparent. Problems only occur when there is a major difference in grid spacing on the individual component grids. In the case of hyperbolic equations, it is necessary that both interpolation and extrapolation be applied at the grid boundaries. When interpolated values are used at outflow boundary points, excessive oscillations in the numerical solution may be the result. The same conclusions would be valid for more complicated systems of hyperbolic equations such as the Euler equations for inviscid flow. Some of the solution values would be extrapolated at the overlap boundary, the exact number depending on the number of characteristics pointing out of the overlap region. It is also possible that similar boundary conditions may be needed for some parabolic equations such as high Reynolds number viscous flow equations. Efforts were also expended on the development of three-dimensional conservative interpolation procedures. Finally, the investigation of grid smoothing procedures were initiated during this reporting period. It was decided that the first grid smoothing algorithms will be based on the concepts of variational grid generation.

Mastin, C. W.↗

Hi-alpha forebody design. Part 1: Methodology base and initial parametrics

The use of Computational Fluid Dynamics (CFD) has been investigated for the analysis and design of aircraft forebodies at high angle of attack combined with sideslip. The results of the investigation show that CFD has reached a level of development where computational methods can be used for high angle of attack aerodynamic design. The classic wind tunnel experiment for the F-5A forebody directional stability has been reproduced computationally over an angle of attack range from 10 degrees to 45 degrees, and good agreement with experimental data was obtained. Computations have also been made at combined angle of attack and sideslip over a chine forebody, demonstrating the qualitative features of the flow, although not producing good agreement with measured experimental pressure distributions. The computations were performed using the code known as cfl3D for both the Euler equations and the Reynolds equations using a form of the Baldwin-Lomax turbulence model. To study the relation between forebody shape and directional stability characteristics, a generic parametric forebody model has been defined which provides a simple analytic math model with flexibility to capture the key shape characteristics of the entire range of forebodies of interest, including chines.

Mason, William H.↗

A Multi-Level Parallelization Concept for High-Fidelity Multi-Block Solvers

The integration of high-fidelity Computational Fluid Dynamics (CFD) analysis tools with the industrial design process benefits greatly from the robust implementations that are transportable across a wide range of computer architectures. In the present work, a hybrid domain-decomposition and parallelization concept was developed and implemented into the widely-used NASA multi-block Computational Fluid Dynamics (CFD) packages implemented in ENSAERO and OVERFLOW. The new parallel solver concept, PENS (Parallel Euler Navier-Stokes Solver), employs both fine and coarse granularity in data partitioning as well as data coalescing to obtain the desired load-balance characteristics on the available computer platforms. This multi-level parallelism implementation itself introduces no changes to the numerical results, hence the original fidelity of the packages are identically preserved. The present implementation uses the Message Passing Interface (MPI) library for interprocessor message passing and memory accessing. By choosing an appropriate combination of the available partitioning and coalescing capabilities only during the execution stage, the PENS solver becomes adaptable to different computer architectures from shared-memory to distributed-memory platforms with varying degrees of parallelism. The PENS implementation on the IBM SP2 distributed memory environment at the NASA Ames Research Center obtains 85 percent scalable parallel performance using fine-grain partitioning of single-block CFD domains using up to 128 wide computational nodes. Multi-block CFD simulations of complete aircraft simulations achieve 75 percent perfect load-balanced executions using data coalescing and the two levels of parallelism. SGI PowerChallenge, SGI Origin 2000, and a cluster of workstations are the other platforms where the robustness of the implementation is tested. The performance behavior on the other computer platforms with a variety of realistic problems will be included as this on-going study progresses.

Hatay, Ferhat F.↗

Foundations for the numerical solution of the Euler equations

The Navier-Stokes equations represent an extremely good model of the physical phenomena encountered in most aeronautical problems. However, the computational resource needed to solve the Navier-Stokes equations are so large that even with today's supercomputers, it is necessary to make use of simpler models. A large number of external aerodynamic problems can be accurately described by a simpler model. This model consists of an outer inviscid flow plus a boundary-layer thickness correction for the vehicle shape. The outer inviscid model may be represented by the potential equation or by the Euler equation. The present paper provides the foundations for the numerical solution of the Euler equations. The governing equations are considered, taking into account conservation laws, the medium, the differential form of the conservation laws, generalized solutions, shock-fitting, and characteristics. Attention is also given to initial and boundary conditions, existence and uniqueness, and rotational phenomena.

Salas, M. D.↗

Compressive behavior of unidirectional fibrous composites

The longitudinal compressive behavior of unidirectional fiber composites was investigated by using the Illinois Institute of Technology Research Institute (IITRI) test method with thick and thin test specimens. The test data obtained are interpreted by means of stress/strain curves from back-to-back strain gages, examination of fracture surfaces by scanning electron microscope, and predictive equations for distinct failure modes including fiber compression failure. Euler buckling, delamination, and flexure. The results show that longitudinal compressive fracture is induced by a combination of delamination, flexure, and fiber tier breaks. No distinct fracture surface characteristics can be associated with unique failure modes. An equation is described that can be used to extract the longitudinal compressive strength from the longitudinal tensile and flexural strengths of the same composite system.

Sinclair, J. H.↗

Compression behavior of unidirectional fibrous composite

The longitudinal compression behavior of unidirectional fiber composites is investigated using a modified Celanese test method with thick and thin test specimens. The test data obtained are interpreted using the stress/strain curves from back-to-back strain gages, examination of fracture surfaces by scanning electron microscope, and predictive equations for distinct failure modes including fiber compression failure, Euler buckling, delamination, and flexure. The results show that the longitudinal compression fracture is induced by a combination of delamination, flexure, and fiber tier breaks. No distinct fracture surface characteristics can be associated with unique failure modes. An equation is described which can be used to extract the longitudinal compression strength knowing the longitudinal tensile and flexural strengths of the same composite system.

Sinclair, J. H.↗

A zonal approach to solution of the Euler equations

A technique for the solution of the one- and two-dimensional Euler equations in a partitioned flow field is presented. The field is divided into distinct 'zones', each of which is computed separately. An implicit boundary procedure based on the characteristic propagaion of information and flux splitting methods is applied at the zonal interfaces. Numerical results are presented for both quasi-one-dimensional nozzle flows with shock waves and the unsteady shock-tube problem. These calculations demonstrate the capability of shock propagation through arbitrarily located zonal boundaries in a stable, conservative, and accurate manner. Two-dimensional results include the zonal computation of flows over blunt bodies and airfoils.

Hessenius, K. A.↗

Conical Euler analysis and active roll suppression for unsteady vortical flows about rolling delta wings

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling, highly swept delta wings undergoing either forced motions or free-to-roll motions that include active roll suppression. The flow solver of the code involves a multistage, Runge-Kutta time-stepping scheme that uses a cell-centered, 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 simultaneously integrating in time the rigid-body equation of motion with the governing flow equations. Results are presented for a delta wing with a 75 deg swept, sharp leading edge at a free-stream Mach number of 1.2 and at 10 deg, 20 deg, and 30 deg angle of attack alpha. At the lower angles of attack (10 and 20 deg), forced-harmonic analyses indicate that the rolling-moment coefficients provide a positive damping, which is verified by free-to-roll calculations. In contrast, at the higher angle of attack (30 deg), a forced-harmonic analysis indicates that the rolling-moment coefficient provides negative damping at the small roll amplitudes. A free-to-roll calculation for this 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, which is characteristic of highly swept delta wings. This limit cycle oscillation may be actively suppressed through the use of a rate-feedback control law and antisymmetrically deflected leading-edge flaps. Descriptions of the conical Euler flow solver and the free-to roll analysis are included in this report. Results are presented that demonstrate how the systematic analysis of the forced response of the delta wing can be used to predict the stable, neutrally stable, and unstable free response of the delta wing. These results also give insight into the flow physics associated with unsteady vortical flows about delta wings undergoing forced motions and free-to-roll motions, including the active suppression of the wing-rock type phenomenon. The conical Euler methodology developed is directly extend able to three-dimensional calculations.

Lee-Rausch, Elizabeth M.↗

X-33 Aerodynamic and Aeroheating Computations for Wind Tunnel and Flight Conditions

This report provides an overview of hypersonic Computational Fluid Dynamics research conducted at the NASA Langley Research Center to support the Phase II development of the X-33 vehicle. The X-33, which is being developed by Lockheed-Martin in partnership with NASA, is an experimental Single-Stage-to-Orbit demonstrator that is intended to validate critical technologies for a full-scale Reusable Launch Vehicle. As part of the development of the X-33, CFD codes have been used to predict the aerodynamic and aeroheating characteristics of the vehicle. Laminar and turbulent predictions were generated for the X 33 vehicle using two finite- volume, Navier-Stokes solvers. Inviscid solutions were also generated with an Euler code. Computations were performed for Mach numbers of 4.0 to 10.0 at angles-of-attack from 10 deg to 48 deg with body flap deflections of 0, 10 and 20 deg. Comparisons between predictions and wind tunnel aerodynamic and aeroheating data are presented in this paper. Aeroheating and aerodynamic predictions for flight conditions are also presented.

Hollis, Brian R.↗

Assessment of analytic methods for the prediction of aerodynamic characteristics of arbitrary bodies at supersonic speeds

Trends toward the automation of the design process for airplanes and missiles accentuate the need for analytic techniques for the prediction of aerodynamic characteristics. A number of computer codes have been developed or are under development which show promise of significantly improving the estimation of aerodynamic characteristics for arbitrarily-shaped bodies at supersonic speeds. The programs considered range in complexity from a simple linearized solution employing slender body theory to an exact finite difference solution of the Euler equations. The results from five computer codes are compared with experimental data to determine the accuracy, range of applicability, ease of use, and computer time and cost of the programs. The results provide a useful guide for selecting the appropriate method for treating bodies at the various levels of an automated design process.

Landrum, E. J.↗

Steady state computations for wave propagation problems

The behavior of difference approximations of hyperbolic partial differential equations as time t goes to infinity is studied. The rate of convergence to steady state is analyzed theoretically and experimentally for the advection equation and the linearized Euler equations. The choice of difference formulas and boundary conditions strongly influences the rate of convergence in practical steady-state calculations. In particular it is shown that upwind difference methods and characteristic boundary conditions have very attractive convergence properties.

Engquist, Bjorn↗

Multidisciplinary aeroelastic analysis of a generic hypersonic vehicle

This paper presents details of a flutter and stability analysis of aerospace structures such as hypersonic vehicles. Both structural and aerodynamic domains are discretized by the common finite element technique. A vibration analysis is first performed by the STARS code employing a block Lanczos solution scheme. This is followed by the generation of a linear aerodynamic grid for subsequent linear flutter analysis within subsonic and supersonic regimes of the flight envelope; the doublet lattice and constant pressure techniques are employed to generate the unsteady aerodynamic forces. Flutter analysis is then performed for several representative flight points. The nonlinear flutter solution is effected by first implementing a CFD solution of the entire vehicle. Thus, a 3-D unstructured grid for the entire flow domain is generated by a moving front technique. A finite element Euler solution is then implemented employing a quasi-implicit as well as an explicit solution scheme. A novel multidisciplinary analysis is next effected that employs modal and aerodynamic data to yield aerodynamic damping characteristics. Such analyses are performed for a number of flight points to yield a large set of pertinent data that define flight flutter characteristics of the vehicle. This paper outlines the finite-element-based integrated analysis procedures in detail, which is followed by the results of numerical analyses of flight flutter simulation.

Gupta, K. K.↗