A Vertical Two-Dimensional Model to Simulate Tidal Hydrodynamics in A Branched Estuary

A vertical (laterally averaged) two-dimensional hydrodynamic model is developed for tides, tidal current, and salinity in a branched estuarine system. The goveming equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. An explicit scheme is employed to solve the continuity equations. The momentum and mass balance equations are solved implicitly in the Cartesian coordinate system. The tributaries are govemed by the same dynamic equations. A control volume at the junctions is designed to conserve mass and volume transport in the finite difference schemes, based on the physical principle of continuum medium of fluid. Predictions by the developed model are compared with the analytic solutions of steady wind-driven circulatory flow and tidal flow. The model results for the velocities and water surface elevations coincide with analytic results. The model is then applied to the Tanshui River estuarine system. Detailed model calibration and verification have been conducted with measured water surface elevations,tidal current, and salinity distributions. The overall performance of the model is in qualitative agreement with the available field data. The calibrated and verified numerical model has been used to quantify the tidal prism and flushing rate in the Tanshui River-Tahan Stream, Hsintien Stream, and Keelung River.     

The numerical study of wave-induced pore water pressure response in highly permeable seabed

The coupling numerical model of wave interaction with porous medium is used to study waveinduced pore water pressure in high permeability seabed.In the model,the wave field solver is based on the two dimensional Reynolds-averaged Navier-Stokes(RANS) equations with a k-ε closure,and Forchheimer equations are adopted for flow within the porous media.By introducing a Velocity-Pressure Correction equation for the wave flow and porous flow,a highly efficient coupling between the two flows is implemented.The numerical tests are conducted to study the effects of seabed thickness,porosity,particle size and intrinsic permeability coefficient on regular wave and solitary wave-induced pore water pressure response.The results indicate that,as compared with regular wave-induced,solitary wave-induced pore water pressure has larger values and stronger action on seabed with different parameters.The results also clearly show the flow characteristics of pore water flow within seabed and water wave flow on seabed.The maximum pore water flow velocities within seabed under solitary wave action are higher than those under regular wave action.     

Wave-Induced Loads on Very Large FPSOs at Restricted Water Depth

The effects of water depth on the wave-induced vertical bending moment and shearing force on a very large FPSO are studied by experiments and computations for regular and irregular waves. The restricted water depth composite Green function is employed to develop a program for the computation of the hydrodynamic coefficients of the very large FPSO at shallow water. A thrce-segment model with 1 : 100 scale is tested in the State Key Laboratory of Ocean Engineering at Shanghai Jiao Tong University for the verification of the nmnerical method. The experimental and computational results show that the water depth has a substantial effect on wave-induced loads. The wave-induced vertical loads increase with the decrease of water depth for shallow water. Especially, for ultra-shallow water these loads increase very evidently with the decrease of water depth. The long-term prediction values of wave-induced vertical loads increase with the decrease of the ratio of water depth to draught. The long-term prediction values of wave-induced vertical loads are about 8% larger than those for deep water when the ratio of water depth to draught is 3.0. However, water depth hardly affects the longterm prediction values of wave-induced loads when the ratio of water depth to draught is larger than 5.0.     

Quasi-3D Numerical Simulation of Tidal Hydrodynamic Field

Based on the 2D horizontal plane numerical model,a quasi-3D numerical model is establishedfor coastal regions of shallow water.The characteristics of this model are that the velocity profiles can be ob-tained at the same time when the equations of the value of difference between the horizontal current velocityand its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent withthe results of the 2D model.The circulating flow in the rectangular area induced by wind is simulated and ap-plied to the tidal flow field of the radial sandbanks in the South Yellow Sea.The computational results fromthis quasi-3D model are in good agreement with analytical results and observed data.The solution of the finitedifference equations has been found to be stable,and the model is simple,effective and practical.     

Numerical Modelling of Standing Waves with Three-Dimensional Non-Linear Wave Propagation Model

Based on the theoretical high-order model with a dissipative term for non-linear and dispersive wave in water of varying depth, a 3-D mathematical model of non-linear wave propagation is presented. The model, which can be used to calculate the wave particle velocity and wave pressure, is suitable to the complicated topography whose relative depth ratio of the characteristic water depth to the characteristic wavelength in deep-water) is equal to or smaller than one. The governing equations are discretized with the improved 2-D Crank-Nicolson method in which the first-order derivatives are corrected by Taylor series expansion, .and the general boundary conditions with an arbitrary reflection coefficient and phase shift are adopted in the model. The surface elevation, horizontal and vertical velocity components and wave pressure of standing waves are numerically calculated. The results show that the numerical model can effectively simulate the complicated standing waves, and the general boundary conditions     

A finite element model of circulation in shallow water with moving boundary on tidal-flat

To examine the circulation in shallow water with tidal flat, a finite element model for the numerical solution of the shallow water equations was developed by means of standard Galerkin's method. The domain computed was covered with triangular finite elements, and water elevation and velocity were approximated by linear interpolation functions, and the lumped coefficients were used to substitute for solving the high order algebraic equation system. The time-dependent land-water boundary changes are treated mathematically by interrelating the location of the land-water boundary with the instantaneous tidal level. The implicit scheme was adopted for the terms of the bottom friction and the Coriolis effect in the motion equation so that the numerical stability of the model has been improved.The model was applied to the tidal current on shoaling water with large tidal flat off Pikou, and a comparison between observed and calculated values showed good agreement, the flow pattern being reproduced. The result     

A New Wetting and Drying Method for Moving Boundary in Shallow Water Flow Models

To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models,a new wetting and drying (WD) method is proposed.In the new method,a "predicted water depth" is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow.Compared with previous WD method,besides the water elevation,more factors,such as the flow velocity and the surface shear stress,are taken into account in the new method to determine the moving boundary.In addition,a formula is deduced to determine the threshold,as critical water depth,which needs to be preset before simulations.The new WD method is tested with five cases including three 1D ones and two 2D ones.The results show that the new WD method can simulate the wetting and drying process,in both typical and practical cases,with smooth manner and achieves effective estimation of the retention volume at shallow water body.     

A New Wetting and Drying Method for Moving Boundary in Shallow Water Flow Models

To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a “predicted water depth” is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in both typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.     

Nonlinear Effect of Wave Propagation in Shallow Water

—In this paper,a nonlinear model is presented to describe wave transformation in shallow wat-er with the zero-vorticity equation of wave-number vector and energy conservation equation.Thenonlinear effect due to an empirical dispersion relation(by Hedges)is compared with that of Dalrymple＇sdispersion relation.The model is tested against the laboratory measurements for the case of a submergedelliptical shoal on a slope beach,where both refraction and diffraction are significant.The computation re-sults,compared with those obtained through linear dispersion relation.show that the nonlinear effect ofwave transformation in shallow water is important.And the empirical dispersion relation is suitable for re-searching the nonlinearity of wave in shallow water.     

An Explicit High Resolution Scheme for Nonlinear Shallow Water Equations

The present study develops a numerical model of the two-dimensional fully nonlinear shallow water equations （NSWE） for the wave run-up on a beach. The finite volume method （FVM） is used to solve the equations, and a second-order explicit scheme is developed to improve the computation efficiency. The numerical fluxes are obtained by the two dimensional Roe＇ s flux function to overcome the errors caused by the use of one dimensional fluxes in dimension splitting methods. The high-resolution Godunov-type TVD upwind scheme is employed and a second-order accuracy is achieved based on monotonic upstream schemes for conservation laws （MUSCL） variable extrapolation; a nonlinear limiter is applied to prevent unwanted spurious oscillation. A simple but efficient technique is adopted to deal with the moving shoreline boundary. The verification of the solution technique is carried out by comparing the model output with documented results and it shows that the solution technique is robust.     

Numerical simulation for the two-dimensional nonlinear shallow water waves

This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equations. Numerical method for modelling the short waves is investigated in detail. The forces, such as Coriolis forces, wind stress, atmosphere and bottom friction, are considered. A two-dimensional implicit difference scheme of Boussinesq equations is proposed. The low-reflection outflow open boundary is suggested. By means of this model,both velocity fields of circulation current in a channel with step expansion and the wave diffraction behind a semi-infinite breakwater are computed, and the results are satisfactory.     

A Higher-Efficient Non-Hydrostatic Finite Volume Model for Strong Three-Dimensional Free Surface Flows and Sediment Transport

In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3-D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection–diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional (2-D) grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.     

Numerical simulation of solitary and randomwave propagation through vegetation based on VOFmethod

A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.     

COMPUTATION OF DILUTION OF RADIOACTIVE WASTE AND HEAT DISCHARGED FROM NUCLEAR POWER PLANT

- In this paper a two-dimensional unsteady shallow water equations and convection-diffusion equations have been considered to describe the diluting process of the radioactive and heat discharged from a nuclear power plant. The theory of characteristic and eccentric difference scheme has been used.In order to obtain the distribution of the concentration in far-field and near-field, three different siges of mesh have been used. The flow field has been verified with the field data, and the computed temperature in the near-field agrees with the measurements in the normal physical model test.     

Investigation of Scouring Mechanism Around A Circular Pier

In the present study,the flow field around a circular pier is investigated with experimental measurements and numerical simulations.The transient flow patterns during erosion are studied in detail.The results show that the traditional equations of particle motion are not perfect for the calculation of the sand motion under this complex flow situation.The scouring agents,such as turbulent intensity,the fluctuating pressure and the vertical pressure gradient,having many effects on the sand motion with the increasing scouring depth,need to be considered in modifying the traditional model.     

Effects of Topography and Current on Horizontal Irrotational Waves in Shallow Water

Based on the Boussinesq assumption,derived are couple equations of free surface elevationand horizontal velocities for horizontal irrotational flow,and analytical expressions of the correspondingpressure and vertical velocity.After the free surface elevation and horizontal velocity at a certain depth areobtained by numerical method,the pressure and vertical velocity distributions can be obtained by simplecalculation.The dispersion at different depths is the same at the O(ε)approximation.The waveamplitude will decrease with increasing time due to viscosity,but it will increase due to the matching ofviscosity and the bed slope.thus,flow is unstable.Numerical or analytical results show that the waveamplitude.velocity and length will increase as the current increases along the wave direction.but theamplitude will increase.and the wave velocity and length will decrease as the water depth decreases.     

Three Dimensional Baroclinic Numerical Model for Simulating Fresh and Salt Water Mixing in the Yangtze Estuary

For simulating fresh and salt water mixing in estuaries, a three dimensional nonlinear baroclinic numerical model is developed, in which the gradients of horizontal pressure contain the gradient of barotropic pressure arising from the gradi-ent of tidal level and the gradient of baroclinic pressure due to the gradient of salinity. The Eulerian-Lagrangian method is employed to descretize both the momentum equations of tidal motion and the equation of salt water diffusion so as to im-prove the computational stability and accuracy. The methods to provide the boundary conditions and the initial conditions are proposed, and the criterion for computational stability of the salinity fields is presented. The present model is used for modeling fresh and salt water mixing in the Yangtze Estuary. Computations show that the salinity distribution has the characteristics of partial mixing pattern, and that the present model is suitable for simulalion of fresh and salt waler mixing in ihe Yanglze Esluary.     

Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.     

A New Approach to High-Order Boussinesq-Type Equations with Ambient Currents

A new approach to high-order Boussinesq-type equations with ambient currents is presented. The current velocity is assumed to be uniform over depth and of the same magnitude as the shallow water wave celerity. The wave velocity field is expressed in terms of the horizontal and vertical wave velocity components at an arbitrary water depth level. Linear operators are introduced to improve the accuracy of the kinematic condition at the sea bottom. The dynamic and kinematic conditions at the free surface are expressed in terms of wave velocity variables defined directly on the free surface. The new equations provide high accuracy of linear properties as well as nonlinear properties from shallow to deep water, and extend the applicable range of relative water depth in the case of opposing currents.     

A finite difference approximation for dynamic calculation of vertical free hanging slender risers in re-entry application

The dynamic calculations of slender marine risers, such as Finite Element Method (FEM) or Modal Expansion Solution Method (MESM), are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation (FDA) method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bernoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser’s free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.