A variationally bounded scheme for delayed detached eddy simulation: Application to vortex-induced vibration of offshore riser

Abstract We present a bounded and positivity preserving variational (PPV) method for the turbulence transport equation of Spalart–Allmaras based delayed detached eddy simulation (DDES). We employ the developed solver to simulate the vortex-induced vibration of a slender flexible riser immersed in a turbulent flow. The fluid-structure interface problem is solved by recently proposed partitioned iterative scheme [1], which consists of nonlinear force correction for a stable coupling of the Navier–Stokes equations with the low-mass structure subjected to strong inertial effects from the surrounding incompressible flow. In our fluid-structure model, we consider the flexible cylindrical riser as a long tensioned beam via linear modal analysis. In the present contribution, we first validate the PPV-based DDES solver for the flow past a three-dimensional stationary cylinder at two subcritical Reynolds numbers of R e = 3900 and R e = 140 , 000 based on the diameter of the cylinder. We next simulate a three-dimensional flexible riser model under a pinned-pinned condition at R e = 4000 and compare our results with that of the experimental measurements. We assess the response characteristics at various locations along the span of the riser and discuss the orbital trajectories at those locations. In particular, the numerical study emphasizes on (i) The effectiveness of the modal analysis in modeling the riser dynamics and a positivity-preserving variational method in DDES for turbulence modeling, (ii) validation of response amplitudes and motion trajectories, and (iii) new insights on the trajectories and vortical structures of the flexible riser undergoing vortex-induced vibrations (VIV). We confirm a standing wave response and complex chaotic response of riser VIV observed in the recent experiments.

[1]  Vaibhav Joshi,et al.  A positivity preserving variational method for multi-dimensional convection-diffusion-reaction equation , 2017, J. Comput. Phys..

[2]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

[3]  R. Govindarajan,et al.  Vortex shedding patterns, their competition, and chaos in flow past inline oscillating rectangular cylinders , 2010, 1012.2369.

[4]  Yuri Bazilevs,et al.  Free-Surface Flow and Fluid-Object Interaction Modeling With Emphasis on Ship Hydrodynamics , 2012 .

[5]  M. M. Zdravkovich,et al.  Flow induced oscillations of two interfering circular cylinders , 1985 .

[6]  Rainald Löhner,et al.  Improved error and work estimates for high‐order elements , 2013 .

[7]  T. P. Miyanawala,et al.  Partitioned iterative and dynamic subgrid-scale methods for freely vibrating square-section structures at subcritical Reynolds number , 2016 .

[8]  Enhao Wang,et al.  Numerical simulation of vortex-induced vibration of a vertical riser in uniform and linearly sheared currents , 2016 .

[9]  G. Karniadakis,et al.  Vortex-induced vibrations of a long flexible cylinder in shear flow , 2011, Journal of Fluid Mechanics.

[10]  Jason Dahl,et al.  Vortex-induced vibration of a circular cylinder with combined in-line and cross-flow motion , 2008 .

[11]  Rajeev K. Jaiman,et al.  A stable second-order partitioned iterative scheme for freely vibrating low-mass bluff bodies in a uniform flow , 2016 .

[12]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[13]  Victor M. Calo,et al.  Improving stability of stabilized and multiscale formulations in flow simulations at small time steps , 2010 .

[14]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[15]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[16]  G. Karniadakis,et al.  Resonant vibrations of bluff bodies cause multivortex shedding and high frequency forces. , 2007, Physical review letters.

[17]  Chao Yan,et al.  Detailed Investigation of Detached-Eddy Simulation for the Flow Past a Circular Cylinder at Re=3900 , 2012 .

[18]  Mike Campbell,et al.  Actual VIV Fatigue Response of Full Scale Drilling Risers: With and Without Suppression Devices , 2008 .

[19]  Stefan Turek,et al.  A Monolithic FEM/Multigrid Solver for an ALE Formulation of Fluid-Structure Interaction with Applications in Biomechanics , 2006 .

[20]  E. Lamballais,et al.  Experimental and numerical studies of the flow over a circular cylinder at Reynolds number 3900 , 2008 .

[21]  Thomas J. R. Hughes,et al.  Conservation properties for the Galerkin and stabilised forms of the advection–diffusion and incompressible Navier–Stokes equations , 2005 .

[22]  P. Spalart.  Detached-Eddy Simulation , 2009 .

[23]  Miguel R. Visbal,et al.  A time‐implicit high‐order compact differencing and filtering scheme for large‐eddy simulation , 2003 .

[24]  Donald Rockwell,et al.  Flow structure from an oscillating cylinder Part 2. Mode competition in the near wake , 1988, Journal of Fluid Mechanics.

[25]  Thomas J. R. Hughes,et al.  What are C and h ?: inequalities for the analysis and design of finite element methods , 1992 .

[26]  Rajeev K. Jaiman,et al.  Industrial application of RANS modelling: capabilities and needs , 2009 .

[27]  Vivek Jaiswal,et al.  The Effectiveness of Helical Strakes in the Suppression of High-Mode-Number VIV , 2006 .

[28]  Dominique Pelletier,et al.  Perspective on the geometric conservation law and finite element methods for ALE simulations of incompressible flow , 2009, J. Comput. Phys..

[29]  Philippe R. Spalart,et al.  Detached-Eddy Simulations Past a Circular Cylinder , 2000 .

[30]  Halvor Lie,et al.  Experimental investigation of vortex-induced vibration of long marine risers , 2005 .

[31]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[32]  Jaime Peraire,et al.  Navier-Stokes Solution Using Hybridizable Discontinuous Galerkin methods , 2011 .

[33]  J. Fröhlich,et al.  Hybrid LES/RANS methods for the simulation of turbulent flows , 2008 .

[34]  Rajeev K. Jaiman,et al.  Fully Coupled Fluid-Structure Interaction for Offshore Applications , 2009 .

[35]  Charbel Farhat,et al.  The discrete geometric conservation law and the nonlinear stability of ALE schemes for the solution of flow problems on moving grids , 2001 .

[36]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[37]  Hamn-Ching Chen,et al.  Vertical Riser VIV Simulation in Uniform Current , 2010 .

[38]  R Spalart Philippe,et al.  Young-Person''s Guide to Detached-Eddy Simulation Grids , 2001 .

[39]  J. W. G. van Nunen.  Pressures and Forces on a Circular Cylinder in a Cross Flow at High Reynolds Numbers , 1974 .

[40]  Joris Degroote,et al.  Partitioned Simulation of Fluid-Structure Interaction , 2013 .

[41]  L. Franca,et al.  Stabilized finite element methods. II: The incompressible Navier-Stokes equations , 1992 .

[42]  Samuel Holmes,et al.  Simulation of Riser VIV Using Fully Three Dimensional CFD Simulations , 2006 .

[43]  J. Kim Vandiver,et al.  Drag Coefficients of Long Flexible Cylinders , 1983 .

[44]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[45]  Yiannis Constantinides,et al.  NUMERICAL PREDICTION OF BARE AND STRAKED CYLINDER VIV , 2006 .

[46]  Vivek Jaiswal,et al.  Insights on vortex-induced, traveling waves on long risers , 2009 .

[47]  Yiannis Constantinides,et al.  Numerical prediction of VIV and comparison with field experiments , 2008 .

[48]  C. Autermann,et al.  崩壊Bs0→Ds(*)Ds(*) , 2007 .

[49]  Tayfun E. Tezduyar,et al.  Finite elements in fluids: Stabilized formulations and moving boundaries and interfaces , 2007 .

[50]  Mark Bull,et al.  Development of mixed mode MPI / OpenMP applications , 2001, Sci. Program..

[51]  J. Wallace,et al.  The velocity field of the turbulent very near wake of a circular cylinder , 1996 .

[52]  P. Moin,et al.  Numerical studies of flow over a circular cylinder at ReD=3900 , 2000 .

[53]  Amine Ben El Haj Ali,et al.  A positivity preserving finite element–finite volume solver for the Spalart–Allmaras turbulence model , 2007 .

[54]  Rajeev K. Jaiman,et al.  Assessment of conservative load transfer for fluid–solid interface with non‐matching meshes , 2005 .

[55]  A. Roshko.  Experiments on the flow past a circular cylinder at very high Reynolds number , 1961, Journal of Fluid Mechanics.

[56]  Rajeev K. Jaiman,et al.  Conservative load transfer along curved fluid-solid interface with non-matching meshes , 2006, J. Comput. Phys..

[57]  Yiqing Shen,et al.  Large eddy simulation using a new set of sixth order schemes for compressible viscous terms , 2010, J. Comput. Phys..

[58]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuscka-Brezzi condition: A stable Petrov-Galerkin formulation of , 1986 .

[59]  James Forsythe,et al.  Large and Detached Eddy Simulations of a Circular Cylinder Using Unstructured Grids , 2003 .

[60]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations , 1991 .

[61]  Michael S. Triantafyllou,et al.  Chaotic response is a generic feature of vortex-induced vibrations of flexible risers , 2011 .

[62]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[63]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .