{'fraction_non_alphanumeric': 0.032086979581012996, 'fraction_numerical': 0.007425086184036065, 'mean_word_length': 4.824710424710425, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Numerical Simulation 3D of Local Scour around Subsea Caisson Structures in Steady Currents and Waves

This paper studied the local scour processes at subsea structures using using the computational fluid dynamics codes (OPENSOURCE CODE). The model solved the fully threedimensional, Reynolds-averaged Navier–Stokes equation closed with the LES model. Stokes and Cnoidal Waves were generated using Fenton’s Fourier series method. Two empirical models must be used to compute the boundary layer at the interface and bed-load transport. We used the types of structured Cartesian grids. The velocity amplification at the corners is the dominant flow mechanism that controls the maximum scour depth. The results were agreement with the experimental study. Keywords-local scour; horseshoe vortex; velocity amplification


I. INTRODUCTION
When a structure is placed in the flow pattern resulting in one or more of the phenomena: the flow accelerations at the corners, the horseshoe vortex at the upstream boundary and the vortex shedding behind the caisson. The investigations of obstruction have been presented by [1]; [2]. Several formulas for estimating local scour depth have been proposed by dimensional analyses by [2] or the analytical approaches [3]. More extensive studies which occur around a bridge pier have been classified: (1) scouring parameters concerning flow or sediment characteristics [4]; (2) effect of pier characteristics on scour depth [5]; (3) estimation of maximum scour depth [4] and most of the studies and (4) temporal variation of scour depth or geometry [6]. Two-dimensional numerical analyses by [7] or three-dimensional models with the hydrostatic assumption by [8] cannot reliably compute flow fields near the structure. Stokes and Cnoidal Waves were generated using Fenton's Fourier series method. Because it is not possible to compute the flow dynamics about each individual grain of sediment and it is oftentimes difficult to compute the boundary layer at the interface, an empirical model must be used. This study is to answer the following questions: (1) What is the dominant flow mechanism that controls the maximum scour depth? Does the horseshoe vortex play a dominant role in the scour process?(2) How does the process of local scour occur around Caissons?

A. Flow Model
The governing equations for flow calculation are the Reynolds -Averaged Navier-Stokes equations

1) Momentum equations
The equations of motion for the fluid velocity components (u, v, w) in the three coordinate directions are the Navier-Stokes equations with some additional terms: In these equations:(G x , G y , G z ) are body accelerations; (f x , f y , f z ) are viscous accelerations; (bx, by, bz) are flow losses in porous media, and the final terms account for the injection of mass at a source represented by a geometry component. (4) Where: ρ is the fluid density, R SOR is a mass source. The velocity components (u, v, w) are in the coordinate directions (x, y, z). A x ; A y and A z are area fractions for flow in the x, y and z directions.

B. Sediment Scour Model
The sediment scour model estimates the motion of sediment flow by predicting the erosion, advection and deposition of sediment.
Where c s,i is the concentration of the suspended sediment, in units of mass per unit volume and u is the mean velocity of the fluid/sediment mixture. Therefore, one can write momentum balances for each sediment species and the fluidsediment mixture: International Conference on Modelling, Simulation and Applied Mathematics (MSAM 2015) Here:u s,i is the velocity of sediment species i; ρ s,i is the density of the sediment material; f s,i is the volume fraction of sediment species i; P is the pressure; K is the drag function.
The drift velocity is computed from the relative velocity using the definition of the drift and relative velocities:

C. Turbulence Transport Models
Large Eddy Simulation Model. The filter provides a formal definition of the averaging process and separates the resolvable scales from the sub grid scales. We use filtering to derive the resolvable-scale equation (Wilcox, D.C -2006). For incompressible flow, the continuity and Navier-Stokes equations assume the following form. The Large Eddy Simulation (LES) model of turbulence a rose from atmospheric modeling efforts. The basic idea is that one should directly compute all turbulent flow structures that can be resolved by the computational grid and only approximate those features that are too small to be resolved (Smagorinsky-1963). When using the LES model, it is important to remember that the model is inherently three dimensional and time-dependent. Furthermore, fluctuations must be initialized and/or input at inflow boundaries.
The filter provides a formal definition of the averaging process and separates the resolvable scales from the sub grid scales.
(11) And scales velocity fluctuations by the magnitude of L times the mean shear stress. These quantities are combined into the LES kinematic eddy viscosity (12) where c is the range of 0.1 to 0.2 and eij denotes the strain rate tensor. This kinematic eddy viscosity is incorporated into the dynamic viscosity. (13)

D. Stokes and Cnoidal Waves Generator
Cnoidal Waves is a nonlinear oscillatory Wave in shallow water, which has sharper crests and flatter troughs than a Stokes Waves. It is generated using Fenton's Fourier series method (Fenton-1999). Let d and H represent the undisturbed water depth and the Waves height (measured vertically from Waves trough to Waves crest), respectively. A stationary reference system (x, z) is established in such a way that its origin is fixed at the bottom, +x goes horizontally in the Waves propagation direction, and +z directs vertically upward (when using the code, we can place the coordinate origin at any location but +z must be in the upward direction). A current is assumed to exist and its x-component u is a constant. The Waves speed c and the angular Waves frequency ω are related to other parameters as where: λ is the waves length; k is the waves number and T is the waves period. The nonlinear equation system is iteratively solved using Newton's method. Once the unknown variables are obtained, the transient water velocity in the stationary system (x, z) is calculated using the following equations:     -Group 1: P11 Velocity (without Waves).

IV. 3D NUMERICAL SIMULATION RESULTS
The determination of stable scour depth: Stable scour depth is identified when mean kinetic energy is constant (See FIGURE. 2h). Scour process started at four upstream and downstream corners. After 3 s, two distinct scour pits around the two upstream corners were observed clearly and their depths were 3cm (FIGURE. 5a). The flow mechanism contribute to scour around a caisson is the flow accelerations at the corners (Upstream and Downstream). Local scour process around the subsea caisson was studied by 3D Numerical simulation. The mechanism of the scour was studied and the development of the stable scour depth and stable scour time was investigated. The results can be summarized as follows: (1) The velocity amplification at the corners of the caisson is the dominant flow mechanism that controls the maximum depth of scour around the structure. The waves, the horseshoe vortex plays a far less role than the velocity amplification.
(2) The results of scour process indicated that although the input parameters of the simulations and the experiments are a little different but scour progress and the intensity scour are quite similar.