Search NASASearch

Engineering topics

Cebeci, T.

Publications and source records attributed to Cebeci, T..

At least 19 records

Aeroelastic Analysis of Aircraft: Wing and Wing/Fuselage Configurations

A previously developed interface method for coupling aerodynamics and structures is used to evaluate the aeroelastic effects for an advanced transport wing at cruise and under-cruise conditions. The calculated results are compared with wind tunnel test data. The capability of the interface method is also investigated for an MD-90 wing/fuselage configuration. In addition, an aircraft trim analysis is described and applied to wing configurations. The accuracy of turbulence models based on the algebraic eddy viscosity formulation of Cebeci and Smith is studied for airfoil flows at low Mach numbers by using methods based on the solutions of the boundary-layer and Navier-Stokes equations.

Chen, H. H.

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.

Interactive-Boundary-Layer Computations For Oscillating Airfoil

Interactive-boundary-layer method developed for computations of steady flow, extended under assumption of quasi-steady flow, to computations of evolution of two-dimensional flow about oscillating airfoil under light-dynamic-stall conditions. Represents advance toward ability to compute unsteady flows at even greater angles of attack with solutions of equations normally used for description of boundary-layer flows on airfoils prior to stall. Important in practical studies of flow on blades of helicopter rotors, axial compressors, and turbines.

Carr, L. W.

Prediction of stall and post-stall behavior of airfoils at low and high Reynolds numbers

An interactive boundary-layer method, together with the e(super n)-approach to the calculation of transition, has been used to predict the stall and post-stall behavior of airfoils at low and high Reynolds numbers. The turbulence model is based on the Cebeci-Smith algebraic eddy-viscosity formulation with improvements for strong pressure gradient effects and transitional flows at low Reynolds numbers. Comparison of calculated results for incompressible flows indicate good agreement with experiment for a wide range of Reynolds numbers. Preliminary calculations for low Mach number flows with this interactive method with compressibility corrections to the panel method indicate that, at a Mach number of 0.3, the compressibility effect on (C sub Q)max is not negligible.

Cebeci, T.

Analysis of iced wings

A method for computing ice shapes along the leading edge of a wing and a method for predicting its aerodynamic performance degradation due to icing is described. Ice shapes are computed using an extension of the LEWICE code which was developed for airfoils. The aerodynamic properties of the iced wing are determined with an interactive scheme in which the solutions of the inviscid flow equations are obtained from a panel method and the solutions of the viscous flow equations are obtained from an inverse three-dimensional finite-difference boundary-layer method. A new interaction law is used to couple the inviscid and viscous flow solutions. The application of the LEWICE wing code to the calculation of ice shapes on a MS-317 swept wing shows good agreement with measurements. The interactive boundary-layer method is applied to a tapered ice wing in order to study the effect of icing on the aerodynamic properties of the wing at several angles of attack.

Cebeci, T.

Three-dimensional compressible stability-transition calculations using the spatial theory

The e(exp n)-method is employed with the spatial amplification theory to compute the onset of transition on a swept wing tested in transonic cryogenic flow conditions. Two separate eigenvalue formulations are used. One uses the saddle-point method and the other assumes that the amplification vector is normal to the leading edge. Comparisons of calculated results with experimental data show that both formulations give similar results and indicate that the wall temperature has a rather strong effect on the value of the n factor.

Niethammer, R.

Viscous/inviscid interaction procedure for high-amplitude oscillating airfoils

A coupling procedure between an inviscid potential theoretical panel method and an unsteady interactive boundary layer calculation procedure has been developed for calculating unsteady airloads on oscillating airfoils in the high-incidence regime. A coordinate transformation into a time-dependent surface coordinate has been incorporated into the boundary layer code to avoid a quasi-steady treatment of the viscous flow equations in the leading edge region. The upper and lower sides of the airfoil are calculated with the complete set of unsteady equations. The interactive region of the boundary layer calculation can be straightforwardly extended over the complete upper surface of the airfoil while accounting for interaction areas in the region of strong adverse pressure gradients with incipient separation. Comparisons of sample calculations with experimental data show reasonably good agreement even for unsteady drag coefficients.

Geissler, W.

Progress towards the development of an inviscid-viscous interaction method for unsteady flows in turbomachinery cascades

The velocity and heat transfer characteristics of blade passage flows subject to onset velocities that vary in space and time are calculated by means of a general method involving the solution of inviscid and boundary-layer equations coupled with the Hilbert integral. The performance of this boundary layer method for flows with prescribed pressure distribution is evaluated in terms of two model problems that give attention to the stagnation region. One model corresponds to an oscillating airfoil with moving stagnation point; the other is a simulation of a blade which is subject to a uniform onset velocity.

Cebeci, T.

Progress in the prediction of unsteady heat transfer on turbines blades

Progress toward developing a general method for predicting unsteady heat transfer on turbine blades subject to blade-passing frequencies and Reynolds numbers relevant to the Space Shuttle Main Engine (SSME) is discussed. The method employs an invisid/viscous interactive procedure which has been tested extensively for steady subsonic and transonic external airfoil problems. One such example is shown. The agreement with experimental data and with Navier-Stokes calculations yields confidence in the method. The technique is extended to account for wake generated unsteadiness. The flow reversals around the stagnation point caused by the nonuniform onset velocity are accounted for by using the Characteristic Box scheme developed by Cebeci and Stewartson. The coupling between the inviscid and viscous methods is achieved by using a special procedure, which, with a novel inverse finite-difference boundary-layer method, allows the calculations to be performed for a wide range of flow conditions, including separation. Preliminary results are presented for the stagnation region of turbine blades for both laminar and turbulent flows. A laminar model problem corresponding to a flow on a circular cylinder which experiences the periodic passing of wakes from turbine blades is presented to demonstrate the ability of the method to calculate flow reversals around the stagnation region.

Cebeci, T.

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.

A comparison of interactive boundary-layer and thin-layer Navier-Stokes procedures

Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Mehta, U.

Relative advantages of thin-layer Navier-Stokes and interactive boundary-layer procedures

Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Mehta, U.

Computational aspects of unsteady flows

The calculation of unsteady flows and the development of numerical methods for solving unsteady boundary layer equations and their application to the flows around important configurations such as oscillating airfoils are presented. A brief review of recent work is provided with emphasis on the need for numerical methods which can overcome possible problems associated with flow reversal and separation. The zig-zag and characteristic box schemes are described in this context, and when embodied in a method which permits interaction between solutions of inviscid and viscous equations, the characteristic box scheme is shown to avoid the singularity associated with boundary layer equations and prescribed pressure gradient. Calculations were performed for a cylinder started impulsively from rest and oscillating airfoils. The results are presented and discussed. It is conlcuded that turbulence models based on an algebraic specification of eddy viscosity can be adequate, that location of translation is important to the calculation of the location of flow separation and, therefore, to the overall lift of an oscillating airfoil.

Cebeci, T.

Boundary layers on oscillating airfoils

A two-point finite-difference unsteady laminar and turbulent boundary-layer computational method has been used to investigate the properties of the flow around an airfoil (NACA 0012) oscillating through angles of attack up to 18 degrees, and for reduced frequencies of 0.01 and 0.20. The unsteady potential flow was determined using the method of Geissler. The influence of transition location on stall behavior has been investigated, using both experimentally determined transition information, and transition located at the pressure peak; the results show the need for viscous-inviscid interaction in future computations of such flows.

Carr, L. W.

Calculation of boundary layers of oscillating airfoils

A two-point finite difference unsteady laminar and turbulent boundary layer computational method was used to investigate the properties of the flow around an airfoil (NACA 0012) oscillating through angles of attack up to 18 degrees, for reduced frequencies of 0.01 and 0.20. The unsteady potential flow was determined using the unsteady potential flow method of Geissler. The influence of transition location on stal behavior was investigated, using both experimentally determined transition information, and transition located at the pressure peak; the results show the need for viscous-inviscid interaction in future computation of such flows.

Cebeci, T.

A time-dependent approach for calculating steady inverse boundary-layer flows with separation

An unsteady inverse boundary-layer method is developed which can be used to calculate steady flows with separation. Two versions of Keller's box method with the Mechul function formulation developed by Cebeci (1976) are employed, depending on the complexity of the flow. The regular box is employed in regions of positive streamwise velocity component u, whereas the zigzag box is employed in regions where u becomes negative (t greater than 0). The regular box with the FLARE approximation is employed when t = 0 and u becomes negative in some region across the layer. Results of calculations show that the use of a time-dependent inverse boundary-layer method in which time is used as an iteration parameter provides a good approach in improving the accuracy of the solutions obtained from the FLARE approximation.

Cebeci, T.

Calculation of boundary layers near the stagnation point of an oscillating airfoil

The results of an investigation of boundary layers close to the stagnation point of an oscillating airfoil are reported. Two procedures for generating initial conditions, the characteristics box scheme and a quasi-static approach, were investigated, and the quasi-static approach was shown to be appropriate provided the initial region was far from any flow separation. With initial conditions generated in this way, the unsteady boundary layer equations were solved for the flow in the leading edge region of a NACA 0012 airfoil oscillating from 0 to 5 deg. Results were obtained for both laminar and turbulent flow, and, in the latter case, the effect of transition was assessed by specifying its occurrence at different locations. The results demonstrate the validity of the numerical scheme and suggest that the procedures should be applied to calculation of the entire flow around oscillating airfoils.

Cebeci, T.

Prediction of boundary-layer characteristics of an oscillating airfoil

The evolution of unsteady boundary layers on oscillating airfoils is investigated by solving the governing equations by the Characteristic Box scheme. The difficulties associated with computing the first profile on a given time line, and the velocity profiles with partial flow reversal are solved. A sample calculation is performed for an external velocity distribution typical of those found near the leading edge of thin airfoils. The viability of the calculation procedure is demonstrated.

Cebeci, T.