An Experimental and Computational Study on Inverted Flag Dynamics for Simultaneous Wind – Solar Energy Harvesting

This paper presents results from experiments and simplified numerical simulations on the flow-induced dynamics and power generation of inverted flags that combine flexible piezoelectric strips with photovoltaic cells to simultaneously harvest kinetic wind energy and solar radiant energy. Experiments were conducted in a wind tunnel under controlled wind excitation and light exposure, focusing in particular on the dynamics and power generation of the inverted flag harvester. Numerical simulations were carried out using a lattice-Boltzmann fluid solver coupled with a finite element structural solver via the immersed-boundary method, focusing in particular on minimizing the simulation run time. The power generated during the tests shows that the proposed inverted flag harvester is a promising concept, capable of producing enough power (on the order of 1 mW) to supply low-power electronic devices in a range of applications where distributed power generation is needed. Notwithstanding key simplifications implemented in the numerical model to achieve a fast execution, simulations and measurements are in good agreement, confirming that the lattice-Boltzmann method is a viable and time-effective alternative to classic Navier–Stokes-based solvers when dealing with strongly coupled fluid–structure interaction problems characterized by large structural displacements.


Introduction
Energy harvesting is the process by which ambient background energy (e.g., kinetic energy from wind/water flow or motion, solar radiation, structural vibration, temperature or salinity gradients, electromagnetic energy) is recovered to supply small-scale low-power (microwatt to milliwatt power range) electronic devices such as sensors, data-loggers, and data-transmitters for use in distributed sensing, equipment/process monitoring, smart city, and Internet of things applications [1][2][3][4][5][6].Piezoelectric materials, in particular, were quite extensively investigated for harvesting energy from structural and flow-induced movement or vibration [7][8][9][10][11][12][13][14][15][16][17][18][19].Piezoelectric harvesters typically comprise a structural element realized with a piezoelectric material (or bonded with a piezoelectric patch) that is either connected to a vibrating or moving structure, or exposed to fluid flow in such a way that flow-induced vibration or motion takes place.The structural or flow-induced vibration/motion periodically bends the piezoelectric material, which in turn converts mechanical

Piezo-Solar Inverted Flag Harvester
The piezo-solar inverted flag was realized with four PVDF elements (TE Connectivity Ltd, Schaffhausen, Switzerland (www.te.com); model LDT2-028K/L; length/width/thickness: 73/16/0.205mm; density: 1780 kg/m 3 ; Young modulus: 3 GPa; Poisson ratio: 0.34 [33]) and two flexible solar panels (Powerfilm Solar Inc., Ames, IA, USA (www.powerfilmsolar.com);model SP3-37; length/width/thickness: 36.5/64/0.2mm; mass: 0.7 g), one attached at each side of the flag as shown in Figure 1d.As can be noted in Figure 1d, the flexible solar panels are located close to the upstream leading edge of the inverted flag, such that their inclusion does not constrain the deformation of the flexible piezoelectric PVDF strips close to the fixed trailing edge of the flag where most of the mechanical strain is concentrated.It is evident that, in the present inverted flag design, the structural support is provided by the PVDF elements, while the solar panels constitute added masses.As is well known, tip masses can have huge influence on the dynamics of cantilevers, such that the inclusion of the solar panels onto the inverted flag can be expected to significantly affect the flag dynamics.In order to study this effect systematically, we realized two more inverted flags: one flag without solar panels and another with smaller solar panels, as shown in Figure 1b,c.As schematically indicated in Figure 1a, the configuration of the inverted flag can be characterized with three geometrical variables: (1) the overhang length Lp of the flag from the fixed support to the free tip; (2) the flag aspect ratio AR = Lp /Wp (ratio of flag length to flag width); and (3) the flag coverage ratio CR = Ls/Lp (ratio of solar panel length to flag length).The coverage ratio values for the three flags in Figure 1 were 0 (Figure 1b, referred to as no-solar flag), 0.32 (Figure 1c, referred to as half-solar flag), and 0.64 (Figure 1d, referred to as full-solar flag), while the overhang length Lp (57 mm) and aspect ratio AR (0.9) were not varied.

Experimental Set-Up and Methodology
The experimental set-up used to characterize the dynamics and power generation of the inverted flags is schematically shown in Figure 2. The downstream edge of the flags was connected with a purpose-built fixation clamp (see Figure 1b-d) to a vertical pole, which was placed in the middle of a wind tunnel (Armfield Limited, Ringwood, UK (armfieldonline.com))with octagonal cross-section and height/width of 350 mm.The fixation clamp was manufactured from laser cut Perspex (Hobarts Laser Supplies, Leybourne, UK (https://hobarts.com);model: AC.CLR0000.03.3020) and was As can be noted in Figure 1d, the flexible solar panels are located close to the upstream leading edge of the inverted flag, such that their inclusion does not constrain the deformation of the flexible piezoelectric PVDF strips close to the fixed trailing edge of the flag where most of the mechanical strain is concentrated.It is evident that, in the present inverted flag design, the structural support is provided by the PVDF elements, while the solar panels constitute added masses.As is well known, tip masses can have huge influence on the dynamics of cantilevers, such that the inclusion of the solar panels onto the inverted flag can be expected to significantly affect the flag dynamics.In order to study this effect systematically, we realized two more inverted flags: one flag without solar panels and another with smaller solar panels, as shown in Figure 1b,c.As schematically indicated in Figure 1a, the configuration of the inverted flag can be characterized with three geometrical variables: (1) the overhang length L p of the flag from the fixed support to the free tip; (2) the flag aspect ratio AR = L p /W p (ratio of flag length to flag width); and (3) the flag coverage ratio CR = L s /L p (ratio of solar panel length to flag length).The coverage ratio values for the three flags in Figure 1 were 0 (Figure 1b, referred to as no-solar flag), 0.32 (Figure 1c, referred to as half-solar flag), and 0.64 (Figure 1d, referred to as full-solar flag), while the overhang length L p (57 mm) and aspect ratio AR (0.9) were not varied.

Experimental Set-Up and Methodology
The experimental set-up used to characterize the dynamics and power generation of the inverted flags is schematically shown in Figure 2. The downstream edge of the flags was connected with a purpose-built fixation clamp (see Figure 1b-d) to a vertical pole, which was placed in the middle of a wind tunnel (Armfield Limited, Ringwood, UK (armfieldonline.com))with octagonal cross-section and height/width of 350 mm.The fixation clamp was manufactured from laser cut Perspex (Hobarts Laser Supplies, Leybourne, UK (https://hobarts.com);model: AC.CLR0000.03.3020) and was designed to hold the flag in a sandwich arrangement, realizing a cantilever boundary condition at the fixed edge of the flags.Construction details for the fixation clamp are provided in Appendix A.
designed to hold the flag in a sandwich arrangement, realizing a cantilever boundary condition at the fixed edge of the flags.Construction details for the fixation clamp are provided in Appendix A.
Two light-emitting diode (LED) lights (Vibesta B.V., Helmond, The Netherland (www.vibesta.com);model: Capra 12 Daylight) were located on both sides of the wind tunnel and oriented perpendicularly to the rest position of the flag.The light intensity the flags were exposed to while at rest was measured with a portable light meter at 1.8 klx, and was not changed during the tests.Daylight intensity varies from about 10 2 lx at sunset/sundown to about 10 5 lx in bright sunlight, with mid-range light intensities of 1-2 klx corresponding to a typical overcast day.The present results for solar power generation, therefore, can be regarded as indicative of the solar energy harvesting potential of the proposed flag harvester during an overcast day.Prior to the tests, the flow in the wind tunnel was characterized with a calibrated single-probe hot-wire anemometer (Dantec Dynamics, Bristol, UK (www.dantecdynamics.com);probe type 55P15; 5-μm-diameter tungsten wire of 2 mm length operated in constant temperature mode).The boundary layer extended less than 5 mm from the walls of the wind tunnel at the lowest wind speed setting employed (1.5 m/s), and the velocity profile across the wind tunnel (excluding the boundary layer) was uniform to within 1%, i.e., velocity variations were on the order of the measuring error of the hot-wire probe.This assures that the inverted flags were always exposed to a fully developed velocity profile during the tests.
As shown in Figure 3, the streamwise turbulence intensity measured in the wind tunnel flow was around 0.5-0.8%within the range of free stream flow velocities from 4 m/s to 10 m/s of interest here, corresponding to mild turbulence.During the tests with the inverted flags, the wind speed in the wind tunnel was deduced (error within ±5%) from static/dynamic pressure measurements, using a calibrated pressure transducer (Sensirion AG, Staefa, Switzerland (www.sensirion.com);model: SDP816) connected to the wind tunnel with four static pressure ports located upstream of the inverted flag.Ambient temperature and pressure during the tests were 293 ± 1 K and 101 ± 1 kPa.Two light-emitting diode (LED) lights (Vibesta B.V., Helmond, The Netherland (www.vibesta.com);model: Capra 12 Daylight) were located on both sides of the wind tunnel and oriented perpendicularly to the rest position of the flag.The light intensity the flags were exposed to while at rest was measured with a portable light meter at 1.8 klx, and was not changed during the tests.Daylight intensity varies from about 10 2 lx at sunset/sundown to about 10 5 lx in bright sunlight, with mid-range light intensities of 1-2 klx corresponding to a typical overcast day.The present results for solar power generation, therefore, can be regarded as indicative of the solar energy harvesting potential of the proposed flag harvester during an overcast day.
Prior to the tests, the flow in the wind tunnel was characterized with a calibrated single-probe hot-wire anemometer (Dantec Dynamics, Bristol, UK (www.dantecdynamics.com);probe type 55P15; 5-µm-diameter tungsten wire of 2 mm length operated in constant temperature mode).The boundary layer extended less than 5 mm from the walls of the wind tunnel at the lowest wind speed setting employed (1.5 m/s), and the velocity profile across the wind tunnel (excluding the boundary layer) was uniform to within 1%, i.e., velocity variations were on the order of the measuring error of the hot-wire probe.This assures that the inverted flags were always exposed to a fully developed velocity profile during the tests.
As shown in Figure 3, the streamwise turbulence intensity measured in the wind tunnel flow was around 0.5-0.8%within the range of free stream flow velocities from 4 m/s to 10 m/s of interest here, corresponding to mild turbulence.During the tests with the inverted flags, the wind speed in the wind tunnel was deduced (error within ±5%) from static/dynamic pressure measurements, using a calibrated pressure transducer (Sensirion AG, Staefa, Switzerland (www.sensirion.com);model: SDP816) connected to the wind tunnel with four static pressure ports located upstream of the inverted flag.Ambient temperature and pressure during the tests were 293 ± 1 K and 101 ± 1 kPa.The motion of the flag was recorded with a digital camera (Panasonic Co., Osaka, Japan (www.panasonic.com);model: LUMIX DMC-FZ200; recording frequency: 200 fps; resolution: 480 × 640 pixels) located on top of the wind tunnel.In order to facilitate the optical tracking, a thin red paper tape was added to the upper border of the flags, as can be noted in Figure 1.The videos were post-processed with the Image Processing Toolbox of MATLAB (version R2015a), employing a tracking methodology previously used for flow-induced vibration and flexible fluid−structure interaction studies [34][35][36].The motion of the flag was characterized with the instantaneous angular position () (measured to within ±2-3°), which is the angle between the tip of the flag at time  and the flag rest position; and with the angular span of motion, which is the angle described between the two extreme positions spanned by the flag when observed from the top (see insert in Figure 4a).Locally, the flag possesses a variable curvature when deflected, such that the instantaneous angular position and the angular span of motion used here give a simple representative measure of the displacement for the whole flag.The power spectral density (PSD) was obtained by applying a fast Fourier transform (FFT) to the time series of the instantaneous angular position (using MATLAB built-in functions), and the frequency of flag flapping was identified as the dominant frequency (when present) in the power spectral density (uncertainty within 0.2-0.3Hz, estimated from the full width at half maximum of the dominant peak in the power spectrum).
In addition, we reconstructed the trajectory of the system in phase-space (i.e., the system attractor) using the delayed vector method proposed by Bradley and Kantz [37]: a methodology developed for nonlinear time-series analysis that is particularly useful for fluid−structure interaction problems [34][35][36].The topology of the reconstructed attractor of the inverted flag, in fact, gives a useful qualitative characterization of the flag dynamics, which completes and corroborates the quantitative information provided by the angular position and flapping frequency.The concept of phase-space representation, rather than a classic analysis in time or frequency domain, is the key point in nonlinear time-series analysis.The topology of the trajectory in phase-space of a nonlinear system, in fact, provides important qualitative information regarding the fundamental dynamics of the systems being investigated.The problem is that, in experimental studies, one normally observes a time series of scalar measurements of some quantity that depends on the current state of the system, and not the trajectory of the system in phase-space.The delayed vectors method [37] allows reconstructing the trajectory of the system in phase-space from a time series of scalar measurements, and can, therefore, be used in experimental studies such as the present one where only time series of scalar measurements of some quantity (the instantaneous angular position in the present case) are available.In the present case, the trajectory in phase-space was reconstructed by plotting the timedelayed instantaneous angular position ( + ) versus ().If the delay  is chosen properly (on the order of 40-60 ms in the present case), then the topology of the reconstructed attractor is representative of the underlying system dynamics [37].
The power generated by the inverted flag harvesters was collected and processed through LabVIEW 2017, while the data were gathered with a National Instruments (Austin, TX, USA (www.ni.com))NI-USB-6225 external DAQ device.Following common practice in energy-harvesting studies, optimal load resistances that allowed peak power generation were empirically identified and employed.In order to collect the data from the individual solar panels and PVDF elements, separate circuits were built using prototyping breadboards that allowed circuit changes to be made quickly.The motion of the flag was recorded with a digital camera (Panasonic Co., Osaka, Japan (www.panasonic.com);model: LUMIX DMC-FZ200; recording frequency: 200 fps; resolution: 480 × 640 pixels) located on top of the wind tunnel.In order to facilitate the optical tracking, a thin red paper tape was added to the upper border of the flags, as can be noted in Figure 1.The videos were post-processed with the Image Processing Toolbox of MATLAB (version R2015a), employing a tracking methodology previously used for flow-induced vibration and flexible fluid-structure interaction studies [34][35][36].The motion of the flag was characterized with the instantaneous angular position θ(t) (measured to within ±2-3 • ), which is the angle between the tip of the flag at time t and the flag rest position; and with the angular span of motion, which is the angle described between the two extreme positions spanned by the flag when observed from the top (see insert in Figure 4a).Locally, the flag possesses a variable curvature when deflected, such that the instantaneous angular position and the angular span of motion used here give a simple representative measure of the displacement for the whole flag.The power spectral density (PSD) was obtained by applying a fast Fourier transform (FFT) to the time series of the instantaneous angular position (using MATLAB built-in functions), and the frequency of flag flapping was identified as the dominant frequency (when present) in the power spectral density (uncertainty within 0.2-0.3Hz, estimated from the full width at half maximum of the dominant peak in the power spectrum).
approximation.Being located downstream of the flag, the support clamp that holds the flag in the wind tunnel does not significantly affect the flow development around the flag.Therefore, the support clamp was not modeled in the simulations, although the fixed edge of the flag was subjected to a local cantilever-type boundary condition that properly accounted for the support clamp.
The numerical simulations were set up to match the experimental conditions.However, to reduce the simulation run-time, the fluid viscosity was increased by three orders of magnitude, thus decreasing the Reynolds number proportionally.This follows the rationale successfully proposed for simulating flexible filaments [49] and inverted flags [50], and is here extended to composite inverted flags realized with flexible PVDF elements and flexible solar panels.In fluid−structure interaction problems characterized by large displacements, such as flexible filaments and the inverted flags of interest here, pressure forces control the interaction between the flow and the structure, whilst viscous effects play a comparatively minor role.This explains why numerical models simplified by artificially increasing the fluid viscosity can still capture the dynamics with reasonable accuracy [49,50], while achieving a much faster execution.Although not critical for the single-flag problem analyzed here, a fast execution will be instrumental in future studies to analyze the interaction among several flags and tailor the inverted flag harvesters to the intended final application.It is worth noting that more complex collision operators [51] could potentially allow stable solutions at higher Reynolds numbers.However, in light of the above discussion, this was deemed unnecessary for the present work.The regularized boundary condition [52] was used to impose the velocity at the inlet and far field boundaries, and a convective condition was imposed at the outlet [53].Preliminary testing showed that a relaxation time scale of  = 0.57 was sufficient to ensure a stable solution while, at the same time, limiting the lattice velocity to reduce compressibility effects.At the highest tested flow speed (12 m/s), the lattice velocity was found to be  = 0.035, which results in a Mach number of  = 0.061 .This is well within the weakly compressible regime such that errors due to compressibility effects are limited.
All calculations were performed on a custom-built computer with a 3-GHz Intel Core i7-5960X processor (eight cores and two threads/core), and with 64 GB of random access memory (RAM).

Inverted Flag Dynamics
The angular span of motion and the frequency of oscillation of the inverted flags, measured while gradually and stepwise increasing the wind speed from 0 m/s (no flow) to a maximum of about 10 m/s, are displayed in Figure 4.In particular, the absence of a dominant frequency of motion is indicated in Figure 4b with a null value of the frequency.As can be noted in Figure 4, at low wind speed (below about 4-5 m/s), the angular span of motion and oscillating frequency are both zero, indicating that the inverted flags maintain their rest equilibrium configuration aligned with the incoming flow.As soon as the wind speed is large enough, the rest configuration becomes unstable and the flags start oscillating.During these In addition, we reconstructed the trajectory of the system in phase-space (i.e., the system attractor) using the delayed vector method proposed by Bradley and Kantz [37]: a methodology developed for nonlinear time-series analysis that is particularly useful for fluid-structure interaction problems [34][35][36].The topology of the reconstructed attractor of the inverted flag, in fact, gives a useful qualitative characterization of the flag dynamics, which completes and corroborates the quantitative information provided by the angular position and flapping frequency.The concept of phase-space representation, rather than a classic analysis in time or frequency domain, is the key point in nonlinear time-series analysis.The topology of the trajectory in phase-space of a nonlinear system, in fact, provides important qualitative information regarding the fundamental dynamics of the systems being investigated.The problem is that, in experimental studies, one normally observes a time series of scalar Fluids 2019, 4, 87 6 of 20 measurements of some quantity that depends on the current state of the system, and not the trajectory of the system in phase-space.The delayed vectors method [37] allows reconstructing the trajectory of the system in phase-space from a time series of scalar measurements, and can, therefore, be used in experimental studies such as the present one where only time series of scalar measurements of some quantity (the instantaneous angular position in the present case) are available.In the present case, the trajectory in phase-space was reconstructed by plotting the time-delayed instantaneous angular position θ(t + τ) versus θ(t).If the delay τ is chosen properly (on the order of 40-60 ms in the present case), then the topology of the reconstructed attractor is representative of the underlying system dynamics [37].
The power generated by the inverted flag harvesters was collected and processed through LabVIEW 2017, while the data were gathered with a National Instruments (Austin, TX, USA (www.ni.com))NI-USB-6225 external DAQ device.Following common practice in energy-harvesting studies, optimal load resistances that allowed peak power generation were empirically identified and employed.In order to collect the data from the individual solar panels and PVDF elements, separate circuits were built using prototyping breadboards that allowed circuit changes to be made quickly.The data acquisition program was written as a VI (Virtual Instrument) in LabVIEW 2017 using the standard DAQ-mx library.This program gathered, saved, and displayed the data in real time to allow an immediate impression of the power generated to be seen and to allow the load resistance to be set to give the highest power output.The sampling rate was set at 1 kHz to allow sufficient resolution of data through a flapping cycle for any particular flapping characteristics to be seen.

Numerical Model
Numerical simulations of the inverted flag dynamics were conducted by adopting a partitioned approach whereby the fluid dynamics and the structural dynamics are handled via their own separate field solvers, a rather common approach for fluid-structure interaction problems characterized by large structural displacements.The lattice-Boltzmann method [38], an alternative to traditional Navier-Stokes-based computational fluid dynamics solvers, was employed to solve the fluid dynamics.The governing equation for the lattice-Boltzmann method is the Boltzmann transport equation, given in discrete form by where x represents the spatial coordinates, t is time, c i is the ith component of the lattice velocity vector, and F i (x, t) is a term accounting for the presence of an external force f (x, t) (e.g., gravity) [39].The probability distribution function f i (x, t) describes the proportion of fluid molecules within an elemental volume located at x and time t moving with velocity c i .Equation ( 1) clearly shows a relaxation toward a local equilibrium.The equilibrium function f eq i (x, t) is obtained via a Taylor series [40] or Hermite polynomial [41] expansion of the Maxwell-Boltzmann distribution function and the relaxation time scale τ is directly related to the fluid viscosity [38].Finally, the macroscopic fluid quantities can be recovered by taking moments of the probability function as follows: The main advantage of the lattice-Boltzmann method is its amenability to parallel implementation, as well as the regular square lattice that aids with the solver coupling.Although not critical for the single-flag problem analyzed here, the parallel implementation capability will be instrumental in future studies to analyze the interaction among several flags.In fact, we envision that, in a real application, a harvester could combine several inverted flags properly tuned (via optimization of their geometry and/or mechanical properties) to yield the desired performance over the whole range of wind speeds and light conditions of interest for that particular application.The computational burden associated with analyzing and optimizing a multi-flag configuration will necessarily require a parallel implementation for time-efficient computation.
A co-rotational finite element method was adopted to solve the structural dynamics.Since large deformations are expected in the present work, a nonlinear Newton-Raphson solver was utilized.The governing equations for the system are given by [42] M Ün+1,k where n and k are the time step and iteration counters, M is the mass matrix, K is the tangent stiffness matrix, F ext is the external load vector (including fluid forces), F int represents the internal forces within the structure due to the stresses, Ü represents the nodal accelerations, and ∆U represents the incremental nodal displacements.The system is iterated over until the incremental nodal displacements are sufficiently small, after which the solution is advanced to the next time step [43].Second-order time integration is achieved via the implicit Newmark scheme.Since the inverted flag deformation is mainly two-dimensional (2D), 2D simulations were deemed appropriate.In particular, Euler-Bernoulli beam elements were used to model the flags, which is appropriate for structures with large length-to-thickness ratios such as the inverted flags of interest here.The immersed boundary method [44] was used to couple the separate fluid dynamics and structural dynamics field solvers.This approach allows the fluid and structure to be solved on their own separate grids, and facilitates the transfer of information between them, thus simplifying the meshing and helping in capturing large deformations of the structure.However, specialized interpolation/spreading operators are required to transfer information between the two grids.Using lowercase notation for quantities on the fluid grid (Eulerian frame) and uppercase for the boundary grid (Lagrangian frame), these operators are given by where x = (x, y, z), X = (q, r, s), ϕ is a quantity defined on the fluid grid, Φ is the same quantity defined on the structure grid, δ is the discrete form of the Dirac delta function, and is a scaling factor which ensures reciprocity between the interpolation and spreading steps [45].Note that, for the regular Cartesian grid used here, ∆x = ∆y = ∆z = 1 in lattice units.Furthermore, for the one-dimensional (1D) line boundary used in this work, ∆q = ∆r = 1.The spacing between immersed boundary points was chosen to match the lattice spacing such that ∆s = ∆x = 1.
The coupled lattice-Boltzmann immersed boundary algorithm decomposes the velocity into a predicted and a force-corrected term [46].Firstly, the lattice-Boltzmann step is computed without the immersed boundary to obtain a predicted velocity field v * .From Equation (3), the corrected fluid velocity is, thus, given by Using the interpolation operator to convert this to the Lagrangian frame gives Since the velocity of the boundary V(X, t) = I[v(x, t)] is known, Equation ( 9) can be rearranged to solve for the corrective force term F(X, t).The corrective force is converted back to the Eulerian frame using the spreading operator and finally the velocity field is updated by adding the corrective force to the predicted velocity field via Equation (8).With regard to the structural calculation, F(X, t) is passed to the finite element method (FEM) solver such that F ext = F(X, t) in the absence of any additional external forcing (e.g., gravity).
To ensure stability, a block Gauss-Seidel coupling scheme was adopted [47], whereby the field solvers are iterated over, within the current time step, until the interface conditions are satisfied.To accelerate the convergence of these iterations, a dynamic relaxation factor was used to relax the displacement of the structure after each iteration [47].The resulting strongly coupled fluid-structure interaction numerical model is capable of handling the large nonlinear interactions associated with the inverted flag configuration.For more details regarding the model, the reader is directed to previous work [48].
To promote the onset of any instability, the flag was initialized with a small transverse displacement (1% of its length), i.e., the simulated flag at rest was slightly deflected on one side.After the initial transient, the simulations were run for approximately 20 flapping periods.The variation in thickness across the flag, associated with the presence of the solar panel, was properly accounted for in the structural solver.However, the immersed boundary, which is the geometry that the fluid "sees", was represented as a line.This means that the variation in thickness across the flag due to the solar panel was not accounted for in the fluid solver.Although critical for the mechanical response of the flag, the variation in thickness across the flag was small enough to neglect its effect on the resulting flow field.For the scope of the present study, this was considered an acceptable approximation.Being located downstream of the flag, the support clamp that holds the flag in the wind tunnel does not significantly affect the flow development around the flag.Therefore, the support clamp was not modeled in the simulations, although the fixed edge of the flag was subjected to a local cantilever-type boundary condition that properly accounted for the support clamp.
The numerical simulations were set up to match the experimental conditions.However, to reduce the simulation run-time, the fluid viscosity was increased by three orders of magnitude, thus decreasing the Reynolds number proportionally.This follows the rationale successfully proposed for simulating flexible filaments [49] and inverted flags [50], and is here extended to composite inverted flags realized with flexible PVDF elements and flexible solar panels.In fluid-structure interaction problems characterized by large displacements, such as flexible filaments and the inverted flags of interest here, pressure forces control the interaction between the flow and the structure, whilst viscous effects play a comparatively minor role.This explains why numerical models simplified by artificially increasing the fluid viscosity can still capture the dynamics with reasonable accuracy [49,50], while achieving a much faster execution.Although not critical for the single-flag problem analyzed here, a fast execution will be instrumental in future studies to analyze the interaction among several flags and tailor the inverted flag harvesters to the intended final application.It is worth noting that more complex collision operators [51] could potentially allow stable solutions at higher Reynolds numbers.However, in light of the above discussion, this was deemed unnecessary for the present work.The regularized boundary condition [52] was used to impose the velocity at the inlet and far field boundaries, and a convective condition was imposed at the outlet [53].Preliminary testing showed that a relaxation time scale of τ = 0.57 was sufficient to ensure a stable solution while, at the same time, limiting the lattice velocity to reduce compressibility effects.At the highest tested flow speed (12 m/s), the lattice velocity was found to be v = 0.035, which results in a Mach number of Ma = 0.061.This is well within the weakly compressible regime such that errors due to compressibility effects are limited.
All calculations were performed on a custom-built computer with a 3-GHz Intel Core i7-5960X processor (eight cores and two threads/core), and with 64 GB of random access memory (RAM).

Inverted Flag Dynamics
The angular span of motion and the frequency of oscillation of the inverted flags, measured while gradually and stepwise increasing the wind speed from 0 m/s (no flow) to a maximum of about 10 m/s, are displayed in Figure 4.In particular, the absence of a dominant frequency of motion is indicated in Figure 4b with a null value of the frequency.
As can be noted in Figure 4, at low wind speed (below about 4-5 m/s), the angular span of motion and oscillating frequency are both zero, indicating that the inverted flags maintain their rest equilibrium configuration aligned with the incoming flow.As soon as the wind speed is large enough, the rest configuration becomes unstable and the flags start oscillating.During these oscillations, the angular span of motion gradually increases, as the wind speed is gradually increased, reaching a maximum of about 190 • at a wind speed of 9 m/s.The frequency of oscillation, on the other hand, gradually decreases as the wind speed is gradually increased, although its variation is rather limited.The large-amplitude flapping oscillation persists until the wind speed is large enough to deflect the flag onto one side; in Figure 4, this corresponds to a wind velocity of about 10 m/s where the angular span of motion and the oscillating frequency are both zero.
The dynamic response of the full-solar flag at selected wind velocities is further analyzed in Figures 5-7, which include (for the three representative wind speeds of 4.7 m/s, 9.0 m/s, and 9.7 m/s) (a) the instantaneous angular position time series (5-s sample); (b) the instantaneous angular position power spectral density (PSD); (c) the instantaneous angular position autocorrelation function; (d) the instantaneous angular position histogram (5-s sample; to help visualize the trend in the data, the histogram includes a red line corresponding to a normal density function fitted to the data); and (e) the reconstructed attractor in phase-space.As can be seen in Figure 5, at low wind speed (4.7 m/s), the flag is aligned with the incoming flow but appears to be not completely static.Rather, it undergoes small-amplitude oscillations (±1 • ) with no discernible dominant frequency, with a very fast-decaying autocorrelation function, with a bell-shaped histogram, and with a reconstructed attractor that covers the whole space.
At low wind speed, therefore, the flag dynamics appears to be a small-amplitude random vibration around the rest equilibrium configuration.The instantaneous angular position, however, is comparable in magnitude with the measuring error, such that the flag could also be actually static, and the measured motion could be the consequence of random noise pick-up from the instrumentation.The numerical simulations (described later) suggest that the flag is actually static, although, due to the limitations of the numerical model previously discussed, the possibility of a small random vibration cannot be ruled out.It is worth highlighting that the optical tracking system and the numerical model were conceived to capture the large flag displacements during flapping; thus, a relatively poor resolution of very small displacements is, to some extent, unavoidable and has no consequence on the conclusions of the present study.
At higher wind speed (9.0 m/s) in Figure 6, on the other hand, the instantaneous angular position time series is periodic and its power spectral density includes one dominant frequency with a few harmonics, the autocorrelation function is periodic and slowly decaying, the histogram is bridge-shaped, and the reconstructed attractor has a ring-like topology.This indicates that the flag dynamics at this wind speed involves a large-amplitude limit-cycle flapping oscillation.
Finally, at the highest wind speed tested (9.7 m/s) in Figure 7, the flag dynamics is qualitatively similar to what observed at low wind speed (4.7 m/s) in Figure 5; the flag apparently undergoes small-amplitude oscillations (±2-3 • ) with no discernible dominant frequency, the autocorrelation function decays fast, the histogram is bell-shaped, and the reconstructed attractor covers the whole space.Therefore, when the flag is fully deflected onto one side at very high wind speed, its dynamics appears to be a small-amplitude random vibration around the deflected configuration.Similarly to that previously discussed for the low wind speed of 4.7 m/s, the instantaneous angular position is comparable in magnitude with the measuring error; thus, the small-amplitude vibration cannot be duly resolved with the present set-up.As previously noted, this limitation has no consequence on the conclusions of the present study.
The dynamics analysis carried out for the no-solar and the half-solar flags yielded results qualitatively analogous to those obtained for the full-solar flag presented in Figures 5-7.For the sake of brevity, these are not documented here.In conclusion, all three main dynamic responses of inverted flags identified in previous research and previously discussed (i.e., static-aligned, flapping, and fully deflected) are, therefore, clearly observable in the present results, and this qualitatively validates the present approach.
reaching a maximum of about 190° at a wind speed of 9 m/s.The frequency of oscillation, on the other hand, gradually decreases as the wind speed is gradually increased, although its variation is rather limited.The large-amplitude flapping oscillation persists until the wind speed is large enough to deflect the flag onto one side; in Figure 4, this corresponds to a wind velocity of about 10 m/s where the angular span of motion and the oscillating frequency are both zero.
The dynamic response of the full-solar flag at selected wind velocities is further analyzed in Figures 5-7, which include (for the three representative wind speeds of 4.7 m/s, 9.0 m/s, and 9.7 m/s) (a) the instantaneous angular position time series (5-s sample); (b) the instantaneous angular position power spectral density (PSD); (c) the instantaneous angular position autocorrelation function; (d) the instantaneous angular position histogram (5-s sample; to help visualize the trend in the data, the histogram includes a red line corresponding to a normal density function fitted to the data); and (e) the reconstructed attractor in phase-space.As can be seen in Figure 5, at low wind speed (4.7 m/s), the flag is aligned with the incoming flow but appears to be not completely static.Rather, it undergoes small-amplitude oscillations (±1°) with no discernible dominant frequency, with a very fast-decaying autocorrelation function, with a bell-shaped histogram, and with a reconstructed attractor that covers the whole space.As evident in Figure 4, the inclusion of the solar panels has a destabilizing effect on the dynamics of the inverted flag.The added mass of the solar panels, in fact, reduces the wind speed required to trigger the large-amplitude limit-cycle oscillation, widens the wind speed window where limit-cycle oscillations are sustained, and increases the angular span of motion at any given wind speed.These observations are consistent with the expected destabilizing effect of the solar panels that constitute a tip added mass onto the inverted flag, which can be regarded as a cantilever.The frequency of oscillation decreases when adding the solar panels, consistently indicating that a heavier flag flaps at a Fluids 2019, 4, 87 11 of 20 lower frequency with respect to a lighter one.The transition from the large-amplitude flapping to the small-amplitude vibration around the fully deflected configuration, however, is not significantly affected by the inclusion of the solar panels.It is evident from the results presented that the dynamics of the inverted flags remains qualitatively the same when the solar panels are included, because all three main dynamics modes of inverted flags (static-aligned, flapping, and fully deflected) are still observable.Importantly, the inclusion of the solar panels widens the wind speed window where limit-cycle oscillations are sustained, which is the mode of operation of interest for energy harvesting with PVDF elements.
vibration around the rest equilibrium configuration.The instantaneous angular position, however, is comparable in magnitude with the measuring error, such that the flag could also be actually static, and the measured motion could be the consequence of random noise pick-up from the instrumentation.The numerical simulations (described later) suggest that the flag is actually static, although, due to the limitations of the numerical model previously discussed, the possibility of a small random vibration cannot be ruled out.It is worth highlighting that the optical tracking system and the numerical model were conceived to capture the large flag displacements during flapping; thus, a relatively poor resolution of very small displacements is, to some extent, unavoidable and has no consequence on the conclusions of the present study.
At higher wind speed (9.0 m/s) in Figure 6, on the other hand, the instantaneous angular position time series is periodic and its power spectral density includes one dominant frequency with a few harmonics, the autocorrelation function is periodic and slowly decaying, the histogram is bridgeshaped, and the reconstructed attractor has a ring-like topology.This indicates that the flag dynamics at this wind speed involves a large-amplitude limit-cycle flapping oscillation.Finally, at the highest wind speed tested (9.7 m/s) in Figure 7, the flag dynamics is qualitatively similar to what observed at low wind speed (4.7 m/s) in Figure 5; the flag apparently undergoes small- Noteworthy, the inverted flags exhibited a rather pronounced hysteresis in their dynamics, as illustrated for the full-solar flag in Figure 8, where the angular span of motion is plotted for increasing and then decreasing wind speed.As can be noted, once the wind speed gradually increased to a value large enough to fully deflect the flag, then the wind velocity was reduced considerably before the large-amplitude limit-cycle flapping resumed.This hysteresis can be traced back to the mechanical properties of the PVDF elements that make up the flags.In particular, the high flexibility of the PVDF elements comes at the expense of a limited elasticity, which prevents the flags from resuming the flapping motion once they reach fully deflected configuration at high wind speed, unless the wind speed is reduced substantially.
amplitude oscillations (±2-3°) with no discernible dominant frequency, the autocorrelation function decays fast, the histogram is bell-shaped, and the reconstructed attractor covers the whole space.Therefore, when the flag is fully deflected onto one side at very high wind speed, its dynamics appears to be a small-amplitude random vibration around the deflected configuration.Similarly to that previously discussed for the low wind speed of 4.7 m/s, the instantaneous angular position is comparable in magnitude with the measuring error; thus, the small-amplitude vibration cannot be duly resolved with the present set-up.As previously noted, this limitation has no consequence on the conclusions of the present study.The dynamics analysis carried out for the no-solar and the half-solar flags yielded results qualitatively analogous to those obtained for the full-solar flag presented in Figures 5-7.For the sake of brevity, these are not documented here.In conclusion, all three main dynamic responses of inverted flags identified in previous research and previously discussed (i.e., static-aligned, flapping, and fully deflected) are, therefore, clearly observable in the present results, and this qualitatively validates the present approach.
As evident in Figure 4, the inclusion of the solar panels has a destabilizing effect on the dynamics of the inverted flag.The added mass of the solar panels, in fact, reduces the wind speed required to The hysteresis in the inverted flag response, clearly not ideal for final applications, will be alleviated in future optimization studies using custom-made PVDF elements with mechanical properties tailored to the intended application, instead of the off-the-shelf PVDF elements used here.A different approach to alleviate the inverted flag hysteresis was explored by Silva-Leon et al. [54], who included a metal shim in their inverted flags.In their design, the structural support to the flag was provided by the metal shim, whilst the PVDF elements and the solar panels were added masses attached onto the metal shim.The higher elasticity of the metal shim helped in reducing the hysteresis, although it also made their flags more rigid and, therefore, more stable, such that the wind speed required to sustain large-amplitude flapping increased substantially.large to fully deflect the flag, then the wind velocity was reduced considerably before the large-amplitude limit-cycle flapping resumed.This hysteresis can be traced back to the mechanical properties of the PVDF elements that make up the flags.In particular, the high flexibility of the PVDF elements comes at the expense of a limited elasticity, which prevents the flags from resuming the flapping motion once they reach the fully deflected configuration at high wind speed, unless the wind speed is reduced substantially.The hysteresis in the inverted flag response, clearly not ideal for final applications, will be alleviated in future optimization studies using custom-made PVDF elements with mechanical properties tailored to the intended application, instead of the off-the-shelf PVDF elements used here.A different approach to alleviate the inverted flag hysteresis was explored by Silva-Leon et al. [54], who included a metal shim in their inverted flags.In their design, the structural support to the flag was provided by the metal shim, whilst the PVDF elements and the solar panels were added masses attached onto the metal shim.The higher elasticity of the metal shim helped in reducing the hysteresis, although it also made their flags more rigid and, therefore, more stable, such that the wind speed required to sustain large-amplitude flapping increased substantially.

Numerical Simulations
The comparison between measurements and numerical predictions is provided in Figure 9a,b, where the predicted angular span of motion and frequency of oscillation are compared with the corresponding measurements for the no-solar and the full-solar flags.Overall, the agreement between measurements and predictions is encouraging.This indicates that, despite the simplifications previously discussed, the present numerical modeling approach is capable of reproducing the inverted flag dynamics.In particular, the location and extent of the flapping is predicted reasonably well by the numerical model, and the predicted values of the angular span of

Numerical Simulations
The comparison between measurements and numerical predictions is provided in Figure 9a,b, where the predicted angular span of motion and frequency of oscillation are compared with the corresponding measurements for the no-solar and the full-solar flags.Overall, the agreement between measurements and predictions is encouraging.This indicates that, despite the simplifications previously discussed, the present numerical modeling approach is capable of reproducing the inverted flag dynamics.In particular, the location and extent of the flapping is predicted reasonably well by the numerical model, and the predicted values of the angular span of motion and frequency of oscillation compare well with the measurements.However, there is a noticeable difference for the no-solar flag that is predicted to flap at lower wind speeds with respect to the measurements.This can be explained by three-dimensional (3D) effects in the deformation of the flags, which were shown to increase stability in the conventional flag configuration [55] but are not captured by the present 2D numerical model.Interestingly, the agreement between predictions and measurements is better for the full-solar flag, even though the wind velocity for the onset of flapping is slightly over predicted.As better discussed in the next paragraph, this indicates that the stiffening of the flag due to the inclusion of the solar panels may reduce 3D effects, as compared with the no-solar flag, thus yielding a better agreement between measurements and 2D simulations.
It is important to note that the assumption of high aspect ratio implicit in a 2D simulation is expected to have an impact on the point of onset of flapping, where it was shown that, as the aspect ratio reduces, the velocity required to initiate flapping is increased.Aspect ratio effects are beyond the scope of the present work; instead, the reader is directed to the works of Sader et al. [56] and Tavallaeinejad et al. [57] for further details.
Predicted envelopes of motion for the no-solar and the full-solar flags are provided in Figure 9c, while accompanying videos are included as Supplementary Materials.Whilst the critical flow velocity for each state is different between the two cases, both exhibit a transition from a static zero-deflection state, through small-amplitude and large-amplitude flapping states, and finally a large-deflection static state.More importantly, in the full-solar case, the bending deformation is localized toward the root of the flag, with very little bending toward the tip due to the stiffening effect of the solar panels.As previously noted, this stiffening reduces 3D effects in the deformation of the flag, which is a beneficial outcome with respect to the applicability of 2D numerical simulations.On the other hand, the stiffening of the flag reduces the local deformation of the PVDF elements, thus reducing the mechanical strain and, therefore, the piezoelectric power generation, even though the additional power generation from the solar panels compensates for this (this is quantitatively assessed in Section 3.3 later on).Therefore, in the present inverted flag design, the solar panels are not only a destabilizing added mass; they also affect the local rigidity of the flag along its length, possibly reducing 3D effects.
Finally, the results of a sensitivity study carried out with the present numerical model for the full-solar flag are presented in Figure 10, where the predicted angular span and frequency of oscillation are displayed as functions of the wind speed for (a) varying flag length L p , (b) varying Young's modulus E (hence, stiffness of the flag), and (c) varying coverage ratio (CR).It is evident that the flag length and the Young's modulus have a significant on the dynamics of the flag, whilst the coverage ratio has a comparatively lower influence.In particular, shorter/longer flags flap at higher/lower wind velocities with higher/lower frequencies, consistently indicating that a longer cantilever is more prone to instability and, being heavier, flaps at a lower frequency with respect to a shorter and lighter one.Furthermore, increasing the Young's modulus makes the flag more rigid, thus increasing the wind velocity required to sustain flapping and the corresponding flapping frequency.This confirms that, as noted previously, having the PVDF elements custom-made to the desired mechanical properties will be a viable strategy to tailor the inverted flag harvester to the intended application, and the present numerical model would be instrumental in informing the harvester optimization.
The run-time of each simulation included in Figure 10 was on the order of 10 h (single thread): fast enough to allow time-efficient exploration of the (highly dimensional) parameter space of inverted flags.In addition to informing the future optimization of the harvester, the numerical tool will also be instrumental in informing future experiments on the fundamental coupled dynamics of multiple flag arrangements.In fact, having a preliminary indication on the dynamics of a multi-flag arrangement, for given flag geometry and constitutive material properties, would provide guidance on how to actually realize the flags, and would also inform the design of the experimental apparatus to track the flag dynamics and measure the flow field.motion and frequency of oscillation compare well with the measurements.However, there is a noticeable difference for the no-solar flag that is predicted to flap at lower wind speeds with respect to the measurements.This can be explained by three-dimensional (3D) effects in the deformation of the flags, which were shown to increase stability in the conventional flag configuration [55] but are not captured by the present 2D numerical model.Interestingly, the agreement between predictions and measurements is better for the full-solar flag, even though the wind velocity for the onset of flapping is slightly over predicted.As better discussed in the next paragraph, this indicates that the stiffening of the flag due to the inclusion of the solar panels may reduce 3D effects, as compared with the no-solar flag, thus yielding a better agreement between measurements and 2D simulations.It is important to note that the assumption of high aspect ratio implicit in a 2D simulation is expected to have an impact on the point of onset of flapping, where it was shown that, as the aspect ratio reduces, the velocity required to initiate flapping is increased.Aspect ratio effects are beyond the scope of the present work; instead, the reader is directed to the works of Sader et al. [56] and Tavallaeinejad et al. [57] for further details.other hand, the stiffening of the flag reduces the local deformation of the PVDF elements, thus reducing the mechanical strain and, therefore, the piezoelectric power generation, even though the additional power generation from the solar panels compensates for this (this is quantitatively assessed in Section 3.3 later on).Therefore, in the present inverted flag design, the solar panels are not only a destabilizing added mass; they also affect the local rigidity of the flag along its length, possibly reducing 3D effects.Finally, the results of a sensitivity study carried out with the present numerical model for the full-solar flag are presented in Figure 10, where the predicted angular span and frequency of oscillation are displayed as functions of the wind speed for (a) varying flag length Lp, (b) varying Young's modulus E (hence, stiffness of the flag), and (c) varying coverage ratio (CR).It is evident that the flag length and the Young's modulus have a significant impact on the dynamics of the flag, whilst the coverage ratio has a comparatively lower influence.In particular, shorter/longer flags flap at higher/lower wind velocities with higher/lower frequencies, consistently indicating that a longer cantilever is more prone to instability and, being heavier, flaps at a lower frequency with respect to a shorter and lighter one.Furthermore, increasing the Young's modulus makes the flag more rigid, thus increasing the wind velocity required to sustain flapping and the corresponding flapping frequency.This confirms that, as noted previously, having the PVDF elements custom-made to the desired mechanical properties will be a viable strategy to tailor the inverted flag harvester to the

Power Generation
The piezoelectric power, calculated from the root-mean-square (RMS) value of the measured voltage generated by the inverted flags for increasing wind speed, is provided in Figure 11a.As can be noted, up to a wind speed of about 7.5 m/s, the piezoelectric power generation is congruent with the angular span of motion in Figure 4a, consistently indicating that a larger angular span of motion yields a larger strain and, therefore, a larger power generation from the PVDF elements.This is no longer the case, however, at wind speeds above about 7.5 m/s, when the half-solar and the full-solar flags follow similar trends, while the no-solar flag yields more power.This confirms the simulation results previously discussed, indicating that the stiffening of the flag due to the inclusion of the solar panels reduces the local deformation of the PVDF elements, thus reducing the mechanical strain and, consequently, the piezoelectric power generation.
The piezoelectric, solar, and total power outputs generated by the full-solar flag for increasing wind speed are presented in Figure 11b.In particular, the solar power levels shown are the average values of power generated, while the total power is the sum of the piezoelectric and solar contributions.As can be noted in Figure 11b, the solar power output is higher when the flag is at rest, while it decreases during flapping.This is due to the lighting configuration adopted here, which is indicative in magnitude of an overcast day, although the lighting is more directional than what would be observed in a real ambient application.It is nonetheless evident in Figure 11b that the proposed piezo/solar inverted flag harvester is a promising concept, capable of harvesting kinetic wind energy and solar radiant energy at the same time and generating enough power (about 1 mW) to supply low-power electronic devices.In particular, the additional power generated by the solar panels more than compensates for the lower piezoelectric power output due to the stiffening of the flag from the inclusion of the solar panels.Perhaps even more importantly, the possibility to harvest wind and solar energies at the same time will accommodate for the inherent intermittency of these energy sources in real ambient applications.Wind and solar energies, in fact, have intermittencies that tend to mutually compensate, because the wind typically blows during stormy conditions when the sun does not shine, whilst, conversely, the sun typically shines on calm days with little wind.
to track the flag dynamics and measure the flow field.

Power Generation
The piezoelectric power, calculated from the root-mean-square (RMS) value of the measured voltage generated by the inverted flags for increasing wind speed, is provided in Figure 11a.As can be noted, up to a wind speed of about 7.5 m/s, the piezoelectric power generation is congruent with the angular span of motion in Figure 4a, consistently indicating that a larger angular span of motion yields a larger strain and, therefore, a larger power generation from the PVDF elements.This is no longer the case, however, at wind speeds above about 7.5 m/s, when the half-solar and the full-solar flags follow similar trends, while the no-solar flag yields more power.This confirms the simulation results previously discussed, indicating that the stiffening of the flag due to the inclusion of the solar panels reduces the local deformation of the PVDF elements, thus reducing the mechanical strain and, consequently, the piezoelectric power generation.The piezoelectric, solar, and total power outputs generated by the full-solar flag for increasing wind speed are presented in Figure 11b.In particular, the solar power levels shown are the average values of power generated, while the total power is the sum of the piezoelectric and solar contributions.As can be noted in Figure 11b, the solar power output is higher when the flag is at rest, while it decreases during flapping.This is due to the lighting configuration adopted here, which is indicative in magnitude of an overcast day, although the lighting is more directional than what would be observed in a real ambient application.It is nonetheless evident in Figure 11b that the proposed piezo/solar inverted flag harvester is a promising concept, capable of harvesting kinetic wind energy and solar radiant energy at the same time and generating enough power (about 1 mW) to supply lowpower electronic devices.In particular, the additional power generated by the solar panels more than compensates for the lower piezoelectric power output due to the stiffening of the flag from the inclusion of the solar panels.Perhaps even more importantly, the possibility to harvest wind and solar energies at the same time will accommodate for the inherent intermittency of these energy sources in real ambient applications.Wind and solar energies, in fact, have intermittencies that tend To better assess the performance of the present inverted flag harvester, Table 1 provides a comparison against other concepts available in the literature.Note that, for consistency, the comparison here is only limited to harvesters based on PVDF elements and tested in controlled wind conditions.Following common practice in energy harvesting, the power densities included in Table 1 were calculated based on the total volume of the harvesters.As can be noted, the power density of the present inverted flag harvester is rather good, and the inclusion of the solar panels almost doubles the power density.This is encouraging, particularly when considering that the present work was essentially a proof of concept, which motivates further study to improve and optimize the proposed harvester.

Conclusions
The flow-induced dynamics and power generation of an energy harvester based on the inverted flag design were investigated with experiments and numerical simulations.The energy harvester was realized using off-the-shelf flexible piezoelectric PVDF strips and flexible mini solar panels to simultaneously harvest wind kinetic energy and solar radiant energy, with a view at compensating the inherent intermittency of ambient energy sources.In the proposed design, the PVDF strips made up the flag harvester, thus providing the structural support, while the solar panels were incorporated as tip masses close to the flag free edge.The performance of the energy harvester was experimentally investigated under controlled wind (up to 10 m/s wind speed) and light exposure (1.8 klx), and simulations based on a lattice-Boltzmann/finite-element/immersed-boundary numerical method were used to corroborate and integrate the measurements.The measurements indicate that the proposed harvester is a promising concept, capable of producing enough power (on the order of 1 mW) main mounting plates, as shown in Figure A1.This way, contact of the upper and lower plates to the electrical terminals of the PVDFs is avoided, allowing a strong restraint of the flag when experiencing high levels of vibration.The small holes in the lower plate shown in Figure A1 are for clamping the upper and lower plates together through matching holes in the upper plate.The larger holes are for attaching the whole unit to the wind tunnel pole.The smaller holes are spaced so that the trailing wires from the PVDFs can pass between them.

Figure 1 .
Figure 1.Configuration of the piezo-solar flags used in this study: a flag is made of four flexible polyvinylidene fluoride (PVDF) elements and two mini solar panels attached to both sides at the tip.(a) The dimensions used to describe the geometry of a flag; (b) no-solar flag (coverage ratio (CR) = 0); (c) half-solar flag (CR = 0.32); (d) full-solar flag (CR = 0.64) (the flow direction is left to right).

Figure 1 .
Figure 1.Configuration of the piezo-solar flags used in this study: a flag is made of four flexible polyvinylidene fluoride (PVDF) elements and two mini solar panels attached to both sides at the tip.(a) The dimensions used to describe the geometry of a flag; (b) no-solar flag (coverage ratio (CR) = 0); (c) half-solar flag (CR = 0.32); (d) full-solar flag (CR = 0.64) (the flow direction is left to right).

Figure 2 .
Figure 2. Schematics of the experimental set-up used to characterize the dynamics and power generation of the flags subjected to controlled wind and light excitation.

Figure 2 .
Figure 2. Schematics of the experimental set-up used to characterize the dynamics and power generation of the flags subjected to controlled wind and light excitation.

Figure 4 .
Figure 4. Experimental measurements of the dynamics of the inverted flags: (a) angular span of motion vs. wind speed; (b) dominant frequency of motion vs. wind speed.

Figure 4 .
Figure 4. Experimental measurements of the dynamics of the inverted flags: (a) angular span of motion vs. wind speed; (b) dominant frequency of motion vs. wind speed.

Figure 5 .
Figure 5. Representative dynamics of aligned-vibrating regime for the full-solar flag at wind speed of 4.7 m/s: (a) angular position time series (5 s); (b) angular position power spectral density (PSD); (c) angular position autocorrelation function (ACF); (d) angular position histogram; (e) system trajectory in phase-space (reconstructed attractor).

Figure 5 .
Figure 5. Representative dynamics of aligned-vibrating regime for the full-solar flag at wind speed of 4.7 m/s: (a) angular position time series (5 s); (b) angular position power spectral density (PSD); (c) angular position autocorrelation function (ACF); (d) angular position histogram; (e) system trajectory in phase-space (reconstructed attractor).

Figure 6 .
Figure 6.Representative dynamics of flapping regime for the full-solar flag at wind speed of 9.0 m/s: (a) angular position time series (5 s); (b) angular position power spectral density (PSD); (c) angular position autocorrelation function (ACF); (d) angular position histogram; (e) system trajectory in phase-space (reconstructed attractor).

6 .
Representative dynamics of flapping regime for the full-solar flag at wind speed of 9.0 m/s: (a) angular position time series (5 s); (b) angular position power spectral density (PSD); (c) angular position autocorrelation function (ACF); (d) angular position histogram; (e) system trajectory in phase-space (reconstructed attractor).

Figure 7 .
Figure 7. Representative dynamics of deflected regime for the full-solar flag at wind speed of 9.7 m/s: (a) angular position time series (5 s); (b) angular position power spectral density (PSD); (c) angular position autocorrelation function (ACF); (d) angular position histogram; (e) system trajectory in phase-space (reconstructed attractor).

Figure 7 .
Figure 7. Representative dynamics of deflected regime for the full-solar flag at wind speed of 9.7 m/s: (a) angular position time series (5 s); (b) angular position power spectral density (PSD); (c) angular position autocorrelation function (ACF); (d) angular position histogram; (e) system trajectory in phase-space (reconstructed attractor).

Figure 8 .
Figure 8. Representative hysteresis cycle for the full-solar flag.

Figure 8 .
Figure 8. Representative hysteresis cycle for the full-solar flag.

Figure 9 .
Figure 9.Comparison of the numerical solver simulation results with experimental measurements: (a) no-solar flag; (b) full-solar flag.(i) Angular span of motion vs. wind speed; (ii) dominant frequency vs. wind speed; (c) envelopes of motion from the numerical simulations at various wind speeds.

Figure 9 .
Figure 9.Comparison of the numerical solver simulation results with experimental measurements: (a) no-solar flag; (b) full-solar flag.(i) Angular span of motion vs. wind speed; (ii) dominant frequency vs. wind speed; (c) envelopes of motion from the numerical simulations at various wind speeds.

Figure 10 .
Figure 10.Sensitivity analysis on the dynamics of the full-solar flag.Effect of varying (a) the flag length, (b) the Young's modulus, and (c) the coverage ratio.(i) Angular span of motion vs. wind speed; (ii) flapping frequency vs. wind speed.

10 .
Sensitivity analysis on the dynamics of the full-solar flag.Effect of varying (a) the flag length, (b) the Young's modulus, and (c) the coverage ratio.(i) Angular span of motion vs. wind speed; (ii) flapping frequency vs. wind speed.

Figure 11 .
Figure 11.Measured power generation for the inverted flags: (a) root-mean-square (RMS) power generated from the four PVDF elements at different wind speeds; (b) total extracted power from the full-solar flag including the contributions from the four PVDF elements and the two full-solar panels (at 1.8 klx constant illumination).

Figure 11 .
Figure 11.Measured power generation for the inverted flags: (a) root-mean-square (RMS) power generated from the four PVDF elements at different wind speeds; (b) total extracted power from the full-solar flag including the contributions from the four PVDF elements and the two full-solar panels (at 1.8 klx constant illumination).

Figure A1 .
Figure A1.Computer-aided design (CAD) drawing of the support clamp (the upper plate was removed for clarity).

Figure A1 .
Figure A1.Computer-aided design (CAD) drawing of the support clamp (the upper plate was removed for clarity).

Table 1 .
Power density comparison of different harvesters based on polyvinylidene fluoride (PVDF) elements.