Search NASA⌕ Search

SEARCH · Search NASA

Results for “modified 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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

CFD Predictions of Sonic-Boom Characteristics for Unmodified and Modified SR-71 Configurations

Shaped sonic-boom signatures refer to signatures that look something other than the typical N-waves. Shaped sonic-boom signatures such as "flat-top," "ramp-type," or "hybrid-type" waveforms have been shown to reduce the subjective loudness without requiring reductions in overpressure peaks. The shaping of sonic-boom signatures requires increasing the shock rise time and changes in frequency spectra. So far, a flat-top waveform was shown to be achievable in wind tunnels; however, the influence of long propagation distance and real atmosphere on shaped signatures should be addressed using flight tests. Two different approaches have been proposed for sonic-boom minimization flight tests. The first approach, proposed by Eagle Aerospace, is for a flight test using a modified BQM-34 "FIREBEE" remotely piloted vehicle. The 30-foot long FIREBEE has a steady state flight condition at the Mach number and altitude of interest, and it can be recovered by helicopter from the water. As an alternative approach, a modified SR-71 vehicle has been proposed by the McDonnell Douglas Corporation. Benefits of the SR-71 include its variable geometry supersonic inlets, small cockpit bulge, higher Mach number capabilities, slender design, and longer length (105 foot). The present investigation addresses the sonic-boom analysis for the second vehicle.The objective of the current investigation is to assess the feasibility of a modified SR-71 configuration, with McDonnell Douglas-designed fuselage modifications, intended to produce shaped sonic-boom signatures on the ground. The present study describes the use of a higher-order computational fluid dynamics (CFD) method to predict the sonic-boom characteristics for both unmodified and modified SR-71 configurations. An Euler unstructured grid methodology is used to predict the near-field, three-dimensional pressure patterns generated by both SR-71 models. The computed near-field pressure signatures are extrapolated to specified distances below the aircraft down to impingement on the ground using the code MDBOOM. Comparisons of the near-field pressure signatures with available flight-test data are presented in the current paper.

Fouladi, Kamran↗
NASA NTRS1997

Euler/Navier-Stokes Solvers Applied to Ducted Fan Configurations

Due to noise considerations, ultra high bypass ducted fans have become a more viable design. These ducted fans typically consist of a rotor stage containing a wide chord fan and a stator stage. One of the concerns for this design is the classical flutter that keeps occurring in various unducted fan blade designs. These flutter are catastrophic and are to be avoided in the flight envelope of the engine. Some numerical investigations by Williams, Cho and Dalton, have suggested that a duct around a propeller makes it more unstable. This needs to be further investigated. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading be available. Aerodynamic solvers based on unsteady three-dimensional analysis will provide accurate and fast solutions and are best suited for aeroelastic analysis. The Euler solvers capture significant physics of the flowfield and are reasonably fast. An aerodynamic solver Ref. based on Euler equations had been developed under a separate grant from NASA Lewis in the past. Under the current grant, this solver has been modified to calculate the aeroelastic characteristics of unducted and ducted rotors. Even though, the aeroelastic solver based on three-dimensional Euler equations is computationally efficient, it is still very expensive to investigate the effects of multiple stages on the aeroelastic characteristics. In order to investigate the effects of multiple stages, a two-dimensional multi stage aeroelastic solver was also developed under this task, in collaboration with Dr. T. S. R. Reddy of the University of Toledo. Both of these solvers were applied to several test cases and validated against experimental data, where available.

Keith, Theo G., Jr.↗
NASA NTRS1993

Calculation of AGARD Wing 445.6 Flutter Using Navier-Stokes Aerodynamics

The flutter characteristics of the first AGARD standard aeroelastic configuration for dynamic response, Wing 445.6, are studied using an unsteady Navier-Stokes algorithm in order to investigate a previously noted discrepancy between Euler flutter characteristics and the experimental data. The algorithm, which is a three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1), was previously modified for the time-marching, aeroelastic analysis of wings using the unsteady Euler equations. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. In this paper, the aeroelastic method is extended and evaluated for applications that use the Navier- Stokes aerodynamics. The paper presents a brief description of the aeroelastic method and presents unsteady calculations which verify this method for Navier-Stokes calculations. A linear stability analysis and a time-marching aeroelastic analysis are used to determine the flutter characteristics of the isolated 45 deg. swept-back wing. Effects of fluid viscosity, structural damping, and number of modes in the structural model are investigated. For the linear stability analysis, the unsteady generalized aerodynamic forces of the wing are computed for a range of reduced frequencies using the pulse transfer-function approach. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. This stability analysis is used to determine the flutter characteristics of the wing at free-stream Mach numbers of 0.96 and 1.141 using the generalized aerodynamic forces generated by solving the Euler equations and the Navier-Stokes equations. Time-marching aeroelastic calculations are performed at a free-stream Mach number of 1.141 using the Euler and Navier-Stokes equations to compare with the linear V-g flutter analysis method. The V-g analysis, which is used in conjunction with the time-marching analysis, indicates that the fluid viscosity has a significant effect on the supersonic flutter boundary for this wing while the structural damping and number of modes in the structural model have a lesser effect.

Lee-Rausch, Elizabeth M.↗
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 NTRS1992

Computational fluid dynamics application to hypersonic flow over a Martian entry vehicle (MEV) - A correlation with experiment

A class of vehicles for a mission to Mars are analyzed for aerodynamic characteristics using advanced Computational Fluid Dynamics (CFD). The general configuration is a modified cone-conical-frustum geometry where the nose radius has a large influence on the flowfield. Inviscid-compressible flow using the Euler equations and viscous-compressible flow using the thin-layer Navier-Stokes equations is applied to the configuration. The surface modeling, grid generation and application of state-of-the-art CFD software are described. The effects of nose radius, angle of attack, and hypersonic velocity on the flight characteristics of the vehicle are discussed. The numerical simulations demonstrate the merits of the inviscid and viscous software. Results are compared with wind tunnel experiments.

Stewart, John E.↗
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 NTRS1992

A computational examination of directional stability for smooth and chined forebodies at high-alpha

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low-speed, angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows over the experimentally investigated F-5A forebody and chine type configuration, previously computed by the authors, were recomputed with better grid topology and resolution. The results were obtained using a modified version of CFL3D (developed at NASA Langley) to solve either the Euler equations or the Reynolds equations employing the Baldwin-Lomax turbulence model with the Degani-Schiff modification to account for massive crossflow separation. Based on the results, it is concluded 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 parametric study of various forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. An unconventional approach for 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, Ramakrishnan↗
NASA NTRS1989

A performance comparison of integration algorithms in simulating flexible structures

Asymptotic formulas for the characteristic root errors as well as transfer function gain and phase errors are presented for a number of traditional and new integration methods. Normalized stability regions in the lambda h plane are compared for the various methods. In particular, it is shown that a modified form of Euler integration with root matching is an especially efficient method for simulating lightly-damped structural modes. The method has been used successfully for structural bending modes in the real-time simulation of missiles. Performance of this algorithm is compared with other special algorithms, including the state-transition method. A predictor-corrector version of the modified Euler algorithm permits it to be extended to the simulation of nonlinear models of the type likely to be obtained when using the discretized structure approach. Performance of the different integration methods is also compared for integration step sizes larger than those for which the asymptotic formulas are valid. It is concluded that many traditional integration methods, such as RD-4, are not competitive in the simulation of lightly damped structures.

Howe, R. M.↗
NASA NTRS1985

Assessment of preliminary prediction techniques for wing leading-edge vortex flows at supersonic speeds

A theoretical investigation of the aerodynamics of sharp leading-edge delta wings at supersonic speeds has been conducted. The primary objective of this was to determine the applicability of existing theoretical methods to predict wing leading-edge separated-flow characteristics at conditions conductive to high-lift supersonic flight. Predicted results from two modified linear-theory methods (LTSTAR and VORCAM) are compared with experimental data. Comparison of the two methods for uncambered wings revealed that the LTSTAR code is in much better agreement with experimentally measured vortex strength, vortex position, and total lifting characteristics than the VORCAM code. Selected analysis was also performed with an Euler code, SWINT. The results of this study indicated that the SWINT code was not well suited to the analysis of wings with separated flow at high lift and low supersonic speeds.

Wood, R. M.↗
NASA NTRS1991

A staggered mesh finite difference scheme for the computation of hypersonic Euler flows

A shock capturing finite difference method for systems of hyperbolic conservation laws is presented which avoids the need to solve Riemann problems while being competitive in performance with other current methods. A staggered spatial mesh is employed, so that complicated nonlinear waves generated at cell interfaces are averaged over cell interiors at the next time level. The full method combines to form a conservative version of the modified method of characteristics. The advantages of the method are discussed, and numerical results are presented for the two-dimensional double ellipse problem.

Sanders, Richard↗
NASA NTRS1999

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗
NASA NTRS1982

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.↗
NASA NTRS1990

Grid generation and flow computation about a Martian entry vehicle

A number of vehicles are currently being proposed for a manned mission to Mars. One of these vehicles has a modified blunt-nosed cone configuration. Experimental results were obtained for this vehicle in 1968. They show lift-over-drag ratios comparable to those needed for Mars entry. Computations are performed to verify the earlier results and to further describe the flight characteristics of this vehicle. An analytical method is used to define the surface of this vehicle. A single-block volume grid is generated around the vehicle using the algebraic Two-Boundary Grid Generation algorithm (TBGG) and transfinite interpolation. Euler solutions are then obtained from a Langley Aerodynamic Upward Relaxation Algorithm (LAURA) at Mach 6.0 and angles of attack of 0, 6, and 12 deg. The lift coefficient determined from the LAURA code agree very well with the experimental results. The drag and pitching moment coefficients, however, are underestimated by the code since viscous effects are not considered. Contour plots of the flowfield show no evidence of separation for angles of attack up to 12 deg.

Stewart, J. E.↗
NASA NTRS1999

Supersonic Aftbody Closure Wind-Tunnel Testing, Data Analysis, and Computational Results

This paper reports on the model, test, and results from the Langley Supersonic Aftbody Closure wind tunnel test. This project is an experimental evaluation of the 1.5% Technology Concept Aircraft (TCA) aftbody closure model (Model 23) in the Langley Unitary Plan Wind Tunnel. The baseline TCA design is the result of a multidisciplinary, multipoint optimization process and was developed using linear design and analysis methods, supplemented with Euler and Navier-Stokes numerical methods. After a thorough design review, it was decided to use an upswept blade attached to the forebody as the mounting system. Structural concerns dictated that a wingtip support system would not be feasible. Only the aftbody part of the model is metric. The metric break was chosen to be at the fuselage station where prior aft-sting supported models had been truncated. Model 23 is thus a modified version of Model 20. The wing strongback, flap parts, and nacelles from Model 20 were used, whereas new aftbodies, a common forebody, and some new tails were fabricated. In summary, significant differences in longitudinal and direction stability and control characteristics between the ABF and ABB aftbody geometries were measured. Correcting the experimental data obtained for the TCA configuration with the flared aftbody to the representative of the baseline TCA closed aftbody will result in a significant reduction in longitudinal stability, a moderate reduction in stabilizer effectiveness and directional stability, and a moderate to significant reduction in rudder effectiveness. These reductions in the stability and control effectiveness levels of the baseline TCA closed aftbody are attributed to the reduction in carry-over area.

Allen, Jerry↗
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 NTRS1992

The natural flow wing-design concept

A wing-design study was conducted on a 65 degree swept leading-edge delta wing in which the wing geometry was modified to take advantage of the naturally occurring flow that forms over a slender wing in a supersonic flow field. Three-dimensional nonlinear analysis methods were used in the study which was divided into three parts: preliminary design, initial design, and final design. In the preliminary design, the wing planform, the design conditions, and the near-conical wing-design 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 flow solver was employed to determine the aerodynamic characteristics of the baseline standard delta wing and to investigate modifications to the airfoil thickness, leading-edge radius, airfoil maximum-thickness position, and wing upper to lower surface asymmetry on the baseline near-conical wing. The final design employed an Euler solver to analyze the best wing configurations found in the initial design and to extend the study of wing asymmetry to develop a more refined wing. Benefits resulting from each modification are discussed, and a final 'natural flow' wing geometry was designed that provides an improvement in aerodynamic performance compared with that of a baseline conventional uncambered wing, linear-theory cambered wing, and near-conical wing.

Wood, Richard M.↗
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 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↗