Engineering topics

Ehlers, F. E.

Publications and source records attributed to Ehlers, F. E..

At least 19 records

NASA NTRS1984

Development and applications of algorithms for calculating the transonic flow about harmonically oscillating wings

A finite difference method to solve the unsteady transonic flow about harmonically oscillating wings was investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The differential equation for the unsteady velocity potential is linear with spatially varying coefficients and with the time variable eliminated by assuming harmonic motion. An alternating direction implicit procedure was investigated, and a pilot program was developed for both two and three dimensional wings. This program provides a relatively efficient relaxation solution without previously encountered solution instability problems. Pressure distributions for two rectangular wings are calculated. Conjugate gradient techniques were developed for the asymmetric, indefinite problem. The conjugate gradient procedure is evaluated for applications to the unsteady transonic problem. Different equations for the alternating direction procedure are derived using a coordinate transformation for swept and tapered wing planforms. Pressure distributions for swept, untaped wings of vanishing thickness are correlated with linear results for sweep angles up to 45 degrees.

Ehlers, F. E.↗

NASA NTRS1984

A method for computing the leading-edge suction in a higher-order panel method

Experimental data show that the phenomenon of a separation induced leading edge vortex is influenced by the wing thickness and the shape of the leading edge. Both thickness and leading edge shape (rounded rather than point) delay the formation of a vortex. Existing computer programs used to predict the effect of a leading edge vortex do not include a procedure for determining whether or not a vortex actually exists. Studies under NASA Contract NAS1-15678 have shown that the vortex development can be predicted by using the relationship between the leading edge suction coefficient and the parabolic nose drag. The linear theory FLEXSTAB was used to calculate the leading edge suction coefficient. This report describes the development of a method for calculating leading edge suction using the capabilities of the higher order panel methods (exact boundary conditions). For a two dimensional case, numerical methods were developed using the double strength and downwash distribution along the chord. A Gaussian quadrature formula that directly incorporates the logarithmic singularity in the downwash distribution, at all panel edges, was found to be the best method.

Ehlers, F. E.↗

NASA NTRS1983

A three degree-of-freedom, typical section flutter analysis using harmonic transonic air forces

Typical section flutter solutions are obtained using generalized air forces calculated with OPTRAN2, a program based on the transonic small perturbation potential equation, which uses potentials from a nonlinear solution to the steady flow to obtain linear unsteady air forces. Flutter results are calculated with a conventional, linear V-g type solution. Results are presented for a NACA 64A010 airfoil oscillating in plunge and pitch, and with a quarter-chord control surface in the Mach number range of 0.80 to 0.90. The flutter boundary shows a significant 'transonic' bucket, together with changes in flutter mode as, with increasing Mach numbers, a shock appears on the airfoil and then moves aft. Single degree-of-freedom instabilities associated with both the pitch and the control surface motions are identified.

Weatherill, W. H.↗

NASA NTRS1982

A harmonic analysis method for unsteady transonic flow and its application to the flutter of airfoils

A finite difference method for solving the unsteady transonic flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The differential equation for the unsteady velocity potential is linear with spatially varying coefficients and with the time variable eliminated by assuming harmonic motion. A study is presented of the shock motion associated with an oscillating airfoil and its representation by the harmonic procedure. The effects of the shock motion and the resulting pressure pulse are shown to be included in the harmonic pressure distributions and the corresponding generalized forces. Analytical and experimental pressure distributions for the NACA 64A010 airfoil are compared for Mach numbers of 0.75, 0.80 and 0.842. A typical section, two-degree-of-freedom flutter analysis of a NACA 64A010 airfoil is performed. The results show a sharp transonic bucket in one case and abrupt changes in instability modes.

Ehlers, F. E.↗

NASA NTRS1980

Further investigation of a finite difference procedure for analyzing the transonic flow about harmonically oscillating airfoils and wings

Analytical and empirical studies of a finite difference method for the solution of the transonic flow about harmonically oscillating wings and airfoils are presented. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances. The steady velocity potential is obtained first from the well-known nonlinear equation for steady transonic flow. The unsteady velocity potential is then obtained from a linear differential equation in complex form with spatially varying coefficients. Since sinusoidal motion is assumed, the unsteady equation is independent of time. An out-of-core direct solution procedure was developed and applied to two-dimensional sections. Results are presented for a section of vanishing thickness in subsonic flow and an NACA 64A006 airfoil in supersonic flow. Good correlation is obtained in the first case at values of Mach number and reduced frequency of direct interest in flutter analyses. Reasonable results are obtained in the second case. Comparisons of two-dimensional finite difference solutions with exact analytic solutions indicate that the accuracy of the difference solution is dependent on the boundary conditions used on the outer boundaries. Homogeneous boundary conditions on the mesh edges that yield complex eigenvalues give the most accurate finite difference solutions. The plane outgoing wave boundary conditions meet these requirements.

Weatherill, W. H.↗

NASA NTRS1980

Analysis of transonic flow about harmonically oscillating airfoils and wings

A finite difference method for analyzing the unsteady transonic flow about harmonically oscillating wings is discussed. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting partial differential equations for small disturbances. Initial solutions were obtained using relaxation procedures, but the solution range proved to be limited in terms of Mach number and reduced frequency. Recent two-dimensional results are presented which have been obtained with direct solution procedures in which the difference equations are solved 'all at once' and these provide reasonable correlation for practical values of Mach number and reduced frequency.

Weatherill, W. H.↗

NASA NTRS1979

A users guide for A344: A program using a finite difference method to analyze transonic flow over oscillating airfoils

The design and usage of a pilot program for calculating the pressure distributions over harmonically oscillating airfoils in transonic flow are described. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equations for small disturbances. The steady velocity potential which must be obtained from some other program, was required for input. The unsteady equation, as solved, is linear with spatially varying coefficients. Since sinusoidal motion was assumed, time was not a variable. The numerical solution was obtained through a finite difference formulation and either a line relaxation or an out of core direct solution method.

Weatherill, W. H.↗

NASA NTRS1979

An investigation of several factors involved in a finite difference procedure for analyzing the transonic flow about harmonically oscillating airfoils and wings

Analytical and empirical studies of a finite difference method for the solution of the transonic flow about harmonically oscillating wings and airfoils are presented. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances. Since sinusoidal motion is assumed, the unsteady equation is independent of time. Three finite difference investigations are discussed including a new operator for mesh points with supersonic flow, the effects on relaxation solution convergence of adding a viscosity term to the original differential equation, and an alternate and relatively simple downstream boundary condition. A method is developed which uses a finite difference procedure over a limited inner region and an approximate analytical procedure for the remaining outer region. Two investigations concerned with three-dimensional flow are presented. The first is the development of an oblique coordinate system for swept and tapered wings. The second derives the additional terms required to make row relaxation solutions converge when mixed flow is present. A finite span flutter analysis procedure is described using the two-dimensional unsteady transonic program with a full three-dimensional steady velocity potential.

Ehlers, F. E.↗

NASA NTRS1979

A higher order panel method for linearized supersonic flow

The basic integral equations of linearized supersonic theory for an advanced supersonic panel method are derived. Methods using only linear varying source strength over each panel or only quadratic doublet strength over each panel gave good agreement with analytic solutions over cones and zero thickness cambered wings. For three dimensional bodies and wings of general shape, combined source and doublet panels with interior boundary conditions to eliminate the internal perturbations lead to a stable method providing good agreement experiment. A panel system with all edges contiguous resulted from dividing the basic four point non-planar panel into eight triangular subpanels, and the doublet strength was made continuous at all edges by a quadratic distribution over each subpanel. Superinclined panels were developed and tested on s simple nacelle and on an airplane model having engine inlets, with excellent results.

Ehlers, F. E.↗

NASA NTRS1978

A Mach line panel method for computing the linearized supersonic flow over planar wings

A method is described for solving the linearized supersonic flow over planar wings using panels bounded by two families of Mach lines. Polynomial distributions of source and doublet strength lead to simple, closed form solutions for the aerodynamic influence coefficients, and a nearly triangular matrix yields rapid solutions for the singularity parameters. The source method was found to be accurate and stable both for analysis and design boundary conditions. Similar results were obtained with the doublet method for analysis boundary conditions on the portion of the wing downstream of the supersonic leading edge, but instabilities in the solution occurred for the region containing a portion of the subsonic leading edge. Research on the method was discontinued before this difficulty was resolved.

Ehlers, F. E.↗

NASA NTRS1978

The practical application of a finite difference method for analyzing transonic flow over oscillating airfoils and wings

Separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances was performed. The steady velocity potential was obtained first from the well known nonlinear equation for steady transonic flow. The unsteady velocity potential was then obtained from a linear differential equation in complex form with spatially varying coefficients. Since sinusoidal motion is assumed, the unsteady equation is independent of time. The results of an investigation into the relaxation-solution-instability problem was discussed. Concepts examined include variations in outer boundary conditions, a coordinate transformation so that the boundary condition at infinity may be applied to the outer boundaries of the finite difference region, and overlapping subregions. The general conclusion was that only a full direct solution in which all unknowns are obtained at the same time will avoid the solution instabilities of relaxation. An analysis of the one-dimensional form of the unsteady transonic equation was studied to evaluate errors between exact and finite difference solutions. Pressure distributions were presented for a low-aspect-ratio clipped delta wing at Mach number of 0.9 and for a moderate-aspect-ratio rectangular wing at a Mach number of 0.875.

Weatherill, W. H.↗

NASA NTRS1978

An improved higher order panel method for linearized supersonic flow

An improved higher order panel method for linearized supersonic flow is described. Each panel, defined by four points on the surface, is divided into eight subpanels in such a way that all subpanel and panel edges are contiguous. By prescribing a quadratic distribution of the doublet on each subpanel, the doublet strength is made strictly continuous on the paneled surface. A linear source distribution is also used. Numerical results are smoother and in better agreement with experiment than the previous method with less strict continuity. A brief discussion of superinclined panels used to eliminate interior interference in nacelles is included.

Ehlers, F. E.↗

NASA NTRS1978

A procedure for analyzing transonic flow over harmonically oscillating airfoils

Finite difference procedures were successfully used to solve the steady transonic flow about airfoils and appear to provide a practical means for calculating the corresponding unsteady flow. The purpose of the paper is to describe a finite difference procedure derived from the equations for the potential flow by assuming small perturbations and harmonic motion. The velocity potential is divided into steady and unsteady parts, and the resulting unsteady equation is linearized on the basis of small amplitudes of oscillation. The steady velocity potential, which must be calculated first, is described by the classical nonlinear transonic differential equation.

Weatherill, W. H.↗

NASA NTRS1977

A user's guide for V174, a program using a finite difference method to analyze transonic flow over oscillating wings

The design and usage of a pilot program using a finite difference method for calculating the pressure distributions over harmonically oscillating wings in transonic flow are discussed. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The steady velocity potential which must be obtained from some other program, is required for input. The unsteady differential equation is linear, complex in form with spatially varying coefficients. Because sinusoidal motion is assumed, time is not a variable. The numerical solution is obtained through a finite difference formulation and a line relaxation solution method.

Butler, T. D.↗

NASA NTRS1977

Application of a finite difference method to the analysis of transonic flow over oscillating airfoils and wings

A finite difference method for solving the unsteady flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearing the resulting unsteady differential equation for small disturbances. Solutions are obtained using relaxation procedures. The means for improving the solution stability characteristics of the relaxation process are explored. A direct procedure is formulated which permits obtaining solutions for combinations of Mach number and reduced frequency for which the relaxation process has proved unstable. The pressure distribution for an aspect ratio 5 rectangular wing oscillating in pitch is presented.

Weatherill, W. H.↗

NASA NTRS1976

A higher order panel method for linearized supersonic flow

A higher order panel method using linearly varying source and quadratically varying doublet distributions for computing linearized supersonic flow over arbitrary wings and bodies is described. With tangential mass flux boundary conditions on circular cones at angles of attack, the source panel method gave pressure distributions which are in excellent agreement with exact nonlinear theory. The numerical source panel method was demonstrated to be very stable by randomly panelling an axially symmetric body. For wings with thickness and camber, both the source and doublet panel methods yield excellent agreement with exact linearized solutions. Computational experiments indicate that the use of exact surface paneling in lieu of linearized boundary conditions for wings provides better agreement with higher order theory and with experiment.

Ehlers, F. E.↗

NASA NTRS1976

On the computation of the transonic perturbation flow field around two- and three-dimensional oscillating wings

A finite difference method for solving the unsteady flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. Solutions are obtained using relaxation procedures. It is determined that there is a limit on reduced frequency, which is a function of Mach number and size of mesh region, above which the relaxation procedures will not converge. It is found that row line relaxation is more efficient than column relaxation and results are presented for a rectangular wing in harmonic pitch.

Weatherill, W. H.↗

NASA NTRS1976

Advanced surface paneling method for subsonic and supersonic flow

Numerical results illustrating the capabilities of an advanced aerodynamic surface paneling method are presented. The method is applicable to both subsonic and supersonic flow, as represented by linearized potential flow theory. The method is based on linearly varying sources and quadratically varying doublets which are distributed over flat or curved panels. These panels are applied to the true surface geometry of arbitrarily shaped three dimensional aerodynamic configurations.

Erickson, L. L.↗