Flow-Induced Vibration of Cantilever Type Elastic Material in Straight Tricylinder

Featured Application: The findings of this study can be utilized to establish a rational foundation for the arrangement of fluid–structure coupling low-frequency energy harvesting devices, particularly in aquatic environments. Abstract: The hydrodynamic performance and wake evaluation results of a cantilever flexible cylindrical triangular array at Re = 9000 and with L/D spacing ratios ranging from 3 to 7 are presented in this paper. The improved delayed separation eddy current simulation (IDDES) SST k-ω turbulence model and tetrahedral mesh are employed to solve the three-dimensional instantaneous flow field. The objective is to find the L/D ratio that the upstream cylinder has the most pronounced effect on the downstream cylinders, to quantify whether a flexible upstream cylinder is more favorable than the rigid upstream cylinder, and to know the flexible cylinder wake distribution. The findings indicate that, at Re = 9000, when the spacing ratio is 3~7 and the spacing ratio L/D is 6, the amplitude of the downstream cylinder is most affected by the wake of the upstream cylinder; the amplitude of downstream cylinders of the rigid upstream cylinder exerted a 1% or 2% greater effect than the flexible cylinder, with vortices predominantly concentrated in the middle and upper parts of the flexible cylinder rather than at its top.


Introduction
The fluid-structure interaction of slender flexible cylinders in a cross-flow is pivotal in various engineering applications, particularly in nuclear reactors, low-frequency energy harvesting, and offshore engineering involving cantilever slender flexible structures [1,2].The elongated offshore structures are prone to VIV when exposed to ocean currents.Sarpkaya [3] highlights that VIV is characterized by nonlinearity, self-excitation, self-limitation, and multiple degrees of freedom.Even the simplest form of VIV, such as the single-cylinder vortex-induced vibration in two-dimensional uniform flow, remains incompletely understood.In the case of slender and flexible cylinders, it has been observed that the dynamics of VIV response are chaotic, complex, and consist of unpredictable motions [4].A vibrating flexible cylinder subjected to cross-flow generates a complex vortex wake that significantly influences the trajectory of response motion, phase synchronization between fluid force structural motion, as well as the interaction between the in-line (IL) and cross-flow (CF) responses.For cantilevered flexible cylinders with fewer constraints imposed on them, predicting the characteristics of wake flow becomes more challenging.
At present, the study of vortex-induced vibration of flexible cylinders has attracted the attention of many scholars [5,6].Chaplin et al. [7] checked various early empirical models (VIVA, VIVANA, VI-Co Mo, SHEAR7, and ABAVIV) according to the experimental results of flexible cylinder vortex-induced vibration and found that there were some errors between the model prediction results and the experiment.He et al. [8] studied three relatively new models, SST-PANS, SST-SAS, and SST-IDDES, and the results showed that SST-IDDES had the best simulation results.The numerical results from SST-IDDES show its comparable capabilities for the simulation of massively separated hydrodynamic flows and its potential application in the prediction of industrial turbulent flows for vortex-induced motions (VIM).Wu [9] studied the coupling mechanism of the cross-flow direction and forward-flow direction of flexible cylindrical VIV and found that the fluid force coefficient was affected by the phase between CF and IL displacement.In addition, when the rigid cylinder VIV occurs, the structure does not deform, and the fluid force at different section locations is roughly the same, while the flexible cylinder has a large length-diameter ratio (the ratio of length to diameter), and the form of vortex shedding and the fluid force change irregularly along the axis.VIV has typical three-dimensional characteristics, and the spatial distribution of fluid force needs further study.
In research on three cylinders, Yang [10] studied the flow past three circular cylinders in equilateral triangular arrangements at Re = 50-200 from the viewpoint of vortex-shedding suppression.Their results showed that for one upstream cylinder and the other two side-by-side downstream, the vortex shedding from the upstream cylinder is suppressed for Re = 100-200.However, for two upstream cylinders side-by-side and one downstream cylinder, the suppression occurs on the downstream cylinder at Re ≤ 175.At a high Re = 1.1 × 10 4 ~6 × 10 4 , Dahl [11] studied the vortex-induced vibration of a cylinder with two degrees of freedom via experiments and found that the amplitude of the lateral displacement response could reach a maximum of 1.35 times the diameter.Gao et al. [12] simulated the flow past three circular cylinders with the incidence angle of 30° for 200 ≤ Re ≤ 3900 and 1.25 ≤ D/L ≤ 4.0, and five flow patterns were identified.The boundary spacing ratios between the five regimes are affected by the Reynolds number.With the increase in the spacing ratio or the Reynolds number, the three-dimensionality of the wake flow becomes stronger.Zhang [13] used three-dimensional numerical simulations to investigate the flow around three cylinders in this arrangement at the super-critical Reynolds number Re = 3 × 10 6 , concentrating on the influence of the spacing ratio among cylinders.When L/D is greater than 3.5, the separation point of the cylinder in the upper stream is close to that of a single cylinder, indicating that the minimum L/D for negligible interaction among the cylinders is 3.5.
Learning about flexible slender risers, Wang [14] studied the VIV of a vertical riser subject to uniform and linearly sheared currents.A standing wave response is observed for the single-mode IL and CF vibrations.Dual resonance is found to occur at most of the locations along the riser.Han [15] developed a bidirectional fluid-structure coupling method.By comparing and analyzing the vortex forms of riser wake under different flow rates, it is found that the riser has a slight 3D effect at low flow rates, while the riser has a strong 3D effect at high flow rates.Lin [16] used a Q3D FSI numerical method that couples the strip theory-based DVM and the FEM of structural dynamics developed to simulate the FIV of the flexible cylinder.The study answered two questions: What is the difference between the flexible cylinder FIV response with an upstream wake interference and the isolated flexible cylinder VIV response?How does the difference between the FIV response with the upstream wake interference and the isolated cylinder VIV response come about?Deng [17] used numerical investigations into the VIV of a flexible cylinder with different length-to-diameter ratios, named the aspect ratio.Comparisons among simulations with different aspect ratios were analyzed in detail.The maximum vibration amplitude, the dominant vibration mode, and the non-dimensional dominant vibration frequency in both directions were apparently enlarged with the increase in the cylinder aspect ratio.The typical vibration trajectory changes from the butterfly type to the stripe type with L/D increases.Wang [18] also used the two-dimensional method to study the flexible slender rod and made a comparison with the experiment to prove the feasibility of this method, but these methods were not a true three-dimensional simulation.
This study focuses on the numerical development and analysis of the VIV response of cantilever flexible cylindrical triangular arrays at subcritical Reynolds numbers.To enhance the investigation into the VIV characteristics of such arrays, we analyze the distribution of flexible cylindrical structures with triangular rigid connections using a combination of the SST-IDDES model and the fluid-structure coupling model.

The Governing Equations
When calculating the flow field, it is essential to ensure the fulfillment of both the continuity equation and momentum equation as a prerequisite, with the governing equation presented as follows.
The fluid will cause the vibration of the flexible cylinder, and its governing equation is as follows: In order to ensure that the fluid-structure coupling follows the most basic conservation principle at the fluid-structure coupling interface, the following requirements should be met: To solve the fluid-structure coupling problem, it is necessary to list the equations of motion in parallel.Based on the continuity equation and momentum equation of the fluid, the motion equation of the fluid is as follows: When the flexible cylinder is subjected to the pressure exerted by the fluid, its vibration equation is as follows: ⋅  ̈ +    ̇ +   ⋅ ⅇ −   ⋅   =   (6) By connecting Equations ( 5) and ( 6), the equation required to solve the bidirectional fluid-structure coupling is obtained as follows: On the fluid-solid coupling surface, each node has the same degree of freedom of displacement and pressure.Then, combined with the known boundary conditions, initial conditions, and given relevant calculation parameters, the coupling equation is solved, and the solution of the equation is obtained.After determining the solution of the contact point of the fluid-solid coupling surface, the solution vector on the coupling surface can be obtained via Equation (7) to solve the fluid-solid coupling dynamics problem.

Numerical Simulation
The computational fluid dynamics solver ANSYS Fluent version 2021 R1 is based on the finite volume method and is utilized for solving the three-dimensional unsteady incompressible Navier-Stokes equation with IDDES turbulent closure.
The mechanical solver ANSYS Transient Structural and the computational fluid dynamics solver ANSYS Fluent via ANSYS System Coupling to solve the problem of the interaction between the pressure generated by the fluid and the displacement generated by the solid, as well as the deformation of the flexible cylinder.Figure 1 is the flow chart of the coupling algorithm.

Computational Domain
The top view of the flow field is constructed based on the cylindrical array structure, with the origin of coordinates designated as O. Figure 2 illustrates the arrangement of cylinders in a cylindrical configuration.The flow field model consists primarily of two components: the near-wall flow of the cylinder and the outflow field on the opposite wall.
The overall dimensions of the flow field are 100 D in length, 40 D in width, and height.The near-wall flow of the cylinder and the outflow field on the opposite wall, with a distance of 2 D separating it from adjacent walls along its length and width directions.The wake direction encryption area is up to 50 D, which is convenient to study.For boundary conditions, pressure inlet and outlet are employed, while standard non-slip wall boundary conditions are applied to both the flexible cylinder and solid wall.To enhance computational efficiency and precision in capturing the vibration and complex turbulence characteristics of the cylinder, a grid-based encryption technique is employed for the fluid-structure coupling surface of the flexible cylinder.In order to minimize DES gray space error, this study utilizes near-wall encrypted flow field to transition the fluid-structure coupling surface into the outflow field grid.The outflow field grid represents a turbulent region dominated by large cavities, with LES cutoff scale primarily dependent on local grid scale and vortex system structure size, determining appropriate grid size selection.

The Grid Dependence Study
The investigation of grid dependence is a pivotal step in numerical computation as the accuracy of the solution can be influenced by the resolution of the grid.To conduct a prescribed computational fluid dynamics (CFD) simulation, it is imperative to assess both the accuracy and time required for obtaining a solution.Therefore, four distinct types of grids will be generated with increasing grid resolution.For Re = 10000, flow simulations were performed at these different mesh densities and compared against Fujarra's (1997) [19] experiment and C.T. Yamamoto's [20].
Table 1 1, resistance coefficient is less than 1%, while the Strouhal number remains within a 1% range, even with gradual grid refinement.This indicates that the grid independence study has reached convergence.The fine grid is deemed sufficient for subsequent calculations, as shown in Figure 3.

Lift Drag Coefficient
The drag coefficient   and lift coefficient   are more effective in reflecting the primary stress on the cylinder for VIV.The lift and drag coefficient quantify the resistance experienced by the flexible cylinder in the direction of incoming flow (x-axis) and the lift force exerted on it in the lateral direction (z-axis).Figure 4 illustrates the time history curve of t   and   at an L/D distance ratio ranging from 3 to 7. At L/D ratios of 3, 4, 5, and 6, upstream cylinder one exhibits a lower to   by approximately 16-27% compared to downstream cylinders two and three.Moreover, there is a drag difference between downstream cylinders two and three, amounting to around 4-8%.This phenomenon occurs due to the influence of a wake generated during VIV migration on downstream cylinders as they deform towards the middle O point.Consequently, upstream cylinder one becomes locked with a downstream cylinder, resulting in continuous deviation towards this particular cylinder while experiencing reduced resistance.This effect is particularly evident when considering spacing ratios L/D of either five or six.The coefficient   generally ranges between −0.1 and +0.1 without significant amplitude differences observed between upstream and downstream cylinders; however, occasional interference from downstream cylinders may cause fluctuations around zero for upstream cylinder motion primarily influenced by wake effects from upstream cylinder, it causing negative deformation along Z-axis direction by negative forces for downstream cylinder two, whereas an opposite trend can be observed for downstream cylinder three instead.When considering a spacing ratio L/D equal to seven, these interference phenomena between cylinders are significantly reduced, as clearly depicted in Figure 5.Note that some profiles exhibit a noticeable oscillation in Figure 4.In current transient structure-fluid coupling simulations, we use a relatively large time step yet still achieve convergent solutions for each case due to the computational resource required for such three-dimensional bidirectional interactions.
The surface pressure coefficient of a cylinder with a spacing ratio of D/L = 3−7 is illustrated in Figure 5.It can be observed that the differences in pressure coefficients among different spacing ratios are relatively small; Figure 5 shows the pressure at the top section of the cylinder.After its integration, the component in the X direction is the cross-section resistance   , and the component in the Z direction is the cross-section lift   .The resistance is mainly obtained by subtracting the pressure near 0° from the pressure near ±180°.The lift force is mainly obtained by subtracting the pressure around −90° from the pressure around 90°.Of course, this is just the force on the cross-section, and the resistance to lifting on the cylinder has to be integrated into the Y direction.We can briefly judge the force on the section according to the difference in pressure at different angles.

Vorticity Field
Figure 6 illustrates typical vorticity results.One can see that the shedding vortex in the lower region of the cylinder, where the height ratio H/D is smaller, remains relatively unchanged irrespective of the period.Only a minor vortex shedding occurs in this area; the presence of the cylinder significantly influences its trajectory by inducing a persistent skewing effect.The phenomenon of vortex shedding becomes more pronounced in the middle section of the cylinder and exhibits noticeable variations with changes in the period.The influence exerted by the cylinder on these vortices gradually diminishes.Towards the upper part of the cylinder, there is a weakening in vortex shedding phenomena, and overall minimal alterations are observed, regarding the full development of vortices.Additionally, from flow diagrams, it becomes apparent that a distinct vortex structure exists when H/D = 15 at mid-cylinder, which greatly affects convection-induced vibration phenomena.

Streamline Field
Based on the observations from Figures 7a and 5, it can be inferred that the wake of the upstream cylinder has an approximate influence range of 2 D along the Z−axis.In Figure 7b, it is evident that for L/D = 3 and L/D = 4, the wake exhibits a highly disordered pattern, indicating its impact on the downstream cylinder's wake.For L/D = 5 and L/D = 6, vortices are still discernible; however, some flow lines remain relatively straight and unaffected by vortices.The vortex from upstream and downstream mixed to produce a bigger scale vortex, which indicates that the upstream cylinder wake has a continuous influence on the downstream cylinder wake.When L/D = 7, there is a significant reduction in vortex activity within the wake region; numerous flow lines align closely with straight trajectories, while only minimal mixing occurs between vortices generated by both cylinders.Notably absent are any large-scale vortices affecting all three cylinders simultaneously, thus indicating a sharp decrease in the influence exerted by an open wake originating from an upstream source.It is evident in Figure 7 that the longitudinal vorticity primarily is concentrated in the middle and lower regions.Combining with the vorticity shown in Figure 6 further reveals that the intensity of vorticity is higher near H/D = 15 and H/D = 22.5 in the middle region.The development of eddy currents at H/D = 30 on top is predominantly horizontal due to its location above the cylinder's top.As a result, it remains relatively unaffected by the upstream cylinder, leading to a stable flow with reduced Y−axis direction vorticity and dominance in the horizontal direction.In the area near the middle of the column, the influence of the upstream column is greater so that the vorticity formed by it changes obviously in the vertical direction.With an increase in spacing ratio L/D, the wake tends to flatten both horizontally and vertically.
The Z−axis flow velocity diagram of different cylinders along the X-axis section when D/L = 6 is presented in Figure 8.The Z−axis flow diagram allows for a quick visualization of the distribution of the flexible cylinder in the Y−axis vortex.For instance, when H/D = 7.5, numerous regions with opposite flow directions can be observed near cylinder 1 in Figure 8a, indicating the formation of vortices and indirectly reflecting their intensity.In Figure 8a, concentrated regions opposing the flow direction can be seen around H/D = 7.5, with lengths reaching up to 10 D, suggesting significant vorticity in these areas.In Figure 8b, parts within the same region exhibiting opposite flow directions appear more dispersed, indicating lower vorticity levels here compared to Figure 8a.In Figure 8c, parts within the same region displaying opposite flow directions are concentrated; however, the length of this region is only 5 D, implying smaller vorticity than that observed in Figure 8a but larger than that in Figure 8b.From Figure 5, it can be observed that at this point in time, wake streams from cylinder 1 and cylinder 2 as well as from cylinder 3 mutually influence each other; furthermore, when H/D = 7.5, it is evident that there is relatively high vortex intensity associated with cylinder 3, while relatively low vortex intensity corresponds to cylinder 2 at this point in time, which supports this conclusion.

Amplitude
For a slender and flexible cylinder, it has been observed that the response dynamics of VIV exhibit chaotic and complex behavior, characterized by non-repeatable motions [4].Consequently, unlike a rigid cylinder, the vibration of a flexible cylinder does not display obvious periodic changes.However, it does exhibit distinct fluctuation intervals which can provide insights into its movement trends.Figure 9 illustrates the Z−axis vibration curve of three cylinders with an L/D ratio of 3 and an H/D ratio of 30.The Z−axis vibration of the three cylinders tends to remain generally stable within an amplitude interval of approximately 0.1 D. In particular, for the upstream cylinder 1, its Z−axis motion is symmetric along Z = 0; whereas for cylinder 2, under the influence of flow field effects, its motion path inclines towards the positive direction on the Z−axis with a deviation around 0.3 D. On the other hand, inflow from cylinder 3 causes its motion trajectory to deviate in the negative direction on the Z−axis with an approximate deviation magnitude also at about −0.3 D. Overall, both cylinders' deviation trajectories are symmetrical along Z = 0.

Figure 1 .
Figure 1.Flow chart of the coupling algorithm.

Figure 2 .
Figure 2. The schematics of the computational domain for three cylinders in an equilateral-triangle arrangement: (a) 3D view of the computation domain; (b) top view of the computation domain.

A
time step  of 0.5 is chosen to maintain Courant number Co, ||/  < 1,   in the LES area of the hybrid grid, where   represents maximum mesh element size.The model employs liquid water as jet viscosity of 1.003 × 10 −3 Pa•s, density of 998.2 kg/m 3 , inlet flow rate of 1.5 m/s, Re = 9000, Re = UD/μ, and performs 20 iteration steps per unit time.Flexible cylinder parameters D = 6 mm, density = 950 kg/m 3 , Young's modulus = 1100 MPa, Poisson's ratio = 0.42, tensile yield strength = 25 MPa, and compressive yield strength = 33 MPa.

Figure 7 .
Figure 7. Main view (a) and right view (b) of a flow chart with L/D = 3−7.The purple bead flow line originates from the cylinder, and the green streamer flow line passes between the cylinders.

Figure 9 .
Figure 9.The Z−axis vibration curve of the cylinder with L/D ratio of 3 and H/D ratio of 30.The deformation of flexible cylinder A at an L/D ratio of 3 and an H/D ratio of 30 is illustrated in Figure 10.In the figure, cylinder 2 and cylinder 3 indicate that the upstream cylinder 1 exhibits flexibility in the flow field.Cylinder 2* and cylinder 3* denote that in the flow field, the upstream cylinder 1 is rigid.The simulation results demonstrate that, in comparison with the rigid cylinder, the upstream flexible cylinder exhibits an overall decrease of approximately 0.02 D in terms of downstream stream-induced vibration amplitude A. Moreover, the amplitude range of cylinder 2 and cylinder 3 is roughly 0.01 D smaller than that of cylinder 2* and cylinder 3*.Regardless of whether it is cylinder 2, cylinder 3, or their counterparts denoted as cylinders 2* and 3*, all exhibit maximum amplitudes exceeding 1 D. Considering a spacing ratio L/D = 2 may lead to collision risks; hence, this study commences from a spacing ratio L/D = 3.

Figure 10 .
Figure 10.Where the L/D ratio is 3 and the H/D ratio is 30, the deformation of flexible cylinder is A. The variation in amplitude A of L/D with different spacing ratios at an H/D ratio of 30 is illustrated in Figure 11.It can be observed from the figure that, within the range of the L/D ratio between 3 and 6, both the average values of amplitude A/D and the amplitude interval increase as the spacing ratio L/D increases.The maximum occurs at H/D = 6, where the maximum amplitude A/D reaches 1.23, and the amplitude interval reaches 0.1D.When H/D exceeds 7, although the average relative amplitude A/D approaches it, with a spacing ratio of L/D = 3 exists a larger fluctuation interval of 0.1 D compared to when L/D = 3.By combining Figures 9-11 together, it can be concluded that the

Table 1 .
reports the results of the grid convergence study.Strouhal number   (dimensionless desvorticity frequency), drag coefficient   and dimensionless number Y + Mesh parameter.
The grid dependence testing involved the use of four different mesh, with the total time average resistance coefficient   and Strouhal number   parameters.As shown in Table