Nonlinear Reduced-Order Models for Aerodynamic Lift of Oscillating Airfoils

Author(s):  
Muhammad Saif Ullah Khalid ◽  
Imran Akhtar

For the present study, setting Strouhal Number as the control parameter, we perform numerical simulations for the flow over oscillating NACA-0012 airfoil at Reynolds Number 103. Temporal profiles of the unsteady lift, and their respective spectra clearly indicate the solution to be a period-1 attractor for low Strouhal numbers. This study reveals that aerodynamic forces produced by the oscillating airfoils are independent of the initial kinematic conditions that proves existence of the limit cycle. Frequencies present in the oscillating lift force are composed of the fundamental harmonics, and its odd harmonics. Using these numerical simulations, we observe the shedding frequencies nearly equal to the excitation frequencies in all the cases. Hence, considering it as a primary resonance case, we model the unsteady lift force through a modified van der Pol oscillator. Using the method of multiple scales and spectral analysis of the steady-state CFD solutions, we estimate the frequencies and the damping terms in the reduced-order model. We prove the applicability of this model to all the planar motions of airfoil; heaving, pitching and flapping. With the increasing Strouhal number, the nonlinear damping terms for all types of motion approach similar magnitudes. Another important aspect in one of the currently-proposed model is capturing the time-averaged value of the aero-dynamic lift force.

Author(s):  
Muhammad Saif Ullah Khalid ◽  
Imran Akhtar

For the present study, setting Strouhal number as the control parameter, we perform numerical simulations for the flow over oscillating NACA-0012 airfoil at a Reynolds number of 1000. This study reveals that aerodynamic forces produced by oscillating airfoils are independent of the initial kinematic conditions suggesting the existence of limit cycle. Frequencies present in the oscillating lift force are composed of the fundamental harmonics and its odd harmonics. Using these numerical simulations, we analyze the shedding frequencies close to the excitation frequencies. Hence, considering it as a primary resonance case, we model the unsteady lift force with a modified van der Pol oscillator. Using the method of multiple scales and spectral analysis of the steady-state computational fluid dynamics (CFD) solutions, we estimate the frequencies and the damping terms in the reduced-order model (ROM). We show the applicability of this model to all planar motions of the airfoil; heaving, pitching, and flapping. With increasing the Strouhal number, the nonlinear damping terms for all types of motion approach similar magnitudes. Another important aspect in one of the proposed model is its ability to capture the time-averaged value of the aerodynamic lift force. We also notice that increase in the magnitude of the lift force is due to the effect of destabilizing linear damping parameter.


Author(s):  
Hassen M. Ouakad ◽  
Mohammad I. Younis ◽  
Fadi M. Alsaleem ◽  
Ronald Miles ◽  
Weili Cui

In this paper, we investigate theoretically and experimentally the static and dynamic behaviors of electrostatically actuated clamped-clamped micromachined arches when excited by a DC load superimposed to an AC harmonic load. A Galerkin based reduced-order model is used to discretize the distributed-parameter model of the considered shallow arch. The natural frequencies of the arch are calculated for various values of DC voltages and initial rises of the arch. The forced vibration response of the arch to a combined DC and AC harmonic load is determined when excited near its fundamental natural frequency. For small DC and AC loads, a perturbation technique (the method of multiple scales) is also used. For large DC and AC, the reduced-order model equations are integrated numerically with time to get the arch dynamic response. The results show various nonlinear scenarios of transitions to snap-through and dynamic pull-in. The effect of rise is shown to have significant effect on the dynamical behavior of the MEMS arch. Experimental work is conducted to test polysilicon curved microbeam when excited by DC and AC loads. Experimental results on primary resonance and dynamic pull-in are shown and compared with the theoretical results.


Author(s):  
Dumitru I. Caruntu ◽  
Ezequiel Juarez

Abstract This paper deals with the frequency-amplitude response of primary resonance of electrostatically actuated Double-Walled Carbon Nanotubes (DWCNT) and Single-Walled Carbon Nanotubes (SWCNT) cantilever resonators. Their responses are compared. Both the DWCNT and SWCNT are modeled as Euler-Bernoulli cantilever beams. Electrostatic and damping forces are applied on both types of resonators. An AC voltage provides a soft electrostatic actuation. For the DWCNT, intertube van der Waals forces are present between the carbon nanotubes, coupling the deflections of the tubes and acting as a nonlinear spring between the two carbon nanotubes. Electrostatic (for SWCNT and DWCNT) and intertube van der Waals (for DWCNT) forces are nonlinear. Both resonators undergo nonlinear parametric excitation. The Method of Multiple Scales (MMS) is used to investigate the systems under soft excitations and weak nonlinearities. A 2-Term Reduced-Order-Model (ROM) is numerically solved for stability analysis using AUTO-07P, a continuation and bifurcation software. The coaxial vibrations of DWCNT are considered in this work, in order to draw comparisons between DWCNT and SWCNT. Effects of damping and voltage of the frequency-amplitude response are reported.


2017 ◽  
Vol 139 (4) ◽  
Author(s):  
Karin Mora ◽  
Oded Gottlieb

The dynamic motion of a parametrically excited microbeam-string affected by nonlinear damping is considered asymptotically and numerically. It is assumed that the geometrically nonlinear beam-string, subject to only modulated alternating current voltage, is closer to one of the electrodes, thus resulting in an asymmetric dual gap configuration. A consequence of these novel assumptions is a combined parametric and hard excitation in the derived continuum-based model that incorporates both linear viscous and nonlinear viscoelastic damping terms. To understand how these assumptions influence the beam's performance, the conditions that lead to both principal parametric resonance and a three-to-one internal resonance are investigated. Such conditions are derived analytically from a reduced-order nonlinear model for the first three modes of the microbeam-string using the asymptotic multiple-scales method which requires reconstitution of the slow-scale evolution equations to deduce an approximate spatio-temporal solution. The response is investigated analytically and numerically and reveals a bifurcation structure that includes coexisting in-phase and out-of-phase solutions, Hopf bifurcations, and conditions for the loss of orbital stability culminating with nonstationary quasi-periodic solutions and chaotic strange attractors.


Author(s):  
Dumitru I. Caruntu ◽  
Le Luo

This paper investigates electrostatically actuated Carbon Nano-Tubes (CNT) cantilevers biosensors using the Reduced Order Model (ROM) method. Forces acting on the CNT are electrostatic, damping, and van der Waals. The electrostatic actuation is given by soft AC voltage. Van der Waals forces are significant for gaps between the CNT and a ground plate lower than 100 nm. Both forces electrostatic and van der Waals are nonlinear. CNT undergoes nonlinear parametric dynamics. ROM is used to investigate the system under soft excitations and/or weak nonlinearities. The frequency-amplitude response is found in the case of primary resonance and compared to the Method of Multiple Scales (MMS). The CNT biosensor is to be used for mass detection applications.


Author(s):  
Ezequiel Juarez ◽  
Dumitru I. Caruntu

In this paper, the Reduced Order Method (ROM) and the Method of Multiple Scales (MMS) are used to investigate the influences of dimensionless damping and voltage parameters on the amplitude-frequency response of an electrostatically actuated double-walled carbon nanotube (DWCNT). The forces responsible for the nonlinearities in the vibrational behavior are intertube van der Waals and electrostatic forces. Soft AC excitation and small viscous damping forces are assumed. Herein, the coaxial case is investigated. In this mode of vibration, the outer and inner carbon nanotubes move synchronously (in-phase) with the same maximum tip deflection. The DWCNT structure is modelled as a cantilever beam with Euler-Bernoulli beam assumptions since the DWCNT is characterized with high length-diameter ratio. The results shown assume steady-state solutions in the first-order MMS solution. The analytical approximate solutions provided by MMS are validated numerically by two-term (2T) Time Reponses and AUTO-07P. The two methods in this paper are found to be in excellent agreement at lower amplitudes. Additionally, the two methods are assessed for their advantages and limitations. The importance of the results in this paper are the effect of damping and voltage on the stability of the DWCNT vibration.


2002 ◽  
Vol 124 (4) ◽  
pp. 988-993 ◽  
Author(s):  
V. Esfahanian ◽  
M. Behbahani-nejad

An approach to developing a general technique for constructing reduced-order models of unsteady flows about three-dimensional complex geometries is presented. The boundary element method along with the potential flow is used to analyze unsteady flows over two-dimensional airfoils, three-dimensional wings, and wing-body configurations. Eigenanalysis of unsteady flows over a NACA 0012 airfoil, a three-dimensional wing with the NACA 0012 section and a wing-body configuration is performed in time domain based on the unsteady boundary element formulation. Reduced-order models are constructed with and without the static correction. The numerical results demonstrate the accuracy and efficiency of the present method in reduced-order modeling of unsteady flows over complex configurations.


2017 ◽  
Vol 9 (6) ◽  
pp. 1485-1505
Author(s):  
Lingchang Meng ◽  
Fengming Li

AbstractThe nonlinear transverse vibrations of ordered and disordered two-dimensional (2D) two-span composite laminated plates are studied. Based on the von Karman's large deformation theory, the equations of motion of each-span composite laminated plate are formulated using Hamilton's principle, and the partial differential equations are discretized into nonlinear ordinary ones through the Galerkin's method. The primary resonance and 1/3 sub-harmonic resonance are investigated by using the method of multiple scales. The amplitude-frequency relations of the steady-state responses and their stability analyses in each kind of resonance are carried out. The effects of the disorder ratio and ply angle on the two different resonances are analyzed. From the numerical results, it can be concluded that disorder in the length of the two-span 2D composite laminated plate will cause the nonlinear vibration localization phenomenon, and with the increase of the disorder ratio, the vibration localization phenomenon will become more obvious. Moreover, the amplitude-frequency curves for both primary resonance and 1/3 sub-harmonic resonance obtained by the present analytical method are compared with those by the numerical integration, and satisfactory precision can be obtained for engineering applications and the results certify the correctness of the present approximately analytical solutions.


Author(s):  
Bowen Yan ◽  
Yangjin Yuan ◽  
Dalong Li ◽  
Ke Li ◽  
Qingshan Yang ◽  
...  

The semi-periodic vortex-shedding phenomenon caused by flow separation at the windward corners of a rectangular cylinder would result in significant vortex-induced vibrations (VIVs). Based on the aeroelastic experiment of a rectangular cylinder with side ratio of 1.5:1, 2-dimensional (2D) and 2.5-dimensional (2.5D) numerical simulations of the VIV of a rectangular cylinder were comprehensively validated. The mechanism of VIV of the rectangular cylinder was in detail discussed in terms of vortex-induced forces, aeroelastic response, work analysis, aerodynamic damping ratio and flow visualization. The outcomes showed that the numerical results of aeroelastic displacement in the cross-wind direction and the vortex-shedding procedure around the rectangular cylinder were in general consistence with the experimental results by 2.5D numerical simulation. In both simulations, the phase difference between the lift and displacement response increased with the reduced wind speed and the vortex-induced resonance (VIR) disappeared at the phase difference of approximately 180∘. The work done by lift force shows a close relationship with vibration amplitudes at different reduced wind speeds. In 2.5D simulations, the lift force of the rectangular cylinder under different wind speeds would be affected by the presence of small-scale vortices in the turbulence flow field. Similarly, the phase difference between lift force and displacement response was not a constant with the same upstream wind speed. Aerodynamic damping identified from the VIV was mainly dependent on the reduced wind speed and negative damping ratios were revealed at the lock-in regime, which also greatly influenced the probability density function (PDF) of wind-induced displacement.


2021 ◽  
Author(s):  
Richard Hollenbach ◽  
Robert Kielb ◽  
Kenneth Hall

Abstract This paper expands upon a multi-degree-of-freedom, Van der Pol oscillator used to model buffet and Nonsynchronous Vibrations (NSV) in turbines. Two degrees-of-freedom are used, a fluid tracking variable incorporating a Van der Pol oscillator and a classic spring, mass, damper mounted cylinder variable; thus, this model is one of fluid-structure interaction. This model has been previously shown to exhibit the two main aspects of NSV. The first is the lock-in or entrainment phenomenon of the fluid shedding frequency jumping onto the natural frequency of the oscillator, while the second is a stable limit cycle oscillation (LCO) once the transient solution disappears. Improvements are made to the previous model to better understand this aeroelastic phenomenon. First, an error minimizing technique through a system identification method is used to tune the coefficients in the Reduced Order Model (ROM) to improve the accuracy in comparison to experimental data. Secondly, a cubic stiffness term is added to the fluid equation; this term is often seen in the Duffing Oscillator equation, which allows this ROM to capture the experimental behavior more accurately, seen in previous literature. The finalized model captures the experimental cylinder data found in literature much better than the previous model. These improvements also open the door for future models, such as that of a pitching airfoil or a turbomachinery blade, to create a preliminary design tool for studying NSV in turbomachinery.


Sign in / Sign up

Export Citation Format

Share Document