Aerodynamic stability derivative calculations using the compressible source and doublet panel method
Dimitriadis, Grigorios
2024 • In Journal of Aircraft
This work presents the development of closed form expressions for the aerodynamic stability derivatives of wing and aircraft flying at subsonic compressible conditions using the unsteady subsonic 3D Source and Doublet Panel Method (SDPM). The stability derivatives are obtained directly from the frequency domain version of the nonlinear unsteady subsonic Bernoulli equation and can be calculated for negligible additional computational cost around the true geometry of the wing or aircraft, for each frequency of interest. The methodology is validated by applying it to an experimental test case of a swept wing oscillating in yaw and sideslip in the wind tunnel, demonstrating that the values of the aerodynamic stability derivatives obtained from the SDPM are in good agreement with the experimental data at low and moderate angles of attack. Then, the methodology is applied to a Blended Wing Body Unmanned Aerial Vehicle configuration and the SDPM stability derivative predictions are compared to estimates obtained from both steady CFD calculations and the USAF DATCOM methodology. It is shown that the SDPM approach is accurate and can yield more information about the stability of the aircraft than the DATCOM, given that the latter was developed for conventional aircraft configurations.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Continuation of Higher Order Harmonic Balance Solutions for Nonlinear Aeroelastic Systems
Dimitriadis, Grigorios
2008 • In Journal of Aircraft, 45 (2), p. 523-537
The Harmonic Balance method is a very useful tool for characterizing and predicting the response of nonlinear dynamic systems undergoing periodic oscillations, either self-excited or due to harmonic excitation. The method and several of its variants have been applied to nonlinear aeroelastic systems over the last two decades.
This paper presents a detailed description of several Harmonic Balance methods and a continuation framework allowing the methods to follow the response of dynamic systems from the bifurcation point to any desired parameter value, while successfully negotiating further fold bifurcations. The continuation framework is described for systems undergoing sub-critical and super-critical Hopf bifurcations as well as a particular type of explosive bifurcation.
The methods investigated in this work are applied to a nonlinear aeroelastic model of a Generic Transport Aircraft featuring polynomial or freeplay stiffness nonlinearity in the control surface. It is shown that high order Harmonic Balance solutions will capture accurately the complete bifurcation behavior of this system for both types of nonlinearity. Low order solutions can become inaccurate in the presence of numerous folds in the Limit Cycle Oscillation branch but can still yield practical engineering information at a fraction of the cost of higher order solutions. Time domain Harmonic Balance schemes are shown to be more computationally expensive than the standard Harmonic Balance approach.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Bifurcation analysis of aircraft with structural nonlinearity and freeplay using numerical continuation
Dimitriadis, Grigorios
2008 • In Journal of Aircraft, 45 (3), p. 893-905
In recent years the aeroelastic research community has carried out substantial work on the characterization and prediction of nonlinear aeroelastic phenomena. Of particular interest is the calculation of Limit Cycle Oscillations (LCO), which cannot be accomplished using traditional linear methods. In this paper, the prediction of the bifurcation and post-bifurcation behavior of nonlinear subsonic aircraft is carried out using Numerical Continuation. The analysis does not make use of continuation packages such as AUTO or MatCont. Two different continuation techniques are detailed, specifically adapted for realistic aeroelastic models. The approaches are demonstrated on model of a simple pitch plunge airfoil with cubic stiffness and an aeroelastic model of a transport aircraft with two different types of nonlinearity in the control surface. It is shown that one of the techniques yields highly accurate predictions for LCO amplitudes and periods while the second method trades off some accuracy for computational efficiency.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Experimental Validation Of The Constant Level Method For Identification Of Non-Linear Multi-Degree-Of-Freedom Systems
Dimitriadis, Grigorios
2002 • In Journal of Sound and Vibration, 258 (5), p. 829-845
System identification for non-linear dynamical systems could find use in many applications such as condition monitoring, finite element model validation and determination of stability. The effectiveness of existing non-linear system identification techniques is limited by various factors such as the complexity of the system under investigation and the type of non-linearities present. In this work, the constant level identification approach, which can identify multi-degree-of-freedom systems featuring any type of non-linear function, including discontinuous functions, is validated experimentally. The method is shown to identify accurately an experimental dynamical system featuring two types of stiffness non-linearity. The full equations of motion are also extracted accurately, even in the presence of a discontinuous non-linearity.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Shooting-Based Complete Bifurcation Prediction for Aeroelastic Systems with Freeplay
Dimitriadis, Grigorios
2011 • In Journal of Aircraft, 48 (6), p. 1864-1877
In recent years there have been several applications of numerical continuation approaches to aeroelastic systems with freeplay. While some of these have been successful, the general application of the method to such systems remains problematic. Numerical continuation can fail in the presence of complex bifurcations, numerous nearby periodic solution branches and other factors. In this paper, a three-part procedure for applying numerical continuation to aeroelastic systems with freeplay is proposed, designed to ensure that the complete periodic behavior is identified, even for systems with very complex bifurcation diagrams. First, the equivalent linearization approach is used to determine approximations to the periodic solutions of the nonlinear system. Then, a shooting-based technique is applied separately to each linearized approximation in order to pinpoint the nearest exact periodic solution. This process results in a cloud of periodic solutions, representing points on all the solution branches and sub-branches. Finally, a branch-following shooting procedure is applied to this cloud of points in order to obtain a complete description of every branch of periodic solutions. The methodology is applied to a simple aeroelastic system with three degrees of freedom and freeplay in the control surface. This system has been often studied but never fully characterised. It is shown that the proposed method succeeds in describing the complete bifurcation behaviour of the system and explaining its limit cycle response.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A combined Multiple Time Scales and Harmonic Balance approach for the transient and steady-state response of nonlinear aeroelastic systems
Berci, Marco; Dimitriadis, Grigorios
2018 • In Journal of Fluids and Structures, 80, p. 132-144
The majority of methods for calculating the dynamic response of nonlinear aeroelastic systems considers only steady-state periodic behaviour. An exception is the Multiple Time Scales method, which can estimate both the transient and steady-state solutions of such systems; nevertheless, this approach is only accurate close to the Hopf bifurcation. This paper proposes a novel combined approach whereby the transient response is obtained from the Multiple Time Scales method and the asymptotic periodic behaviour is corrected using the Harmonic Balance method. This consistent and efficient framework mutually empowers both techniques and accounts for large parameter variations around the critical condition. The effect of cubic aero-structural nonlinearities on the dynamic response of a generic aeroelastic system is then investigated. Both the Multiple Time Scales and Harmonic Balance methods are adopted and perfect agreement of the explicit results is demonstrated, albeit near the system instability. In contrast, the proposed combined solution is valid for a wider range of perturbations, is analytical and has negligible computational cost while retaining accuracy. The role of key parameters and terms on the core mechanism of the dynamic behaviour is rigorously identified and discussed, from both physical and mathematical points of view. Galloping is finally considered as the simplest but complete application to a fundamental yet practical problem, featuring full conceptual complexity while exploiting the solid synthesis capability of the newly proposed analytical approach. Excellent agreement was found in all cases with results from the numerical integration of the nonlinear equations of motion in the time domain.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Experimental and numerical investigations of the torsional flutter oscillations of a 4:1 rectangular cylinder
Andrianne, Thomas; Dimitriadis, Grigorios
2013 • In Journal of Fluids and Structures, 41, p. 64-88
The torsional flutter oscillations of a 4:1 rectangular cylinder around its
pitching axis are investigated through wind tunnel experiments and numerical
simulations. The rectangle’s responses to different initial conditions and
turbulence excitations at various wind tunnel airspeeds are recorded. Timeresolved
Particle Image Velocimetry measurements are taken at two different
airspeeds, when the rectangle undergoes Limit Cycle Oscillations. Aeroelastic
simulations are carried out using the Discrete Vortex Method and the
resulting responses are compared to the experimental measurements.
The Common-base Proper Orthogonal Decomposition method is used to
analyse and compare the measured and simulated unsteady flow fields around
the rectangle. A discussion of the participation of each mode in the different
states of the flow-field is presented, at two different amplitudes of oscillation.
The Motion Induced Vortex (MIV) is identified as the fundamental cause
of the torsional flutter phenomenon and its role over a complete cycle is
studied. MIV-induced oscillations can be started either by a suitable initial
disturbance or by a second, nearly linear self-excited instability that causes negative aerodynamic damping. The combination of these two instabilities
results in a complete description of the torsional flutter of the rectangle.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Diagnosis of Process Faults in Chemical Systems Using a Local Partial Least Squares Approach
Kruger, Uwe; Dimitriadis, Grigorios
2008 • In AIChE Journal, 54 (10), p. 2581-2596
This article discusses the application of partial least squares (PLS) for monitoring
complex chemical systems. In relation to existing work, this article proposes the integration
of the statistical local approach into the PLS framework to monitor changes in
the underlying model rather than analyzing the recorded input/output data directly. As
discussed in the literature, monitoring changes in model parameters addresses the
problems of nonstationary behavior and presents an analogy to model-based
approaches. The benefits of the proposed technique are that (i) a detailed mechanistic
plant model is not required, (ii) nonstationary process behavior does not produce false
alarms, (iii) parameter changes can be non-Gaussian, (iv) Gaussian monitoring statistics
can be established to simplify the monitoring task, and (v) fault magnitude and signatures
can be estimated. This is demonstrated by a simulation example and the analysis
of recorded data from two chemical processes.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Damping identification of lightly damped linear dynamic systems using Common-base Proper Orthogonal Decomposition
Andrianne, Thomas; Dimitriadis, Grigorios
2012 • In Mechanical Systems and Signal Processing, 28, p. 492-506
This paper presents a new technique to identify the damping of linear systems. It is developed from the Proper Orthogonal Decomposition (POD) of the free response of the system and extended to the recently proposed Common-base POD (CPOD). The present application of CPOD considers simultaneously several free responses of the system to different initial conditions. The eigen-decomposition of the co-variance matrix leads to a unique vector basis which is likely to contain more information about the dynamics of the system than a vector basis obtained by the classic POD technique.
The ability of the technique to estimate the mode shapes and the modal damping is demonstrated on a simulated mass-spring-damper system. Two different distributions of masses are considered in order to confront the CPOD analysis to the intrinsic limitation of POD, i.e. that the mode shapes are identified exactly only if the mass matrix is proportional to the identity matrix. It is shown that the identification of the damping is still possible when the modes are not orthonormal. The robustness of the technique is demonstrated in the presence of noise in the responses of the system and through an experimental application with comparison with other identifications techniques.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Unsteady Lifting Line Theory Using theWagner Function for the Aerodynamic and Aeroelastic Modeling of 3D Wings
Boutet, Johan; Dimitriadis, Grigorios
2018 • In Aerospace, 5 (3), p. 92
A method is presented to model the incompressible, attached, unsteady lift and pitching
moment acting on a thin three-dimensional wing in the time domain. The model is based on
the combination of Wagner theory and lifting line theory through the unsteady Kutta–Joukowski
theorem. The results are a set of closed-form linear ordinary differential equations that can be solved
analytically or using a Runge–Kutta–Fehlberg algorithm. The method is validated against numerical
predictions from an unsteady vortex lattice method for rectangular and tapered wings undergoing
step or oscillatory changes in plunge or pitch. Further validation is demonstrated on an aeroelastic
test case of a rigid rectangular finite wing with pitch and plunge degrees of freedom.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Empirical modelling of the bifurcation behaviour of a bridge deck undergoing across-wind galloping
Andrianne, Thomas; Dimitriadis, Grigorios
2014 • In Journal of Wind Engineering and Industrial Aerodynamics, 135, p. 129-135
This work presents an empirical model capable to describe the galloping bifurcation behaviour of a bridge deck. It is based on a general polynomial form proposed by Novak, which we limit to the 5th order. The advantage of choosing this function for modelling the vertical force coefficient is that asymmetry of the even terms is enforced in order to reproduce the sub-critical aeroelastic behaviour of the bridge deck. The coefficients of the polynomial are identified from several pairs of displacement amplitudes and the corresponding airspeeds, measured in a wind tunnel during dynamic tests on the sectional bridge model. The identification is carried out using a first order harmonic balance technique. A stability analysis is presented in order to highlight the need of such a model to catch the complete bifurcation behaviour of the system. The resulting force coefficient of this full order model is compared to the well known models of Parkinson and Novak. Finally, the concept universal curve is used in order to discuss the galloping responses of square and rectangular cylinders, in comparison with the one of the bridge deck.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Induced Drag Calculations with the Unsteady Vortex Lattice Method for Cambered Wings
Lambert, Thomas; Dimitriadis, Grigorios
2017 • In AIAA Journal, 55 (2), p. 668-672
The Unsteady Vortex Lattice Method (UVLM) is an approach widely used to estimate the aerodynamic loads in unsteady subsonic flows. It is based on modeling the camber surface of a lifting body by means of bound vortex rings. Even though this method has been known and used for several decades, there is little discussion of the modeling of the leading-edge suction in the literature. To address this concern, Simpson et al. [1] presented a comparison of two different ways to model this effect for the case of uncambered airfoils and wings in harmonic pitch or plunge motions. They concluded that the Joukowski method converges significantly faster than the Katz technique as the number of chorwise panels is increased. The present paper is an extension of the study by Simpson et al. to cambered lifting surfaces. It shows that the presence of camber can change radically the convergence performance of the two methods. For cambered wings, the Katz approach converges significantly faster than the Joukowski technique.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Bifurcation Behavior of Airfoil Undergoing Stall Flutter Oscillations in Low-Speed Wind Tunnel
Dimitriadis, Grigorios; Li, Jing
2009 • In AIAA Journal, 47 (11), p. 2577-2596
Stall flutter is a nonlinear aeroelastic phenomenon that can affect several types of aeroelastic systems such as helicopter rotor blades, wind turbine blades, and highly flexible wings. Although the related aerodynamic
phenomenon of dynamic stall has been the subject of many experimental studies, stall flutter itself has rarely been investigated. This paper presents a set of experiments conducted on a NACA0012 airfoil undergoing stall flutter oscillations in a low-speed wind tunnel. The aeroelastic responses are analyzed with the objective of characterizing the local bifurcation behavior of the system. It is shown that symmetric stall flutter oscillations are encountered as a result of a subcritical Hopf bifurcation, followed by a fold bifurcation. The cause of these bifurcations is the occurrence of dynamic stall, which allows the transfer of energy from the freestream to the wing. A second bifurcation occurs at the system’s static divergence airspeed. As a consequence, the wing starts to undergo asymmetric stall flutter bifurcations at only positive (or only negative) pitch angles. The dynamic stall mechanism itself does not change but the flow only separates on one side of the wing.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Application of a 3D unsteady surface panel method with flow separation model to horizontal axis wind turbines
Prasad, Chandra Shekhar; Dimitriadis, Grigorios
2017 • In Journal of Wind Engineering and Industrial Aerodynamics, 166, p. 74-89
This work describes the development and application of a 3D unsteady surface panel method with a separation model to the problem of simulating the flow around the blades of a horizontal axis wind turbine. The present method is intended as a design tool to capture the 3D time-dependent characteristics of both attached and separated flow conditions and is an extension of previous 2D approaches. Flow separation is modelled using a loose coupling procedure between the inviscid panel method and a quasi-3D viscous boundary layer solution. A separated wake is shed at the predicted separation points and propagated at the local flow velocity, just like the trailing edge wake. The methodology is demonstrated on the NREL phase-VI wind turbine test case and the model predictions are compared to experimental measurements.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Demonstrating the identification of nonlinear vibrating systems to undergraduate students
Dimitriadis, Grigorios; Vio, Gareth Arthur
2007 • In International Journal of Mechanical Engineering Education, 35 (4), p. 336-360
The identification of nonlinear dynamic systems is increasingly becoming a necessary part of vibration testing and there is significant research effort devoted to it. However, as the current methodologies are still not suitable for the identification of general nonlinear systems the subject is very rarely introduced to undergraduate students. In this paper, recent progress in developing an expert approach to nonlinear system identification is used in order to demonstrate the subject within the context of an undergraduate course or as an introductory tool for postgraduate students. The demonstration is based around a software package of an Expert System designed to apply systematically a wide range of identification approaches to the system under investigation. It is shown that the software can be used to demonstrate the need for nonlinear system identification, the complexity of the procedure, the possibility of failure and the good chances of success when enough physical information about the system is available.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A method for identification of non-linear multi-degree-of-freedom systems
Dimitriadis, Grigorios; Cooper, Jonathan E
1998 • In Proceedings of the Institution of Mechanical Engineers. Part G, Journal of Aerospace Engineering, 212 (4), p. 287-298
System identification methods for non-linear aeroelastic systems could find uses in many aeroelastic applications such as validating finite element models and tracking the stability of aircraft during flight flutter testing. The effectiveness of existing non-linear system identification techniques is limited by various factors such as the complexity of the system under investigation and the type of non-linearities present. In this work, a new approach is introduced which can identify multi-degree-of-freedom systems featuring any type of non-linear function, including discontinuous functions. The method is shown to yield accurate identification of three mathematical models of aeroelastic systems containing a wide range of structural non-linearities.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Flutter Prediction from Flight Flutter Test Data
Dimitriadis, Grigorios; Cooper, Jonathan E
2001 • In Journal of Aircraft, 38 (2), p. 355-367
The most common approach to flight flutter testing is to track estimated modal damping ratios of an aircraft over a number of flight conditions. These damping trends are then extrapolated to predict whether it is safe to move to the next test point and also to determine the utter speed. In the quest for more reliable and efficient flight flutter testing procedures, a number of alternative data analysis methods have been proposed. Five of these approaches are compared on two simulated aeroelastic models. The comparison is based on both the accuracy of prediction and the efficiency of each method. It is found that, for simple aeroelastic systems, the Nissim and Gilyard method (Nissim, E., and Gilyard, G. B., “Method for Experimental Determination of Flutter Speed by Parameter Identification,” AIAA Paper 89-1324, 1989) yields the best flutter predictions and is also the least computationally expensive approach.However, for larger systems, simpler approaches such as the damping fit and envelope function methods are found to be most reliable.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Characterization of the Behaviour of a Simple Aeroservoelastic System with Control Nonlinearities
Dimitriadis, Grigorios; Cooper, Jonathan E
2000 • In Journal of Fluids and Structures, 14 (8), p. 1173-1193
The characterization of the behaviour of nonlinear aeroelastic systems has become a very important research topic in the Aerospace Industry. However, most work carried to-date has concentrated upon systems containing structural or aerodynamic nonlinearities. The purpose of this paper is to study the stability of a simple aeroservoelastic system with nonlinearities in the control system and power control unit. The work considers both structural and control law nonlinearities and assesses the stability of the system response using bifurcation diagrams. It is shown that simple feedback systems designed to increase the stability of the linearized system also stabilize the nonlinear system, although their effects can be less pronounced. Additionally, a nonlinear control law designed to limit the control surface pitch response was found to increase the flutter speed considerably by forcing the system to undergo limit cycle oscillations instead of fluttering. Finally, friction was found to affect the damping of the system but not its stability, as long as the amplitude of the frictional force is low enough not to cause stoppages in the motion.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Limit Cycle Oscillation Control and Suppression
Dimitriadis, Grigorios; Cooper, Jonathan E
1999 • In Aeronautical Journal, 103 (1023), p. 257-263
The prediction and characterization of the Limit Cycle Oscillation (LCO) behaviour of nonlinear aeroelastic systems has become of great interest recently. However, much of this work has concentrated on determining the existence of LCOs. This paper concentrates on LCO stability. By considering the energy present in di®erent limit cycles, and also using the Harmonic Balance Method, it is shown how the stability of limit cycles can be determined. The analysis is then extended to show that limit cycles can be controlled, or even suppressed,
by the use of suitable excitation signals. A basic control scheme is developed to achieve this, and is demonstrated on a simple simulated nonlinear aeroelastic system.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Experimental study of wings undergoing active root flapping and pitching
Abdul Razak, Norizham; Dimitriadis, Grigorios
2014 • In Journal of Fluids and Structures, 49, p. 687-704
This paper presents the results of experiments carried out on mechanical wings undergoing active root flapping and pitching in the wind tunnel. The objective of the work is to investigate the effect of the pitch angle oscillations and wing profile on the aerodynamic forces generated by the wings. The experiments were repeated for different reduced frequency, airspeed, flapping and pitching kinematics, geometric angle of attack and wing sections (one symmetric and two cambered airfoils). A specially designed mechanical flapper was used, modelled on large migrating birds. It is shown that, under pitch leading conditions, good thrust generation can be obtained at a wide range of Strouhal numbers if the pitch angle oscillation is adjusted accordingly. Consequently, high thrust was measured at both the lowest and the highest tested Strouhal numbers. Furthermore, the work demonstrates that the aerodynamic forces can be sensitive to the Reynolds number, depending on the camber of the wings. Under pitch lagging conditions, where the effective angle of attack amplitude is highest, the symmetric wing was affected by the Reynolds number, generating less thrust at the lowest tested Reynolds value. In contrast, under pure flapping conditions, where the effective angle of attack amplitude was lower but still significant, it was the cambered wings that demonstrated Reynolds sensitivity.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Subsonic source and doublet panel methods
Sanchez Martinez, Mariano; Dimitriadis, Grigorios
2022 • In Journal of Fluids and Structures
This work presents a source and doublet surface panel method for solving unsteady subsonic flows around wing and bodies. The approach is based on Morino's theory but features three modifications. The impermeability boundary condition is defined as zero mass flux normal to the surface instead of zero normal velocity, a new methodology is introduced to calculate the aerodynamic influence coefficient matrices and the pressure distribution around the surface is evaluated using a nonlinear unsteady Bernoulli equation. The method is demonstrated on three experimental test cases from the literature; it is shown that it can predict more accurately the steady and oscillating pressure distribution around the surface than Morino's original linear approach. Crucially, the modified approach can represent asymmetries in the pressure distribution due to both asymmetric motion and shape. While higher fidelity alternatives exist for solving such flows, the present method has very low computational cost and would be suitable for optimization procedures during preliminary design, as long as the flow remains subcritical and attached at all times.
Peer Reviewed verified by ORBi
Short communication (Scientific journals)
Comment on Flutter Prediction from Flight Flutter Test Data
Dimitriadis, Grigorios; Cooper, Jonathan E
2006 • In Journal of Aircraft, 43 (3), p. 862-863
In a previous paper entitled “Flutter Predictions from Flight Flutter Test Data” the authors applied a number of different flutter prediction methods to data from two simulated aeroelastic aircraft models and compared the resulting flutter predictions. The two simulated models were a simple three-degree-of-freedom Hancock wing
model and the Sim-2 model of a generic four-engined civil transport. One of the methods examined in the paper was the Nissim and Gilyard method (NGM). Because of difficulties encountered with the Sim-2 model, the authors failed to apply the NGM successfully to it, and only results for the Hancock model were presented
in the paper. With the aid of Eli Nissim, the authors have now succeeded in applying the method to the Sim-2 model and to obtain quality flutter predictions from it. In this short Comment, the initial problems encountered will be described and then the solutions will be outlined. Finally, flutter predictions for the method will be
presented and compared to the flutter predictions obtained from the other methods by Dimitriadis and Cooper.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A time–frequency technique for the stability analysis of impulse responses from nonlinear aeroelastic systems
Dimitriadis, Grigorios; Cooper, Jonathan E
2003 • In Journal of Fluids and Structures, 17 (8), p. 1181–1201
A time–frequency method is proposed for the analysis of response time histories from nonlinear aeroelastic systems. The approach is based on a time-varying curve-fit of the short time Fourier transform of the impulse response. It is shown that the method can be used in order to obtain a clear picture of the sub-critical stability of a number of aeroelastic systems with a variety of structural and aerodynamic nonlinearities. Additionally, frequency and amplitude information can be obtained for both the linear and nonlinear signatures of the response signals in the sub- and postcritical regions. Finally, it is shown that, given certain types of nonlinear functions, sub-critical damping trends can be extrapolated to predict bifurcation airspeeds.
Peer Reviewed verified by ORBi
Article (Scientific journals)
The impact of circularity defects on bridge stay cable dry galloping stability
2018 • In Journal of Wind Engineering and Industrial Aerodynamics, 181, p. 14-26
The present work studies the influence of circularity defects, on the aerodynamic behaviour of stay cables of cable-stayed bridges. It focuses on wind tunnel tests on High-Density Polyethylene cable covers with and without helical fillets, in a range of Reynolds numbers from the sub-critical to the critical regime. The paper considers the impact of circularity defects on the aerodynamic stability of cable sheaths by testing various amplitudes of imposed ovalization. The defects are artificially applied on real sheaths whose original cross-sections are close to circular. The experiment consists in measuring surface pressures to investigate how the amplitude of ovalization influences the flow around the sheaths, especially in the critical Reynolds number regime when transition in the boundary layer occurs. The analysis is based on bifurcation diagrams and Proper Orthogonal Decomposition. The investigation demonstrates that important circularity defects can significantly increase the bi-stable nature of the flow around a sheath at the critical regime. Nevertheless, the introduction of a helical fillet de-correlates the flow around the sheath, causing jumps in lift that have different sign along its length.
Peer Reviewed verified by ORBi
Article (Scientific journals)
On the aerodynamic characteristics of inclined stay cable sheaths in the critical flow regime
2020 • In Journal of Fluids and Structures, 99, p. 103147
This paper studies the presence of the bistable flow activity around inclined cable models. It presents results from wind tunnel test on static original High-Density Polyethylene cable covers in a range of Reynolds numbers from the subcritical to the critical regime, inclined at angles of 90° (vertical), 60° and 45°. It has already been shown that, in the critical regime, turbulent transition in the boundary layers around a circular cylinder exhibits bistable behavior at zero inclination, i.e. the boundary layer on either side is intermittently turbulent or laminar, leading to intermittent asymmetry of the flow and resulting aerodynamic loads. The present wind tunnel tests consist in measuring the pressure patterns around inclined circular cylinders using numerous pressure taps. Bifurcation diagrams are created in order to quantify the occurrence of bistability and Proper Orthogonal Decomposition is used to identify the geometrical flow modes that are governed by the bistability. The results show that the bistable phenomenon is only identified in the Reynolds number range corresponding to the first transition between the TrBL0 and TrBL1 regimes. Furthermore, the inclination angle seems to have a significant impact on the energy of the flow mode related to unsteady asymmetric reattachment. This paper also treats the potential sensitivity of the aerodynamic drag and lift forces to the circularity defect of the covers and presents the positive impact of the helical fillet on controlling the bistable activity at a critical inclination angle of 60°.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A continuation method for determining the speed dependent modal properties of large MDOF aeroelastic systems
Heremans, Julien; Dimitriadis, Grigorios; Denoël, Vincent
2024 • In Journal of Wind Engineering and Industrial Aerodynamics, 252, p. 105804
Peer Reviewed verified by ORBi
Article (Scientific journals)
A modified Leishman-Beddoes model for airfoil sections undergoing dynamic stall at low Reynolds numbers
2020 • In Journal of Fluids and Structures, 93, p. 102852
A modified Leishman-Beddoes[1] model is proposed to predict the aerodynamic load responses of airfoils undergoing dynamic stall at low Reynolds, low Mach numbers and low to very high equivalent reduced pitch rates. The modifications include using Wagner function aerodynamics for the attached flow model, defining and using a time-varying reduced pitch rate, equating the time delays applied to the normal force, the separation point location and the vortex shedding and using Sheng and Galbraith's dynamic stall onset criteria [2]. The steady and unsteady force and pitching moment measurements of three airfoils with distinct stall mechanisms (a flat plate, a NACA0012, and a NACA0018) are used to demonstrate the validity of this new model. Dynamic tests consist of oscillating the airfoils in pitch around the quarter chord with a mean angle A0 = 10 deg and with different prescribed reduced frequencies from k = 0.015, to k = 0.16, and amplitudes from A = 5 deg to 20 deg, at Reynolds numbers of the order of Re=1.8x104. The resulting equivalent reduced pitch rates range from r' = 0.001 to 0.06. The computations of the values of the different model parameters are demonstrated and it is observed that the values of the time delay parameters change continuously and smoothly with time instead of jumping discontinuously at discrete time instances as proposed by Leishman and Beddoes [1]. It is shown that the modified Leishman-Beddoes model is in good agreement with experiments for low to moderate equivalent reduced pitch rates (r'<=0.04). It also significantly improves the dynamic load prediction in this range of pitch rates, in comparison to the original Leishman-Beddoes model. Nevertheless, as for the original model, the modified model presented here becomes inaccurate at the highest equivalent reduced pitch rates of r' > 0.05.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Vortex Lattice simulations of attached and separated flows around flapping wings
Lambert, Thomas; Abdul Razak, Norizham; Dimitriadis, Grigorios
2017 • In Aerospace, 4 (2), p. 22
Flapping flight is an increasingly popular area of research, with applications to micro-unmanned air vehicles and animal flight biomechanics. Fast but accurate methods for predicting the aerodynamic loads acting on flapping wings are of interest for designing such aircraft and optimising thrust production. In this work, the unsteady Vortex Lattice method is used in conjunction with three load estimation techniques in order to predict the aerodynamic lift and drag time histories produced by flapping rectangular wings. The load estimation approaches are the Katz, Joukowski and simplified Leishman-Beddoes techniques. The simulations' predictions are compared to experimental measurements from a flapping and pitching wing presented by Razak and Dimitriadis [1]. Three types of kinematics are investigated, pitch-leading, pure flapping and pitch lagging. It is found that pitch-leading tests can be simulated quite accurately using either the Katz or Joukowski approaches as no measurable flow separation occurs. For the pure flapping tests, the Katz and Joukowski techniques are accurate as long as the static pitch angle is greater than zero. For zero or negative static pitch angles these methods underestimate the amplitude of the drag. The Leishman-Beddoes approach yields better drag amplitudes but can introduce a constant negative drag offset. Finally, for the pitch-lagging tests the Leishman-Beddoes technique is again more representative of the experimental results, as long as flow separation is not too extensive. Considering the complexity of the phenomena involved, in the vast majority of cases the lift time history is predicted with reasonable accuracy. The drag (or thrust) time history is more challenging.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Flutter and stall flutter of a rectangular wing in a wind tunnel
Norizham, Abdul Razak; Andrianne, Thomas; Dimitriadis, Grigorios
2011 • In AIAA Journal, 49 (10), p. 2258-2271
The aeroelastic behavior of a rectangular wing with pitch and plunge degrees of freedom was observed experimentally using pressure, acceleration and PIV measurements. The wing was set at different static angles of attack and wind tunnel airspeeds. The wing's dynamic behavior was governed by a two-parameter bifurcation from steady to Limit Cycle Oscillations (LCO), the two parameters being the airspeed and the static angle of attack. At the lowest static angle, the wing underwent a classical flutter phenomenon that was transformed into a supercritical Hopf bifurcation at higher angles. The latter was combined with a fold bifurcation at intermediate angles of attack.
All LCOs observed were either low amplitude oscillations with time-varying amplitude or high amplitude oscillations with nearly steady amplitude. They were caused by two different types of dynamic stall phenomena. During low amplitude LCOs the periodically stalled flow covered only the rear part of the wing. During high amplitude LCOs, trailing edge and leading edge separation occured. Trailing edge separation was characterized by a significant amount of unsteadiness, varying visibly from cycle to cycle. The occurrence of leading edge separation was much more regular and had the tendency to stabilize the amplitude of the LCO motion.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Experimental passive flutter suppression using a linear tuned vibration absorber
The current drive for increased efficiency in aeronautic structures such as aircraft, wind turbine blades and helicopter blades often leads to weight reduction. A con- sequence of this tendency can be increased flexibility, which in turn can lead to un- favourable aeroelastic phenomena involving large amplitude oscillations and non- linear effects such as geometric hardening and stall flutter. Vibration mitigation is one of the approaches currently under study for avoiding these phenomena. In the present work, passive vibration mitigation is applied to a nonlinear experimental aeroelastic system by means of a linear tuned vibration absorber. The aeroelastic apparatus is a pitch and flap wing that features a continuously hardening restoring torque in pitch and a linear restoring torque in flap. Extensive analysis of the sys- tem with and without absorber at pre-critical and post-critical airspeeds showed an improvement in flutter speed of around 36%, a suppression of a jump due to stall flutter, and a reduction in LCO amplitude. Mathematical modelling of the exper- imental system is used to demonstrate that optimal flutter delay is achieved when two of the system modes flutter at the same flight condition. Nevertheless, even this optimal absorber quickly loses effectiveness as it is detuned. The wind tunnel mea- surements showed that the tested absorbers were much slower to lose effectiveness than those of the mathematical predictions.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A Modal Frequency-Domain Generalised Force Matrix for the Unsteady Vortex Lattice Method
Dimitriadis, Grigorios; Giannelis, N. F.; Vio, G. A.
2018 • In Journal of Fluids and Structures, 76, p. 216-228
The unsteady Vortex Lattice method is becoming an increasingly popular aerodynamic modelling method for incompressible aeroelastic problems, such as flexible low-speed aircraft, wind turbines and flapping flight. It leads to discrete time aeroelastic state space equations, which must be solved in a time-marching framework. Eigenvalue or singular value decompositions of the discrete time equations can be used in order to perform stability analysis but such procedures must be accompanied by model order reduction because the size of the equations is large. This work proposes a modal frequency domain implementation of the Vortex Lattice method, resulting in a modal generalised force matrix. Model order reduction is implicit in the modal approach and stability analysis can be carried out using industry-standard flutter analysis techniques, such as the p-k method. The approach is validated by comparison to wind tunnel flutter data obtained from rectangular cantilever flat plate wings of different aspect ratios and sweep angles. It is found that the aeroelastic model predictions follow the experimental trends for both flutter speed and frequency but tend to be moderately conservative.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Impact of roughness and circularity-defect on bridge cables stability
2015 • In Journal of Wind Engineering and Industrial Aerodynamics, 137, p. 1-13
The stay cables on cable-stayed bridges have been the subject of considerable research in order to understand the origins of the dry galloping phenomenon. Surface irregularity is one of the last parameters that have not been thoroughly examined. This paper focuses on static wind tunnel tests on original High Density Polyethylene cable covers in a range of Reynolds numbers from the sub-critical regime to the critical, corresponding to values ranging from Re = 9.6x104 to Re = 3.3x105. The experiment consists in testing cable covers of various diameters in order to investigate the effect of surface irregularity (roughness and circularity defect) on the mechanism of dry galloping excitation. Previous studies reported that dry galloping is caused by the appearance of a negative pressure bubble on one side of the circular cylinder at the critical Reynolds number range, leading to a rapid drop in the drag coefficient and the appearance of a non-negligible lift force. The results of the present investigation demonstrate that there is a clear correlation between the single bubble pressure pattern and the circularity defect along the tube. They further show that surface roughness has little effect on the location of the bubble. The paper also treats the spatial and temporal correlation of the instantaneous pressure pattern along the tube with respect to the circularity defect.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Unsteady aerodynamic modeling methodology based on dynamic mode interpolation for transonic flutter calculations
2019 • In Journal of Fluids and Structures, 84, p. 218-232
A new unsteady aerodynamic modeling methodology for calculating transonic flutter characteristics is presented. The main idea of the methodology is to obtain the unsteady flow response to small amplitude periodic deformations of a structure over a large range of oscillation frequencies through the interpolation of the most dominant fluid dynamic modes obtained from Dynamic Mode Decomposition (DMD) of a few reference unsteady simulations at different oscillation frequencies. These simulations can be carried out by solving the Euler or RANS equations. The methodology can then be used to obtain a frequency-domain generalized aerodynamic force matrix, and stability analysis can be performed using standard flutter analysis methods such as the p-k method. The proposed methodology provides a very good estimate of the flutter boundary for the 2D Isogai airfoil and 3D AGARD 445.6 wing models, but at a lower computational cost than the traditional higher-fidelity Fluid-Structure Interaction (FSI) simulations.
Peer Reviewed verified by ORBi
Article (Scientific journals)
The influence of flight style on the aerodynamic properties of avian wings as fixed lifting surfaces
Lees, John J.; Dimitriadis, Grigorios; Nudds, Robert L.
2016 • In PeerJ, 4:e2495
The diversity of wing morphologies in birds reflects their variety of flight styles and the associated aerodynamic and inertial requirements. Although the aerodynamics underlying wing morphology can be informed by aeronautical research, important differences exist between planes and birds. In particular, birds operate at lower, transitional Reynolds numbers than do most aircraft. To date, few quantitative studies have investigated the aerodynamic performance of avian wings as fixed lifting surfaces and none have focused upon the differences between wings from different flight style groups. Dried wings from 10 bird species representing 3 distinct flight style groups were mounted on a force/torque sensor within a wind tunnel in order to test the hypothesis that wing morphologies associated with different flight styles exhibit different aerodynamic properties. Morphological differences manifested primarily as differences in drag rather than lift. Maximum lift coefficients did not differ between groups, whereas minimum drag coefficients were lowest in undulating flyers (Corvids). The lift to drag ratios were lower than in conventional aerofoils and data from free-flying soaring species; particularly in high frequency, flapping flyers (Anseriformes), which do not rely heavily on glide performance. The results illustrate important aerodynamic differences between the wings of different flight style groups that cannot be explained solely by simple wing-shape measures. Taken at face value, the results also suggest that wing-shape is linked principally to changes in aerodynamic drag, but, of course, it is aerodynamics during flapping and not gliding that is likely to be the primary driver.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Numerical and experimental study of the flow around a 4:1 rectangular cylinder at moderate Reynolds number
2019 • In Journal of Wind Engineering and Industrial Aerodynamics, 189, p. 289-303
This paper presents the results of investigations into the flow around a rectangular cylinder with a chord-to-depth ratio equal to 4. The studies are performed through wind tunnel dynamic pressure measurements along a cross section combined with two-dimensional Unsteady Reynolds-Averaged Navier-Stokes (URANS) and three-dimensional Delayed-Detached Eddy Simulation (DDES). These experimental and numerical studies are complementary and combining them allows a better understanding of the unsteady dynamics of the flow. These studies aim mainly at determining the effects of the rectangle incidence and freestream velocity on the variation of the flow topology and the aerodynamic loads, and at assessing the capability of the industrially affordable URANS and DDES approaches to provide a sufficiently accurate estimation of the flow for different incidences. The comparison of experimental and numerical data is performed using statistics and Dynamic Mode Decomposition. It is shown that the rectangular cylinder involves complex separation-reattachment phenomena that are highly sensitive to the Reynolds number. In particular, the time-averaged lift slope increases rapidly with the Reynolds number in the range 7.8e3 < Re < 1.9e4 due to the modification of the time-averaged vortex strength, thickness and distance from the surface. Additionally, it is shown that both URANS and DDES simulations fail to accurately predict the flow at all the different incidence angles considered. The URANS approach is able to qualitatively estimate the spatio-temporal variations of vortices for incidences below the stall angle alpha = 4°. Nonetheless, URANS does not predict stall, while DDES correctly identifies the stall angle observed experimentally.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Space Debris Interaction Across a Two-dimensional Oblique Shock Wave
This work presents, compares and discusses results obtained with two indirect methods for the cal
culation of aerodynamic forces and pitching moment from 2D Particle Image Velocimetry (PIV) measurements. Both methodologies are based on formulations of the momentum balance: the integral Navier-Stokes equations and the “flux equation” proposed by Noca et al. (1999), which has been extended to the computation of moments. The indirect methods are applied to spatio-temporal data for different separated flows around a plate with a 16:1 chord-to-thickness ratio. Experimental data are obtained in a water channel for both a plate undergoing a large amplitude imposed pitching motion and a static plate at high angle of attack. In addition to PIV data, direct measurements of aerodynamic loads are carried out to assess the quality of the indirect calculations. It is found that indirect methods are able to compute the mean and the temporal evolution of the loads for two-dimensional flows with a reasonable accuracy. Nonetheless, both methodologies are noise sensitive and, the parameters impacting the computation should thus be chosen carefully. It is also shown that results can be improved through
the use of Dynamic Mode Decomposition (DMD) as a pre-processing step.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A Class of Methods for the Analysis of Blade Tip Timing Data from Bladed Assemblies Undergoing Simultaneous Resonances—Part I: Theoretical Development
Gallego-Garrido, Jon; Dimitriadis, Grigorios; Wright, Jan Robert
2007 • In International Journal of Rotating Machinery, 2007, p. 11
Blade tip timing is a technique for the measurement of vibrations in rotating bladed assemblies. Although the fundamentals of the technique are simple, the analysis of data obtained in the presence of simultaneously occurring synchronous resonances is problematic. A class of autoregressive-based methods for the analysis of blade tip timing data from assemblies undergoing two simultaneous resonances has been developed. It includes approaches that assume both sinusoidal and general blade tip responses. The methods can handle both synchronous and asynchronous resonances. An exhaustive evaluation of the approaches was performed on simulated data in order to determine their accuracy and sensitivity. One of the techniques was found to perform best on asynchronous resonances and one on synchronous resonances. Both methods yielded very accurate vibration frequency estimates under all conditions of interest.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Aerodynamic instability of circular slender structures due to bistable flow activity
2021 • In Structural Engineering International, 32 (1), p. 71-77
Slender structures with circular cylindrical cross-sections are very common in civil engineering, industrial equipment and bridge engineering. Such long cylinders subjected to wind are columns, piles, long cables, chimney, pylons and others. They usually have smooth surfaces and their shape is supposed to be circular, at first glance. It has been observed that if long cables are subjected to wind load, an aerodynamic phenomenon that generates instability can be experienced. This phenomenon (called dry galloping) has been recently theorized to be linked to an unpredictable jump state occurring once the boundary layers of the circular cylinder changes from laminar to turbulent. With a non-stationary character, the force generated at the critical flow regime is named bi-stable activity. In this work, the birth of bi-stable flow activity around a circular cylinder made of High-Density Polyethylene is studied experimentally. The experiment consists of testing a real cable sheath equipped with a large number of pressure taps in a wind tunnel. The pressure pattern is linearly decomposed and geometric modes governing the bi-stable activity as a function of circularity defect are created. The results show that the energy of the non-stationary load observed at the critical flow regime depends strongly on the circularity defect. Furthermore, the inclination angle of 60° confers a significant increase in the energy of the bi-stable flow mode when an unsteady asymmetric reattachment event occurs. The paper also highlights the variability of the mean drag coefficient as a function of the macroscopic defect.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Bifurcation analysis and limit cycle oscillation amplitude prediction methods applied to the aeroelastic galloping problem
Vio, Gareth Arthur; Dimitriadis, Grigorios; Cooper, Jonathan E
2007 • In Journal of Fluids and Structures, 23 (7)
A global stability and bifurcation analysis of the transverse galloping of a square section beam in a normal steady flow has been implemented. The model is an ordinary differential equation with polynomial damping nonlinearity. Six methods are used to predict bifurcation, the amplitudes and periods of the ensuing Limit Cycle Oscillations: (i) Cell mapping, {ii} Harmonic Balance, (iii) Higher Order Harmonic Balance,(iv) Centre Manifold linearization, (v) Normal Form and (vi) Numerical Continuation. The resulting stability predictions are compared with each other and with results obtained from numerical integration. The advantages and disadvantages of each technique are discussed. It is shown that, despite the simplicity of the system, only two of the methods succeed in predicting its full response spectrum. These are Higher Order Harmonic Balance and Numerical Continuation.
Peer Reviewed verified by ORBi
Article (Scientific journals)
The aerodynamics of big ears in the brown long-eared bat Plecotus auritus
2008 • In Acta Chiropterologica, 10 (2), p. 313-321
Wings are the most obvious adaptation bats have for powered flight and differences in wing morphology are known to correlate with flight behaviour. However, the function(s) of ancillary structures such as the ears and tail, which may also play an important role during flight, are less well understood. Here we constructed a simplified model of a bat body with ears based upon morphological measurements of a brown long-eared bat (Plecotus auritus) to examine the aerodynamic implications of flying with large ears. The forces and moments produced by the model were measured using a sensitive 6-component force and torque balance during wind tunnel testing. The large ears of the model bat produced positive lift as well as positive drag of the same order of magnitude. At small ears angles (0° to 10°), increasing the angle of the ears resulted in an increase of the lift-to-drag ratio. At higher ear angles (> 10°) separation of the flow occurred which caused a large decrease in the lift-to-drag ratio produced. To maximise the benefit from the ears (i.e., lift-to-drag ratio) our model predicts that a horizontal free flying P. auritus should hold its ears at an approximate angle of 10°. The results of the pitching moment coefficient are inconclusive in determining if the large ears are important as flight control structures. The additional drag produced by the ears has consequences for the foraging behaviour of P. auritus with reductions in its flight speed and foraging range.
Peer reviewed
Article (Scientific journals)
Experimental and Numerical Study of Mini-UAV Propeller Performance in Oblique Flow
2017 • In Journal of Aircraft, 54 (3), p. 1076-1084
This paper presents the modelling of the performance of small propellers used for Vertical Take Off and Landing Micro Aerial Vehicles (VTOL MAVs) operating at low Reynolds numbers and in oblique flow. Blade Element Momentum Theory (BEMT), Vortex Lattice Method (VLM) and momentum theory for oblique flow are used to predict propeller performance. For validation, the predictions for a commonly used propeller for VTOL MAVs are compared to a set of wind tunnel experiments. Both BEMT and VLM succeed in predicting correct trends of the forces and moments acting upon the propeller shaft, although accuracy decreases significantly in oblique flow. For the dataset analysed here, combining the available data of the propeller in purely axial flow with momentum theory for oblique flow and applying a correction factor for the wake skew angle results in more accurate performance estimates at all elevation angles.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Flutter calculations using the unsteady source and doublet panel method
The Source and Doublet Panel Method (SDPM) developed by Morino in the 1970s can model unsteady compressible ideal flow around wings and bodies. In this work, the SDPM is adapted to the calculation of aeroelastic solutions for wings. A second order nonlinear version of Bernoulli's equation is transformed to the frequency domain and written in terms of the generalized mode shapes and displacements. It is shown that the pressure component at the oscillating frequency is a linear function of the generalized displacements, velocities and accelerations and can therefore be used to formulate a linear flutter problem. The proposed approach has several advantages over the usual Doublet Lattice Method (DLM) approach: the exact geometry is modelled, including thickness, camber and twist effects, the motion of the surface can be represented using all six degrees of freedom, the pressure calculation is of higher order and the aerodynamic mass, damping and stiffness terms are calculated explicitly. The complete procedure is validated using the experimental data from the weakened AGARD 445.6 wing, a NACA 0012 rectangular wing with pitch and plunge degrees of freedom and an experimental model of a T-tail, yielding flutter predictions that lie closer to the experimental observations than those obtained from the DLM, particularly in the case of the T-tail.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A Class of Methods for the Analysis of Blade Tip Timing Data from Bladed Assemblies Undergoing Simultaneous Resonances—Part II: Experimental Validation
Gallego-Garrido, Jon; Dimitriadis, Grigorios; Carrington, Ian B.et al.
2007 • In International Journal of Rotating Machinery, 2007, p. 10
Blade tip timing is a technique for the measurement of vibrations in rotating bladed assemblies. In Part I of this work a class of methods for the analysis of blade tip timing data from bladed assemblies undergoing two simultaneous synchronous resonances was developed. The approaches were demonstrated using data from a mathematical simulation of tip timing data. In Part II the methods are validated on an experimental test rig. First, the construction and characteristics of the rig will be discussed. Then, the performance of the analysis techniques when applied to data from the rig will be compared and analysed. It is shown that accurate frequency estimates are obtained by all the methods for both single and double resonances. Furthermore, the recovered frequencies are used to calculate the amplitudes of the blade tip responses. The presence of mistuning in the bladed assembly does not affect the performance of the new techniques.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Reduced Order Analysis of Aeroelastic Systems with Freeplay using an Augmented Modal Basis
Silva, Gustavo H. C.; Dal Ben Rossetto, Gustavo; Dimitriadis, Grigorios
2015 • In Journal of Aircraft, 52 (4), p. 1312-1325
Freeplay in aircraft control surface actuators is a common source of nonlinearity that can cause undesirable aeroelastic phenomena and, as a consequence, certification authorities place strict limits on the size of the deadband gaps of actuators. In recent years many authors have worked on determining the effects of freeplay on the dynamic behavior of aeroelastic systems; these efforts have yielded significant improvements in the understanding of the phenomena, but have concentrated on simple models. Investigations of industry-relevant aeroelastic models of complete aircraft with freeplay in the actuators have been much rarer, the main difficulty being the selection of an appropriate modal basis that can represent the full range of possible dynamic phenomena. This work presents a novel implementation of the residual vectors approach to create consistent and numerically stable reduced order models suitable for industrial-standard aeroelastic models of aircraft. The methodology relies on the piecewise linear nature of freeplay to create a single reduction basis that represents the dynamics of the system both inside and outside the freeplay deadband gap. The method is demonstrated on an Embraer generic test bench aircraft, showing that the resulting reduced order model is efficient and effective. Nonlinear analysis is carried out using equivalent linearization as well as time integrations of the full nonlinear system. It is shown that, while equivalent linearization is an indispensable tool for a preliminary mapping of the main aeroelastic instabilities, the time integration-based nonlinear analysis is an essential complementary tool to confirm the characteristics of the system behavior.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Experimental and numerical study of the flight of geese
2015 • In Aeronautical Journal, 119 (1217), p. 1-30
The flight of barnacle geese at airspeeds representing high-speed migrating flight is investigated using experiments and simulations. The experimental part of the work involved the filming of three barnacle geese (Branta Leucopsis) flying at different airspeeds in a wind tunnel. The video footage was analysed in order to extract the wing kinematics. Additional information, such as wing geometry and camber was obtained from a 3D scan of a dried wing. An unsteady vortex lattice method was used to simulate the aerodynamics of the measured flapping motion. The simulations were used in order to successfully reproduce the measured body motion and thus obtain estimates of the aerodynamic forces acting on the wings. It was found that the mean of the wing pitch angle variation with time has the most significant effect on lift while the difference in the durations of the upstroke and downstroke has the major effect on thrust. The power consumed by the aerodynamic forces was also estimated; it was found that increases in aerodynamic power correspond very closely to climbing motion and vice versa. Root-mean-square values of the power range from 100 W to 240 W. Finally, it was observed that tandem flying can be very expensive for the trailing bird.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Numerical investigation into discontinuity-induced bifurcations in an aeroelastic system with coupled non-smooth nonlinearities
Vishal, Sai; Raaj, Ashwad; Bose, Chandanet al.
2022 • In Nonlinear Dynamics
The present study focuses on investigating the bifurcation characteristics of a pitch–plunge aeroelastic system possessing coupled non-smooth nonlinearities, both in structural and aerodynamic fronts. To this end, a freeplay nonlinearity is considered in the stiffness of the pitch degree-of-freedom. The effects of dynamic stall arising due to large instantaneous angles-of-attackare incorporated using the semi-empirical Leishman–Beddoes aerodynamic model. A systematic response analysis is carried out to discern the bifurcation characteristics of the aeroelastic system considering the airspeed as the system parameter. At low airspeeds, a series of dynamical transitions, including aperiodic responses, occur predominantly due to the structural freeplay nonlinearity while the flow remains attached to the surface of the wing. However, beyond a critical value of airspeed, the system response is dominated by high amplitude pitch-dominated limit-cycle oscillations, which can be attributed to stall flutter. It is demonstrated that the freeplay gap plays a key role in combining the effects of structural and aerodynamic nonlinearities. At higher values of the freeplay gap, interesting discontinuity-induced bifurcation scenarios, such as grazing and boundary equilibrium bifurcations arise due to coupled nonlinear interactions, which can significantly impact the safety of the aeroelastic system.
Peer reviewed
Article (Scientific journals)
On the use of control surface excitation in flutter testing
Wright, Jan R; Wong, Jerry; Cooper, Jonathan Eet al.
2003 • In Proceedings of the Institution of Mechanical Engineers. Part G, Journal of Aerospace Engineering, 217 (6), p. 317-332
Flutter testing is aimed at demonstrating that the aircraft flight envelope is flutter free.
Response measurements from deliberate excitation of the structure are used to identify and track
frequency and damping values against velocity. In this paper, the common approach of using a flight
control surface to provide the excitation is examined using a mathematical model of a wing and
control surface whose rotation is restrained by a simple actuator. In particular, it is shown that it is
essential to use the demand signal to the actuator as a reference signal for data processing. Use of the
actuator force (or strain) or control angle (or actuator displacement) as a reference signal is bad
practice because these signals contain response information. It may also be dangerous in that the
onset of flutter may not be seen in the test results.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Two-domain and three-domain limit cycles in a typical aeroelastic system with freeplay in pitch
Verstraelen, Edouard; Dimitriadis, Grigorios; Dal Ben Rossetto, Gustavoet al.
2017 • In Journal of Fluids and Structures, 69, p. 89-107
Freeplay is a significant source of nonlinearity in aeroelastic systems and is strictly regulated by airworthiness authorities. It splits the phase plane of such systems into three piecewise linear subdomains. Depending on the location of the freeplay, limit cycle oscillations can result that span either two or three of these subdomains. The purpose of this work is to demonstrate the existence of two-domain cycles both theoretically and experimentally. A simple aeroelastic system with pitch, plunge and control deflection degrees of freedom is investigated in the presence of freeplay in pitch. It is shown that two-domain and three-domain cycles can result from a grazing bifurcation and propagate in the decreasing airspeed direction. Close to the bifurcation, the two limit cycle branches interact with each other and aperiodic oscillations ensue. Equivalent linearization is used to derive the conditions of existence of each type of limit cycle and to predict their amplitudes and frequencies. Comparisons with measurements from wind tunnel experiments demonstrate that the theory describes these phenomena with accuracy.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A potential role for bat tail membranes in flight control
Gardiner, James D; Dimitriadis, Grigorios; Codd, Jonathan Ret al.
2011 • In PLoS ONE, 6 (3), p. 18214
Wind tunnel tests conducted on a model based on the long-eared bat Plecotus auritus,
indicated that the position of the tail membrane (uropatagium) can significantly influence flight control. Adjusting tail position by increasing the angle of the legs relative to the body, has a two-fold effect; increasing leg-induced wing camber (i.e. locally increased camber and angle of attack of the inner wing surface) and increasing the angle of attack of the tail membrane. We also used our model to examine the effects of flying with and without a tail membrane. For the bat model with a tail membrane increasing leg angle increased the lift, drag and nose-down pitching moment produced. However, removing the tail membrane significantly reduced the change in pitching moment with increasing leg angle, but it had a much smaller effect on the level of lift and drag produced. The tail membrane, therefore, is potentially important for controlling the level of pitching moment produced by bats and an aid to flight control, specifically improving agility and manoeuvrability. Although the tail of bats is different from that of birds, in that it is only divided from the wings by the legs, it nonetheless, may, in addition to its prey capturing function, fulfil a similar role in aiding flight control.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Blade-Tip Timing Measurement Of Synchronous Vibrations Of Rotating Bladed
Dimitriadis, Grigorios; Carrington, Ian B.; Wright, Jan Robertet al.
2002 • In Mechanical Systems and Signal Processing, 16 (4), p. 599-622
Blade-tip timing (BTT) is a promising method for the detection, measurement and analysis of blade vibrations in rotating bladed assemblies. However, the intricacies of the method when applied to real rotating structures undergoing synchronous (Engine Ordered) vibrations are not yet fully understood. In this paper, a mathematical model is developed to simulate data from typical BTT tests of rotating assemblies. The simulator is then used in order to provide a qualitative analysis of several phenomena that can be associated with the synchronous vibrations of rotating assemblies, including mistuning, coupling, excitation at multiple Engine Orders and simultaneous synchronous and asynchronous responses. It is concluded that none of these phenomena on its own will render the identification of the frequency and amplitude of blade vibrations impossible. However, there is no single BTT data analysis method that is able to deal with all of these phenomena.
Peer Reviewed verified by ORBi
Article (Scientific journals)
The aerodynamic cost of head morphology in bats: maybe not as bad as it seems
Vanderelst, Dieter; Peremans, Herbert; Abdul Razak, Norizhamet al.
2015 • In PLoS ONE, 10 (3), p. 0118545
At first sight, echolocating bats face a difficult trade-off. As flying animals, they would benefit from a streamlined geometric shape to reduce aerodynamic drag and increase flight efficiency. However, as echolocating animals, their pinnae generate the acoustic cues necessary for navigation and foraging. Moreover, species emitting sound through their nostrils often feature elaborate noseleaves that help in focussing the emitted echolocation pulses. Both pinnae and noseleaves reduce the streamlined character of a bat’s morphology. It is generally assumed that by compromising the streamlined charactered of the geometry, the head morphology generates substantial drag, thereby reducing flight efficiency. In contrast, it has also been suggested that the pinnae of bats generate lift forces counteracting the detrimental effect of the increased drag. However, very little data exist on the aerodynamic properties of bat pinnae and noseleaves. In this work, the aerodynamic forces generated by the heads of seven species of bats, including noseleaved bats, are measured by testing detailed 3D models in a wind tunnel. Models of Myotis daubentonii, Macrophyllum macrophyllum, Micronycteris microtis, Eptesicus fuscus, Rhinolophus formosae, Rhinolophus rouxi and Phyllostomus discolor are tested. The results confirm that non-streamlined facial morphologies yield considerable drag forces but also generate substantial lift. The net effect is a slight increase in the lift-to-drag ratio. Therefore, there is no evidence of high aerodynamic costs associated with the morphology of bat heads
Peer Reviewed verified by ORBi
Article (Scientific journals)
Subcritical, nontypical and period doubling bifurcations of a Delta Wing in a low speed wind tunnel
2011 • In Journal of Fluids and Structures, 27 (3), p. 408-426
Limit Cycle Oscillations (LCOs) involving Delta wings are an important area of research in modern aeroelasticity. Such phenomena can be the result of geometric or aerodynamic nonlinearity. In this paper, a flexible half-span Delta wing is tested in a low speed wind tunnel in order to investigate its dynamic response. The wing is designed to be more flexible than the models used in previous research on the subject in order to expand the airspeed range in which LCOs occur. The experiments reveal that this wing features a very rich bifurcation behavior. Three types of bifurcation are observed for the first time for such an aeroelastic system: subcritical bifurcations, period doubling/period halving and nontypical bifurcations. They give rise to a great variety of LCOs, even at very low angles of attack.The LCOs resulting from the nontypical bifurcation display Hopf-type behavior, i.e. have fundamental frequencies equal to one of the linear modal frequencies. All of the other LCOs have fundamental frequencies that are unrelated to the underlying linear system modes.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Identification of multi-degree of freedom non-linear systems using an extended modal space model
Platten, Michael F; Wright, Jan Robert; Dimitriadis, Grigorioset al.
2009 • In Mechanical Systems and Signal Processing, 23 (1), p. 8-29
The identification of non-linear dynamic systems is an increasingly important area of research, with potential application in many industries. Current non-linear identification methodologies are, in general, mostly suited to small systems with few degrees of freedom and few non-linearities. In order to develop a practical identification approach for real engineering structures, the capability of such methods must be significantly extended. In this paper, it is shown that such an extension can be achieved using multi-exciter techniques in order to excited specific modes or degrees of freedom of the system under investigation.
A novel identification method for large non-linear systems is presented, based on the use of
a multi-exciter arrangement using appropriated excitation applied in bursts. This proposed Non-linear Resonant Decay Method is applied to a simulated system with 5 degrees of freedom and an experimental clamped panel structure. The technique is essentially a derivative of the Restoring Force Surface method and involves a non-linear curve fit performed in modal space. The effectiveness of the resulting reduced order model in representing the non-linear characteristics of the system is demonstrated. The potential of the approach for the identification of large continuous non-linear systems is also discussed.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A comparison of blade tip timing data analysis methods
Carrington, Ian B; Wright, Jan Robert; Cooper, Jonathan Eet al.
2001 • In Proceedings of the Institution of Mechanical Engineers. Part G, Journal of Aerospace Engineering, 215 (6), p. 301-312
The experimental determination of the vibration characteristics of rotating engine blades is very important for fatigue failure considerations. One of the most promising techniques for measuring the frequency of blade vibrations is blade tip timing. In this paper, three vibration analysis methods were specifically formulated and applied to the tip timing problem for the first time, using data obtained from a simple mathematical blade tip timing simulation. The results from the methods were compared statistically in order to determine which of the techniques is more suitable. One of the methods, the global autoregressive instrumental variables approach, produced satisfactory results at realistic noise levels. However, all of the techniques produced biased results under certain circumstances.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Identification of a Nonlinear Wing Structure Using an Extended Modal Model
Platten, M. F.; Wright, J. R.; Cooper, J. E.et al.
2009 • In Journal of Aircraft, 46 (5), p. 1614-1626
The nonlinear resonant decay method identifies a nonlinear dynamic system using a model based in linear modal space comprising the underlying linear system and a small number of additional terms that represent the nonlinear behavior. In this work, the method is applied to an aircraftlike wing/store/pylon experimental structure that consists of a rectangular wing with two stores suspended beneath it by means of nonlinear pylons with a nominally hardening characteristic in the store rotation degree of freedom. The nonlinear resonant decay method is applied to the system using multishaker excitation. The resulting identified mathematical model features five modes, two of which are strongly nonlinear, one is mildly nonlinear, and two are completely linear. The restoring force surfaces obtained from the mathematical model are in close agreement with those measured from the system. This experimental application of the nonlinear resonant decay method indicates that the method could be suitable for the identification of nonlinear models of aircraft in ground vibration testing.
Peer Reviewed verified by ORBi
Article (Scientific journals)
On the mechanism of divergent cable movement: Contribution of the bistable flow activity in the critical flow regime
2022 • In Journal of Wind Engineering and Industrial Aerodynamics, 228, p. 105103
This paper studies the contribution of the bistable flow activity on the installation of incipient conditions to provoke divergent stayed-cable movement. The experiment consists of performing dynamic tests on a long (6.1 m) and genuine inclined and smooth-surfaced cable in a wind tunnel. The cable has an outer diameter of 0.219 m and is equipped with 128 pressure taps distributed on seven (7) rings. The range of the Reynolds numbers is covering the critical flow regime around a circular cylinder. The evolution of mean and instantaneous aerodynamic coefficients is discussed and Proper Orthogonal Modes of the pressure pattern are created. The results indicate that both transitory regimes in the critical flow regime could contribute to install cable instability incipient conditions. Seemingly, only the unsteady states in the critical flow regime appear to participate on imposing large cable displacement. The influence of the circularity defect on disturbing flow fields around the circular cylinder has also been confirmed. Furthermore, the bistable flow activity has been clearly identified to contribute strongly in the installation of large-amplitude wind-induced vibrations.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Aeroservoelastic Simulations for Horizontal Axis Wind Turbines
2017 • In Proceedings of the Institution of Mechanical Engineers. Part A, Journal of Power and Energy, 231 (2), p. 103-117
This paper describes the development of a complete methodology for the aeroservoelastic modeling of horizontal axis wind turbines at the conceptual design stage. The methodology is based on the implementation of unsteady aerodynamic modeling, advanced description of the control system and nonlinear finite element calculations in the SWT wind turbine design package. The aerodynamic modeling is carried out by means of fast techniques, such as the Blade Element Method and the unsteady Vortex Lattice Method, including a free wake model. The complete model also includes a description of a doubly fed induction generator and its control system for variable speed operation. The SWT software features a non-linear finite element solver with multi-body dynamics capability. The full methodology is used to perform complete aeroservoelastic simulations of a realistic 2MW wind turbine model. The interaction between the three components of the approach is carefully analyzed and presented here.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Non-Linear Identification in Modal Space Using a Genetic Algorithm Approach for Model Selection
Platten, Michael F; Wright, Jan Robert; Worden, Keithet al.
2007 • In International Journal of Applied Mathematics and Mechanics, 3 (1), p. 72-89
The Non-Linear Resonant Decay Method is an approach for the identification of non-linear
systems with large numbers of degrees of freedom. The identified non-linear model is
expressed in linear modal space and comprises the modal model of the underlying linear
system with additional terms representing the non-linear behaviour. Potentially, a large
number of these non-linear terms will exist but not all of them will be significant. The
problem of deciding which and how many terms are required for an accurate identification has
previously been addressed using the Forward Selection and Backward Elimination techniques.
In this paper, a Genetic Algorithm optimisation is proposed as an alternative to those methods.
A simulated lumped parameter non-linear dynamic system is used to demonstrate the
proposed optimisation. The use of separate data sets for the identification and validation of the
modal model is also investigated. It is found that the Genetic Algorithm approach yields
significantly better results than the Backward Elimination and Forward Selection algorithms
in many cases.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Transonic corrections for the unsteady compressible Source and Doublet Panel Method
2024 • In Journal of Aircraft, Articles in advance
The steady and unsteady aerodynamic predictions of the Doublet Lattice Method (DLM) are routinely corrected for camber, twist, dihedral and transonic effects using steady reference data obtained from experiments or higher fidelity calculations. In this work, a novel correction technique is presented for the unsteady compressible Source and Doublet Panel Method (SDPM) in order to correct for transonic effects. Although the SDPM is more computationally expensive than the DLM, it can predict complete steady and oscillatory pressure distributions around the entire surface of the exact geometry of a wing or body. Consequently, it does not require corrections for geometric properties such as camber, twist and dihedral, only for transonic flow phenomena. The corrected predictions are compared to unsteady pressure distributions measured experimentally for the BSCW and LANN wings. It is shown that that the predictions of the SDPM correction method are accurate as long as the reference pressure data are of adequate quality. Finally, the SDPM with transonic correction is applied to the prediction of the flutter boundary of the AGARD 445.6 wing and it is demonstrated that the proposed approach is capable of predicting the transonic dip.
Peer Reviewed verified by ORBi
Article (Scientific journals)
A discrete adjoint full potential formulation for fast aerostructural optimization in preliminary aircraft design
Crovato, Adrien; Prado, Alex Prado; Cabral, Pedro Higinoet al.
2023 • In Aerospace Science and Technology, 138
Preliminary aircraft design is often carried out with the help of multidisciplinary optimization processes. Because there is a strong coupling between the flow and the structure of the aircraft, and because new composite materials have a unique capability of being designed for anisotropic directional stiffness and strength, these processes must model the aeroelastic behavior of the aircraft in order to be effective. Since many design variables are involved during the preliminary design stage, the optimization problems, formulated using high or medium fidelity models, can be solved using the adjoint method. Moreover, the level of fidelity of the fluid model and the associated simulation technique must also be selected with care, as they tend to be the main contributors to the computational cost. The novel contribution of the present work is twofold. Firstly, the discretized gradients of the full potential flow equation, which is a medium-fidelity model, are derived analytically so that they can be used in adjoint optimization problems. Moreover, the full potential flow solution and the computation of the gradients are implemented in an open-source and readily available finite element code. Secondly, aerodynamic shape and aerostructural optimization calculations are carried out on example wings to demonstrate the effectiveness and the computational efficiency of the proposed method. Overall, the results show that the newly implemented discrete adjoint nonlinear potential flow formulation is able to quickly optimize both the shape and the structural parameters of a typical wing. More specifically, the twist distribution along the wingspan is adapted to reduce the induced drag, the airfoils become more supercritical so that the shock strength and the associated wave drag are reduced, and the thickness of the structural elements is tailored to the loads to reduce the internal stresses. The present methodology is therefore able to deliver results at a low computational cost which are sufficiently accurate for the early design stages. Furthermore, the results obtained using the proposed methodology could be used as a starting point for the optimization calculations performed in later design stages.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Integrating Experimental and Computational Fluid Dynamics approaches using Proper Orthogonal Decomposition Techniques
The concept of Proper Orthogonal Decomposition (POD) is used to integrate Experimental Fluid Dynamics (EFD) and Computational Fluid Dynamics (CFD) approaches. The key idea is to take advantage of the optimality of the POD technique and its capability to extract the most energetic patterns of complex aerodynamic flow fields.
First, the concept of Modal Assurance Criterion (MAC) is used to obtain a simple quantitative criterion to compare EFD measurements to CFD results. The comparison is based on the POD modes, extracted from each set of data. The analysis of the energy content of the modes allows to draw important conclusions about the role of the latter. The method is applied in the study of the flow field around a rectangular cylinder, which is either static or oscillating in a low-speed flow field.
The second EFD/CFD integration technique deals with the reconstruction of a flow field from measured data, making use of CFD simulation results. The POD modes are first extracted from several CFD data sets, using a snapshot POD approach. Then the entire flow field of measured data can be reconstructed using a gappy POD method. The technique is applied to the transonic flow around a civil aircraft type wind tunnel model. The EFD measurements consist in pressure coefficient data from pressure ports or pressure-sensitive paint. It is shown that the complete flow field can be reconstructed from the pressure data with satisfactory accuracy and at relatively low computational cost.
The work demonstrates the potential of the POD technique to integrate EFD and CFD data, leading to a combined, validated and complete analysis of the flow under consideration.
Peer Reviewed verified by ORBi
Article (Scientific journals)
CUPyDO - An integrated Python environment for coupled fluid-structure simulations
Thomas, David; Cerquaglia, Marco Lucio; Boman, Romainet al.
2019 • In Advances in Engineering Software
CUPyDO, a fluid-structure interaction (FSI) tool that couples existing independent fluid and solid solvers into a single synchronization and communication framework based on the Python language is presented. Each coupled solver has to be wrapped in a Python layer in order to embed their functionalities (usually written in a compiled language) into a Python object, that is called and used by the coupler. Thus a staggered strong coupling can be achieved for time-dependent FSI problems such as aeroelastic flutter, vortex-induced vibrations (VIV) or conjugate heat transfer (CHT). The synchronization between the solvers is performed with the predictive block-Gauss-Seidel algorithm with dynamic under-relaxation. The tool is capable of treating non-matching meshes between the fluid and structure domains and is optimized to work in parallel using Message Passing Interface (MPI). The implementation of CUPyDO is described and its capabilities are demonstrated on typical validation cases. The open-source code SU2 is used to solve the fluid equations while the solid equations are solved either by a simple rigid body integrator or by in-house linear/nonlinear Finite Element codes (GetDP/Metafor). First, the modularity of the coupling as well as its ease of use is highlighted and then the accuracy of the results is demonstrated.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Aeroelastic system identification using transonic CFD data for a wing/store configuration
Vio, Gareth Arthur; Dimitriadis, Grigorios; Cooper, Jonathan E.et al.
2007 • In Aerospace Science and Technology, 11 (2-3), p. 146-154
This paper is part of a study investigating the prediction of the aeroelastic behaviour of aircraft subjected to non-linear aerodynamic forces. The main objective of the work is the characterization of the dynamic response of aeroelastic models resulting from coupled Computational Fluid Dynamic and Finite Element calculations. Of interest here is the identification of the flight condition at which the response bifurcates to limited or divergent amplitude self-sustained oscillations without carrying out a comprehensive set of full, computationally expensive, time-marching calculations. The model treated in this work is a three-dimensional wing in a transonic flowfield. Short datasets of pre-bifurcation behaviour are analysed to determine the system’s stability and degree of non-linearity. It is found that the calculated responses on the run-up to a transonic Limit Cycle Oscillation show little or no evidence of non-linearity. The non-linearity appears abruptly at the bifurcation flight condition. The variation of the local Mach number over the wing’s surface in the steady-state case is used to demonstrate that the non-linearity is due to a shock wave that can move along the surface. At Mach numbers where this is not possible the system behaves in a linear manner and its stability can be analysed using linear methods.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Quantitative Identification of Functional Connectivity Disturbances in Neuropsychiatric Lupus Based on Resting-State fMRI: A Robust Machine Learning Approach.
Neuropsychiatric systemic lupus erythematosus (NPSLE) is an autoimmune entity comprised of heterogenous syndromes affecting both the peripheral and central nervous system. Research on the pathophysiological substrate of NPSLE manifestations, including functional neuroimaging studies, is extremely limited. The present study examined person-specific patterns of whole-brain functional connectivity in NPSLE patients (n = 44) and age-matched healthy control participants (n = 39). Static functional connectivity graphs were calculated comprised of connection strengths between 90 brain regions. These connections were subsequently filtered through rigorous surrogate analysis, a technique borrowed from physics, novel to neuroimaging. Next, global as well as nodal network metrics were estimated for each individual functional brain network and were input to a robust machine learning algorithm consisting of a random forest feature selection and nested cross-validation strategy. The proposed pipeline is data-driven in its entirety, and several tests were performed in order to ensure model robustness. The best-fitting model utilizing nodal graph metrics for 11 brain regions was associated with 73.5% accuracy (74.5% sensitivity and 73% specificity) in discriminating NPSLE from healthy individuals with adequate statistical power. Closer inspection of graph metric values suggested an increased role within the functional brain network in NSPLE (indicated by higher nodal degree, local efficiency, betweenness centrality, or eigenvalue efficiency) as compared to healthy controls for seven brain regions and a reduced role for four areas. These findings corroborate earlier work regarding hemodynamic disturbances in these brain regions in NPSLE. The validity of the results is further supported by significant associations of certain selected graph metrics with accumulated organ damage incurred by lupus, with visuomotor performance and mental flexibility scores obtained independently from NPSLE patients.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Effect of Levels of Fidelity on Steady Aerodynamic and Static Aeroelastic Computations
Crovato, Adrien; Almeida, Hugo; Vio, Garethet al.
2020 • In Aerospace, II
Static aeroelastic deformations are nowadays considered as early as in the preliminary aircraft design stage, where low-fidelity linear aerodynamic modeling is favored because of its low computational cost. However, transonic flows are essentially nonlinear. The present work aims at assessing the impact of the aerodynamic level of fidelity used in preliminary aircraft design. Several fluid models ranging from the linear potential to the Navier–Stokes formulations were used to solve transonic flows for steady rigid aerodynamic and static aeroelastic computations on two benchmark wings: the Onera M6 and a generic airliner wing. The lift and moment loading distributions, as well as the bending and twisting deformations predicted by the different models, were examined, along with the computational costs of the various solutions. The results illustrate that a nonlinear method is required to reliably perform steady aerodynamic computations on rigid wings. For such computations, the best tradeoff between accuracy and computational cost is achieved by the full potential formulation. On the other hand, static aeroelastic computations are usually performed on optimized wings for which transonic effects are weak. In such cases, linear potential methods were found to yield sufficiently reliable results. If the linear method of choice is the doublet lattice approach, it must be corrected using a nonlinear solution.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Intracardiac or transesophageal echocardiography for left atrial appendage occlusion: an updated systematic review and meta-analysis.
2025 • In International Journal of Cardiovascular Imaging, 41 (3), p. 489 - 505
[en] BACKGROUND: Intracardiac echocardiography (ICE) appears to be a potential alternative for percutaneous left atrial appendage occlusion (LAAO) to transesophageal echocardiography (TEE). Thus, a meta-analysis was performed comparing ICE vs. TEE for LAAO guidance.
METHODS: A comprehensive literature search was performed using MEDLINE, Scopus and Web of Science electronic databases from their inception to November 2023.
RESULTS: 18 studies (124,230 patients) were included. Technical success was higher in ICE- compared to TEE-guidance (OR: 1.36, 95% CI 1.14 to 1.63, p = 0.006) and fewer devices employed (SMD: -0.22, 95% CI -0.43 to -0.01, p = 0.04, I2 = 62%). ICE guidance related with more pericardial effusion/tamponade and iatrogenic residual shunts (logRR: 0.62, 95% CI 0.36 to 0.89, p < 0.001 and RR: 1.53, 95% CI 1.12 to 2.09, p = 0.02, I2 = 1%, respectively). More vascular complications were noted in ICE group (logRR: 0.45, 95% CI 0.11 to 0.78, p = 0.009).
CONCLUSION: ICE-guided imaging is an effective alternative to TEE in LAAO, as it shows better efficacy than TEE, considering technical success. However, the higher rates of adverse events should be carefully considered.
Peer Reviewed verified by ORBi
Article (Scientific journals)
Non-linear aeroelastic prediction for aircraft applications
Henshaw, M.J. de C.; Badcock, Ken J.; Vio, G. A.et al.
2007 • In Progress in Aerospace Sciences, 434 (4-6), p. 65-137
Current industrial practice for the prediction and analysis of flutter relies heavily on linear methods and this has led to overly conservative design and envelope restrictions for aircraft. Although the methods have served the industry well, it is clear that for a number of reasons the inclusion of non-linearity in the mathematical and computational aeroelastic prediction tools is highly desirable. The increase in available and affordable computational resources, together with major advances in algorithms, mean that non-linear aeroelastic tools are now viable within the aircraft design and qualification environment. The Partnership for Unsteady Methods in Aerodynamics (PUMA) Defence and Aerospace Research Partnership (DARP) was sponsored in 2002 to conduct research into non-linear aeroelastic prediction methods and an academic, industry, and government consortium collaborated to address the following objectives:
(1) To develop useable methodologies to model and predict non-linear aeroelastic behaviour of complete aircraft.
(2) To evaluate the methodologies on real aircraft problems.
(3) To investigate the effect of non-linearities on aeroelastic behaviour and to determine which have the greatest effect on the flutter qualification process.
These aims have been very effectively met during the course of the programme and the research outputs include:
(a) New methods available to industry for use in the flutter prediction process, together with the appropriate coaching of industry engineers.
(b) Interesting results in both linear and non-linear aeroelastics, with comprehensive comparison of methods and approaches for challenging problems.
(c) Additional embryonic techniques that, with further research, will further improve aeroelastics capability.
This paper describes the methods that have been developed and how they are deployable within the industrial environment. We present a thorough review of the PUMA aeroelastics programme together with a comprehensive review of the relevant research in this domain. This is set within the context of a generic industrial process and the requirements of UK and US aeroelastic qualification. A range of test cases, from simple small DOF cases to full aircraft, have been used to evaluate and validate the non-linear methods developed and to make comparison with the linear methods in everyday use. These have focused mainly on aerodynamic non-linearity, although some results for structural non-linearity are also presented. The challenges associated with time domain (coupled computational fluid dynamics–computational structural model (CFD–CSM)) methods have been addressed through the development of grid movement, fluid–structure coupling, and control surface movement technologies. Conclusions regarding the accuracy and computational cost of these are presented. The computational cost of time-domain methods, despite substantial improvements in efficiency, remains high. However, significant advances have been made in reduced order methods, that allow non-linear behaviour to be modelled, but at a cost comparable with that of the regular linear methods. Of particular note is a method based on Hopf bifurcation that has reached an appropriate maturity for deployment on real aircraft configurations, though only limited results are presented herein. Results are also presented for dynamically linearised CFD approaches that hold out the possibility of non-linear results at a fraction of the cost of time coupled CFD–CSM methods. Local linearisation approaches (higher order harmonic balance and continuation method) are also presented; these have the advantage that no prior assumption of the nature of the aeroelastic instability is required, but currently these methods are limited to low DOF problems and it is thought that these will not reach a level of maturity appropriate to real aircraft problems for some years to come. Nevertheless, guidance on the most likely approaches has been derived and this forms the basis for ongoing research. It is important to recognise that the aeroelastic design and qualification requires a variety of methods applicable at different stages of the process. The methods reported herein are mapped to the process, so that their applicability and complementarity may be understood. Overall, the programme has provided a suite of methods that allow realistic consideration of non-linearity in the aeroelastic design and qualification of aircraft. Deployment of these methods is underway in the industrial environment, but full realisation of the benefit of these approaches will require appropriate engagement with the standards community so that safety standards may take proper account of the inclusion of non-linearity.