Search

We present the development and benchmarking of an in-house fluid-structure interaction (FSI) solver. An implicit partitioned approach is utilized to couple a sharp-interface immersed boundary (IB) method based flow solver and a finite-element method based structural solver. In the present work, the coupling is accelerated using a dynamic under-relaxation scheme. The revised coupling is around two to three times faster and numerically stable, as compared to the one that uses a constant under-relaxation parameter. The solver is validated against two FSI benchmarks in which a thin, finite thickness, elastic splitter plate is attached to the lee side of a circular or square rigid cylinder, subjected to laminar flow. In these two-dimensional benchmarks, the flow induces a wave-like deformation in the plate, and it attains a periodic self-sustained oscillation. We employ the FSI solver to analyze the flow-induced vibration (FIV) of the plate in a uniform laminar free-stream flow for a wide range of mass ratio and bending stiffness at Reynolds number ($Re$) of 100, based on the diameter of the cylinder. At the given $Re$, two-dimensional numerical simulations show that the FIV of the plate effectively depends only on the mass ratio and bending stiffness. The largest displacement of the plate vibration is found to occur in the lock-in region, where the vortex shedding frequency of the coupled fluid-structure system is close to the natural frequency of the splitter plate. We briefly discuss wake structures and phase plots for different cases of mass ratio and bending stiffness.

The generation of electrical power from Vortex-Induced Vibration (VIV) of a cylinder is investigated numerically. The cylinder is free to oscillate in the direction transverse to the incoming flow. The cylinder is attached to a magnet that can move along the axis of a coil made from conducting wire. The magnet and the coil together constitute a basic electrical generator. When the cylinder undergoes VIV, the motion of the magnet creates a voltage across the coil, which is connected to a resistive load. By Lenz's law, induced current in the coil applies a retarding force to the magnet. Effectively, the electrical generator applies a damping force on the cylinder with a spatially varying damping coefficient. For the initial investigation reported here, the Reynolds number is restricted to Re < 200, so that the flow is laminar and two-dimensional (2D). The incompressible 2D Navier-Stokes equations are solved using an extensively validated spectral-element based solver. The effects of the electromagnetic (EM) damping constant xi_m, coil dimensions (radius a, length L), and mass ratio on the electrical power extracted are quantified. It is found that there is an optimal value of xi_m (xi_opt) at which maximum electrical power is generated. As the radius or length of the coil is increased, the value of xi_opt is observed to increase. Although the maximum average power remains the same, a larger coil radius or length results in a more robust system in the sense that a relatively large amount of power can be extracted when xi_m is far from xi_opt, unlike the constant damping ratio case. The average power output is also a function of Reynolds number, primarily through the increased maximum oscillation amplitude that occurs with increased Reynolds number at least within the laminar range, although the general qualitative findings seem likely to carry across to high Reynolds number VIV.

The flow-induced deformation of a thin, elastic splitter plate attached to the rear of a circular cylinder and subjected to laminar pulsatile inflow is investigated. The cylinder and elastic splitter plate are contained within a narrow channel and the Reynolds number is mostly restricted to Re = 100, primarily covering the two-dimensional flow regime. An in-house fluid-structure interaction code is employed for simulations, which couples a sharp-interface immersed boundary method for the fluid dynamics with a finite-element method to treat the structural dynamics. The structural solver is implicitly (two-way) coupled with the flow solver using a partitioned approach. This implicit coupling ensures numerical stability at low structure-fluid density ratios. A power spectrum analysis of the time-varying plate displacement shows that the plate oscillates at more than a single frequency for pulsatile inflow, compared to a single frequency observed for steady inflow. The multiple frequencies obtained for the former case can be explained by beating between the applied and plate oscillatory signals. The plate attains a self-sustained time-periodic oscillation with a plateau amplitude in the case of steady flow, while the superimposition of pulsatile inflow with induced plate oscillation affects the plateau amplitude. Lock-in of the plate oscillation with the pulsatile inflow occurs at a forcing frequency that is twice of the plate natural frequency in a particular mode and this mode depends on the plate length. The plate displacement as well as pressure drag increases at the lock-in condition. The percentage change in the maximum plate displacement, and skin-friction and pressure drag coefficients on the plate, due to pulsatile inflow is quantified. The non-linear dynamics of the plate and its coupling with the pulsatile flow are briefly discussed.