Investigation of Fluid-Structure Interaction Induced Bending for Elastic Flaps in a Cross Flow

: With the recent increase in the design of light and ﬂexible structures, numerical investigations of ﬂuid and structure together play a signiﬁcant role in most engineering applications. Therefore, the current study presents an examination of ﬂuid-structure interaction involving ﬂexible structures. The problem is numerically solved by a commercial software ANSYS-Workbench. Two-way coupled three-dimensional transient simulations are carried out for the ﬂexible ﬂaps of di ﬀ erent thicknesses in glycerin for a laminar ﬂow and Reynolds number ranging from 3 < Re < 12. The bending line of the ﬂaps is compared with experimental data for di ﬀ erent alignments of the ﬂaps relative to the ﬂuid ﬂow. The study reports the computation of the maximum tip-deﬂection and deformation of ﬂaps ﬁxed at the bottom and mounted normal to the ﬂow. Additionally, drag coe ﬃ cients for ﬂexible ﬂaps are computed and ﬂow regimes in the wake of the ﬂaps are presented. As well, the study gives an understanding on how the ﬂuid response changes as the structure deforms and the model is appropriate to predict the behavior of thick and comparatively thinner ﬂaps. The results are su ﬃ ciently encouraging to consider the present model for analyzing turbulent ﬂow processes against ﬂexible objects.


Introduction
Fluid-structure interaction (FSI) states the mutual interaction of a structure with a flowing fluid and develops the fluid and solid motions and the dependency of these motions on each other. Such kind of problems involve, the velocity and pressure of fluid effects and are affected by the structure deformation. In recent years, it has become a prominent characteristic of various engineering projects because design involves lighter and more flexible structures. The conduct of these structures with a fluid generates multi physical situations which are usually complex to handle. Numerical calculations play an important role in most engineering areas and have become more challenging due to the prediction of different flow regimes which cannot be seen in experiments (EXP) but it is still a challenge to have accurate numerical simulation of a FSI problem because of strong fluid-structure coupling. FSI problems are ever present in nature for example flapping wings, flags [1][2][3][4][5][6][7][8][9], aquatic animals [10], flexible vegetation [11][12][13], in hemodynamics [1,14] and human hair and hair sensors [15]. Thin structures are of great interest for microfluidic procedures in engineering and biomedical applications [16][17][18] and structural dynamics [18][19][20].
For investigating fluid motion together with the structure, extensive theoretical, experimental and numerical studies across flexible plates are carried out in the context of plate flapping, vibration, presented [33] can be used as a benchmark for the formulation of numerical problems involving 2D FSI with immersed structures that have large displacements. While 2D benchmark cases, free-surface flows and test cases of an elastic obstacle and plate have been presented for monolithic approaches and validated with other numerical methods [34,35]. Similarly, for the analysis of hemodynamic problems, a new and more sophisticated computational formulation named immersed structural potential method (ISPM), based on the original IBM, is presented. The novel method accounts for the computation of FSI force field at the fluid mesh and description of cardiovascular tissue scheme by viscoelastic fiber-reinforced constitutive models. A new time-integration method is also introduced for the computation of a deformation gradient tensor. Three advancement cases of deformable cylinders and flapping membranes are presented. The study is applicable for highly complex FSI problems incorporating large deformations under pulsating flows and sinking and floating rigid bodies [36]. To describe the deflection of flaps and moving boundaries, the fluid-solid interface-tracking/interface-capturing technique (FSITICT) is applied to the FSI problem with ALE interface tracking and Eulerian interface tracking. The method is demonstrated by test FSI cases. First, test computations are performed for two flapping membranes and the second test case is composed of one flapping membrane and a movable, elastic arterial-structure. The flapping membranes and fluid are modelled in fully Eulerian coordinates in both test cases, while the elastic structure in the second case is computed in Lagrangian coordinates. The overall objective is to reduce the computational cost and test the applicability of the method for large deformations in order to analyze different cardiac cycles [37].
In the light of the above survey, it is concluded that the previous research efforts are made to encounter the dynamics of flapping structures in uniform flows. Prediction of the bending line and tip-displacement of a flexible flap fixed at the bottom in 3D is rarely seen in the literature. With the simple geometry of a flap in a rectangular fluid domain, the problem is fundamental in nature and can be used in applications that coincide with the flow, forcing functions and the solid boundary conditions. Therefore, this study emphasizes on the deformation characteristics of the flaps over time, flap bending positions and different flap orientations. The focus of this study is to explore the motion of an elastic flap immersed in a rectangular fluid domain and to analyze the bending of flexible structures in 3D under laminar flow conditions. This paper is organized as follows-first, an insight into experimental procedure is presented then a computational model is demonstrated in terms of different flap thicknesses and positions. Later, the response of flexible flaps is analyzed in terms of bending and deformation in addition, the results of flap bending are offered in comparison with the experimental work. The effect of the different clamping angles is discussed as well. Similarly, the computation of drag coefficient at various clamping positions is a part of the current research. Flow regimes are revealed at various time steps and velocity contours in the wake are shown. Finally, some discussions are presented.

Material and Method
A numerical two-way FSI model is presented for validating the towing-tank experimental data for flexible flaps. For the validation purposes, transient structure (flap) and unsteady flow are modeled using a commercial software ANSYS-Workbench. Each field is demonstrated separately by defining its geometry, properties, mesh, boundary conditions and respective solvers. Coupling is introduced between individual domains to establish a two-way FSI.
The experimental and numerical setups, settings and conditions are presented in the following subsets.

Experiments
A summary of the experimental data is given in the present section, whereas detail can be found in the reference [32]. Flaps of different thicknesses are fixed in the center of a plate and immersed into the working fluid.
Flaps are designed with length of 100 mm, width of 20 mm and thicknesses of 5 and 10 mm. Several measurements for the deformable flaps are taken by varying the flow conditions (Reynolds number, Re = 3 to 12) and flap arrangements that is, at clamping angles from 0 • to 90 • and inclined angle of 45 • . The bending of the flaps was recorded for each case. As the effective Reynolds number is small and viscous forces are dominating, only bending response is observed in the experiments instead of vortex shedding, which would cause periodic vibrations.

Numerical Model
FSI modeling mainly deals with the information about how the fluid domain acts as a function of structure and how structural domain behaves as a function of fluid domain, at the interface, information is shared between the structural and fluid solver thus founding FSI [38]. In the present case, an elastic structure was surrounded by incompressible Newtonian fluid flow. Ω f (t) is defined as a time dependent fluid domain and Ω s (t) is defined as a time dependent structural domain. It is assumed that for all the time t, Ω f (t) is entirely occupied by the fluid and Ω s (t) is entirely occupied by the structure. Besides, is an interface where the elastic structure interacts with the fluid. A schematic of FSI in terms of structural domain, fluid domain and interface for a present case of an elastic flap immersed in a rectangular fluid domain is shown in Figure 1. Flaps are designed with length of 100 mm, width of 20 mm and thicknesses of 5 and 10 mm. Several measurements for the deformable flaps are taken by varying the flow conditions (Reynolds number, Re = 3 to 12) and flap arrangements that is, at clamping angles from 0 o to 90 o and inclined angle of 45 o . The bending of the flaps was recorded for each case. As the effective Reynolds number is small and viscous forces are dominating, only bending response is observed in the experiments instead of vortex shedding, which would cause periodic vibrations.

Numerical Model
FSI modeling mainly deals with the information about how the fluid domain acts as a function of structure and how structural domain behaves as a function of fluid domain, at the interface, information is shared between the structural and fluid solver thus founding FSI [38]. In the present case, an elastic structure was surrounded by incompressible Newtonian fluid flow. Ωf(t) is defined as a time dependent fluid domain and Ωs(t) is defined as a time dependent structural domain. It is assumed that for all the time t, Ωf(t) is entirely occupied by the fluid and Ωs(t) is entirely occupied by the structure. Besides, Γo(t) = Ωf(t) ∩ Ωs(t) is an interface where the elastic structure interacts with the fluid. A schematic of FSI in terms of structural domain, fluid domain and interface for a present case of an elastic flap immersed in a rectangular fluid domain is shown in Figure 1. For the established FSI solution, the kinematic, dynamic and geometric conditions at the common interface of fluid and structure should be satisfied, numerically at the interface Γo(t), the following conditions are assumed for all the time t: Here, and represent the structure and fluid velocity constituting the kinematic stability.
Dynamic stability is attained by the continuous normal stresses of solid and fluid at the interface. being the unit vector normal at the interface and the subscript i indicates its direction.
According to the geometric condition stability fluid-solid domain always match and do not overlap [38].
Equations of motion for the structure and fluid are expressed as follows:

Structure equations
The flexible flap is modeled as an elastic continuum therefore the balance of linear momentum is given by: For linear elastic material, the linear Hook's law is followed by the structural stress as: For the established FSI solution, the kinematic, dynamic and geometric conditions at the common interface of fluid and structure should be satisfied, numerically at the interface Γ o (t), the following conditions are assumed for all the time t: Here, v s and v f represent the structure and fluid velocity constituting the kinematic stability. Dynamic stability is attained by the continuous normal stresses of solid σ s and fluid σ f at the interface. n being the unit vector normal at the interface and the subscript i indicates its direction. According to the geometric condition stability fluid-solid domain always match and do not overlap [38]. Equations of motion for the structure and fluid are expressed as follows:

Structure equations
The flexible flap is modeled as an elastic continuum therefore the balance of linear momentum is given by: For linear elastic material, the linear Hook's law is followed by the structural stress as: where ρ s is the structural density and σ s ij.j is the structural stress. The subscript ij of the stresses in Equations (3), (4), (6) and (7) represents the components of stress tensor whereas the subscript j, indicates the differential term involved in the linear elastic model.
The structural stress is a function of strain ε and the Lamé parameters H and λ, which are defined as: where u represents the displacement, E and are Young's modulus and Poisson ratio, respectively [38][39][40].

Fluid equations
The governing equations for time-dependent fluid regime are modeled by incompressible continuity and Navier-Stokes equations in terms of velocity v f , static pressure P and fluid density ρ f as follows [38][39][40]: div The fluid stress for an incompressible fluid is simplified to: The dynamic mesh motion for an FSI problem is defined in terms of arbitrary movement of the grid points with respect to their frame of reference by considering the convection of the material points. The material derivative in terms of the movement of fluid mesh can be expressed as [40]: where U is the mesh movement of fluid, also U = K U sur representing K as a stiffness matrix of an elastic media and U sur being the fluid dynamics surface mesh movement.
Modeling of structure and fluid in ANSYS-Workbench are explained in the following sections:

Structural Model
The solid model is designed in 'Transient Structural' (ANSYS-Mechanical) by a (FE) analysis, afterwards the geometry is defined in the ANSYS Design Modeler. To estimate the flow as close as possible to the experiments, all dimensions and parameters are needed to be precisely defined according to the experimental setup and can be seen in Table 1.
Analogous to the experiments, the structural model is considered as a flexible flap that is, the bottom wall is fixed, while all other walls are free. The fixed wall is assigned as 'Fixed support' whereas the remaining part of the flap is given the 'Fluid Solid Interface.' The flap is meshed by hexahedral elements and are depicted in Figure 2a

Computational Model
The fluid domain is modeled by ANSYS-FLUENT for flaps immersed in a glycerin environment as depicted in Figure 2 (c). Because of the low Reynolds number range, a laminar fluid model is chosen. For an effective simulation in comparison to the experiments, a fixed coordinate system is used. The constant flow velocity as 'velocity inlet' is given for the input conditions, plate and flap are assigned as walls with a 'no slip' boundary condition while the top walls, remaining bottom parts and sides are assigned as 'slip walls' and 'pressure outlet' is defined for 0 Pa at the outlet conditions. A non-dimensional Reynolds number Re is defined as Re = w/υ based on the flow velocity , flap width w and fluid kinematic viscosity υ. The inlet flow velocities, , ranging from 0.125-0.5 m/sec are given as input values based on the identified Reynolds number. The engineering parameters used for the computations are listed in Table 1. Due to subsonic incompressible flow field, pressure based solver with coupled scheme is used and second order implicit transient formulation is employed for modeling the unsteady flow.

Computational Model
The fluid domain is modeled by ANSYS-FLUENT for flaps immersed in a glycerin environment as depicted in Figure 2c. Because of the low Reynolds number range, a laminar fluid model is chosen. For an effective simulation in comparison to the experiments, a fixed coordinate system is used. The constant flow velocity as 'velocity inlet' is given for the input conditions, plate and flap are assigned as walls with a 'no slip' boundary condition while the top walls, remaining bottom parts and sides are assigned as 'slip walls' and 'pressure outlet' is defined for 0 Pa at the outlet conditions. A non-dimensional Reynolds number Re is defined as Re = vw/υ based on the flow velocity v, flap width w and fluid kinematic viscosity υ. The inlet flow velocities, v, ranging from 0.125-0.5 m/sec are given as input values based on the identified Reynolds number. The engineering parameters used for the computations are listed in Table 1. Due to subsonic incompressible flow field, pressure based solver with coupled scheme is used and second order implicit transient formulation is employed for modeling the unsteady flow.

Structure-Flap
Length Bending convergence is shown in Figure 3 whereas different meshes before and after flap deformations are depicted in Figure 4. Moreover, the difference between results obtained with mesh 2 and 3 is small, showing that grid convergence is attained. The difference between mesh 2 and 3 in terms of bending x and height y is 0.53% and 1.1% respectively.  Initially, simulations were performed for a flap of 5 mm thickness. To consider mesh independence, three different meshes were tested to explore the effect of mesh resolution on a flap bending. This mesh resolution was needed to compensate for the effect of computational time such as the computational time increasing after coupling the computational fluid dynamics (CFD) and FEA models [41]. Bending line convergence for different meshes was attained with coarser to finer: mesh 1: 52,486 elements, mesh 2: 101,641 elements and mesh 3: 137,693 elements. So, mesh 3 was chosen for the final solution, although this number of elements may vary since flap deforms and dynamic meshing takes place. Bending convergence is shown in Figure 3 whereas different meshes before and after flap deformations are depicted in Figure 4. Moreover, the difference between results obtained with mesh 2 and 3 is small, showing that grid convergence is attained. The difference between mesh 2 and 3 in terms of bending x and height y is 0.53% and 1.1% respectively. As the large fluid domain is present, it is only desirable to have high mesh resolution around the flap and the dynamic mesh region in order to have influence of pressure and deformation gradient. For a 3D dynamic mesh domain in Fluent, a structured mesh cannot be applied [41]. For that reason, a tetrahedron unstructured mesh was chosen around the flap, though the rest of the fluid domain was designed as a structured mesh. To simulate the flap movement, both smoothing and remeshing methods were used in the dynamic mesh region such that the displacement of the flap was large compared to the grid size. Dynamic meshing is needed for both surface and volume mesh in 3D cases and it is necessary to retain element quality for both surface and volume meshes otherwise simulations may result in errors and computational differences [41]. Considering remeshing parameters, mesh scale information was followed and a minimum cell size of 6.10 −4 m and a maximum cell size of 0.04 m were used as input values. In addition, the interface is allocated as 'system coupling' in dynamic mesh zones under dynamic mesh settings with a cell thickness of 0.6 As the large fluid domain is present, it is only desirable to have high mesh resolution around the flap and the dynamic mesh region in order to have influence of pressure and deformation gradient. For a 3D dynamic mesh domain in Fluent, a structured mesh cannot be applied [41]. For that reason, a tetrahedron unstructured mesh was chosen around the flap, though the rest of the fluid domain was designed as a structured mesh. To simulate the flap movement, both smoothing and remeshing methods were used in the dynamic mesh region such that the displacement of the flap was large compared to the grid size. Dynamic meshing is needed for both surface and volume mesh in 3D cases and it is necessary to retain element quality for both surface and volume meshes otherwise simulations may result in errors and computational differences [41]. Considering remeshing parameters, mesh scale information was followed and a minimum cell size of 6 × 10 −4 m and a maximum cell size of 0.04 m were used as input values. In addition, the interface is allocated as 'system coupling' in dynamic mesh zones under dynamic mesh settings with a cell thickness of 0.6 mm along the wall. To provide a balance among stable mesh movement and total computational time, a time step ∆t of 0.001 s is used for the coupled FSI problem for 2 s of physical flow time where the flap attained the steady-state condition.

Ten mm thick flap:
In order to meet the experimental conditions, FSI analysis for a 10 mm thick flap is performed for different flap orientations, that is, varying the clamped angle from 0° to 90°-for the angle of 0° the flap width is facing the flow, correspondingly angles are changed from 0°-90°. At 90°, the flap is rotated completely and is parallel to the flow direction. Analysis is also performed for an inclined angle of 45 o . The 10 mm flap is modeled by the same structural mesh conditions as a 5 mm thick flap. Also for the fluid domain, the same mesh parameters are employed. Though the total number of elements vary in the range of 137,000-140,000 for different angles. The same boundary conditions and models are specified for 10 mm as described for the 5 mm thick flap. However, minimum and maximum length scales of remeshing differ, depending upon the flap positioning. Flap orientations designed in Design Modeler are shown in Figure 5a

5 mm Thick Flap:
A comparison between FSI and the experimental bending of a flap at different Reynolds numbers is expressed in Figure 6a   Time-dependent total deformation curve of an elastic flap is processed by ANSYS-Transient Structural. Total deformation consists of directional deformations, d, in all three directions that is, X, Y and Z for a given flap. Resultant deformation varies with time from smallest to highest Reynolds number till it reaches its maximum value and then becomes constant. The deformation curves for all four Reynolds numbers at a physical flow time of 2 sec after reaching steady-state can be seen in Figure 7 while Figure 8 indicates the flap visualization versus time till reaching steady-state for a Reynolds number of 9. The deformation is influenced by the magnitude of inflow velocity that is, maximum structural deformation occurs at the tip of the flap and has a least value for Reynolds 3 and increases by 0.0449 m reaching the peak value of 0.0982 m for a Reynolds number of 12. Besides, less computational time is required by the flap for the higher velocity flow in converging to the steady-state value. Thus the FSI model consists of the observations for the maximum deflections offered by the flap.

10 mm Thick Flap:
As described in section 2.2.2, a flap of 10 mm thickness is used to simulate experimental data with different clamping orientations. For the highest Reynolds number of 12, the maximum deformation for a 10 mm thick flap is 0.53 times the 5 mm thickness for angle 0°. However, it continues to decrease with increase in the clamping angle and shows the least at 90° which is 0.37 times the deformation at clamping angle of 0°. Similarly, as the Reynolds number and the clamping angle increases, the x-coordinate of flap bending decreases indicating the minimum flap deflection with highest clamping angle (see Figure 9 b) that is, the maximum x-coordinate for angle 0° is ∆x = 0.039m whereas for 90° it bends to x = 0.011 m. Similarly, a small flap deflection is also observed for inclined angle of 45°. Deformation at each position can be seen in Table 2   maximum structural deformation occurs at the tip of the flap and has a least value for Reynolds 3 and increases by 0.0449 m reaching the peak value of 0.0982 m for a Reynolds number of 12. Besides, less computational time is required by the flap for the higher velocity flow in converging to the steady-state value. Thus the FSI model consists of the observations for the maximum deflections offered by the flap.

10 mm Thick Flap:
As described in section 2.2.2, a flap of 10 mm thickness is used to simulate experimental data with different clamping orientations. For the highest Reynolds number of 12, the maximum deformation for a 10 mm thick flap is 0.53 times the 5 mm thickness for angle 0°. However, it continues to decrease with increase in the clamping angle and shows the least at 90° which is 0.37 times the deformation at clamping angle of 0°. Similarly, as the Reynolds number and the clamping angle increases, the x-coordinate of flap bending decreases indicating the minimum flap deflection with highest clamping angle (see Figure 9 b) that is, the maximum x-coordinate for angle 0° is ∆x = 0.039m whereas for 90° it bends to x = 0.011 m. Similarly, a small flap deflection is also observed for inclined angle of 45°. Deformation at each position can be seen in Table 2

10 mm Thick Flap
As described in Section 2.2.2, a flap of 10 mm thickness is used to simulate experimental data with different clamping orientations. For the highest Reynolds number of 12, the maximum deformation for a 10 mm thick flap is 0.53 times the 5 mm thickness for angle 0 • . However, it continues to decrease with increase in the clamping angle and shows the least at 90 • which is 0.37 times the deformation at clamping angle of 0 • . Similarly, as the Reynolds number and the clamping angle increases, the x-coordinate of flap bending decreases indicating the minimum flap deflection with highest clamping angle (see Figure 9b) that is, the maximum x-coordinate for angle 0 • is ∆x = 0.039m whereas for 90 • it bends to x = 0.011 m. Similarly, a small flap deflection is also observed for inclined angle of 45 • . Deformation at each position can be seen in Table 2     where is the drag force exerted by the fluid on the body and is the fluid density, S is the body area perpendicular to the flow direction and v is the velocity of the flow. For the present simulations, the drag coefficient, CD, for three representative flexible flap cases is computed, that is, for a flap with a thickness of 10 mm and the clamping angles of 0°, 90° and the inclined angle of 45° at a Reynolds number of 12. The frontal area of the flap facing the flow is defined for drag calculations meant for a current 3D case. Figure 11 represents the variation of drag coefficients versus time. The mean CD during t = 1.5-2 sec for clamping angle of 90° is 3.61 and 1.95 for clamping position of 0°. Whereas for an inclined angle of 45° it is smaller than the other two cases and gives a mean of 1.45. It appears that the more bended flap is a smaller obstacle to the flow and therefore reduces the resistance from the mainstream. Moreover, the drag coefficient values are comparable with a benchmark case of a flexible flat plate in a cross flow [28].   where is the drag force exerted by the fluid on the body and is the fluid density, S is the body area perpendicular to the flow direction and v is the velocity of the flow. For the present simulations, the drag coefficient, CD, for three representative flexible flap cases is computed, that is, for a flap with a thickness of 10 mm and the clamping angles of 0°, 90° and the inclined angle of 45° at a Reynolds number of 12. The frontal area of the flap facing the flow is defined for drag calculations meant for a current 3D case. Figure 11 represents the variation of drag coefficients versus time. The mean CD during t = 1.5-2 sec for clamping angle of 90° is 3.61 and 1.95 for clamping position of 0°. Whereas for an inclined angle of 45° it is smaller than the other two cases and gives a mean of 1.45. It appears that the more bended flap is a smaller obstacle to the flow and therefore reduces the resistance from the mainstream. Moreover, the drag coefficient values are comparable with a benchmark case of a flexible flat plate in a cross flow [28].

Flow Regimes
The velocity evaluation along with flap position for thickness 5 and 10 mm intended for clamped angle 0 • and Reynolds 12 is displayed in the present section at various time steps. As discussed above (Section 3.2) and shown in Figure 12, the first two time steps at t = 0.05 s and 0.1 s corresponds to the initial deflection of the flap for both thicknesses. Nevertheless, the maximum deformation reached by the 10 mm thick flap is less for the 5 mm thick flap and attained in less time, whereas, for 5 mm thickness, the flap continues to bend and attain steady-state at approximately t = 1 s.

Flow Regimes
The velocity evaluation along with flap position for thickness 5 and 10 mm intended for clamped angle 0° and Reynolds 12 is displayed in the present section at various time steps. As discussed above (section 3.2) and shown in Figure 12, the first two time steps at t = 0.05 sec and 0.1 sec corresponds to the initial deflection of the flap for both thicknesses. Nevertheless, the maximum deformation reached by the 10 mm thick flap is less for the 5 mm thick flap and attained in less time, whereas, for 5mm thickness, the flap continues to bend and attain steady-state at approximately t = 1 sec. The steady-state 3D configurations of clamping angle 0°, 90° and inclined angle 45° at maximum inlet velocity is shown in Figure 13 Figure 13d,e and also for 5 mm thick flap representing the flow for minimum and maximum Reynolds number in Figure 13f,g, respectively. The higher Reynolds number leads to the higher velocity of the flow at the flap-tip and further downstream in the wake of the flap. Besides, the maximum bending of the flap in Figure 13g bends the streamlines near the flap downstream.
Reynolds number range. Moreover, the flow past the flap is shown and compared for Reynolds numbers of 6 and 12 for clamping angle of 90° in Figure 13d,e and also for 5mm thick flap representing the flow for minimum and maximum Reynolds number in Figure 13f,g, respectively. The higher Reynolds number leads to the higher velocity of the flow at the flap-tip and further downstream in the wake of the flap. Besides, the maximum bending of the flap in Figure 13g bends the streamlines near the flap downstream.

Discussion
The present investigation focuses on the numerical investigation of flexible flaps in a glycerin environment for different thicknesses and flap orientations. The deformation of the flaps was measured to observe the various flap clamping positions under laminar flow and also the maximum tip-deflection of the flap was observed. For both experimental and simulation conditions, the flaps settle down to a steady-state and the bending patterns attained by current simulation model are very close to that observed in the experiments. The focus of the present FSI simulations is to design a 3D interaction between structure and fluid in order to narrow the literature gap for flexible structures.
The results indicate that flap bending increases with the increase in Reynolds number in the range of 3 < Re < 12. However, for a different clamping position of 0 • to 90 • indicates minimum flap deflection with the highest clamping angle. Correspondingly, for an inclined position of 45 • less deflection is perceived as compared to its undeformed shape. The objective is to judge the flap positions and flow regimes at varying computational time and to observe transient deformation response under different flow velocities. Likewise, the drag coefficients are also predicted for various flap alignments. Henceforward, the numerical model is suitable to accurately predict the response of thick and relatively thinner structures.
The flap bending in the current work is observed for low Reynolds number; however, the bending and periodic nature of flap motion is schedule by another 3D FSI model for higher Reynolds number in case of water as a working fluid. For this purpose, experiments will be performed for thinner flaps and the numerical model will be helpful in understanding the flow physics around the flaps. This will be further helpful in designing the flexible structures in hemodynamic, biomedical and structural dynamics.

Conflicts of Interest:
The authors declare no conflict of interest.

Nomenclature
The following symbols are used in the manuscript: