Mohammad Tavallaeinejad

Mathias Legrand

Michael Païdoussis

Nonlinear dynamics of slender inverted flags in uniform steady flows

Keywords: inverted flag, bistability, nonlinear dynamics, fluid-structure interactions Nomenclature

à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

motion. A familiar example is the flapping motion of leaves in the wind. There has been ongoing research on this fluid-structure interaction problem over the past few years, with the aim of understanding the essential dynamics. These studies may be classified into three categories: experimental investigations, computational fluid dynamics (CFD) analyses utilizing numerical solvers for the Navier-Stokes equations, and analytical studies using simplified continuum expressions for the aerodynamic forces.

Starting with analytical and experimental studies, a linear stability analysis was first carried out by Guo & Païdoussis [START_REF] Qin | Stability of rectangular plates with free side-edges in twodimensional inviscid channel flow[END_REF] who employed inviscid potential flow theory to predict the onset of oscillatory instability. They showed that free-clamped plates (inverted flags) flutter at critical velocities inversely proportional to their mass ratios 1 . The flapping motion of inverted flags subjected to uniform air and water flows was explored in the pioneering experimental work of Kim et al. [START_REF] Kim | Flapping dynamics of an inverted flag[END_REF]. Three regimes were reported as the flow velocity was varied: the straight (undeflected) stationary regime, the flapping regime with large amplitude, and a fully deflected static regime. Bistable states were also observed in the transition from straight-to-flapping and flapping-to-deflected regimes. This work was pursued by Sader et al. [START_REF] Sader | Large-amplitude flapping of an inverted flag in a uniform steady flow -a vortex-induced vibration[END_REF][START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF] who conducted experiments on inverted flags of various aspect ratios. Their experimental observations revealed that the aspect ratio of a flag can dramatically alter the dynamics. In particular, slender flags exhibit a buckling behaviour via a static divergence instability. Intermediately wide flags, on the other hand, undergo an abrupt flapping motion prior to a fully deflected regime, with a highly curved shape, via buckling. Moreover, they reported that the undeflected static equilibrium of inverted flags with smaller aspect ratios loses stability at higher critical flow velocities. In their first study, in conjunction with mathematical modelling based on a linear aerodynamic theory and a scaling analysis, Sader et al. [START_REF] Sader | Large-amplitude flapping of an inverted flag in a uniform steady flow -a vortex-induced vibration[END_REF] concluded that flapping motion of inverted flags subjected to heavy fluid loading (light flags) is a vortex-induced vibration (VIV). For light fluid loading (heavy flags), however, the inverted flag cannot undergo a VIV, despite displaying large-amplitude flapping; the comprehension of the mechanism involved requires further work. In their second study, in addition to experimental measurements, Sader et al. [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF] undertook a nonlinear theoretical investigation and provided an analytical framework to explore the stability and dynamics of inverted flags in the asymptotic limit of zero aspect ratio, within the restricted framework of steady fluid forces. Moreover, they formulated a mathematical solution to predict the onset of the static bifurcation for inverted flags with various aspect ratios. More specifically, they employed a Padé approximant to correlate the two solutions for zero-aspect-ratio and two-dimensional flow limits, yielding a formula applicable over the entire range of aspect ratios studied; a very good quantitative agreement with measurements was found. They offered a reasonable qualitative explanation that the deflected shape of the flag, originating from the the static divergence point, presents a sharp leading edge to the impinging flow, and the resulting unsteady vortex shedding causes the flag to flap. This analysis was subsequently extended by the authors [START_REF] Tavallaeinejad | Nonlinear static response of lowaspect-ratio inverted flags subjected to a steady flow[END_REF] to deal with intermediate aspect ratios, as well as asymmetric behaviour of inverted flags in the presence of an initial curvature along the plate.

All analytical models to date are restricted to steady-flow theories, and do not provide any mathematical formulation or numerical framework to explore the effect of time-dependent reactive forces. Such modelling disallows the possibility of capturing flapping motion, which is a feature of inverted flags of sufficiently large aspect ratio. The need of modelling such effects has motivated the present study.

From the computational fluid dynamics perspective, various investigations have contributed to the field. For instance, Tang et al. [START_REF] Tang | Dynamics of an inverted flexible plate in a uniform flow[END_REF] performed a three-dimensional (3D) simulation using an immersed boundary lattice-Boltzmann approach focusing on low Reynolds number flows. They concluded that wider plates oscillate with a larger amplitude and undergo flapping motion at lower critical flow velocities. They also found that the stretchedstraight state of lighter plates is more stable. The effect of an initial angle of attack was also explored and it was shown that an asymmetric flapping regime emerges, which becomes narrower with an increasing angle of attack. In another study, Shoele and Mittal [START_REF] Shoele | Flutter instability of a thin flexible plate in a channel[END_REF] investigated the energy-harvesting performance of inverted piezoelectric flags using an immersed boundary method. They also examined the effect of the initial inclination of an inverted flag to the impinging flow and showed that the dynamic response of the system is relatively insensitive to small deviations in the initial inclination angle. Gurugubelli and Jaiman [START_REF] Pardha | Self-induced flapping dynamics of a flexible inverted foil in a uniform flow[END_REF] employed a second-order scheme combining the finite element procedure with the arbitrary Lagrangian-Eulerian (ALE) approach to conduct a two-dimensional (2D) numerical simulation on the inverted flag problem and found that the mass ratio affects the transition from large-amplitude flapping motion to moderate flapping around the deflected position of the inverted flag. In another study [START_REF] Pardha | Large amplitude flapping of an inverted elastic foil in uniform flow with spanwise periodicity[END_REF], they performed 3D numerical experiments, considering periodic boundary conditions in the spanwise direction while ignoring the side vortices to explore the influence of trailing edge vortices on the large-amplitude flapping of inverted flags. In a recent study, Goza et al. [START_REF] Goza | Global modes and nonlinear analysis of inverted-flag flapping[END_REF] employed an immersed boundary method to study the nonlinear response and stability of inverted flags and explored the effect of Reynolds number as well as the mass ratio parameter on the system dynamics for a wide range of flow velocities. They concluded that the mass ratio has a negligible effect on the onset of flapping. They also concluded that heavy flags are subject to a 'non-classical-VIV' flapping mechanism, in which, as opposed to 'classical VIV', the flapping regime is not the result of frequency synchronisation between flapping and vortex shedding.

The analytical work by Sader et al. [START_REF] Sader | Large-amplitude flapping of an inverted flag in a uniform steady flow -a vortex-induced vibration[END_REF] and the numerical simulations by Goza et al. [START_REF] Goza | Global modes and nonlinear analysis of inverted-flag flapping[END_REF] showed the susceptibility of VIV to the mass ratio parameter. In contrast to light flags, the synchronisation of flapping and vortex shedding is lost in the case of heavy flags. Goza et al. [START_REF] Goza | Global modes and nonlinear analysis of inverted-flag flapping[END_REF] moreover demonstrated that at very low Reynolds numbers, where vortex shedding is essentially absent, large-amplitude flapping still occurs for heavy flags. Hence, they concluded that the mechanism generating the flapping motion of heavy flags is distinct from VIV. They attributed the flapping to the instability of the deformed equilibria of the system.

The non-VIV flapping motion identified by Goza et al. [START_REF] Goza | Global modes and nonlinear analysis of inverted-flag flapping[END_REF] has motivated the current study. The emergence of flapping motion for low mass ratio flags in low Reynolds number flow, where vortex shedding does not occur, impelled us to consider the importance of the effect of unsteady forces arising from changes in the fluid-flow momentum due to motion of the immersed flag. Of particular interest is the investigation of flow-induced oscillation of heavy inverted flags with mass ratios of O.1/, with the aid of well-established inviscid flow theories available in the literature. Indeed, the nonlinear large-amplitude flapping of an inverted flag is a challenging fully coupled fluid-structure interaction problem for which computational techniques are particularly expensive. Despite advances in CFD techniques and computing power, analytical models are still preferable for conducting simplified numerical simulations and interpreting experimental observations in a straightforward fashion. To further support the use of inviscid potential flow theory, it has been observed experimentally that the flow around the leading edge of the inverted flag reattaches at middle stroke, where the angle of attack is relatively large [START_REF] Yu | Vortex dynamics behind a self-oscillating inverted flag placed in a channel flow: Time-resolved particle image velocimetry measurements[END_REF]. In addition, analytical inviscid expressions for the normal force acting on low-aspect-ratio plates have yielded good predictions for angles of attack up to approximately 40 ı , before massive separation at the leading edge and stall occurs [START_REF] Bollay | A non-linear wing theory and its application to rectangular wings of small aspect ratio[END_REF][START_REF] Torres | Aerodynamic characteristics of low aspect ratio wings at low Reynolds numbers, fixed and flapping wing aerodynamics for micro air vehicle applications[END_REF][START_REF] Devoria | A vortex model for forces and moments on low-aspect-ratio wings in side-slip with experimental validation[END_REF]. Accordingly, inviscid forces are dominant, at least at the region with moderate angles of attack, thus justifying the preference for inviscid potential flow theory.

To the best of our knowledge, there is no analytical study in the literature on the nonlinear dynamical response of inverted flags subjected to axial flow accounting for velocity-dependent forces. The present paper aims at producing such a model in the limit case of slender flexible plates, and testing its capability to predict the dynamical behaviour of heavy slender inverted flags in axial flow, such as the onset of instabilities, frequency and amplitude of flapping, as well as elucidating the physical mechanisms involved in the transition between different dynamical states.

To this end, continuum representations of fluid-flow forces are detailed on the basis of the large amplitude elongated body theory [START_REF] James | Large-amplitude elongated-body theory of fish locomotion[END_REF] to formulate the reactive force and Bollay's model to express the nonlinear vortex lift mechanism, with no account for the vortex shedding mechanism from the leading edge. Hamilton's principle is employed to balance the kinetic and potential energies of the inverted flag and the virtual work of flow forces and damping. More specifically, the kinematics is covered in Section 2.1, constitutive equations in Section 2.2, and conditions of equilibrium leading to final equation of motion in Section 2.3. The resultant partial-integro-differential equation is discretized spatially by means of a conventional Galerkin procedure and is recast in a set of nonlinear ordinary differential equations (ODEs). The nonlinear response of the system is explored in detail by utilizing direct time-stepping integration together with bifurcation analysis. In particular, the ODEs are solved in time using Gear's backward differentiation, and bifurcation diagrams are constructed with the aid of the pseudo-arclength continuation technique. The sensitivity of the dynamics to the aspect ratio, mass ratio, incidence angle, and viscous drag is also explored. The results are summarized in the form of bifurcation diagrams to identify transitions between various regimes.

Analytical model

In this section, a nonlinear model is developed for the inverted flag/plate problem. First, the problem is defined mathematically. Expressions for fluid forces are formulated separately, and the nonlinear integro-partial-differential equations governing the rotation of the mid-plane are derived. These equations are non-dimensionalized and then discretized via the Galerkin method.

Problem geometry

The system considered is shown in Figure 1(a) consisting of a vertical cantilevered thin plate subjected to an inviscid axial flow impinging to its free end. The plate has length L, width H , thickness h, transverse moment of inertia I , cross-sectional area A, surface area S , and flexural rigidity D. The mechanical properties of the plate, i.e., its density, Poisson ratio, Young's modulus, and internal damping coefficient are denoted by p , , E, and , respectively; f stands for the density of the fluid flowing with mean flow velocity U . A right-handed Cartesian reference system .OI x; y; z/ is considered, with axes x i and z k being the axial and the transverse directions, respectively. The transverse and longitudinal motions w.x; t/ and u.x; t/ of a generic point at a distance z from the mid-plane on the cross-section of the flag are shown in Figure 1(b); t denotes time. The inextensibility assumption [START_REF] Semler | The non-linear equations of motion of pipes conveying fluid[END_REF] reads .1 C @ x u.x; t// 2 C @ x w.x; t/ 2 D 1, with @ x denoting the spatial derivative; it enforces a constant flag length during deformation. As a consequence, w and u can be expressed in terms of only the slope , which is also the rotation angle of the cross-section. Given that @ x w.x; t/ D cos .x; t/ and 1 C @ x u.x; t/ D sin.x; t/, the corresponding velocities are formulated as

P
Defining the instantaneous angle of attack ˛as the instantaneous inclination between the chord-line and the relative velocity vector, one obtains

˛D atan P w U C P u : (3)

Fluid modelling

The purpose of the present study is to gain insight into the factors affecting the dynamics of small aspect ratio inverted flags. While this could be treated by fully-coupled fluid-structure interaction (FSI) solvers and CFD schemes, it is of interest to develop a model which is predominantly analytical. This provides an opportunity to better understand the underlying aerodynamic effects. Besides, direct numerical simulation of such a complex FSI problem is particularly expensive. Hence, an analytical expression of the fluid flow forces is developed in the present study. The flow is analyzed as if it were inviscid in order to find the pressure-related forces, provided that there is no flow separation at the leading edge, even for relatively large angles of attack. The strong vortex formation along the lateral edges of the flag is approximated by Bollay's nonlinear wing theory [START_REF] Bollay | A non-linear wing theory and its application to rectangular wings of small aspect ratio[END_REF]. Then, semi-analytical expressions are added for viscous forces acting on the plate. It is nevertheless recognized that the applicability of inviscid theory depends on whether the flow remains attached. Hence, the predictions for aspect ratios larger than unity, for which massive separation occurs, are limited to small angles of attack [START_REF] Devoria | A vortex model for forces and moments on low-aspect-ratio wings in side-slip with experimental validation[END_REF].

The proposed model is based on the large amplitude elongated body theory [START_REF] James | Large-amplitude elongated-body theory of fish locomotion[END_REF] to obtain the reactive fluid forces associated with the relative motion of the inverted flag with respect to the incident flow. This model has been successfully employed in a number of studies, showing the reliability of predictions of the flapping motion in comparison with Navier-Stokes simulations [START_REF] Candelier | Three-dimensional extension of Lighthill's largeamplitude elongated-body theory of fish locomotion[END_REF] and a full three-dimensional panel method [START_REF] Yu | Extension of Lighthill's slender-body theory to moderate aspect ratios[END_REF] for aspect ratios of O.1/ [START_REF] Doaré | Piezoelectric coupling in energy-harvesting fluttering flexible plates: linear stability analysis and conversion efficiency[END_REF][START_REF] Yadykin | The flow-induced vibration of a flexible strip hanging vertically in a parallel flow. Part 1: Temporal aeroelastic instability[END_REF][START_REF] Eloy | The origin of hysteresis in the flag instability[END_REF]. According to the large amplitude elongated body theory, the reactive force relates to the balance of the rate of change of solid and fluid momenta and the pressure force acting in the tangential direction corresponding to the quadratic terms in Bernoulli's equation [START_REF] James | Aquatic animal propulsion of high hydromechanical efficiency[END_REF]. As such, the inviscid pressure force is expressed as

F N D 4 f H 2 h @ t .V n n/ C @ x V n V n 1 2 V 2 n i (4)
The n component of F N gives the magnitude of the reactive force (per unit length) acting on the flag, as shown in Figure 2. Inserting the expressions for the normal and tangential components of the relative velocity in Eq. (4) gives

.x; i/ .z; k/ .x;t/ F F N F R n
Figure 2: Schematic of the forces acting on a generic infinitesimal element of the plate.

F N D 4 f H 2 R w cos R u sin C 2 P . P u cos C P w sin / C 3 2 @ x . P w 2 P u 2 / cos 2 C @ x . P u 2 =2 3 P u P w cos sin P w 2 / C U.. P u 3 P u cos 2 3 P w sin cos / @ x C 2 P cos / 1 2 U 2 .3 cos 2 1/ @ x : (5)
The nonlinear vortex lift associated to the circulation of the side-edge vortices and the distribution of spanwise vorticity is devised based on Bollay's nonlinear wing theory [START_REF] Bollay | A non-linear wing theory and its application to rectangular wings of small aspect ratio[END_REF]. Briefly, Bollay obtained the induced velocity originating from the system of vortices at a generic point of the plate and calculated the normal force per unit length associated with the cross-flow component of the flow velocity. The expression for the normal force given in [START_REF] Bollay | A non-linear wing theory and its application to rectangular wings of small aspect ratio[END_REF] is not explicit, though approximate expressions can be sought based on the assumption made by Gersten in [START_REF] Gersten | A nonlinear lifting-surface theory especially for low aspect-Ratio wings[END_REF] 2 or analogously based on Polhamus' work [START_REF] Polhamus | A concept of the vortex lift of sharp-edge delta wings based on a leading-edge-suction analogy[END_REF][START_REF] Tavallaeinejad | Nonlinear static response of lowaspect-ratio inverted flags subjected to a steady flow[END_REF]. Hence, the normal force reads

F R D f H U 2 2 .K p .A/ sin ˛cos ˛C K v .A/ sin 2 ˛/ (6)
where K p and K v are coefficients reflecting the circulation and side-edge vortex flow, and are determined for various aspect ratios, according to the technique implemented in [START_REF] Tavallaeinejad | Nonlinear static response of lowaspect-ratio inverted flags subjected to a steady flow[END_REF]; A denotes the aspect ratio of the flag. The curves obtained for the K p and K v coefficients over a range of aspect ratios is shown in Figure 3; the reader is referred to a previous study [START_REF] Tavallaeinejad | Nonlinear static response of lowaspect-ratio inverted flags subjected to a steady flow[END_REF] by the authors for more details. These values depend only on the aspect ratio of the inverted flag. For instance, for the case of an infinitely long plate (i.e., zero aspect ratio), K v D 2 and K p D 0. The value of K v (K p ) decreases (increases) with an increase in the flag aspect ratio and approaches zero (2) as the aspect ratio approaches infinity. Inserting Eq. ( 3) in Eq. ( 6), one obtains

0 0.5 1 1.5 2 0 1 2 3 4 A Kp, Kv
F R D f H U 2 2. P u C U / 2 C 2 P w 2 K p .. P u C U / 2 P w 2 / sin cos . P w P u C U P w/ cos 2 CK v . P u C U / 2 sin 2 C P w 2 cos 2 2. P w P u C U P w/ sin cos : (7)
The uniform distribution of tangential aerodynamic force acting on the flag due to the viscous effects of the fluid flow simplifies to

F D 1 2 f H U 2 K .Re/ cos ˛(8)
where K is the coefficient of the effective drag caused by viscous stresses acting on the flag and can be determined using the semi-empirical expressions proposed by Taylor [START_REF] Taylor | Analysis of the swimming of long and narrow animals[END_REF] (see, for instance [START_REF] Schouveiler | Flow-induced vibrations of high mass ratio flexible filaments freely hanging in a flow[END_REF]) or utilizing a Blasius-type solution for either turbulent or laminar flows (see, for instance [START_REF] Yadykin | The flow-induced vibration of a flexible strip hanging vertically in a parallel flow. Part 1: Temporal aeroelastic instability[END_REF]). This allows for the approximate treatment of the viscous drag. The effect of K on the system dynamics is discussed through extensive numerical tests in Section 3.6.

Equation of motion

The inverted flag is idealized as a cantilevered elastic beam shown in Figure 1. The conventional nonlinear Euler-Bernoulli beam theory is considered and the equation of motion is obtained via Hamilton's principle. The variation of the kinetic energy of the inverted flag system is given by

Z t 2 t 1 ıT .t/ dt D Z t 2 t 1 p A Z L 0 . P w ı P w C P u ı P u/ dx C p I Z L 0 P ı P dx dt (9)
I D H h3 =12 being the second moment of inertia. The expression for the variation of the strain energy of the system is expressed as

Z t 2 t 1 ıV.t/ dt D Z t 2 t 1 D Z L 0 @ x .x; t/ @ x ı .x; t/ dx dt; (10)
where D D EI denotes the flexural rigidity of a slender inverted flag; for a flag of large aspect ratio, the plane-strain flexural rigidity D D EI =.1 2 / is considered. The virtual work associated with the fluid-dynamic forces can be expressed as

ıW f .t/ D Z L 0 .F N C F R / cos F .x; t/ sin ıw C .F N F R / sin F cos ıu dx: (11)
The virtual work done by the Kelvin-Voigt structural damping and the dissipation due to fluid viscosity is

ıW d .t/ D D Z L 0 @ xx P ı dx c Z L 0 @ t Z x 0 cos .; t/ ds ıu C @ t Z x 0 sin .; t/ ds ıw dx ( 12
)
where is the material viscosity damping and c the fluid damping coefficient [START_REF] Eloy | The origin of hysteresis in the flag instability[END_REF] 3 . Inserting expressions ( 8), ( 7) and ( 5) into Eq. ( 11) and substituting the resultant along with Eqs. ( 9), ( 10) and ( 12) into the expression for the extended Hamilton's principle

Z t 2 t 1 .ıT .t/ ıV.t/ C ıW f .t/ C ıW d .t// dt D 0 (13) yields p I @ tt p A sin @ tt Z L x Z s 0 cos .; t/ d ds C p A cos @ tt Z L x Z s 0 sin .; t/ d ds D @ txx D @ xx C cos Z x L .F N .s; t/ F
Utilizing the dimensionless parameters where q r .t/ are the unknown time-dependent generalized coordinates and M is the number of modes used in the analysis. Plugging Eq. ( 16) into Eq. ( 1) and Eq. ( 3), inserting the resultant in Eq. ( 15), and applying Galerkin's technique, yields a set of M second-order ODEs with trigonometric nonlinearities.

x D x=L, t D t=, ˇD I =.AL 2 /, k D K =2, k p D K p =2, k v D K v =2, s D =, f D f , D f L=. p h/, A D H=L, … D f SL 2 U 2 =D,

System dynamics

The set of ODEs is solved by means of the AUTO software [START_REF] Doedel | Continuation and Bifurcation Software for Ordinary Differential Equations (with HomCont)[END_REF][START_REF] Modarres-Sadeghi | Nonlinear dynamics of a slender flexible cylinder subjected to axial flow[END_REF] which uses the pseudo-arclength continuation technique to construct bifurcation diagrams. This technique performs continuation of both stable and unstable solution branches. Time integration is performed via Gear's backward differentiation formula [START_REF] William | The automatic integration of ordinary differential equations[END_REF]. For the problem at hand, the DIVPAG routine of the IMSL library is used. Time histories of the amplitude of oscillation are used to plot phase-plane portraits, Fast Fourier transforms (FFT), and static as well as flapping flag shapes at points of interest in the parameter space. The FFT is obtained once the transient response has decayed, and the mean amplitude of each harmonic is calculated. Flag shapes are retrieved through Eq. ( 1). In all bifurcation diagrams, continuous and dashed lines represent stable and unstable branches of the solution, respectively. Each colour used therein corresponds to a specific regime in the range of flow velocities considered.

Numerical parameters

In this study, six modes are employed to discretize Eq. ( 15) (i.e. M D 6); this is sufficient to reach convergence. Moreover, the following parameters are utilized [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF]: L D 30 cm, h D 1 mm, p D 1200 kg m 3 , f D 1:2 kg m 3 , E D 2:4 GPa, D 0:38, and s D 0:008. The viscous damping coefficient c d is replaced by 2! 1 where is the modal damping ratio and ! 1 is the first dimensionless natural frequency. The modal damping ratio is set at 0:05. The dynamics of the system depends on five dimensionless parameters: the dimensionless flow velocity …, the flag aspect ratio A, the incidence angle 0 , and the viscous drag coefficient k . Here, … is selected as the bifurcation parameter while other parameters are kept constant. velocity is increased, a sequence of regimes emerges as follows:

Possible regimes

i) The undeflected static equilibrium where the inverted flag is stable. Small disturbances generate a motion that ultimately dies out. ii) The undeflected static equilibrium loses stability at point A where … D 0:91 via a supercritical Hopf bifurcation from which emanates limit-cycle motion 4 . In other words, the inverted flag then undergoes a large-amplitude flapping motion around the undeflected static equilibrium. This regime survives over a range of … 2 OE0:91 8:71 and the amplitude of flapping increases with flow velocity. Figure 5 shows the time history, phase-plane portrait and FFT-plot harmonic content with the frequency of oscillation f D 0:40 at … D 6:00, from which a periodic flapping motion is identified which develops after the transient regime. iv) The flapping motion ceases upon further increase in the flow velocity. The system response switches to a fully-deflected regime via a Hopf bifurcation at point D where … D 9:47. Beyond that, the response of the system is attracted to stable large-amplitude deflected equilibria, and the inverted flag bends over one side.

The tip rotation increases with flow velocity and the tip transverse deflection remains almost constant at its maximum value.

It should be noted that the response of the system is symmetric with respect to the undeflected equilibrium state and the inverted flag may flap around deflected equilibria (buckled state) on either side. In this paper, flapping around one side is plotted only, so that one can properly distinguish flapping around the undeflected static equilibrium and flapping around the deflected static equilibria. Both trivial and nontrivial unstable equilibria of the system, connecting the two Hopf bifurcation points, are shown in the bifurcation diagram as dashed lines. It is seen that, for the system studied here, the nontrivial static solution emerges from the bifurcation point E where … D 10:6 and folds at point F where … D 6:6.

The variation of the dimensionless frequency with … is plotted in Figure 7. With increasing …, the flapping frequency decreases first and then recovers slightly at … D 2:95. The frequency increases sharply and monotonically for … 2 OE8:53 9:47, wherein the flapping around the deflected equilibrium takes place. This change in trend is accompanied by an abrupt increase in frequency. In order to compare the response of the system in different regimes further, the (flapping) shapes of the inverted flag system are plotted in Figure 8. Subfigure 8(a) shows the undeflected flag for velocities below the critical value, 0 that is … 0:91. In subfigure 8(b), the inverted flag displays large-amplitude flapping motion at … D 6:00. It can be noticed that the transverse motion of the system is predominantly in the fundamental mode of a cantilevered beam. The deformed-flapping regime (oscillation around the deflected equilibrium) can be observed in subfigure 8(c) where … D 9:00, and finally, the fully-deflected state is shown in subfigure 8(d) for … D 20:00. It is nevertheless recognized that the validity of the proposed model deteriorates when the flag exhibits highly curved shapes. Indeed, at very large angles of attack, massive separation occurs leading to the normal force stall, which in turn can cause a sudden decrease in the amplitude of oscillation. The normal force stall phenomenon cannot be captured by inviscid flow theory, which naturally, neglects such viscous effects. Furthermore, the model ignores the effects of the unsteady vortex shedding at high flow rate in the fully-deflected regime; the results shown for this regime should be regarded with some caution, considering this limitation.

1 -1 0 1 Z 0 1 0 1 -1 0 1 X Z 0 1 X (a) (b) (c) (d)

Sensitivity to aspect ratio

It is clear from Eqs. ( 5) and ( 7) that the coefficients of fluid forces are functions of the aspect ratio. Accordingly, the contribution of the fluid forces to the system dynamics varies with the aspect ratio. In order to investigate the effect of the aspect ratio on the response, four different cases are considered: flags of H D 3 cm, H deflected equilibria emerging from a saddle-node bifurcation at … D 8:6. This behaviour agrees with other existing analytical studies [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF][START_REF] Tavallaeinejad | Nonlinear static response of lowaspect-ratio inverted flags subjected to a steady flow[END_REF]. One underlying reason for this is that when the aspect ratio is small, the contribution of the reactive component of the fluid forces is weak and the vortex force and viscous damping may dominate. In this case, the steady parts of the fluid forces are responsible for initiating the instability, in the form of a static divergence, also known as buckling, via a pitchfork bifurcation. Moreover, the response of the system is strongly subcritical due to the contribution of the nonlinear term in the normal force acting on the inverted flag (see Eq. ( 6)). More specifically, for small aspect ratios, k v ! 2 while k p ! 0. As such, the solution branches meet at relatively large … (e.g. … D 33:1 for A D 0:1, Figure 9(a)) where new solutions branch off from the trivial solution and fold at a relatively small … (e.g. … D 8:6 for A D 0:1). This results in a wide range of bistability.

As the aspect ratio becomes larger, in Figure 9(b)-9(d), the destabilizing reactive component of the fluid forces comes into play. This materializes as a supercritical Hopf bifurcation prior to the pitchfork bifurcation, thereby establishing a limit-cycle oscillation. More specifically, the Coriolis term in Eq. ( 15) proportional to … generates a reaction on the inverted flag in the direction of motion. Inverted flags of sufficiently large aspect ratios undergo a large-amplitude flapping motion via a supercritical Hopf bifurcation. In particular, a Hopf bifurcation materializes at … D 4:41 preceding the saddle-node bifurcation at … D 7:75 for the flag with A D 0:25 (Figure 9(b)), giving rise to oscillation with relatively small amplitude since max t 2OE0 T .w.L; t// D 0:08 at … D 8:0. A range of bistability can be observed for … 2 OE7:75 12:08, wherein the response of the system could be attracted by either the stable limit cycle around the origin or the deflected equilibrium. The unstable solutions to Eq. ( 15) (shown as dashed lines) separate the basin of attraction of the stable limit cycle around the origin from that about the deflected equilibrium. In other words, systems with initial conditions close to the deflected equilibrium are attracted by the fixed point.

Remarkably, for sufficiently large aspect ratios, the system dynamics is more complicated. In particular, for A D 0:35 (Figure 9(c)), the system displays a deformed-flapping motion for … 2 OE8:25 9:95. Comparing the A D 0:35 and A D 0:5 cases (shown in subfigures 9(c) and 9(d), respectively), it is seen that the deformedflapping regime narrows down and the amplitude of flapping grows faster with increasing aspect ratio. It should be realized, however, that for a true inverted flag with A > 0:2, the large amplitude flapping takes place abruptly, and the amplitude of oscillations varies substantially with …. To improve the results, a nonlinear analytical model accounting for the unsteady vortex shedding is called for, which is not a trivial task and virtually cannot be done in its entirety.

Moreover, the proposed model suggests that the critical flow velocity for the onset of large amplitude flapping approaches zero as A increases, as a result of an increase in the magnitude of the Coriolis force (refer to the discussion in Section 3.5). This signals that the validity of the slender body model deteriorates as the aspect ratio increases.

Influence of mass ratio

As evidenced in Eq. ( 15), the mass ratio parameter is associated with the added mass effect and velocity-related terms. This parameter plays no role in the static solution (fixed points) but makes a significant contribution to the stability of the solution branches. More specifically, ignoring time derivatives in Eq. ( 15) results in a set of algebraic equations yielding the static response of the system only. In the absence of time-derivative terms in Eq. ( 15) or in the presence of a heavy damping, this static solution is stable prior to a pitchfork bifurcation, and it regains stability after the folding point (saddle-node bifurcation). On the other hand, introducing time-derivative terms associated with unsteady fluid forces and performing a stability analysis leads to the detection of a Hopf bifurcation.

Figure 10 compares the theoretical results obtained by the proposed model with large damping and experimental data reported by Sader et al. [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF] with additional damping. The Hopf bifurcation does not occur, and hence the deflected equilibrium becomes stable, which agrees with the statement made just above. Interestingly, a fairly good qualitative and quantitative agreement can be seen in this figure. More specifically, Figure 10 shows a similar trend in the response calculated by the present model and the experimental results. Furthermore, it is seen that the subcritical behaviour of the response becomes weaker as A increases.

Quantitative agreement may be assessed by comparing the critical flow velocity for the saddle-node bifurcation (point F in Figure 4(a)). For instance, the present theory predicts … F D 6:6 for a slender inverted flag with A D 0:5, whereas from the experimental data in [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF], … F ' 6:7.

On the other hand, it is evident that there is a discrepancy between theory and experiment with regard to the deflection amplitude. This is likely caused by unsteady flow forces (arising from vortex shedding) leading to a change in the mean forces on the flag, which in turn leads to a different defected equilibrium. Agreement with experiment could well improve substantially if vortex shedding were included in the model.

Figure 11 shows the effect of on the nonlinear response of the inverted flag. The aim here is to investigate the relation between the length of the flag L and f , the critical flow velocity … c for flapping (which is also a function of L), as well as the amplitude of flapping, while the other physical parameters of the system ( f , p , h, and D) are fixed. As seen in Figure 11( amplitude oscillations. The critical velocity precipitating flapping is inversely proportional to the mass ratio. This endorses the linear stability analysis performed by Guo & Païdoussis [START_REF] Qin | Stability of rectangular plates with free side-edges in twodimensional inviscid channel flow[END_REF]. The transition from the limit-cycle in the vicinity of the origin to the one in the neighbourhood of the buckled state (the deformed-flapping regime) is also affected by the mass ratio: increasing delays the transition between the two regimes. For instance, the saddle-node bifurcation occurs at … D 7:7 for D 0:2, while it materializes at … D 8:9 for D 0:4. The stable limit-cycle around the buckled state (in the range … ' 8 to 9 for D 0:2 0:4) disappears for D 0:1. Interestingly, the bifurcated solution exhibits a subcritical behaviour for D 0:4, suggesting that there is hysteresis in the critical … obtained with decreasing and increasing flow velocity, as shown in the inset of Figure 11(a). The corresponding periodic solution loses stability via a saddle-node bifurcation at … D 1:78. The emerging unstable limit cycle folds at a second saddle-node bifurcation point … D 1:72 which corresponds to a restabilization of the periodic solution to one of slightly higher frequency and larger amplitude.

The mass ratio also affects the frequency of oscillation, as shown in Figure 11(b). For all cases, the flapping frequency decreases with the flow velocity. Furthermore, for a range of flow velocities associated with largeamplitude flapping around the undeflected equilibrium, the frequency of oscillation decreases with decreasing . Cases with 0:2 studied here display an abrupt change in the frequency with increasing …, corresponding to the change in the motion from large-amplitude motion around the origin to flapping around a deformed equilibrium (deformed-flapping regime).

Strengths and limitations of the proposed model

As with any theoretical study, the proposed model has both strengths and limitations. In this section, we address the similarities and differences in dynamical behaviour of small-aspect-ratio inverted flags as obtained via the proposed model and other theoretical/experimental studies in the literature, and we discuss the potential sources of anomalies as much as possible, putting the emphasis on physical understanding.

The measurements in [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF]Figure 8] for A D 0:5 show that the undeflected equilibrium is stable for … . 6. This can also be inferred from figure 5 of that study which gives the critical dimensionless flow velocities for the bifurcation of the zero-deflection equilibrium. It can also be seen therein that the large-amplitude flapping is the global attractor for flow rates beyond this critical flow velocity. In other words, the inverted flag with A D 0:5 undergoes large-amplitude flapping at circa … D 6. Furthermore, the sequence of "straight!flapping!fullydeflected" states for A 0:2 is also available from measurements [START_REF] Sader | Large-amplitude flapping of an inverted flag in a uniform steady flow -a vortex-induced vibration[END_REF].

In this sense, the present results are in fair agreement with the experimental observations on the global dynamical behaviour (i.e. the type of instabilities and their sequence) of low-aspect-ratio inverted flags. Additionally, there is agreement with the 3D simulations of Tang et al. [START_REF] Tang | Dynamics of an inverted flexible plate in a uniform flow[END_REF], predicting a similar qualitative behaviour to that in this paper for inverted flags, namely a strong sensitivity of the onset of flapping to the mass ratio parameter, signaling that the loss of stability of the trivial equilibrium is via a Hopf bifurcation.

While comparisons above provide partial agreement between the sequence of flow regimes shown in Figure 8 of the present paper with [24, Figure 3], as well as the computational study of Tang et al. [START_REF] Tang | Dynamics of an inverted flexible plate in a uniform flow[END_REF], we have noted a strong discrepancy between the structure of the bifurcation diagram predicted by the proposed model and that predicted via 2D simulations by Goza et al. [START_REF] Goza | Global modes and nonlinear analysis of inverted-flag flapping[END_REF] for infinitely wide inverted flags. More specifically, they identify a stable deformed equilibrium associated with a buckled configuration of the flag, emerging from a supercritical pitchfork bifurcation.

Moreover, the strong sensitivity shown to the mass ratio parameter in Figure 11 shows a rather severe disagreement with the experimental observations by Kim et al. [START_REF] Kim | Flapping dynamics of an inverted flag[END_REF], who reported no significant change in the critical flow velocity, while varying the mass ratio -although only two data points were obtained, 0:4 in air and 320 in water, and the mechanism of instability for the two systems may be of two totally distinct types, as shown by Sader et al. [START_REF] Sader | Large-amplitude flapping of an inverted flag in a uniform steady flow -a vortex-induced vibration[END_REF] and Goza et al. [START_REF] Goza | Global modes and nonlinear analysis of inverted-flag flapping[END_REF].

Based on the experimental observations of Kim et al. [START_REF] Kim | Flapping dynamics of an inverted flag[END_REF] on the influence of the mass ratio, Sader et al. [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF][START_REF] Sader | Large-amplitude flapping of an inverted flag in a uniform steady flow -a vortex-induced vibration[END_REF] employed steady-flow models and proposed a mathematical formula, connecting the critical velocity of the static saddle-node bifurcation calculated for a zero-aspect-ratio flag and that of the static divergence obtained for a flag of infinite aspect ratio. The model predicts that a flag with A D 0:5 undergoes a static deflection at … D 5:17, which is in good agreement with measurements. Sader et al. [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF] attributed this instability to the onset of flapping observed in measurements for A D 0:5 at circa … D 6. They explained that, at this point, the flag is deflected, presenting a sharp leading edge to the impinging flow, and thus the resulting unsteady vortex shedding causes the flag to oscillate via a vortex-induced vibration mechanism.

By virtue of this conclusion, the primary limitation of the model proposed in this paper may arise from the inherent absence of unsteady vortex shedding at large angles of attack, and the approximation employed to formulate the unsteady force arising from the vortex formation along the lateral edges of the flag.

One may employ a steady-flow model, omitting the time-dependent terms in Eq. ( 15) (ignoring the effect of the accelerated air as it passes over and under the flag). As a result, the bifurcation diagram in Figure 4(a) gives way to a static subcritical pitchfork bifurcation, with a stable static trivial solution, as shown in Figure 12. Following Sader et al. [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF], the onset of flapping can be attributed to the first instability (i.e. the static saddle-node bifurcation, marked by point F) at … D 6:6. Although a very good quantitative agreement can be found, no explanation can be proposed on how the undeflected equilibrium remains stable for velocities lower than point E; omitting unsteady forces leads to a subcritical pitchfork bifurcation for inverted flags of non-zero aspect ratio. Time-derivative terms, however, are essential in performing stability analyses on the branches of the obtained solution.

In this paper, it was found that taking into account the effect of unsteady forces -arising from the changing fluid flow momentum due to the oscillating inverted flag -causes the undeflected equilibrium to become linearly unstable. The negative damping in the system, induced by the Coriolis force, overcomes the stabilizing forces for the case of sufficiently large aspect ratio flags. The mechanism of this instability is found to be via a Hopf bifurcation with a critical velocity inversely proportional to . Our calculations show that this instability is very marginal and weak for very small aspect-ratio flags (the imaginary part of the complex eigenfrequency of the unstable mode is very small in magnitude and could easily become positive when damping is introduced), in line with observations by Sader et al. [START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF]. For larger aspect ratio flags, however, the Coriolis force becomes more dominant and becomes larger with increasing flow velocity. The first mode loses stability via a Hopf bifurcation, leading to large-amplitude limit cycle oscillations. At sufficiently high flow velocities, the vortex lift force, F R , which opposes the Coriolis force, grows in magnitude, and this leads the non-trivial static solution (deflected equilibrium) to regain stability.

The loss of stability via a Hopf bifurcation for inverted flags with A & 0:1 differs from the speculation of a divergence instability proposed in [START_REF] Kim | Flapping dynamics of an inverted flag[END_REF][START_REF] Sader | Large-amplitude flapping of an inverted flag in a uniform steady flow -a vortex-induced vibration[END_REF][START_REF] Sader | Stability of slender inverted flags and rods in uniform steady flow[END_REF] as outlined above. A reason for that could be the limitation of using an approximation of the dissipative effect of the surrounding fluid, instead of performing a full viscous flow calculation. The approximations used in the present study to model 3D effects may not be accurate enough to fully describe the underlying mechanism of the first instability, but the similarity of predictions with those of Tang et al. [START_REF] Tang | Dynamics of an inverted flexible plate in a uniform flow[END_REF] who used numerical methods to solve 3D Navier-Stokes equations leaves us with no definite conclusion. One possible improvement for analytical modelling of the problem, taking into account the 3D effects, may be to combine reactive and vortex-lift effects by extending Bollay's nonlinear lift theory to a plate with arbitrary curvature along its length, retaining the unsteady effects due to the shed vortices.

Further elucidation of 3D effects requires future studies.

Sensitivity to viscous drag

In this section, an attempt has been made to improve the fluid-forces approximated in the foregoing on viscous drag and to explore the effect on the stability and post-critical response of the system (refer to [START_REF] Tang | On the instability and the post-critical behaviour of two-dimensional cantilevered flexible plates in axial flow[END_REF]). The effect of the viscous drag k on the response is investigated numerically, and the results are shown in Figure 13. The system considered is a system with the mechanical parameters detailed in Figure 4, but with k varied. It should be emphasized that the viscous drag is assumed to be evenly distributed along the plate and is associated with a constant viscous drag coefficient k ; also, its contribution in the drag force (in the longitudinal direction) is considered only [START_REF] Tang | On the instability and the post-critical behaviour of two-dimensional cantilevered flexible plates in axial flow[END_REF].

As a first approximation for the viscous drag coefficient, the Blasius-type solution is used for both laminar and turbulent cases, that is

k D ( 0:054 Re 1 2
for a laminar flow; 1:328 Re for a turbulent flow; [START_REF] James | Large-amplitude elongated-body theory of fish locomotion[END_REF] where Re D UL= f denotes the Reynolds number, f being the fluid kinematic viscosity. For the calculation of the viscous drag modelled using Eq. ( 17), Re D 5 10 4 is considered. Furthermore, according to the measurements by Kim et al. [START_REF] Kim | Flapping dynamics of an inverted flag[END_REF], the value of the total drag (i.e., the summation of the viscous drag and the inviscid force component in the drag direction) is reported as being around 0:08 before flapping is initiated, and it jumps to around 1:0 at the critical point. The effect of this parameter is further investigated for a wider range here. stretched-straight state becomes unstable) is insignificant. However, with increasing k , subsequent bifurcations take place at lower flow velocities. In particular, the first Hopf bifurcation shifts slightly to the left in the bifurcation diagram, and the pitchfork bifurcation as well as the second Hopf bifurcation (where the deflected equilibrium becomes unstable) occur at smaller … when k is increased. Therefore, the drag coefficient significantly affects the interval of existence of the flapping regime. Furthermore, at a given value of … in the flapping regime, the amplitude of oscillation decreases with increasing k . The role of the viscous drag is more significant in the range of deformed-flapping regime. This range narrows down with larger k until it vanishes. Moreover, at a fixed flow velocity, the amplitude of the static equilibrium grows with k .

Sensitivity to angle of incidence

It is of interest to explore the dynamics of the system when the neutral plane of the inverted flag is not perfectly parallel to the direction of the undisturbed flow, as illustrated in Figure 14. Considering a small angle of incidence, the

A sensitivity analysis on the effect of the angle of incidence on the dynamics has been performed by considering four angles 0 . The resulting bifurcation diagrams are plotted in Figure 15. As seen in Figure 15(a), a small angle of incidence ( 0 D 0:1, corresponding to 5:7 ı ) breaks the symmetric response of the system, and the pitchfork bifurcation is replaced by a saddle-node bifurcation. The inverted flag exhibits an asymmetric flapping motion with a slight deviation from symmetric flapping. More interesting behaviour can be observed for larger values of 0 . Figure 15(b) shows the bifurcation diagram for 0 D 0:2. Four jumps associated with four saddle-node bifurcations are observed. More specifically, increasing the flow velocity, the amplitude of oscillation increases reaching to max t 2OE0 T . .L; t// D 7=25 at … D 4:73. At this point, the response of the system exhibits a jump and morphs into a large-amplitude limit-cycle. Further increasing the flow velocity leads to an additional jump at … D 6:55. A reverse scenario is expected when the system is initially flapping with a large amplitude and the flow velocity is decreased. The system then undergoes a jump at … D 2:76 and is attracted to the static equilibrium. The first two saddle-node bifurcations bracket a bistable region wherein the response of the system is attracted either by the stable limit cycles or the stable fixed points. The bistable region becomes narrower for larger values of 0 (for instance, for 0 D 0:3 in Figure 15(c)), and eventually is replaced by a gradual growth of oscillation amplitude (see Figure 15(d) for 0 D 0:4).

In Figure 15, a small 0 does not affect the large-amplitude asymmetric flapping regime much, signifying a relatively minor effect of small initial inclination on the amplitude oscillation. In particular, the peak-to-peak amplitude of the tip is 1:64 for 0 D 0:1 and 1:57 for 0 D 0:2, while it is 1:37 for 0 D 0:3 and 0:76 for 0 D 0:4. Thus, increasing 0 beyond 0 D 0:3 decreases the amplitude of asymmetric flapping significantly. In addition, the range of asymmetric flapping motions shrinks for inverted flags with larger angle of incidence, in line with the conclusions of Tang et al. [START_REF] Tang | Dynamics of an inverted flexible plate in a uniform flow[END_REF].

Figure 16 shows the deflected or flapping shapes of the inverted flag with 0 D 0:4 by superimposing views of the flag motion in one cycle. The trend shown agrees with previous experiments conducted by Cossé et al. [4]: flapping emerges only on one side of the inverted flag for sufficiently large angles of incidence, and the flapping amplitude decreases with 0 .

Concluding remarks

The nonlinear response of an inverted flag (thin elastic plate) was investigated in this paper. An analytical model was developed taking into account the geometric nonlinearities and the inertia of the flag. Fluid-related forces were treated independently by employing the extension of the large-amplitude elongated-body theory accounting for reactive forces acting on the flag, as well as a vortex lift mechanism based on Bollay's model, retaining unsteady (velocity-dependent) forces. The nonlinear partial integro-differential equation was first discretized via a Galerkin technique and then solved numerically using Gear's backward differentiation formula. Bifurcation diagrams were obtained through pseudo-arclength continuation. The results are presented as bifurcation diagrams, time histories, phase portraits, FFT plots, and flapping shapes, so as to explore the rich dynamical behaviour of the system. The overall nonlinear responses predicted by the analytical model exhibit several of qualitative similarities to those observed experimentally and reported in the literature. The responses exposed in this work indicate that the change in the flag aspect ratio substantially alters the dynamics. Although inverted flags of small aspect ratio do not exhibit flapping motion, inverted flags of sufficiently large aspect ratios do undergo a flapping motion. More specifically, in the case of slender inverted flags (i.e. flags of very small aspect ratio), the undeflected equilibrium loses stability via a subcritical pitchfork bifurcation giving rise to multiple additional equilibrium states. The unsteady component of the flow increases in strength with increasing aspect ratio. This results in quite intricate bifurcation mechanisms. It was seen that the undeflected equilibrium loses stability via a supercritical Hopf bifurcation in the first mode, the amplitude of which increases with the flow velocity. In particular, the inverted flag exhibits various dynamical regimes as the flow velocity is varied: (i) it initially resides at the undeflected equilibrium; (ii) it then exhibits a large-amplitude flapping motion around the undeflected equilibrium; (iii) this regime switches to a deformed-flapping motion; and finally, (iv) the system remains at the fully deflected static state with a highly curved shape.

The effect of various dimensionless parameters was investigated. It was shown that an initial angle of incidence affects the onset of flapping and the symmetry of the response. The stretched-straight state is replaced by a bent state, the amplitude of which subsequently increases with flow velocity. Increasing the angle of incidence results in a narrower range for the flapping regime, with smaller amplitude and stronger asymmetric oscillations.

It was found that flags with larger mass ratios display larger amplitude flapping and exhibit flapping at lower flow velocities.

It was also shown that the dynamical response is sensitive to viscous drag. Ignoring the viscous drag affects the response of the system considerably. In particular, it affects the range of flapping significantly and dampens the flapping motion around the deflected equilibrium.

Finally, in the present study, an attempt was made to explore the dominant physical mechanisms in the nonlinear dynamics of an inverted flag. Whereas there are a number of similarities in the results obtained here with measurements and results obtained by fully-coupled fluid-elastic simulations available in the literature, some inexplicable differences abide also. Differences were attributed to the limitations of inviscid flow theory. In particular, since the model is based on the assumption of attached flow, it is not expected to capture the aerodynamic forces beyond the angle where massive separation at the leading edge takes place. This phenomenon can cause a change in the mean forces on the flag, which in turn can lead to a different system response. The omission of unsteady fluid forces due to vortex shedding is also a strong limitation of this model. It can be said that the proposed analytical model is essentially useful for investigating the qualitative dynamics of inverted flags, and valuable for performing fast simulations to estimate the range of flow velocities wherein flapping takes place, but the results should be regarded with some caution considering the approximations employed. Improvements may be desirable in order to account for viscous effects and unsteady vortex shedding.

Figure 1 :

 1 Figure 1: (a) Inverted flag, idealized as a thin Euler-Bernoulli beam in axial flow, free at the upstream end and clamped at the downstream end, (b) Generic infinitesimal element of the beam.

Figure 3 :

 3 Figure 3: Variation of K p [ ] and K v [ ] with the aspect ratio of the flag, A.

Figure 4 :

 4 Figure 4: Bifurcation diagrams for D 0:3 and A D 0:5. (a) Maximum/minimum values of the tip slope of the inverted flag over a cycle of steady-state oscillation of period T . [ ] Stable static solution (undeflected stable equilibrium), [ ] unstable static solution, [ ] stable periodic solution (oscillation around undeflected equilibrium), [ ] stable periodic solution (oscillations around deflected equilibria), and [ ] stable static solution (deflected static equilibria). (b) maximum/minimum values of the tip transverse(black)/longitudinal(red) motion over a cycle of steady-state oscillation of period T . The colour scheme in (a) is used in Figs. 5-8.

Figure 5 :Figure 6 :

 56 Figure 5: Periodic response of the inverted flag with D 0:3 and A D 0:5 at … D 6:0 displaying a periodic motion around the trivial equilibrium; (a) time history, (b) phase-plane portrait, (c) and FFT plot indicating the frequency content of the response. Colour scheme from Figure 4(a).

Figure 7 :

 7 Figure 7: Dimensionless frequency f versus dimensionless flow velocity … for an inverted flag with D 0:3 and A D 0:5. Colour scheme from Figure 4(a).

Figure 8 :

 8 Figure 8: Dynamics obtained by superimposed deformations of an inverted flag with D 0:3 and A D 0:5 at different flow velocities displaying (a) straight (undeflected) static equilibrium at … D 0:80, (b) large-amplitude flapping around the undeflected equilibrium at … D 6:00, (c) flapping around the deflected equilibrium at … D 9:00, and (d) fully-deflected state at … D 20:00. The colour scheme is defined in Figure 4(a).

Figure 9 :

 9 Figure 9: Bifurcation diagrams for inverted flags with D 0:3 and various aspect ratios: (a) A D 0:1 with k p D 0:22 and k v D 0:93, (b) A D 0:25 with k p D 0:41 and k v D 0:93 (c) A D 0:35 with k p D 0:55 and k v D 0:94, and (d) A D 0:5 with k p D 0:75 and k v D 0:89. Maximum/minimum values of the tip slope of the inverted flag over a cycle of steady-state oscillation of period T . [ ] Stable static solution (undeflected stable equilibria), [ ] unstable static solution, [ ] stable periodic solution (oscillation around undeflected equilibria), [ ] stable periodic solution (oscillations around deflected equilibria), and [ ] stable static solution (deflected static equilibria).

Figure 10 :

 10 Figure 10: Comparison between static results obtained by the proposed model and experimental data reported in [25] for various A: [ ] A D 1:0, [ ] A D 0:5 [ ] A D 0:13, and [ ] A D 0:067; [ , ] stable and [ , ] unstable deflected equilibria.

Figure 11 :

 11 Figure11shows the effect of on the nonlinear response of the inverted flag. The aim here is to investigate the relation between the length of the flag L and f , the critical flow velocity … c for flapping (which is also a function of L), as well as the amplitude of flapping, while the other physical parameters of the system ( f , p , h, and D) are fixed. As seen in Figure11(a), a larger has the effect of generating instability at lower flow velocities, with larger

FFigure 12 :

 12 Figure 12: Bifurcation diagram of the inverted flag with D 0:0 and A D 0:5. [ ] Stable static solution (undeflected stable equilibria), [ ] unstable static solution.

Figure 13 depictsFigure 13 :

 1313 Figure 13: Bifurcation diagram of the inverted flag with D 0:3 and A D 0:5 and various k .

Figure 14 :

 14 Figure 14: Inverted flag in axial flow with a small angle of incidence 0

Figure 15 :Figure 16 :

 1516 Figure 15: Bifurcation diagrams for inverted flags with D 0:3, A D 0:5, and various initial angles of incidence: (a) 0 D 0:1, (b) 0 D 0:2, (c) 0 D 0:3, and (d) 0 D 0:4. Maximum/minimum values of the tip slope of the inverted flag over a cycle of steady-state oscillation of period T . [ ] Stable static solution, [ ] unstable static solution and [ ] stable periodic solution (asymmetric oscillations).

  and c d D c d L 4 =.D /, in which f is the frequency of flapping and D p p A=DL 2 , yields the following dimensionless nonlinear equation of motion (in which the asterisk notation has been dropped for simplicity): C 2 P w 2 . Also, it is recalled that the boundary conditions for the clamped-free flag read .0; t/ D 0 and x .L; t/ D 0. The discretization of the system in space is obtained via the Galerkin technique. Hence, as a comparison functions, a family of eigenfunctions r , r D 1; : : : ; M , satisfying the geometric boundary conditions of the problem, is selected to approximate the deflection pattern of the plate [21, Section 2.1.4]. The rotation angle .x; t/ is then expanded as

			Z 1	Z s	Z 1	Z s
	ˇR sin @ tt			cos .; t/ d ds	C cos @ tt	sin .; t/ d ds
					x	0	x	0
						Z 1	Z s	Z 1	Z s
	@ xx C c d		cos @ t	sin .; t/ d ds	sin @ t	cos .; t/ d ds
						x	0	x	0
						Z x
		s @ txx C …	cos	C f .s; t/ cos .s; t/ C .s; t/ sin .s; t/	ds
						1
				Z x
		C sin				C f .s; t/ sin .s; t/ C C .s; t/ cos .s; t/	ds	D 0:	(15)
					1
	Using the notation D	p	=…,
	C f D	A 4	R w cos	R u sin C 2 P . P u cos C P w sin /
		C @ x .3 cos 2 . P w 2 P u 2 1/ C 1 C P u 2 6 P u P w cos sin	2 P w 2 /=2
			C @ x . P u 3 P u cos 2	3 P w sin cos / C 2 P cos
		k p .A/ B 2	.. P u C 1/ 2 P w 2 / sin cos C . P w P u C P w/.sin 2	cos 2 /
		k v .A/ B 2	. P u C 1/ sin 2 C 2 P w 2 cos 2	2. P w P u C P w/ sin cos
	and
	C D	k .Re/ B	h .1 C P u/ cos C P w sin	i ;
	in which B D u/ 2 .x; t/ D p .1 C P

M X rD1 r .x/q r .t/ with r .x/ D sin. r x/Csinh. r x/C cos. r / C cosh. r / sin. r / C sinh. r / .cos. r x/ cosh. r x// (16)

In this paper, including the Introduction where others' work is also discussed, the mass ratio is as defined in the nomenclature (fluid/structure); in some other papers, the inverse definition (structure/fluid) is utilized.

Gersten assumed that the flow trails off the surface of the plate at an angle equal to half of the plate's angle of attack. Hence, the expression for the normal force takes on an explicit form in terms of ˛.

This model of damping in quiescent fluid is an approximation taking into account the dissipative effect due to the surrounding fluid. It has been employed in work involving cylinders oscillating in fluid[START_REF] Michael | Dynamics of a long tubular cantilever conveying fluid downwards, which then flows upwards around the cantilever as a confined annular flow[END_REF][START_REF] Abdelbaki | A nonlinear model for a free-clamped cylinder subjected to confined axial flow[END_REF]. A more comprehensive approach to model the viscous forces is through a full viscous flow calculation using computational techniques, a very challenging task beyond the purposes of this paper.

If the limit cycle is unstable or repelling, the Hopf bifurcation is of subcritical type; otherwise, the bifurcation is supercritical [21, Section 2.3].

Acknowledgments

The authors gratefully acknowledge the financial support by the Natural Sciences and Engineering Research Council of Canada, the Solution Mining Research Institute (SMRI) and Pipeline Research Council International (PRCI).


Christophe Marquesuzaà

Patrick Etcheverry

Pantxika Dagorret

Philippe Lopistéguy

Thierry Nodenot

Marta Toribio Fontenla

Models, Process and Tool to Assist Cooperative Scenarization of Distant Learning Modules

Keywords: Pedagogical scenarization, scenario models, scenarization platform, Moodle

published or not. The documents may come   L'archive ouverte pluridisciplinaire

INTRODUCTION

Digital competences are becoming increasingly important, whether in the daily lives of the general public, in companies or in the educational field from elementary school to University.

For example, the private french Orange telecommunications company aims to train 50% of its 100 000 employees by using MOOC platforms.On March 5, 2018, Muriel Pénicaud, the French Minister of Labour, unveiled the outline of her professional training reform project to facilitate access to any kind of training programs. In this reform, the use of MOOC platforms is promoted and become eligible by employees as part of their Training Personal Account (CPF in French). The designers of these MOOCs would have to provide high-quality educational pathways that would be assessed by a National Agency.

Within the university sector, lecturers have recognized skills and long-standing experience in creating face-to-face courses. Many educational institutions now offer online workspaces allowing at least access to digital resources in the course of traditional classroom (teaching/learning sessions). These resources can be used before classes, for example in the context of "flipped classrooms", during the class to access nonphotocopied documents or after the class, for revision. The hybrid form (face-to-face + online support) aims at making people work better together in order to smooth the training time, not to contract it. Generally, these pedagogical resources have strong dependences with the classroom teaching/learning session in which they are used and they often lose their relevance when one tries to exploit them on learning platforms such as MOOCs where classroom teaching/learning sessions disappear.

If lecturers have no difficulties in elaborating educational pathways for classroom teaching/learning session, their experience is currently quite limited when it is about designing courses that are fully online and not including any classroom session. In this context, they need the support of pedagogical engineers to see to what extent existing or new courses can be adapted to suit e-learner's requirements.

While [START_REF] Bakki | A Model to Assist Pedagogical Scenario Building Process in cMOOCs[END_REF] proposes a model to assist pedagogical scenario building process for MOOCs in a connectivist approach, in this paper we propose a methodology and an associated online platform allowing pedagogical engineers to assist lecturers all along the process aiming at designing online learning modules in a constructivist approach.

The expected spinoffs are numerous. From the learner's point of view, it means improving the learning modes and their arrangement and enhancing their motivation by diversifying available teaching methods. As regards lecturers, the aim is to expand their teaching practice by diversifying the ways to share their knowledge. As to the educational institutions, the main interest is to be equipped with a methodological framework facilitating the development of online courses and also to participate actively in the process of raising the level of knowledge and skills of all the actors involved, from the design to the implementation of an education offer.

This paper is structured as follows. In section 2, we present a feedback on a scenarization process used to set up an online degree course. Considering this experience, in section 3 we present related works dealing with pedagogical scenarios modelling but also design processes aiming at elaborating pedagogical scenarios. Section 4 presents several models to deal with both (1) educational scenarios, (2) scenarization methodologies. In section 5, we describe an online environment implementing the previous models. This platform allows lecturers and pedagogical engineers to design, remotely and cooperatively, online learning modules. Finally, we conclude this paper by presenting several ideas concerning potential future work.

FEEDBACK CONCERNING THE IMPLEMENTATION OF A REAL DESIGN PROCESS

We conducted a design work aiming at creating a university degree that can be obtained through an online degree course. This university degree, called DU ASR in French, focuses on systems and network administration. The target audience consists of adult learners living in Latin America and wishing to acquire knowledge and skills related to networking and operating systems in order to work as high-level technicians.

Trainers involved in the design of the pedagogical modules forming this university degree come from various higher education institutions from France (IUT de Bayonne et du Pays Basque), Spain (Universidad del País Vasco / Euskal Herriko Unibertsitatea UPV/EHU), Argentina (ITU Mendoza), Chile (CEDUC Universidad Católica del Norte) and Colombia (Universidad Tecnológica de Pereira).

The design process carried out to elaborate all these online modules has been developed by a pedagogical engineer who worked in cooperation with 18 trainers. Each elaborated modules corresponding to 30 or 60 hours of learning/training that learners can follow remotely and in total autonomy.

For each pedagogical module, the scenarization process is composed of three main steps: 1. Contacting concerned trainers in order to identify the objective of the module 2. Establishing a synthetic description of the chapters composing the module 3. Designing the scenario describing accurately the contents of each chapter. The first step takes the form of a questionnaire composed of 50 questions to the course author. The aim of this questionnaire is to help the pedagogical engineer get acquainted with the author and identify some of his/her pedagogical practices. This questionnaire also allows the pedagogical engineer to identify the level of the trainer's skills and personal practice related to e-teaching/learning.

The second step deals with the general structure of the course. The pedagogical engineer provides the trainer with a commented and illustrated form aiming at formalizing a full specification of the module to bring online. This form allows the trainer to specify the target audiences, the necessary prerequisites, knowledge and skills to be acquired, equipment and software needed, the planned assessment method and a first list of pedagogical sequences (chapters).

Basing on this list of pedagogical sequences, the third step of the design process details the pedagogical scenario of the course. It consists of:

• Dividing each chapter into activities while checking that each activity objective contributes to the general objective of the chapter; • Defining the pedagogical resources which must be mediatized for each activity; • Defining when and how each resource must be produced. An extract of the documents produced at each step can be obtained at the following URL: https://goo.gl/nYnzay. This scenarization task takes the form of an iterative specification process resulting in several exchanges in which the pedagogical engineer brings his/her expertise and his/her recommendations in order to assist lecturers elaborating a qualitative pedagogical scenario adapted to distance learning to be implemented on the Moodle platform. This platform has been chosen for two main reasons:

• it supports a social constructivism design approach which favours rich interactions between pedagogical engineers and lecturers ; • it provides adaptivity to address users requirements when designing e-education courses (Despotovic-Zrakic, M. et al, 2012) ; • it is based on open source code and may be easily customized.

The scenarization process carried out to elaborate the DU ASR learning modules presents undeniable qualities which allowed the pedagogical engineer to remotely design each module with lecturers located in several countries and with different practices. The first learning modules were put online in September 2015 and since March 2017, the 11 modules composing the whole diploma course are available online (https://elearn.univ-pau.fr/course/index.php?categoryid=6).

Despite these satisfactory results, the work carried out to design this online degree course has allowed us to characterize several improvable points:

• The forms currently used at each design stage (Word and Excel files sent by e-mail) are ill-adapted to a cooperative scenarization task which requires several exchanges between lecturers and pedagogical engineers but also sometimes between lecturers themselves when several lecturers are involved in a same learning module). All these exchanges imply to manage the different versions of each document and it is sometimes tiresome to produce a document satisfying at the same time the pedagogical engineers and all the lecturers implied in the design process. • The scenarization design process is currently formalized through document models which offer little flexibility when it becomes necessary to adapt some stages/questions to take into account the specific features of some modules or the profile and/or language of lecturers implied in the scenarization process. • The online publication of a new learning module is tiresome and time-consuming because it is carried out manually and it implies to gather and to take into account all pedagogical elements specified in the documents describing the scenario.

From these observations, we have elaborated models and a platform support:

• allowing pedagogical engineers to create, adapt and make scenarization processes evolve;

• facilitating cooperative work between pedagogical engineers and lecturers;

• automating the export of a pedagogical scenario towards the Moodle LMS (Learning System Management) to facilitate the publication of the corresponding course.

These different contributions are presented in the following sections

RELATED WORK

As defined in [START_REF] Schneider | Conception and implementation of rich pedagogical scenarios through collaborative portal sites: clear focus and fuzzy edges[END_REF] and [START_REF] Peter | Pedagogical scenario modelling, deployment, execution and evolution[END_REF] a pedagogical scenario is an ordered set of learning activities implying actors who use and produce resources (or "learning objects"). A pedagogical scenario can take several forms (text, audio or narrative video) and aims at describing the sequences, and sometimes the objectives, the actors, the stages, even the instructions, the tools and documents used or to be produced. These concepts are also described in (Alario-Hoyos, C. and Muñoz Cristóbal, J., 2012) and [START_REF] Martel | Modeling collaborative learning activities on e-learning platforms[END_REF].

There are two standards to describe a scenario where each scenario is a pedagogical "design" which can be recorded into an XML file:

• The SCORM (Sharable Content Object Reference Model) standard (http://www. scorm.fr) is a compromise between several proposals aiming at standardizing the contents aggregation and data exchanges between a course and a compatible LMS. However, this standard presents various drawbacks concerning learner tracking or the choice of the teaching methods (mono-actor, simple sequencing) • The IMS-LD (Instructional Management Systems Learning Design) (www. imsglobal.org/learningdesign) is a generalization of SCORM with a multiactor model. This model is more flexible and supports an approach which focuses on activities and collaborations rather than on a succession of contents. However, this standard is relatively technical and also has some limitations concerning the management of actors and their production, informal definition of learning goals or prerequisites). So, the most complete models are not necessarily the best and they must be customized to pedagogical practices and needs while making them accessible to the users.

These pedagogical scenario models are integrated into scenarization tools like SCENARI -OPALE [START_REF] Crozat | Scenari -La chaîne éditoriale libre : Structurer et publier textes, images et son[END_REF], LAMS [START_REF] Dalziel | Proceedings of the 9 th International LAMS and Learning Design Conference[END_REF] or G-MOT/TELOS (http://lice.licef.ca/) but also into an LMS such as Moodle.

As part of the experiment conducted for the DU ASR, and because the diploma had to be available on a Moodle platform, the pedagogical engineer scripted the pedagogical scenario of each learning module on the base of his/her own knowledge of the underlying scenario of Moodle (Loiseau, E. et al, 2014). Moodle's built in scenario model operationalizes a scenario with entities like courses, composed of sections in which learners will find resources (files, directories…) and tools (ambiguously referred to as "activities") such as MCQs, forums, wikis, etc.

From a methodological point of view, the pedagogical engineer elaborated each learning module of the DU ASR degree in five steps, according to the ADDIE model (Analysis, Design, Development, Implementation and Evaluation).

ADDIE model has provided the basis for many methods of pedagogical design (Gagne, R.M., 1992), because it identifies the life cycle phases of a learning system according to a classic approach used in software design methods [START_REF] Serhat | ADDIE Model: Instructional Design, web site[END_REF]. ADDIE model leads to a directive design approach. It differs to SAM, the Successive Approximation Models [START_REF] Allen | Leaving ADDIE for SAM: An Agile Model for Developing the Best Learning Experiences[END_REF], which proposes an iterative design approach close to agile methods and relying on three stages: Evaluation, Design and Development, with a set of iterations on the whole process. Iterations can lead to the reconsideration of each result and may lead to carry out the evaluation step at each milestone.

The particular context of DU ASR project (various institutions and interlocutors, different native languages, different pedagogical approaches) led the pedagogical engineer to choose ADDIE model and its structuring character where each step can be repeated as often as necessary to improve the obtained result but necessarily ends with a validation before taking the next step.

Based on the experience gained on the DU ASR degree, our goal is to design models for developing educational scenarios garanteeing the operationalization of these scenarios on Moodle. We distinguish two types of approaches to tackle this problem:

• Design approaches with ad hoc tools allowing pedagogical engineers to elaborate pedagogical scenarios whose specification must, afterwards, be matched with data expected by Moodle (Alario-Hoyos, C. and Muñoz Cristóbal, J., 2012), [START_REF] Gagne | Principles of Instructional Design[END_REF]. Considering the important dissimilarity between the scenarization language and the data model of the LMS, these approaches often entail modifications and semantic losses during the mapping process. • Design approaches based on IDM to transform a given scenario according to a given model towards a specification compliant with the Moodle model. However, these transformations are complex and often induce a manual intervention and semantic losses.

To avoid these semantic losses, we decided as in (Loiseau, E. et al, 2014), to take Moodle's scenario model directly as a basis and to offer pedagogical engineers facilities to elaborate tailor-made design processes aiming at specifying scenarios in accordance with this model.

The pedagogical scenario is led by the teachers' intentions and in accordance with their values [START_REF] Emin | Model and tool to clarify intentions and strategies in learning scenarios design[END_REF]. It is therefore important to be able to adapt the scenarisation approach in order to take into account teachers pedagogical requirements. To this end, we propose a design approach model that is flexible and can be adapted by the pedagogical engineer according to the needs of the teachers with whom he/she cooperates.

We expect that the elaborated design processes be both rigorous and structuring like ADDIE but also flexible like SAM to allow pedagogical engineers to adapt their manner to assist lecturers during the design process as suggested in [START_REF] Peter | Pedagogical scenario modelling, deployment, execution and evolution[END_REF]. With this purpose, we propose a design process model allowing pedagogical engineers to elaborate design processes producing pedagogical scenario descriptions compliant with the scenario model integrated in Moodle.

PEDAGOGICAL SCENARIZATION MODELS

The scenarization models we propose include two aspects: modelling the description of a pedagogical scenario in terms of components, and modelling the scenarization approach, it means the steps followed by the pedagogical engineer when building a new course.

A scenario model

The definition of the elements making up a pedagogical scenario has been achieved in tight cooperation with the pedagogical engineer, who settled a vocabulary derived from the experience she acquired during the DU ASR project and her knowledge of the Moodle platform where the scenario had to be implemented. The model includes the following terminology (Figure 1): a pedagogical scenario (or course module) is composed of a set of chapters (or sections / pedagogical sequences). Each chapter may be composed of sessions where learners use resources (PDF files, videos, MCQs, forums, chats…) to carry out a given activity.

Then, in order to allow the exportation of all scenarios items defined by the pedagogical engineer and lecturers to a Moodle platform, we built a correlation table between this terminology and the concepts which are specific to the Moodle platform (Table 1).

Pedagogical engineer vocabulary

Corresponding Moodle concepts Module -Pedagogical Scenario Course Chapter -Pedagogical Sequence Section Session Activity -Resource Table 1. Matching pedagogical items with Moodle artefacts

A Scenarization approach model

In order to facilitate the capitalization and reuse of the formerly described pedagogical engineer's scenarization approach, we first formalized it. Then, we generalized the concept of scenarization approach to allow the pedagogical engineers to customize their scenarization approach, or to develop their own new approach(es).

According to our point of view, a scenarization approach has to be developed by a pedagogical engineer. The application of an approach is cooperatively implemented between one or more lecturers and one or more pedagogical engineers assisting them in the scenario specification task. The experiment conducted as part of the definition of the modules on the DU ASR showed that it was important to appoint a lecturer as a supervisor when several lecturers participate in the scenarization of a same module. As a supervisor, this person is the privileged interlocutor of the pedagogical engineer for the synchronisation and coordination of the specification task led by the teaching staff.

Hence, the term "pedagogical engineer", as used hereafter may refer to a single or more person coordinating the specification task. Similarly, the word "lecturer" may refer to a unique lecturer or the lecturer responsible for the coordination of a teaching staff.

Our scenarization approach model is presented in the right hand of Figure 2. In this model, a scenarization approach is composed of a set of ordered steps, each intented to specify part of a scenario. The specification objective of each step is defined by the pedagogical engineer. He/She also, depending on the step of the scripting work, defines whether this step is currently accessible or not to the lecturer. As suggested in ADDIE, this property allows the pedagogical engineer to manage the progress of the specification task by giving access to a new step only if the specification work of the previous one is considered satisfactory and completed.

A scenarization step is composed of a set of questions whose content and order are defined by the pedagogical engineer. These questions will be asked to the lecturer in order to let him/her specify the properties of his/her pedagogical scenario. Since we considered that a scenario (module) was composed of chapters which are composed of sessions (Table 1), the questions defined by the pedagogical engineer will have to lead the lecturer to specify these items. • General questions are not directly related to the scenario (module) to develop. These may be questions that allow the pedagogical engineer to better understand the lecturer with whom he/she is working (scenarization abilities, ICTE skills…) or questions that will prepare the lecturer for the identification of chapters and sessions composing his/her scenario (module); • Sequence-type questions aim at precisely specifying the chapters composing his/her scenario (module); • Session-type questions aim at defining the sessions composing each chapter. As shown in Figure 2, the answers given to sequence-type and session-type questions match the specifications of chapters and sessions of the scenario (module) which will later be implemented on the educational platform.

Our study's final model gathers the two former models and some additional classes, like Lecturer and PedagogicalEngineer, and relationships between them. The lecturer is responsible for coordinating the proposals of his/her colleagues when the module is designed by several lecturers. It is he/she who is in direct relationship with the pedagogical engineer.

A PLATFORM FOR COOPERATIVE DESIGN OF DISTANCE OR MIXED LEARNING MODULES

To evaluate the interest and relevance of our approach, we have integrated our models into a cooperative scenarization platform in order to determine their relevance during the design of new modules. The platform allows pedagogical engineers to design their own scenarization approach and implement it while cooperating with lecturers. The design of scenarization approaches is flexible and may be led according to different ways:

• Top-down: firstly defining the goals of the approach, next producing the steps composing the approach and lastly specifying the questions to be asked to the lecturers at each step; • Bottom-up: firstly defining the relevant questions to specify a scenario (module), next sharing these questions according to each step; • Mixed: creating for example a first step and then defining the associated questions, next creating a second step and its associated questions, and so forth. To support this flexible designing approach, the scenarization platform allows the pedagogical engineers to define in any order:

• The questions to add to a bank of questions and which can later be associated with a step.

• The steps to add to a bank of steps and which can later be used to design a specific approach step.

If the scenarization approach is top-down, then the pedagogical engineer will:

CONCLUSION AND FUTURE WORK

The design models implemented on our scenarization platform are currently being evaluated. The evaluation focuses on the relevance (and the limitations) of the models allowing to design new scenarization approaches or to customize existing approaches.

The above scenario model and the scenarization approaches model are both a first stage allowing pedagogical engineers and lecturers to design the frame of a scenario/module (sessions and chapters) and to export it to the educational platform. The proposed models do not currently allow to specify/create the educational resources involved in each learning module. These resources are therefore currently defined directly on Moodle and not on the scenarization platform. Future works therefore aim at extending the scenario model to support the specification of the educational resources involved in each educational activity. This also implies to extend the scenarization approaches model in order to specify this dimension of the scenario.

In order to cover a wider range of educational practices and to streamline the progression of the activities to be conducted to create the learning modules, the general model should be extended to also take into account the pedagogical roles which are commonly identified in the literature. We may here cite for example the roles identified in [START_REF] Peter | Platform support for pedagogical scenarios[END_REF]: module manager, module author, scriptwriter, media creator, tutor, administrative manager… More broadly, it is necessary to define more precisely the type of scenarization approaches that can be designed with our models to identify the teaching practices that are not covered and make our models evolve accordingly.

Figure

  Figure 1. Scenario model

Figure 2 .

 2 Figure 2. General model Our model distinguishes three types of questions:• General questions are not directly related to the scenario (module) to develop. These may be questions that allow the pedagogical engineer to better understand the lecturer with whom he/she is working (scenarization abilities, ICTE skills…) or questions that will prepare the lecturer for the identification of chapters and sessions composing his/her scenario (module); • Sequence-type questions aim at precisely specifying the chapters composing his/her scenario (module); • Session-type questions aim at defining the sessions composing each chapter. As shown in Figure2, the answers given to sequence-type and session-type questions match the specifications of chapters and sessions of the scenario (module) which will later be implemented on the educational platform.Our study's final model gathers the two former models and some additional classes, like Lecturer and PedagogicalEngineer, and relationships between them. The lecturer is responsible for coordinating the proposals of his/her colleagues when the module is designed by several lecturers. It is he/she who is in direct relationship with the pedagogical engineer.