DEVELOPMENT OF A THREE-DIMENSIONAL ITERATIVE METHODOLOGY USING A COMMERCIAL CFD CODE FOR FLOW SCOURING AROUND BRIDGE PIERS

摘要

One of the major concerns affecting the safety of bridges with foundation supports in river-beds is the scouring of riverbed material from bridge supports during floods. Scour is the engineering term for the erosion caused by water around bridge elements such as piers, monopiles, or abutments. Scour holes around a monopile can jeopardize the stability of the whole structure and will require deeper piling or local armoring of the river-bed. About 500,000 bridges in the National Bridge Registry are over waterways. Many of these are considered as vulnerable to scour, about five percent are classified as scour critical, and over the last 30 years bridge failures caused by foundation scour have averaged about one every two weeks. Therefore it is of great importance to predict the correct scour development for a given bridge and flood conditions. Apart from saving time and money, integrity of bridges are important in ensuring public safety. Recent advances in computing boundary motion in combination with mesh morphing to maintain mesh quality in computational fluid dynamic analysis can be applied to predict the scour hole development, analyze the local scour phenomenon, and predict the scour hole shape and size around a pier.The main objective of the present study was to develop and implement a three dimensional iterative procedure to predict the scour hole formation around a cylindrical pier using the mesh morphing capabilities in the STARCCM+ commercial CFD code. A computational methodology has been developed using Python and Java Macros and implemented using a Bash script on a LINUX high performance computer cluster. An implicit unsteady approach was used to obtain the bed shear stresses. The mesh was iteratively deformed towards the equilibrium scour position based on the excess shear stress above the critical shear stress (supercritical shear stress). The model solves the flow field using Reynolds Averaged Navier-Stokes (RANS) equations, and the standard k-e turbulence model. The iterative process involves stretching (morphing) a meshed domain after every time step, away from the bottom where scouring flow parameters are supercritical, and remeshing the relevant computational domain after a certain number of time steps when the morphed mesh compromises the stability of further simulation. The simulation model was validated by comparing results with limited experimental data available in the literature.