Search NASASearch

Engineering topics

Livne, Eli

Publications and source records attributed to Livne, Eli.

Progress on Gust Load Alleviation Wind Tunnel Experiment and Aeroservoelastic Model Validation for a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper discusses a wind tunnel experiment of active gust load alleviation of a flexible wing which took place at University of Washington (UW) in 2019. The experiment performed under a NASA SBIR contract with Scientific Systems Company, Inc (SSCI). The objective of the experiment is to demonstrate active controls of the Variable Camber Continuous Trailing Edge Flap (VCCTEF) system for gust load alleviation and real-time drag optimization. The wind tunnel model is a 8.2% sub-scale Common Research Model (CRM) wing. The wing structure is designed to provide a substantial degree of flexibility to represent that of a modern high-aspect ratio wing. Eight active control surfaces are employed in the VCCTEF. A new gust generator system was designed and installed by UW under a sub-contract with SSCI. The first test entry started in July 2019 and ended in September 2019. During this test entry, many significant issues were found with the hardware and software. The significant issues with the servos prevented the test objective from being completed. A follow-up second test entry in 2020 is being planned. The wing system is being repaired by SSCI. This paper reports on the progress of this experimental effort and the aeroservoelastic (ASE) model validation which was conducted during the test entry.

Nguyen, Nhan

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization

Experimental Investigation of a Flexible Wing with a Variable Camber Continuous Trailing Edge Flap Design

This paper presents experimental results of a flexible wing wind tunnel model with a variable camber continuous trailing edge flap (VCCTEF) design for drag minimization, tested at the University of Washington Aeronautical Labo- ratory (UWAL). The wind tunnel test was designed to explore the relative merit of the VCCTEF concept for improved cruise efficiency through the use of low-cost aeroelastic model test techniques. The flexible wing model is a 10-scaled model of a typical transport wing and is constructed of woven fabric composites and foam core. The wing structural stiffness in bending is tailored to be half of the stiffness of a Boeing 757-era transport wing, while the torsional stiffness is about the same. This stiffness reduction results in a wing tip deflection of about 10 of the wing semi-span. The VCCTEF is a multi-segment flap design having three chordwise camber segments and five spanwise flap sections for a total of 15 individual flap elements. The three chordwise camber segments can be positioned appropriately to create a desired trailing edge camber. Elastomeric material is used to cover the gaps in between the spanwise flap sections, thereby creating a continuous trailing edge. Wind tunnel data indicate a high degree of data correlation and repeata- bility. The VCCTEF can achieve a drag reduction of up to 6.31 and an improvement in the lift-to-drag ratio (LD) of up to 4.85. The paper also presents two methods for estimating the lift coefficient of a rigid wing using a dynamic pressure correction and an aeroelastic deflection correction. Both methods provide good estimates of the rigid-wing lift coefficient.

Aerodynamics

Design Oriented Structural Modeling for Airplane Conceptual Design Optimization

The main goal for research conducted with the support of this grant was to develop design oriented structural optimization methods for the conceptual design of airplanes. Traditionally in conceptual design airframe weight is estimated based on statistical equations developed over years of fitting airplane weight data in data bases of similar existing air- planes. Utilization of such regression equations for the design of new airplanes can be justified only if the new air-planes use structural technology similar to the technology on the airplanes in those weight data bases. If any new structural technology is to be pursued or any new unconventional configurations designed the statistical weight equations cannot be used. In such cases any structural weight estimation must be based on rigorous "physics based" structural analysis and optimization of the airframes under consideration. Work under this grant progressed to explore airframe design-oriented structural optimization techniques along two lines of research: methods based on "fast" design oriented finite element technology and methods based on equivalent plate / equivalent shell models of airframes, in which the vehicle is modelled as an assembly of plate and shell components, each simulating a lifting surface or nacelle / fuselage pieces. Since response to changes in geometry are essential in conceptual design of airplanes, as well as the capability to optimize the shape itself, research supported by this grant sought to develop efficient techniques for parametrization of airplane shape and sensitivity analysis with respect to shape design variables. Towards the end of the grant period a prototype automated structural analysis code designed to work with the NASA Aircraft Synthesis conceptual design code ACS= was delivered to NASA Ames.

Livne, Eli

Integrated Aeroservoelastic Optimization: Status and Direction

The interactions of lightweight flexible airframe structures, steady and unsteady aerodynamics, and wide-bandwidth active controls on modern airplanes lead to considerable multidisciplinary design challenges. More than 25 years of mathematical and numerical methods' development, numerous basic research studies, simulations and wind-tunnel tests of simple models, wind-tunnel tests of complex models of real airplanes, as well as flight tests of actively controlled airplanes, have all contributed to the accumulation of a substantial body of knowledge in the area of aeroservoelasticity. A number of analysis codes, with the capabilities to model real airplane systems under the assumptions of linearity, have been developed. Many tests have been conducted, and results were correlated with analytical predictions. A selective sample of references covering aeroservoelastic testing programs from the 1960s to the early 1980s, as well as more recent wind-tunnel test programs of real or realistic configurations, are included in the References section of this paper. An examination of references 20-29 will reveal that in the course of development (or later modification), of almost every modern airplane with a high authority active control system, there arose a need to face aeroservoelastic problems and aeroservoelastic design challenges.

Livne, Eli

Panel Flutter Constraints: Analytic Sensitivities and Approximations Including Planform Shape Design Variables

Analytical sensitivities of panel flutter constraints with respect to panel shape as well as thickness and material properties are derived and numerically tested. Cases of fixed in-plane loads and cases in which in-plane loads are variable (depending on panel and overall wing shape as well as material and sizing design variables) are considered. Accuracy of approximations and range of move limits required are studied in preparation for integration with nonlinear programming/approximation concept aeroelastic design synthesis methodology.

Livne, Eli

Analytic Sensitivity and Approximation of Skin Buckling Constraints in Wing-Shape Synthesis

Explicit expressions for terms of the stiffness and geometric stiffness matrices are derived for the buckling analysis of trapezoidal fiber composite wing skin panels. The formulation is based on Ritz analysis using simple polynomials, and leads to explicit expressions for the analytic sensitivities of the stiffness and geometric stiffness matrices with respect to layer thickness, fiber directions, and panel shape. Integration with wing box analysis using either the equivalent plate approach or the finite element method, makes it possible to obtain sensitivities of panel buckling constraints with respect to wing planform shape or locations of internal ribs and spars. The analytic sensitivities are used to construct approximations of panel buckling constraints for integrated wing/panel design synthesis.

Livne, Eli

Analytic Sensitivities for Shape Optimization in Equivalent Plate Structural Wing Models

Equivalent plate modeling techniques based on Ritz analysis with simple polynomials prove to be efficient tools for structural modeling of wings in the preliminary design stage. Accuracy problems are encountered, however, when these models are used to obtain finite difference behavior sensitivities with respect to planform shape. The accuracy problems are associated with the poor numerical conditioning of static and eigenvalue equations. As higher-order polynomials are being used to Improve the analysis itself, the more sensitive is the finite difference derivative to the step size used. This article describes a formulation of wing equivalent plate modeling in which it is simple to obtain analytic, explicit expressions for stiffness and mass matrix elements without the need to perform numerical integration. This formulation leads naturally to analytic expressions for the derivatives of displacements, stresses, and natural frequencies with respect to shape design variables. This article examines the accuracy of finite difference derivatives compared with the analytic derivatives, and shows that In some cases it is impossible to obtain any information of value by finite differences. Analytic sensitivities, in this case, are still sufficiently accurate for design optimization.

Livne, Eli

Equivalent plate structural modeling for wing shape optimization including transverse shear

A new technique for structural modeling of airplanes wings is presented taking transverse shear effects into account. The kinematic assumptions of first-order shear deformation plate theory in combination with numerical analysis, where simple polynomials are used to define geometry, construction, and displacement approximations, lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps, and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of wing modeling techniques based on classical plate theory are discussed, and the improved accuracy of the new equivalent plate technique is demonstrated through comparison with finite element analysis and test results. Expressions for analytical derivatives of stiffness, mass, and load terms with respect to wing shape are given. Based on these, it is possible to obtain analytic sensitivities of displacements, stresses, and natural frequencies with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli

Equivalent Plate Structural Modeling for Wing Shape Optimization Including Transverse Shear

A new technique for structural modeling of airplane wings is presented taking transverse shear effects into account. The kinematic assumptions of first-order shear deformation plate theory In combination with numerical analysis, where simple polynomials are used to define geometry, construction, and displacement approximations, lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps, and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of wing modeling techniques based on classical plate theory are discussed, and the Improved accuracy of the new equivalent plate technique is demonstrated through comparison with finite element analysis and test results. Expressions for analytical derivatives of stiffness, mass, and load terms with respect to wing shape are given. Based on these, it is possible to obtain analytic sensitivities of displacements, stresses, and natural frequencies with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli

The delayed-detonation model of a type Ia supernovae. 1: The deflagration phase

The nature of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is investigated by using two-dimensional numerical hydrodynamics simulations. A new algorithm is used to treat the deflagration front. Assuming that it propagates locally at the laminar flame speed, the deflagration is insufficient to unbind the star. Expansion shuts of the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. The burning front does become wrinkled, but the wavelength of the instability is much larger than the computational grid size and is resolved; this is consistent with previous analysis. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. This paper deals with the deflagration phase, from the onset of burning, through expansion and quenching of the flame, to the first contraction.

Arnett, David

The delayed-detonation model of Type Ia supernovae. 2: The detonation phase

The investigation, by use of two-dimensional numerical hydrodynamics simulations, of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is continued. Previously we found that the deflagration is insufficient to unbind the star. Expansion shuts off the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. The burning was allowed to develop into a detonation in these nonspherical models. The detonation grows toward spherical symmetry at late times. At these densities (rho approx. 10(exp 7) to 10(exp 8) g cm(exp -3)), either Ni-56 or nuclei of the Si-Ca group are the dominant products of the burning. The bulk yields are sensitive to the density of the star when the transition to detonation occurs. The relevance of the abundances, velocities, mixing, and total energy release to the theory and interpretation of Type Ia supernovae is discussed.

Arnett, David

An implicit method for two-dimensional hydrodynamics

An implicit method for compressible multidimensional flows is presented. The method, which is strongly oriented toward astrophysical applications, enables one to simulate very subsonic flows by removing the Courant condition upon time steps. It consists of an implicit purely Lagrangian step, followed by an explicit and second-order accurate (at least in one dimension) remapping step, which is optional. When the remapping step is performed the time step is limited by the 'particle crossing time' and otherwise it is limited only by accuracy considerations. The suggested method, which results from a compromise between accuracy and efficiency, is very efficient relative to other methods. It enables the computation of many multidimensional problems in stellar evolution, such as those governed by very subsonic flows, which were not calculable with existing explicit methods.

Livne, Eli

Numerical simulations of the convective flame in white dwarfs

A first step toward better understanding of the mechanism driving convective flames in exploding white dwarfs is presented. The propagation of the convective flame is examined using a two-dimensional implicit hydrodynamical code. The large scales of the instability are captured by the grid while the scales that are smaller than the grid resolution are approximated by a mixing-length approximation. It is found that largescale perturbations (of order of the pressure scale height) do grow significantly during the expansion, leading to a very nonspherical burning front. The combustion rate is strongly enhanced (compared to the unperturbed case) during the first second, but later the expansion of the star suppresses the flame speed, leading to only partial incineration of the nuclear fuel. Our results imply that large-scale perturbations by themselves are not enough to explain the mechanism by which convective flames are driven, and a study of the whole spectrum of relevant perturbations is needed. The implications of these preliminary results on future simulations, in the context of current models for Type Ia supernovae, are discussed.

Livne, Eli

Recent developments in equivalent plate modeling for wing shape optimization

A new technique for structural modeling of airplane wings is presented taking transverse shear effects into account. The kinematic assumptions of first order shear deformation plate theory in combination with numerical analysis based on simple polynomials which define geometry, construction and displacement approximations lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of current equivalent plate wing modeling techniques based on classical plate theory are discussed, and the improved accuracy of the new equivalent plate technique is demonstrated through comparison to finite element analysis and test results. Analytical derivatives of stiffness, mass and load terms with respect to wing shape lead to analytic sensitivities of displacements, stresses and natural modes with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli

Transonic adaptive flutter suppression using approximate unsteady time domain aerodynamics

A digital adaptive controller is applied to the active flutter suppression problem of a wing under time varying flight conditions in subsonic and transonic flow. Linear quadratic controller gain at each time step is obtained using an iterative Riccati solver. The digital adaptive optimal controller is robust with respect to the unknown external loads. Flutter and divergence instabilities are simultaneously suppressed using a trailing-edge control surface and displacement sensing. A new transonic unsteady aerodynamic approximation methodology is developed which enables one to carry out the rapid calculation required for transonic aeroservoelastic applications. This approximation is based on a combination of unsteady subsonic aerodynamics combined with a transonic correction procedure. Aeroservoelastic transient time response is obtained using Roger's approximation, state transition matrices and an iterative time marching algorithm. The aeroservoelastic system in the time domain is modeled using a deterministic ARMA model together with a parameter estimator. Transonic flutter boundaries of a wing structure are computed, in the time domain, using an estimated aeroelastic system matrix and are in good agreement with experimental data for the low transonic Mach number range.

Pak, Chan-Gi

Accurate calculation of control-augmented structural eigenvalue sensitivities using reduced-order models

A method is presented for generating mode shapes for model order reduction in a way that leads to accurate calculation of eigenvalue derivatives and eigenvalues for a class of control augmented structures. The method is based on treating degrees of freedom where control forces act or masses are changed in a manner analogous to that used for boundary degrees of freedom in component mode synthesis. It is especially suited for structures controlled by a small number of actuators and/or tuned by a small number of concentrated masses whose positions are predetermined. A control augmented multispan beam with closely spaced natural frequencies is used for numerical experimentation. A comparison with reduced-order eigenvalue sensitivity calculations based on the normal modes of the structure shows that the method presented produces significant improvements in accuracy.

Livne, Eli