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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

A Discrete Analysis of Non-reflecting Boundary Conditions for Discontinuous Galerkin Method

We present a discrete analysis of non-reflecting boundary conditions for the discontinuous Galerkin method. The boundary conditions considered in this paper include the recently proposed Perfectly Matched Layer absorbing boundary condition for the linearized Euler equation and two non-reflecting boundary conditions based on the characteristic decomposition of the flux on the boundary. The analyses for the three boundary conditions are carried out in a unifled way. In each case, eigensolutions of the discrete system are obtained and applied to compute the numerical reflection coefficients of a specified out-going wave. The dependencies of the reflections at the boundary on the out-going wave angle and frequency as well as the mesh sizes arc? studied. Comparisons with direct numerical simulation results are also presented.

Hu, Fang Q.↗
NASA NTRS2008

Computational Investigation of the NASA Cascade Cyclonic Separation Device

Devices designed to replace the absent buoyancy separation mechanism within a microgravity environment are of considerable interest to NASA as the functionality of many spacecraft systems are dependent on the proper sequestration of interpenetrating gas and liquid phases. Inasmuch, a full multifluid Euler-Euler computational fluid dynamics investigation has been undertaken to evaluate the performance characteristics of one such device, the Cascade Cyclonic Separator, across a full range of inlet volumetric quality with combined volumetric injection rates varying from 1 L/min to 20 L/min. These simulations have delimited the general modes of operation of this class of devices and have proven able to describe the complicated vortex structure and induced pressure gradients that arise. The computational work has furthermore been utilized to analyze design modifications that enhance the overall performance of these devices. The promising results indicate that proper CFD modeling may be successfully used as a tool for microgravity separator design.

Hoyt, Nathaniel C.↗
NASA NTRS1988

A stiffly-stable implicit Runge-Kutta algorithm for CFD applications

A stiffly-stable implicit Runge-Kutta integration algorithm is derived for CFD applications spanning the range of semidiscrete theories. The algorithm family contains the one-step 'theta' algorithms, including backwards Euler and the trapezoidal rule, and provides a versatile framework to identify expressions governing algorithm stability characteristics. Parameters of a Runge-Kutta optimal implicit algorithm, second-order accurate in time and stiffly-stable, are established. This algorithm is implemented within a weak statement finite element semidiscrete formulation for one- and two-dimensional conservation law systems. Numerical results are compared to theta-algorithm solutions, for unsteady quasi-one-dimensional Euler predictions with shocks, and for a specially derived two-dimensional conservation law system modeling the Euler equations.

Baker, A. J.↗
NASA NTRS1997

Computation of Transonic Flows Using Potential Methods

The proposed paper will describe the state of the art associated with numerical solution of the full or exact velocity potential equation for solving transonic, external-aerodynamic flows. The presentation will begin with a review of the literature emphasizing research activities of the past decade. Next, the various forms of the full or exact velocity potential equation, the equation's corresponding mathematical characteristics, and the derivation assumptions will be presented and described in detail. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, will be presented and discussed relative to the more complete Euler or Navier-Stokes formulations. The technical presentation will continue with a description of recently developed full potential numerical approach characteristics. This description will include governing equation nondimensionalization, physical-to-computational-domain mapping procedures, a limited description of grid generation requirements, the spatial discretization scheme, numerical implementation of boundary conditions, and the iteration scheme. The next portion of the presentation will present and discuss numerical results for several two- and three-dimensional aerodynamic applications. Included in the results section will be a discussion and demonstration of a typical grid refinement analysis for determining spatial convergence of the numerical solution and level of solution accuracy. Computer timings for a variety of full potential applications will be compared and contrasted with similar results for the Euler equation formulation. Finally. the presentation will end with concluding remarks and recommendations for future work.

Hoist, Terry L.↗
NASA NTRS1989

Advances in upwind relaxation methods

Numerical techniques for solving the compressible Euler and Navier-Stokes equations are discussed with an emphasis on characteristic-based schemes. Two popular approaches, flux difference splitting and flux vector splitting, are described in one-dimensional Cartesian coordinates and then extended to three-dimensional generalized coordinates. A technique for increasing the spatial accuracy is presented, followed by a discussion of numerical dissipation mechanisms. An introduction to the use of implicit time integration schemes for accelerating the convergence rate to steady-state solutions including Newton's method, relaxation strategies, and approximate factorization techniques and their implementation on a vector processor concludes the chapter.

Walters, R. W.↗
NASA NTRS1991

Distributed parameter estimation for NASA Mini-Mast truss through displacement measurements

Most methods of system identification of large flexible structures by far are based on the lumped parameter approach. Because of the considerable computational burden due to the large number of unknown parameters, distributed parameter approach, which greatly decreases the number of unknowns, has being investigated. In this paper a distributed parameter model for the estimation of modal characteristics of NASA Mini-Mast truss has been formulated. Both Bernoulli-Euler beam and Timoshenko beam equations are used to characterize the lateral bending vibrations of the truss. The measurement of the lateral displacement at the tip of the truss is provided to the maximum likelihood estimator. Closed-form solutions of the partial differential equations and closed-form expressions of the sensitivity functions are derived so that the estimation algorithm is highly efficient. The resulting estimates from test data by using Timoshenko beam model are found to be comparable to those derived from finite element analysis.

Huang, Jen-Kuang↗
NASA NTRS2005

The Apollo Capsule Optimization for Improved Stability and Computational/Experimental Data Comparisons

Numerical optimization was employed on the Apollo Command Module to modify its external shape. The Apollo Command Module (CM) that was used on all NASA human space flights during the Apollo Space Program is stable and trimmed in an apex forward (alpha of approximately 40 to 80 degrees) position. This poses a safety risk if the CM separates from the launch tower during abort. Optimization was employed on the Apollo CM to remedy the undesirable stability characteristics of the configuration. Geometric shape changes were limited to axisymmetric modifications that altered the radius of the apex (R(sub A)), base radius (R(sub O)), corner radius (R(sub C)), and the cone half angle (theta), while the maximum diameter of the CM was held constant. The results of multipoint optimization on the CM indicated that the cross-range performance can be improved while maintaining robust apex-aft stability with a single trim point. Navier-Stokes computations were performed on the baseline and optimized configurations and confirmed the Euler-based optimization results. Euler Analysis of ten alternative CM vehicles with different values of the above four parameters are compared with the published experimental results of numerous wind tunnel tests during the late 1960's. These comparisons cover a wide Mach number range and a full 180-degree pitch range and show that the Euler methods are capable of fairly accurate force and moment computations and can separate the vehicle characteristics of these ten alternative configurations.

Cliff, Susan E.↗
NASA NTRS1992

Perturbation of a normal shock - A test case for unsteady Euler codes

Shock tracking capabilities of modern nonoscillatory schemes for solving the Euler equations are studied by comparing their solutions to solutions obtained using the method of characteristics on two test problems. Both test problems involve accelerating and decelerating shock motions. The schemes tested include first order upwind, second order SONIC-A, third order ENO, and third order TONIC-A. Results show that the second and higher order schemes give nearly identical shock trajectories that are closer to the characteristic solution than the first order upwind scheme. All schemes approach the characteristic solution upon mesh refinement.

Suresh, Ambady↗
NASA NTRS1993

On the static stability of forward swept propfans

An hybrid Euler solver coupled with a NASTRAN structural analysis is used to investigate the static stability characteristics of two forward swept propfan blades. The designs, having same geometry but different structural properties - one being statically stable and the other unstable, are analyzed. The stable design takes more iterations to converge than previously designed aft swept blades. The unstable design diverges to an unrealistic shape within five iterations.

Srivastava, R.↗
NASA NTRS1986

Numerical experiments with a symmetric high-resolution shock-capturing scheme

Characteristic-based explicit and implicit total variation diminishing (TVD) schemes for the two-dimensional compressible Euler equations have recently been developed. This is a generalization of recent work of Roe and Davis to a wider class of symmetric (non-upwind) TVD schemes other than Lax-Wendroff. The Roe and Davis schemes can be viewed as a subset of the class of explicit methods. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. In a recent paper, a comparison of a linearized form of the present implicit symmetric TVD scheme with an implicit upwind TVD scheme originally developed by Harten and modified by Yee was given. Results favored the symmetric method. It was found that the latter is just as accurate as the upwind method while requiring less computational effort. Currently, more numerical experiments are being conducted on time-accurate calculations and on the effect of grid topology, numerical boundary condition procedures, and different flow conditions on the behavior of the method for steady-state applications. The purpose here is to report experiences with this type of scheme and give guidelines for its use.

Yee, H. C.↗
NASA NTRS1997

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

Baysal, Oktay↗
NASA NTRS1993

Acceleration of convergence and spectrum transformation of implicit finite difference operators associated with Navier-Stokes equations

Implicit difference operator spectra are presently computed by applying eigensystem analysis techniques to finite-difference formulations of 2D Euler and Navier-Stokes equations, and attention is given to these iterative methods' convergence and stability characteristics by taking into account the effects of grid geometry, time-step, numerical viscosity, and boundary conditions. On the basis of the eigenvalue distributions for various flow configurations, the feasibility of applying such convergence-acceleration techniques as eigenvalue annihilation and relaxation is discussed. Spectrum-shifting is applied to NASA-Ames' ARC2D flow code, achieving a 20-33 percent efficiency.

Saleem, M.↗
NASA NTRS1995

Efficient implementation of weighted ENO schemes

In this paper, we further analyze, test, modify and improve the high order WENO (weighted essentially non-oscillatory) finite difference schemes of Liu, Osher and Chan. It was shown by Liu et al. that WENO schemes constructed from the r-th order (in L1 norm) ENO schemes are (r+1)-th order accurate. We propose a new way of measuring the smoothness of a numerical solution, emulating the idea of minimizing the total variation of the approximation, which results in a 5-th order WENO scheme for the case r = 3, instead of the 4-th order with the original smoothness measurement by Liu et al. This 5-th order WENO scheme is as fast as the 4-th order WENO scheme of Liu et al., and both schemes are about twice as fast as the 4-th order ENO schemes on vector supercomputers and as fast on serial and parallel computers. For Euler systems of gas dynamics, we suggest computing the weights from pressure and entropy instead of the characteristic values to simplify the costly characteristic procedure. The resulting WENO schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights, and work well for problems which do not contain strong shocks or strong reflected waves. We also prove that, for conservation laws with smooth solutions, all WENO schemes are convergent. Many numerical tests, including the 1D steady state nozzle flow problem and 2D shock entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the WENO scheme using the new smoothness measurement, in resolving complicated shock and flow structures. We have also applied Yang's artificial compression method to the WENO schemes to sharpen contact discontinuities.

Jiang, Guang-Shan↗
NASA NTRS1991

Hi-alpha forebody design. Part 2: Determination of body shapes for positive directional stability

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low speed high angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows of the F-5A forebody and Erickson forebody were recomputed with better and refined grids. The results were obtained using a modified version of cfl3d to solve either the Euler equations or the Reynolds equations employing a form of the Baldwin-Lomax turbulence model. Based on those results, we conclude that current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody model is described, and a systematic study of forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. A novel way of presenting the results is used to illustrate how the positive contribution arises. Based on the results of this initial parametric study, some guidelines for aerodynamic design to promote positive directional stability are presented.

Ravi, R.↗
NASA NTRS1991

Developing and utilizing an Euler computational method for predicting the airframe/propulsion effects for an aft-mounted turboprop transport. Volume 2: User guide

This manual explains how to use an Euler based computational method for predicting the airframe/propulsion integration effects for an aft-mounted turboprop transport. The propeller power effects are simulated by the actuator disk concept. This method consists of global flow field analysis and the embedded flow solution for predicting the detailed flow characteristics in the local vicinity of an aft-mounted propfan engine. The computational procedure includes the use of several computer programs performing four main functions: grid generation, Euler solution, grid embedding, and streamline tracing. This user's guide provides information for these programs, including input data preparations with sample input decks, output descriptions, and sample Unix scripts for program execution in the UNICOS environment.

Chen, H. C.↗
NASA NTRS1986

Unsteady transonic flows past airfoils using the Euler equations

An explicit finite difference solver for the Euler equations was developed and used to calculate the flow past two of the AGARD standard aeroelastic configurations, an NACA 64A010 and an NLR 7301. The algorithm employed uses a modified four step Runge-Kutta time stepping scheme, boundary conditions determined from the time-dependent theory of characteristics and an unsteady automatic grid generation procedure. In general, the calculated results are in good agreement with available experiment. Moreover, they demonstrate the importance of using the Euler equations for super-critical airfoils and in situations where the transonic small disturbance theory results in poor agreement with experiment.

Smith, G. E.↗
NASA NTRS1988

Application of Navier-Stokes analysis to stall flutter

A solution procedure was developed to investigate the two-dimensional, one- or two-dimensional flutter characteristics of arbitrary airfoils. This procedure requires a simultaneous integration in time of the solid and fluid equations of motion. The fluid equations of motion are the unsteady compressible Navier-Stokes equations, solved in a body-fitted moving coordinate system using an approximate factorization scheme. The solid equations of motion are integrated in time using an Euler implicit scheme. Flutter is said to occur if small disturbances imposed on the airfoil attitude lead to divergent oscillatory motions at subsequent times. The flutter characteristics of airfoils in subsonic speed at high angles of attack and airfoils in high subsonic and transonic speeds at low angles of attack are investigated. The stall flutter characteristics are also predicted using the same procedure.

Wu, J. C.↗
NASA NTRS1978

An old integration scheme for compressible flows revisited, refurbished and put to work

A scheme for integrating the Euler equations of compressible flow in any hyperbolic case is presented. The scheme relies on the concept of characteristics but is strictly a finite difference scheme. Improvements in accuracy and physical consistence due to the scheme are discussed and results of its application to complex flows are shown.

Moretti, G.↗