On the closure of Collar’s triangle by optical diagnostics

An experimental methodology is proposed to study aeroelastic systems with optical diagnostics. The approach locally evaluates the three physical mechanisms that produce the forces involved in Collar’s triangle, namely aerodynamic, elastic, and inertial forces. Flow and object surface tracers are tracked by a volumetric particle image velocimetry (PIV) system based on four high-speed cameras and LED illumination. The images are analysed with Lagrangian particle tracking techniques, and the flow tracers and surface markers are separated based on the different properties of their images. The inertial and elastic forces are obtained solely analysing the motion and the deformation of the solid object, whereas the aerodynamic force distribution is obtained with pressure from PIV techniques. Experiments are conducted on a benchmark problem of fluid–structure interaction, featuring a flexible panel installed at the trailing edge of a cylinder. Data are collected in the resonant regime, where the panel exhibits a two-dimensional motion. The estimation of inertial and elastic forces is obtained enforcing a high-order polynomial fit to the surface motion and deformation. The aerodynamic loads on the panel are challenged by the need to devise adaptive boundary conditions complying with the panel motion. The closure of Collar’s triangle yields overall residuals of about one-half of the inertial force taken as reference. The simultaneous measurement of the three forces paves the way to assessing the equilibrium of forces closing the Collar’s triangle. The latter can be intended for uncertainty evaluation or, when only two forces are measured, for estimation of the remaining Collar element.


Introduction
Aeroelastic phenomena occur frequently in nature as well as in a wide number of engineering applications, most notably in civil (Sarkar et al. 1994), energy (Abdelkefi 2016) and power engineering (Marshall et al., 1996), and in the transport sector (Wright and Cooper 2008). In the aeroelastic regime, a flexible body immersed in a flow is subjected to the aeroelastic system of forces, namely aerodynamic, elastic, and inertial. Their mutual relationship was conceptualised in the so-called Collar's triangle (Collar 1947).
Still, nowadays, fluid-structure interaction (FSI) problems are difficult to study both with experimental as well as computational techniques. Computer simulations of the fluid flow coupled with a time-varying structural boundary remain a challenging task. Computational aeroelasticity reduces nowadays practical problems into simpler geometries (Schuster et al. 2003) for affordability. The numerical coupling of independent fluid and structural solvers, while improving the computational efficiency, poses new challenges at the fluid-structure interface, especially in the presence of structural or aerodynamic non-linearities, such as occurring in regimes with unsteady flow separation. This situation justifies the need to advance experimental approaches as a necessary tool for the validation of FSI simulations.
The first challenge encountered with FSI experiments is given by the problem of scaling. While Mach and Reynolds numbers are the most relevant parameters for aerodynamic scaling, FSI problems require consideration of the ratio of fluid and structure inertia alongside the structure relative stiffness (Friedmann 2004;Wan and Cesnik 2014).
Even when the problem is scaled to a satisfactory extent and the experiment is realised, monitoring the fluid flow and structure behaviour is hindered by the technical complexity of the needed instrumentation. Most of the works carrying out simultaneous measurements of flow and structure have been reported only in the last decade, and by combining optical techniques. Many studies have approached the flow measurements by planar particle image velocimetry; only a few have employed volumetric measurements (Kalmbach and Breuer 2013;Mitrotta et al. 2019). The measurement of the structure motions was performed either by digital image correlation (Bleischwitz et al. 2017;Marimon Giovannetti et al. 2017) or based on multiple-point laser triangulation (Kalmbach and Breuer 2013). Moreover, Timpe et al. (2013) directly measured the inertia forces on the structure with a force balance. Furthermore, Mitrotta et al. (2019) used the Lagrangian particle tracking technique  to determine simultaneously the flow velocity and motion of the structure.
The combined use of PIV techniques (Adrian et al. 2011) and optical diagnostics like DIC or alike (DIC, Chu et al. 1985;ground vibration tests, Peeters et al. 2009; Image Pattern Correlation Technique or IPCT, Boden et al. 2014) to perform both tasks is a topic currently under development. Pressure from 3D-PIV enables the evaluation of the fluid flow pressure both in the fluid domain as well as in close proximity of the solid object surface (Liu and Katz 2006;van Oudheusden 2013). Pioneering work that applies pressure from PIV to FSI problems with generic three-dimensional and flexible structures has been conducted by Percin et al. (2017) for the case of a flexible membrane under mechanical actuation.
Finally, the evaluation of inertia forces, traditionally performed by means of accelerometers (as in-ground vibration tests, Peeters et al. 2009;or under in-flight conditions, Kehoe 1995), has also been surrogated by optical diagnostic methods. For example, the videogrammetric model deformation (VMD, Burner and Liu 2001) measurement technique captures the attitude of wings by tracking rows of markers on the model. Similarly, a new sensing approach (de Figueiredo et al. 2020) measures natural frequencies of aeronautical structures using a computer vision system. From the above, it can be stated that the use of optical techniques, such as PIV, offers the potential to provide a holistic analysis for FSI problems, compared to probe-based measurements.
Once it has been shown that aeroelastic forces can be measured by combining different optical techniques to study FSI phenomena, an element of novelty can be introduced: the principle of force equilibrium applied to Collar's triangle of aeroelasticity, formulated as the Collar criterion. The feasibility of invoking the Collar criterion to relate two or more elements of the Collar triangle is scrutinised in the present study. The potential benefit of such approach is to obtain a global evaluation of the experimental uncertainties when all forces are measured. Alternatively, when only two forces are accessed, the third one may be inferred using such equation.
The present work describes the foundational aspects of the experimental hardware and procedures needed to diagnose the aeroelastic forces non-intrusively by use of PIV techniques. The measurements are intended for the closure of Collar's triangle by a single experiment, that is, simultaneously determining all dynamic forces acting within the FSI problem. The demonstration makes use of a benchmark problem for aeroelasticity, i.e. the oscillation of a flexible thin plate trailing a circular cylinder immersed in a fluid stream.

FSI measurement working principle
Let us consider a flexible, thin panel of length L, width W and thickness T ( T << L, W ) in oscillatory motion induced by aerodynamic forces arising from the air flow as illustrated in Fig. 1. The aerodynamic loads consist of pressure p and friction . The object opposes deformation by its stiffness, and acceleration motions by its inertia. Considering a square tile of such panel, i.e. a small element of side length l (see Fig. 1), the equilibrium of forces acting on such element reads as: where F i , F e and F a stand for the inertia, elastic, and aerodynamic forces, respectively. For the current work, the aforementioned vector identity will be reduced to the dominant direction of motion, z or transverse component, under the assumption of small deformations. As a consequence, the vector Eq. 1 is projected along z and rewritten as a scalar relation. Note that z will be used as coordinate along the thickness of the element, with z = 0 corresponding to the middle of the plate in the undeformed state. In addition, for a simplified treatment of the problem, the quantities will be hypothesized uniform along the spanwise direction y; thus, the forces will be related solely to the streamwise component, x. The accuracy of the latter hypothesis is verified for the experiment presented in the remainder of the work.

Inertia force
A tile subject to non-zero elastic and aerodynamic forces will accelerate in proportion to such forces. Based on Newton third principle, the reaction force F i due to inertia acts in the direction that opposes its acceleration a. For a tile element, the expression of the inertia force reads as: (1) In the above equation, m is the density of the material, and the displacement of the plate mid-point along z is denoted by h, as illustrated in Fig. 1. The measurement of the object displacement and motion is typically performed by means of optical markers applied onto the surface of the panel. The markers' positions in time are recorded, while their velocities and accelerations are computed from time derivatives of the positions. The inertial force is then obtained by knowledge of the tile's specific mass and volume.

Elastic force
Elastic forces arise as a reaction of the object opposing deformation. The present discussion pertains to the linear elastic regime. For small displacements, the force can be considered to act normal to the surface and tensional forces (e.g. arising from aerodynamic skin friction) are neglected. In the above conditions, the whole elastic force reduces to the bending stress only.
Considering a thin deformable body, a 1-D simplification of the Kirchhoff-Love plate theory provides a relation between the deformation of the tile and the elastic forces, given its physical attributes. According to such theory, and invoking the strain-displacement and the stress-strain relations, the stress along the body reads: where E is the Young modulus, and ν is the Poisson ratio. The linear dependence on z indicates tensile/compressive stresses at the opposite faces of the panel. The non-uniform normal stress distribution gives rise to a moment reaction. Finally, considering the tile moment equilibrium, the elastic reaction force, which results normal to the surface, can be expressed as, The complete derivation of the Kirchhoff-Love plate theory for 3D solids can be found in the book of Panc (1975). Once the material properties are known, the experimental evaluation of the elastic force relies on the assessment of the 4th spatial derivative of the instantaneous deformed shape h(x). The latter is regarded as a challenging task for an experiment, where even a small extent of measurement noise propagates and amplifies through the multiple operations of spatial differentiation.

Aerodynamics force
Aerodynamics loads are usually decomposed into normal (pressure) and shear (friction) stresses. Friction is negligible in most aeroelastic regimes (Fung 2008). Moreover, the pressure loads on a thin panel reduce the discussion to the pressure difference Δp from its opposite sides. The aerodynamic force acting on the tile reads as: The pressure field can be inferred from the velocity (and acceleration) distribution invoking the momentum equation. The latter in turn requires the numerical solution of the Poisson equation for pressure (van Oudheusden 2013). The instantaneous aerodynamic force is not straightforward to obtain from fluid flow measurements as it typically requires high-quality time-resolved 3D measurements (Beresh 2021). Moreover, the choice of boundary conditions can largely affect the evaluation of the pressure field (Neeteson et al. 2016). To minimize errors due to low spatial resolution or unwanted light reflections close to solid surfaces, the surface pressure is typically estimated imposing Dirichlet boundary conditions for the pressure values in the external flow region, where the Bernoulli relation between pressure and velocity holds true. In the rotational flow region (shear layers and separated wake), the integration of the pressure can follow the 3D propagation scheme recently proposed by Jux et al. (2020). It should be stressed that the evaluation of the unsteady pressure distribution over a generic surface, curved (5) F a = Δplw and in motion, time-resolved three-dimensional flow field measurements are required that surround such surface.

Collar's triangle of forces
Equation 1 can be rewritten combining the explicit terms obtained in Eqs. 2, 4 and 5. The resulting equilibrium of forces is expressed as: From such equation, the structural parameters (plate thickness T) and material properties (density m , Young modulus E, and Poisson ratio ν) need to be known. The corresponding values applying to the current experiment are gathered in Table 1. The remaining elements in the equation, such as the pressure difference Δp , the local acceleration and the spatial derivatives of the deformed shape need to be measured to close the problem.
It is worth mentioning that although the above identity theoretically holds, the experimental determination of such forces has inherent uncertainties. The latter is mostly associated with the measurement techniques adopted to monitor the flow and object motion. Therefore, one can already predict a residual force when computing the force equilibrium for Collar's triangle, referred to as Collar closure criterion (CC) force CC : In the present study, an optical diagnostic system based on PIV techniques is developed for the purpose of determining the above forces. The fluid flow velocity is obtained over an extended three-dimensional domain by tracking large tracers (sub-millimetre helium-filled soap bubbles, HFSB, Scarano et al. 2015) along their trajectories to determine the flow pressure. The acceleration of the structural element is measured by tracking the motion of surface markers with the same algorithm as for the fluid tracers. Finally,  For the measurement of the foil kinematics, a regular pattern of 210 bright surface markers was applied to the foil. The markers had an average diameter d m = 0.8mm and a pitch p m = 1cm . The overall experimental setup is illustrated and provided in Fig. 2.

Flow measurement
The flow was seeded with neutrally buoyant sub-millimetre (mean diameter d ≈ 350μm ) helium-filled soap bubbles (HFSB, Scarano et al. 2015), released by an in-house 200-generator seeding rake installed in the tunnel settling chamber (similarly to Jux et al. 2020). The nominal production rate of the flow tracers is 6 × 10 6 bubbles/s, resulting in a seeding concentration of ~ 10 bubbles/cm 3 for a free stream velocity U ∞ = 2.25m∕s at the test section (seeded streamtube of 60 ×25cm 2 ). For the current experiment, only 100 generators were operated to reduce the seeding concentration to ~ 5 bubbles/cm 3 , so as to satisfy the condition on the maximum number of particles per pixel (ppp = 0.05) for 4-camera tomographic systems (Elsinga et al. 2006;Scarano, 2013).
A measurement volume of 22 × 40 × 20cm 3 was illuminated vertically (spanwise direction) by two LaVision Flashlight 300 LED banks located under the test section. In the current experiments, flow tracers were measured at both sides of the transparent foil. A mirror was installed at the top plane of the test section to increase the light intensity in the measurement volume, thus allowing to identify the tracer particles also in the region behind the foil. The imaging system consisted of four high-speed cameras ( 1024 × 1024 pixels, 12 bit, 20 μm pixel pitch) mounted in a tomographic configuration 90 cm away from the foil, as shown in Fig. 2. The cameras subtended an aperture of approximately 30 degrees. Details of the optical system used in the experiments are summarized in Table 2.
The image acquisition and processing were conducted with the LaVision Davis 10.0.5 software. Sets of 5000 images were recorded in continuous mode at a rate of 1 kHz. Samples of the latter can be visualized in Fig. 3. For the calibration process, a geometrical calibration (Soloff et al.   Schanz et al. 2012). The surface markers were visually distinguished from the flow tracers based on their distinct optical characteristics. In order to track most fluid tracers (spherical bubbles) and most structural markers (2D flat bright circles), the volume self-calibration and the OTF calibration were performed separately with clean runs containing only one tracer type. Both the flow tracers' velocities and the foil markers motions were evaluated by means of the Shake-the-Box Lagrangian Particle Tracking (LPT) algorithm . Based on the work of Battacharya and Vlachos (2020), particle location uncertainty is expected below 0.3 pixels for IPR-based reconstructions and further decreases to 0.1 pixels when the particles are track over multiple time instants with the STB algorithm (Sciacchitano et al. 2021).

Data reduction for load estimation
The following section describes the data reduction procedure to obtain the distribution of forces along the foil from the analysis of the tracers and markers motion.

Structural forces
As explained in Sect. 2, inertial and elastic forces are determined from the time and space derivatives, respectively, of the measured instantaneous foil's shape. In order to avoid noise propagation through the derivation process, a spatio-temporal least-square regression is performed on the locations of the markers to generate a coherent smooth surface. The regression makes use of a high-order polynomial in space (5th for x-direction and 2nd for y-direction) and a 2nd order in time, which is applied to all the measured markers over time stencil of 17 samples (~ 10% of one oscillation cycle). The resulting foil's shape is calculated on a Cartesian grid with 0.5 cm spacing in both streamwise and spanwise (x and y, respectively) directions. Therefore, the foil can be considered as discretised by adjacent square tiles, each centred at a grid node and having a surface area of 0.25 cm 2 (0.5 cm × 0.5 cm). The derivations in time and space to determine the acceleration and elastic force, respectively, are conducted analytically based on the polynomial description of the foil's instantaneous shape. The uncertainty of the elastic and inertia forces is quantified based on the uncertainty of the regression coefficients of the least-square fit (Montgomery et al. 2021) and linear error propagation. The comparison between the raw markers' measurements and the foil's position and acceleration based on the regression analysis is illustrated in Figs. 4 and 5. Figure 4 compares the kinematics in time of a tracked surface marker with that of its projection onto the reconstructed surface. From this figure, the periodicity of the foil's motion both in terms of displacement and acceleration is evident, with the oscillation period equal to T = 175 ms. The use of the least-square regression allows to reconstruct missing time samples throughout the entire sequence ( Fig. 4-left) and to clearly suppress the noise in the foil's acceleration (Fig. 4right). Figure 5 displays the instantaneous location of the markers (in black circles) and the fitted surface at one of the  The computation of the pressure field relies upon the material derivative of the velocity Du∕Dt , which is readily available from the LPT measurements, and the viscous dissipation ∇ 2 u . Due to the high Reynolds ( Re = 2 × 10 4 ), the latter term is neglected for the pressure computation. As shown by van Oudheusden et al. (2007) for a flow problem at similar Reynolds number, even in regions of high flow curvature, the viscous contribution to the pressure gradient is two orders of magnitude smaller than the other contributions. Dense velocity reconstruction from particle tracks is performed with the VIC# algorithm implemented in DaVis 10.1 (Jeon et al. 2018). The approach is based on the Vortex-in-Cell method for VIC + (Schneiders and Scarano 2016) and optimises the vorticity field that best fits the measured velocity and acceleration by minimizing a cost function J built upon the disparity between the computed and the measured variables. Pressure integration is performed simultaneously to the velocity reconstruction by solving the Poisson equation with Neumann conditions at every boundary and a constant Dirichlet condition at a free stream location.

Fluid load
Because the applicability of the VIC# algorithm is limited to Cartesian domains not containing solid objects, the pressure reconstruction is only performed over two separate regions of the measurement domain that span from its outmost boundaries to the limits of the foil motion (approximately 4 cm away from the neutral position of the foil). Hence, the foil does not enter the fluid domain under consideration for the VIC# reconstruction (black rectangles in Fig. 6). For the determination of the actual surface pressure on the foil, the pressure is then extrapolated towards the foil location (from both sides) via a nearest neighbour approach, neglecting the change of static pressure across the shear layer and the wake (Schlichting and Gersten 2016). In order to determine an upper bound of the uncertainty of the reconstructed pressure, the aforementioned approach is performed from the outer regions towards the mid location ( z = 0 ) at downstream locations of the trailing edge of the foil. The directions of extrapolation and the mid location are shown with black arrows and a red cross in Fig. 6, respectively. The difference in static pressure at the latter location when computing it from opposite directions, p = 1.7 Pa, gives an estimation of the uncertainty of pressure jump across the foil.

Results
The flow velocity fields are analysed first, before discussing the loads of Collar's triangle and the resulting structural motion. Despite the 3D flow measurement, given the oscillations feature a quasi-2D behaviour and for sake of conciseness, only 2D flow data will be presented and discussed. The latter assumption is verified on both structural and flow data. The foil deformation presents a 2.5% deviation over different locations along the span of the foil. For the flow, the streamwise velocity component is taken as reference. The standard deviation along the span reaches at most 10% of the local average velocity, which is ascribed to the inherent turbulence of the flow. Hence, spanwise averaging is considered a valid approach in this problem to enhance the data density without any loss of flow or structural information.
A view in the median x-z plane illustrates velocity vectors and contours of the streamwise component (Fig. 7-top  row). The values shown are averaged along the y-direction, where also the foil location is indicated. The velocity field reveals the oscillatory behaviour of the wake as expected in the regime of Kármán vortex shedding. The foil position also reveals its oscillatory motion. However, the foil does not appear to be in phase with the wake axis oscillations. The position and phase of the vortices shed in the cylinder wake are made more evident by displaying the transverse velocity component, where a negative blob followed by a positive one indicates a counter-clockwise rotating Kármán vortex. The dominant frequency of the flow could be already inferred from the sequence, and it corresponds to that of the flapping foil, f ≈ 5.7Hz . The latter is verified from the time evolution of the transverse flow velocity, displayed in blue in Fig. 8. Note that the computed Δp is also included in Fig. 8 and exhibits a different frequency, in particular double to that of the transverse velocity. This can be explained from the fact that low-pressure regions arise at the vortex core locations regardless of the rotation direction of the vortex, thus resulting in double the pressure dominant frequency with respect to the flow. Nevertheless, the measured flow dominant frequency agrees well with the characteristic frequency for cylinder shedding wakes in the sub-critical regime, St ≈ 0.2 (Williamson 1996), apparently, unaffected by the presence of the flexible foil. Figure 7 (bottom row) illustrates a sequence of the pressure field time evolution. Recall that, due to the limitations of the available pressure solvers, an accurate pressure field could be calculated in the outer regions, i.e. away from the foil. Therefore, for the current visualization, the effect of the foil (imparting a small pressure discontinuity at its opposite faces) is neglected, and the pressure equation is solved in the entire measurement domain. Despite the expected errors at the vicinity of the solid interface, the pressure fields provide insights to the structure and phase of the aerodynamic forcing mechanism. The presence of low-pressure blobs, corresponding to the core of Kármán vortices, is observed, alternating from the opposite sides of the foil and travelling at the same convection speed of the vortices. During the passage of a low-pressure blob, the foil is accelerated in the direction of the vortex and continues after the vortex has passed the foil tip under the effect of inertia.
A time evolution of the computed surface pressure along the streamwise direction x is provided in Fig. 9. Note that as the foil deflects, the projection to the x-axis reduces the range to 0 < x∕L<0.9 depending on the motion phase.
Regarding the forcing, the most effective appears to occur when the low-pressure region transits from the foil midlength ( x∕L = 0.6 ) to its end ( x∕L = 0.9 ). In this situation of maximum forcing, both the pressure load (a pulling force) and the foil velocity vector are aligned, thus resulting in a positive work performed by the flow to the foil. This mechanism is enhanced by the high-pressure regions (red orange in Fig. 7), which build up in the outer region when the foil is at the limit of the motion and aid to push the foil back towards the neutral position. The above effects ultimately lead to an oscillatory aeroelastic regime with pronounced periodic behaviour.
The time evolution of the foil's motion is illustrated in Fig. 10. Contours of position, tiles' acceleration, and elastic force per tile are provided on the deformed shape of the foil. Position (left column) and acceleration of the foil (centre column) are in phase-opposition, as expected from the second derivative of an oscillatory motion. The inertia force (not shown here) is opposite to the acceleration; hence, it exhibits the same trend as the foil's displacement. Similarly, to the acceleration, the elastic force opposes to the foil position, as the deformation promotes internal stresses and moments that act against it to recover the undeformed state. This working mechanism is analogous to that of a massspring system in oscillation. Figure 11 illustrates the foil motion by the time evolution of three markers located at the mid-span of the foil and at 1/4, 1/2 and 3/4, respectively, of the foil's chord length. The in-phase displacement, for all positions, indicates that the foil motion is dominated by its first mode of oscillation, and a stationary wave at the frequency f foil = 5.7Hz is observed. This frequency is very close to the first natural frequency (regardless damping) of the studied configuration, f 1,natural = 5.17Hz . This mode of oscillations is consistent with the theory of oscillating cantilevered beams (Den Hartog 1985), where the displacement increases more than linearly from root to tip.

Discussion on the closure of Collar's triangle
The time evolutions of all measured aeroelastic forces acting on tiles of 0.25 cm 2 area (0.5 cm × 0.5 cm) are presented in Fig. 12. The value of such forces is extracted at the three tile locations x∕L = {0.25, 0.5, 0.75}(y∕W = 0.5 ). An observation can be made of the relative phase of the forces with respect to the motion of the foil. The tile's inertia force features a standing wave in-phase with the displacement h, which is consistent with the operator of second derivative between position and acceleration for an oscillatory signal. Therefore, the tile acceleration signal is in phase-opposition, (i.e. π phase-shift) to the displacement. Based on Eq. 2, the inertia (red line in Fig. 12) has opposite sign to the acceleration, leading to a condition of zero phase shift (in-phase signals) between displacement and inertia force. As expected from the measurement technique, the inertia force has the lowest uncertainty of the triangle, inertia = 6 × 10 −6 N. The elastic force (in yellow, Fig. 12) should exhibit a value that is in phase opposition to the displacement and the inertia force. This behaviour is found back at x∕L = 0.75 ( Fig. 12-right). However, the variations in time do not feature a sinusoidal behaviour and the signal appears to be corrupted by spurious fluctuations at a frequency higher than that of the primary oscillation. At points closer to the root of the foil ( x∕L = 0.25 and x∕L = 0.75 ), it appears that these spurious fluctuations compromise the evaluation of the elastic force as the noise level attains the same level of the expected magnitude of the elastic force. Overall, based on error propagation through the spatial regression onto the fourth derivative of the deformation, the elastic force presents the highest uncertainty elastic = 2 × 10 −4 N.
Finally, the aerodynamic force (shown in blue in Fig. 12) displays the lowest force levels among the three components and for the three locations of the foil. Differently to inertia and elastic forces, the aerodynamic force appears as a travelling wave with a different delay to the foil motion at each chord location. The latter can be also observed from the surface pressure in Fig. 9. This effect is explained by the streamwise convection of Kármán vortices featuring lowpressure blobs at their core. The uncertainty for the aerodynamic force resulted in aero = 4.3 × 10 −5 N.
The condition of maximum forcing for the aerodynamic force occurs when the aerodynamic force and the foil velocity are in phase. In such condition, the alignment between the aerodynamic force and the foil velocity vectors transfers the maximum work from the flow to the foil. Therefore, the aerodynamic force must anticipate the foil's motion (displacement) by a phase shift of about ∕2 (phase quadrature). From the present experiments (Fig. 12), it appears that such condition is best attained around the position x∕L = 0.75.
Once the foil reaches a steady-state oscillatory regime (of approximately constant amplitude and frequency), only a small aerodynamic force is required to maintain the system in such regime, where the sum of elastic and inertia forces conserve most of the system kinetic energy. Instead, the work produced by the aerodynamic force is dissipated by structural damping. The evaluation of structural losses, however, goes beyond the scope of the current study and is not examined further.
The closure of Collar's triangle concludes this work and is achieved by imposing the equilibrium of the aeroelastic forces. The equilibrium addresses the dominant effects along the z-component of the force vector. When uncertainties pertaining to the experimental approach and the subsequent modelling to estimate the forces are accounted for by a term CC , Eq. 1 is rewritten as: This residual was anticipated in Sect. 2 as the Collar closure criterion (CC), CC , and quantifies how accurate the measurement of three aeroelastic forces is altogether. For the present experiments, the closure of Collar's triangle (i.e. the magnitude of the residual force, CC ) is evaluated for the afore-analysed tile locations. In view of the oscillatory nature of the forces, the standard deviation (STD) of the error is reported in Table 3 taking the inertia force as a reference for the amplitude.
While the inertia force increases downstream, CC increases towards the root. This is ascribed to the uncertainty of the elastic force, which appears to be maximum at the foil's root, with the small displacements. The ratio between the Collar force and the inertia force as provided in Table 3 can be seen as an indicator of the accuracy of the closure of Collar's triangle of forces. It is evident that from the foil's root to approximately the half chord (x/L = 0.5), the residual of Collar's force strongly distorts the equilibrium. As previously mentioned, this is ascribed to the large uncertainty of (8) F i + F e + F a = CC the elastic force associated with the fourth spatial derivative of displacements, and the vanishing magnitude of the inertia force close to the foil's root. Instead, towards the foil's tip (x/L ≥ 0.75), the residual tends to become a vanishing component, indicating the feasibility of Collar's triangle closure by non-intrusive experimental means.
Aside from the uncertainty of the elastic force the hypothesis of small displacements will require further scrutiny. When the latter assumption is not satisfied, aerodynamic elastic and inertia forces will need to be evaluated with their full vectors and the local orientation of the foil needs to be accounted for.

Conclusions
A non-intrusive measurement methodology for aeroelastic experiments based on the Lagrangian particle tracking technique is proposed to simultaneously estimate aerodynamic, elastic and inertial loads from a deforming model. The method relies on the tracking of flow tracers to characterize the flow, and of surface markers to describe the dynamics of the structure. An experiment was carried out to show the applicability of the proposed methodology. The measurements regarded a flexible transparent foil attached to a solid cylinder and interacting with the wake of the latter. Both the flow tracers (neutrally-buoyant HFSB) and surface markers applied onto the panel were tracked in time with a single optical tomographic system.
The results show the potential of evaluating the aeroelastic forces from Lagrangian particle tracking measurements of the flow tracers and the markers on the flexible structure. The inertia forces were successfully evaluated via timederivation of the markers' displacements. The pressure on the foil surface was obtained extrapolating the solution of the Poisson problem from the outer flow field to the foil. The elastic forces required evaluation of the fourth spatial derivative of the displacement and turned out to be the most challenging quantity to determine.
Flow analysis results show that pressure blobs associated with the Karman vortices determine the unsteady aerodynamic forcing. A phase shift exists between such external loading and the inertial and elastic response.
The Collar closure criterion is introduced that considers the residual of the force equilibrium equation as a measure of the accuracy of the overall experiment. The criterion yields unfavourable results close to the foil root where displacement and acceleration are vanishing. At the foil trailing edge, a consistent combination of the aeroelastic forces is observed, indicating the feasibility of the present optical diagnostic approach for aeroelastic systems characterisation. Table 3 Collar closure criterion force and inertia force for a tile location at y∕W = 0.5 x/L = 0.25 x/L = 0.5 x/L = 0.75 STD F i (mN) 4.4 × 10 -2 9.4 × 10 -2 13.2 × 10 -2 STD CC (mN) 11.2 × 10 -2 8.7 × 10 -2 6.2 × 10 -2 STD CC /STD F i 2.5 0.9 0.47