Proceedings of the 9th International Conference

on Hydrodynamics (ICHD 2010)

October 11 – 15, 2010,Shanghai, China
Editors: Y. S. Wu, S. Q. Dai, H. Liu, L. D. Zhou, X. C. Yang

YSHE Symposium (32)

I1-1 Three-dimensional hydroelasciticity analysis of acoustic responses of ship structures

Ming-song Zhou, You-sheng Wu , Yong-lin Ye

China Ship Scientific Research Center, Wuxi, China

Abstract :

The three-dimensional hydroelasticity theory of ships[1-3]was extended to include the effect of fluid compressibility. This enables the acoustic responses of a ship structure induced by the machinery or wave excitations being predicted with the inclusion of the free surface and the forward speed effect. A method for eliminating the irregular frequencies in the numerical analysis is proposed. The acoustic radiation of an elastic spherical shell was predicted and compared with analytical solution for validation. Illustrated in this paper are also the acoustic radiations of a submerged ring-stiffened cylindrical hull in the near field and far field.

I2-1 Similarity law of flowrate for hot-gas ventilated supercavity flows

Jiu-xi Chen

China Ship Scientific Research Center, Wuxi, China

Abstract :

Based on the complete sets of similarity criteria for modelling experiments of  exhausted-hot gas ventilated supercavity flows which had been derived by author before, considering the effects of phase changes, a similarity law of flowrate for ventilated hot-gas supercavity flows is found and described in this paper. The effects of phase changes on modelling scheme and ventilated-gas flowrate are studied. The phase changes are the vaporization of liquid water on cavity wall, and condensation of vapor contained in the flows of gaseous mixture within cavity.

I2-2 A summary of water entry problem of a wedge based on the fully nonlinear velocity potential theory

Guo-Xiong Wu, Guo-Dong Xu,Wen-Yang Duan

Department of Mechanical Engineering, University College London, UK

Abstract :

This paper presents a comprehensive study of water entry of a wedge through numerical simulation. The mathematical model is based on the assumptions of inviscid and incompressible liquid and irrotational flow, while the solution procedure is based on the boundary element method for the complex velocity potential. In particular, the problem is solved in the stretched coordinate system. The wedges considered include a symmetrical wedge, an asymmetrical wedge and twin wedges, and the motions considered include vertical water entry and oblique water entry at prescribed velocity, and free fall motion in one and three degrees of freedom.

I2-3 On the gas leakage way of supercavity and vehicle vibration

Jing-jun Zhou, Kai-ping Yu, Ming Yang , Xiao-hui Wan

School of Aeronautics, Harbin Institute of Technology, Harbin, China

Abstract :

Steady and unsteady three-dimension numerical simulations were carried out to ventilated supercavitation by solving Navier Stokes equations with finite volume method based on the two-fluid multiphase flow model and SST (Shear Stress Transport) and DES (Detached Eddy Simulation Model) turbulence model were adopted respectively. The cavity developing process with re-entrant jet flow and two vortex tube gasleakage way were studied and the change curve of pressure at a fixed point in cavity with time was given. On the other hand, the dynamics equation combined with moving technology was used to control the oscillation of the vehicle in the cavity, the change law of drag, lift, velocity as well as the displacement were given and analyzed.

I2-4 Numerical simulation of thewater-entry of body based on the Lattice Boltzmann method

Ke Zhang, Kai Yan, Xue-sen Chu, Guan-yi Chen

China Ship Scientific Research Center, Wuxi, China

Abstract :

In this paper, a lattice Boltzmann single-phase model is implemented to the simulation of the process of body entering water vertically. The initial water-entry process of a two dimensional cylinder and the process of a three dimensional disk entering water vertically at constant speed is simulated. The deformation of free surface and the relationship between the relative closure depth of water-entry cavity of the disk and the Froude number is studied. It is demonstrated that the LBM single-phase method is applicable to the water-entry problems.

I2-5 Numerical simulation of water-exit of a cylinder with cavities

Xue-sen Chu, Kai Yan, Zhi Wang, Ke Zhang, Guang Feng, Wei-qi Chen

China Ship Scientific Research Center, Wuxi, China

Abstract :

The water-exit process of a cylindrical body was simulated in this paper. In the simulation, cavities were formed on the nose and at the tail of the body when the body approaches to the free surface, and then collapsed during the passage of body through the water surface. The interaction of free surface with cavity were simulated, and the high pressure due to collapse of cavities were captured.

I2-6 A short review on the homotopy analysis method in fluid mechanics

Sh-ijun Liao

State Key Lab of Ocean Engineering, Shanghai Jiaotong University, Shanghai , China

Abstract :

We give a short review on the current development of homotopy analysis method (HAM), an analytic technique for strongly nonlinear problems, and its applications in fluid mechanics.

I3-1 From water entry to lock entry

Xue-nong Chen

Institute for Nuclear and Energy Technologies, Karlsruhe Institute of Technology, P.O. Box 3640, D-76021 Karlsruhe, Germany

Abstract :

In this paper, the ship lock entry problem is studied, which is physically similar to, but numerically different from, the water entry problem simulated by the author 25 years ago under the instruction of Prof. He. A one dimensional unsteady hydraulic narrow-channel model for the flow coupled to the ship's motion in surge, heave and pitch is formulated and numerically implemented. The calculated ship motions were validated by comparison with model experiments carried out in the Duisburg shallow water tank. The viscous effects are further taken into account in order to investigate the scale effect in the experimental modeling in a framework of new lock project of Panama Canal. Thus, the experimentally validated model is applied to optimize the lock entry time by changing the trust course in the full-scale case.

I3-2 Numerical study of the cavitating flows over underwater vehicle with large angle of attack

YingChen, Chuan-jing Lu , Jian-hongGuo

Department of Engineering Mechanics,School of Naval Architecture, Ocean and Civil Engineering, Shanghai, Jiaotong University, Shanghai, China

Abstract :

A Homogenous-Equilibrium-Model based cavitation code was developed to numerically analyze the three-dimensional cavitating flows over underwater vehicle navigating with large angle of attack. Rayleigh Plesset equation based cavitation model was used together with non-linear eddy-viscosity turbulence model. The computed cavity shapes and pressure distributions along the vehicle body werefound to generally accord with experimental results at different cavitation numbers and angles of attack. The forces acting on the vehicle body was studied qualitatively to explain why the body may get damaged during navigation. It was also discovered interestingly that, the variation trend of dragforce along with cavitation number at the conditions of large angle was completely opposite to that at zero angle.

I3-3 Numerical study of periodically forced-pitching of a supercavitating vehicle

Zhan-cheng Pan , Chuan-jing Lu , Ying Chen, Shi-liang Hu

Department of Engineering Mechanics, School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiaotong University, Shanghai, China

Abstract :

The unsteady behaviors, such as surging, heaving and pitching motion, which often occur during the advancing of supercavitating vehicle, has significant effect on the stability of supercavitaty and the trajectory of the vehicle. This paper presents a 3-dimentional numerical simulation of periodically forced-pitching of supercavitating vehicle. Based on the finite volume method and the pressure-based segregate algorithm, in the framework of Mixture multiphase model, associated with dynamic mesh method, the Reynolds-Averaged Navier-Stokes equations are solved for the ventilated cavitating flow field in a cavitation tunnel. For both steady-state and dynamic cases, the numerical results agree with the experimental results very well. When the vehicle is pitching periodically, the profile of the supercavity doesn’t vary significantly. The pressure inside the cavity fluctuates slightly during the pitching motion, while the pressure fluctuates significantly at the rear of the down-line of the vehicle. The hydrodynamic forces of the vehicle oscillated periodically, but not linear related to the pitching motion.

I3-4 Simulation of unsteady artifiсial supercavities

Guang Feng, Wei-zheng Chen, Xue-sen Chu, Zhi Wang, Ming-hui Zhang, Wei-qi Chen

China Ship Scientific Research Center, Wuxi, China

Abstract :

The ventilated cavity in acceleration process for a supercavitating vehicle was simulated in this paper based on the independent expansion principle of cavity section. The unsteady artificial super-cavities were simulated under the condition of a straight navigation at the fixed depth. The hydrodynamic schemes of gas leakage from the cavity, the gravity effect and the angle of attack of the cavitator were taken into account in the simulation. By calculation for the supercavity of the acceleration process, a good ventilating rule was summarized as follows: a large ventilating rate lasts about 0.3s, then the ventilating rate should be reduced to the rate corresponding to the cruising state.

I3-5 Application of SPH method on free surface flows on GPU

Ben-long Wang, Hua Liu

Engineering Mechanics, Shanghai Jiao Tong University, Shanghai, China

Abstract :

A parallel SPH solver is developed using GPU for free surface flows. As a kind of Lagrange methods, fluid particles are tracked by integrating the momentum and kinematics equations. Mass conservation is automatically satisfied. Therefore SPH is suitable for simulating the flows with complex free surface. With this 3D flow solver, flows after dam-breaking are simulated with good performance using GPU.

I3-6 Rheological characterizations and extrudate swell predictions of an LDPE melt by two KBKZ-type constitutive equations

Shu-xin Huang, Xin Chen, Chuan-jing Lu, etc.

Department of Engineering Mechanics, Shanghai Jiao Tong University, Shanghai, China

Abstract :

The shear viscoelastic properties of an LDPE (1I2A-1) melt at 160 ℃were described by using two KBKZ integral type constitutive equations, i.e. Wagner model and the three- parameter rational type model. The present characterizations were compared with the previous calculations by using PSM model. The results show that there are differences between the descriptions on the stress relaxation property in step strain and the first normal stress difference for the three models, and however, all the three models show the almost same predictions on the shear-thinning viscosity of the ldpe melt. Both Wagner model and PSM model are good choice in describing the rheological properties of the melt due to theirs simplicity. The Tanner theory combined with the Wagner model or the three-parameter rational type model was used to predict the swell ratio of the LDPE melt through a long capillary die, which was also compared with the experimental results and the previous calculations on swell. The predictions by Tanner theory combined with the three different constitutive models on the swell ratio of the LDPE melt are similar, which are apparently lower than the experiments and numerical simulations over the whole experimental shear rate range. The results further indicate that Tanner theory cannot predict the swell ratio of the LDPE melt precisely.


I4-1 Orientation distribution and rheological properties of fiber suspensions flowing through curved expansion and rotating ducts

Qi-hua Zhang, Jian-zhong Lin

School of Aeronautics and Astronautics, Zhejiang University, Hangzhou, China

Abstract :

In dilute fiber suspensions, the Jeffery equation is employed to predict the fiber orientation distribution through a curved expansion and rotating duct. The results show that the fiber flips more quickly in the inlet region, especially in the region close to the concave wall. In the central region and the downstream region, the fiber orientation distributions are more uniform. In semi-dilute fiber suspension, the fiber interactions can be represented by a statistical model developed by Folgar & Tucker. Then a deduced Fokker-Planck equation can be used to describe the fiber orientation distribution probability. In the present work, the Fokker-Planck equation is directly solved by FVM. Then with the probability, the fiber orientation tensors are obtained. Furthermore, the fiber extra stresses are gained with the Batchelor model. The results show that the shear stress and normal stress difference are concentrated around the inlet close to the concave wall regions. In the central regions, these properties are less obvious. And in the downstream, these properties are negligible.


I4-2 Calculations of massive separationaround landing-gear-like geometries

Zhi-xiang Xiao, Jian Liu, Song Fu

Laboratoire de Mécanique des Fluides Ecole Centrale de Nantes, Nantes, France.

Abstract :

The massive separating flows around landing gear-like configurations, such as Rudimentary landing gear (RLG) and tandem cylinders (TC), are calculated using unsteady Reynolds averaged Navier-Stokes (URANS) and delayed-detached-eddy simulation (DDES) based on k-ω-SST model. A number of numerical schemes and dissipation are applied in an effort to compare the averaged as well as the instantaneous flow-fields with the available measurements. It is shown that high-order and low-dissipation scheme is necessary to calculate the small scale structures. Furthermore, DDES delivered better results than URANS.

I4-3 Verification and validation for RANSsimulation of KCS container ship without/with propeller

Zhi-rong Zhang

China Ship Scientific Research Center, Wuxi, China

Abstract :

The free surface flow of a modern container ship KCS without propeller was firstly simulated using three sets of grids. The computed results including resistance, wave elevation and flow field on propeller disk were compared with the experimental data in detail. Verification and validation of resistance and wave profile were performed using recommended procedures proposed by ITTC. Then the viscous flow around KCS with operating propeller behind was also simulated. Both body force approach and sliding mesh approach were applied to consider for the effect of propeller. The results of these two approaches were compared with the measured data. These numerical investigation shows that accurate prediction of propeller/hull interaction using CFD method is becoming feasible and the huge potential of CFD application in ship hydrodynamics performance prediction is demonstrated.

I4-4 Tunnel hood effects on high speed train-tunnel compression wave

Xin-tao Xiang, Lei-ping Xue

School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai China

Abstract :

This paper reports the numerical research of tunnel hood effects on high speed train-tunnel compression wave. The three-dimensional simulation with real geometry is carried out by the implementation of a commercial computational code. The train speed is 350 km/h. The train/tunnel blockage ratio is 0.115. Nine different types of tunnel hoods were studied. The calculation results showed that the hood length, the hood cross sectional area and the ventilation holes might have significant influence on the first compression wave, and inclined entry or asymmetric distribution of the ventilation holes is not available for alleviating the impulsive wave.


I4-5 Simulation of acoustic scattering by the fast BEM approach

Zai-you Yan

Department of aerodynamics, Nanjing University of Aeronautics & Astronautics, Nanjing, China

Abstract :

Boundary element method accelerated by the precorrected-FFT (pFFT) algorithm is developed and implemented for acoustic scattering problems. Because models with very fine meshes can be computed using this method, high frequency responses can be simulated. In the implementation, constant triangular elements are applied. An example of a plane acoustic wave scattering by a rigid sphere is simulated to validate the in-house fast BEM code for acoustic scattering problems. Numerical results are compared with the corresponding analytical solutions and the solutions by the conventional boundary element method. As a result, the pFFT accelerated boundary element method greatly improve the computational efficiencyand reduce the requirement of storage largely. The accuracy by this method is comparable to that by the conventional boundary element method as the grid-order is no less than 4.


I4-6 Vortex-induced vibration on 2D circular riser using a high resolution numerical scheme

Jia-song Wang, Hua Liu, Shi-quan Jiang , Liang-bin Xu,, Peng-liang Zhao

School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai, China

Abstract :

This paper presents a high resolution numerical method for vortex induced vibration (VIV) simulation on the fluid structure interaction (FSI) of circular cylinder which represents a two dimensional marine riser. For the VIV case, the cylinder is elastically mounted and is modeled as a spring mass oscillation system. Based on a new proposed elemental velocity vector transformation (EVVT) method, a finite-volume total variation diminishing (TVD) approach developed recently for solving unsteady Reynolds Averaged Navier-Stokes (URANS) equation with the RNG turbulence model was used to simulate the key hydrodynamic parameters such as lift coefficients. The four-stage Runge-Kutta method is used to solve the dynamic response equation of the structure. The FSI prediction results are compared with the available experimental data and showed a good agreement in a wide range of Reynolds number, which provide a good picture of real physics of phenomenon including the Karman vortex streets with different vortex modes with regard to the reduced velocities.


I5-1 Dwarf solitary waves and low tsunamis

Sunao Murashige, Theodore Yaotsu Wu

School of Systems Information Science, Future University Hakodate, Hakodate, Hokkaido, Japan

Abstract :

This work applies the regularized solitary wave theory to develop accurate computational method for evaluating the dwarf solitary waves, with amplitude-to-water depth ratio α≤10−2, as a useful model of one-dimensional tsunamis propagating in the open ocean. The algebraic branch singularities of these solitary waves magnifying with diminishing wave amplitude,  making their computations insurmountable by existing methods, are removed by the regularized coordinates given by this new theory. Numerical examples show that this new method can produce accurate results even for α  ~= 10−3 or less.


I5-2 Examination on low-frequency QTF of a platform

Bin Teng, Pei-wen Cong, Ying Gou

State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China

Abstract :

For deep water offshore structures, the second order different frequency force may induce large responses in the horizontal directions. The second order different frequency force is usually computed by the Fourier transformation of the second order quadratic transform function (QTF), which is hard to be obtained due to its complexity, time consuming and the influence of ‘irregular frequency’ in its computation. In practical application Newman’s approximation[5] is widely used, without considering its viability. Based on a higher-order boundary element method, a numerical model is established for complex solution of second-order diffractions and radiations at different frequencies in bichromatic incident waves. The second-order quadratic transform functions are obtained by the integration of second order fluid pressure on the body surface. With application of Teng et al’s method[9], full QTF’s are computed for a typical TLP model ISSC. Computation shows that the second-order quadratic functions have distinct distribution at lower frequency ranges. Comparisons are made on the Newman’s approximation with the full difference-frequency QTF.


I5-3 On the simulation of highly nonlinear wave-breakwater interactions

Chi Yang , Hai-dong Lu, Rainald Löhner

Department of Computational and Data Sciences, George Mason University, Fairfax, Virginia, USA

Abstract :

A numerical time domain simulation model has been developed to study the highly nonlinear interactions between waves and rubble mound breakwaters. In this model, a volume of fluid (VOF) technique isused to capture the violent free surface motion. The incompressible Euler/Navier-Stokes equations, written in an arbitrary Lagrangian-Eulerian (ALE) frame, are solved using projection schemes and a finite element method on unstructured grids. A general advancing front technique for filling space with arbitrary separated objects is developed to model the rubbles that are laid down on the sloped surface of the breakwater in a random way. Three case studies are performed to study the effects of rubbles and rubble types on the wave dissipation and wave overtopping.


I5-4 Numerical simulation of wave transformation incorporating porous media wave absorber

Jie-min Zhan, Zhi Dong, Yi Han, Wei Jiang

Department of Applied Mechanics and Engineering, Sun Yat-sen University, Guangzhou, China

Abstract :

Wave propagation over arbitrary depth and interaction with structures are of fundamental importance in coastal and ocean engineering. Accurate numerical simulations providing more detailed information than physical experiments can enhance the understanding of such problems. Early numerical studies are mostly based on potential flow theory and shallow water equations, in which the fluid viscous is not taken into account. Recently with the development of computer technology and CFD methods, it is desirable to directly solve the Navier-Stokes equation. In this paper, a numerical wave tank is established with the classical Navier-Stokes equations and the VOF method. An effective numerical method for wave absorbing utilizing the energy-dissipating property of porous media is proposed. The applicability of the presented viscous models is assessed in detail through the two test cases of wave diffraction around a semi-infinite breakwater and wave focusing by bottom topography. The computed results agree well with experimental data.


I5-5 Numerical simulation of wave impact on the slab

Ya-mei Lan, Wen-hua Guo, Hua Liu, Qiu-hong Song , Jun-ting Yuan

College of Engineering Science & Technology, Shanghai Ocean University, Shanghai, China

Abstract :

In the paper, the wave generating and absorbing were introduced into the RANS equations as the source terms using the relaxation approach. A new module of the wave generating and absorbing function, which is suitable for FLUENT based on the volume of fluid method (VOF), was established. The computational results of the wave pressures on the bottom of the slab were validated for the different relative clearance by the experimental data. Good agreements were found. It is suggested that the numerical wave flume is capable of simulating the phenomenon of wave impact, especially the process of the negative pressure.


I5-6 Development and application of high order Boussinesq model in tsunami studies

Hua Liu, Xi Zhao , Ben-long Wang

Department of Engineering Mechanics, Shanghai Jiao Tong University,Shanghai, China.

Abstract :

The recent progress in development and application of the High order Boussinesq-type equations for numerical simulation of tsunamis is presented. The moving bottom is implemented to simulate the generation of tsunami induced by submarine earthquakes. The wave patterns of tsunami near the source are obtained for the cases of different magnitude of earthquake and the depth of ocean, which confirms that the leading depression N-wave occurs when the earthquake magnitude is large enough. The numerical model is applied in modeling the India Ocean Tsunami in 2004. The computed value of the runup of tsunami is compared with the data from field survey available.


I6-1 Hydrodynamics of coupled river-lake system

Lin-lin Li, Xi-ping Yu

Earth Observatory of Singapore, Nangyang Technological University, Singapore

Abstract :

A numerical model based on the finite difference solution of Saint-Venant equation for one dimensional unsteady river flows and the shallow water equation for vertically integrated two-dimensional flows in the lake is employed to study the propagation of flood waves in Jingjiang River and Dongting Lake. Three different types of flood occurred in 1996, 1999 and 2002, respectively, are computed. All numerical results are shown in fairly good agreement with observed data.


I6-2 Large-eddy simulation of stratification effects on dispersion in urban environments

Zheng-tong Xie

School of Engineering Sciences, University of Southampton,Southampton, UK

Abstract :

This paper focuses on investigating thermal buoyancy effects on dispersion of approach flows and local heat transfer from/to buildings by taking the DAPPLE site as the test case but for varying wind direction and magnitude. Firstly, based on the BT tower data, only weakly unstable conditions (i.e. Richardson number |Ri|<1) of approach flows were considered, with adiabatic boundary conditions at the ground and building surfaces. It was found that the predicted concentration in the far field was significantly reduced. Secondly, numerical experiments were performed, i.e. using various Richardson numbers. It was found that the modelled mean concentration for bulk Richardson number -0.9 was in reasonable agreement with field data at all stations. Thirdly, to understand the effects of weakly local unstable conditions, nonadiabatic boundary conditions at ground and building surfaces were implemented as well as weakly unstable conditions in the approach flows. No significant difference was found comparing with those using adiabatic boundary conditions. These simulations suggest that stratification effects on dispersion in weakly unstable conditions (e.g. in London) are not negligible.


I6-3 Environmental fluid dynamics-jet flow

Xi-kun Wang, Soon Keat Tan

School of Civil and Environmental Engineering, DHI-NTU Centre

Abstract :

Jet flow is a very important research subject in both fundamental fluid dynamics and engineering applications. Jet flow has the essences of fluid dynamics, such as free and wall-bounded shear flows, turbulent flow, eddy and large vortical structures and their stability and control, and so forth. This article serves as an overview of our past and ongoing research activities on various types of jet flow, with particular reference to their application in the field of environmental fluid dynamics. The research objectives, approach, results and their engineering implications of each topic will be presented.


I6-4 An unstructured grid hydrodynamicand sediment transport model for Changjiang Estuary

Ding-man Qi, Gang-feng Ma , Feng-feng Gu, Lin Mou

Estuarine and Coastal Scientific Research Center, Shanghai 201201, China

Abstract :

An unstructured grid hydrodynamic and sediment transport model for Changjiang Estuary is developed in the current paper. The model employs finite volume method to discretize the governing equations. Semi-implicit method originally developed by Casulli is utilized to remove the stability limitations associated with the surface gravity wave. A wetting and drying (WAD) scheme is proposed to account for the moving boundary at the shoals and tidal flats. The model is used to investigate the hydrodynamics and sediment transport in the Changjiang Estuary. Comparisons with the measured data show that the model can predict water level and tidal current very well. The variations of sediment concentration are also reasonably captured by the model.


I6-5 Research on the hydraulic characteristics of gap closing with large framed cages

Yong-jin Lu, Xiao-tao Du

Shanghai Water Engineering Design & Research Institute,Shanghai, China

Abstract :

Gap closing under high current velocity is a Gordian knot in water conservancy project, jackstone gap closing with large frame cages is a new technology to solve this problem. Associated with practice in the project of Qingcaosha reservoir and analysis on hydraulic numerical simulation and physical experiment, several hydraulic characteristic rules of this technology were discovered, it is of instructional value in practices of similar projects.


I6-6 Unsteady behaviors of two moving bodies through inviscid unbounded liquid

Run Sun, Hong-yu Ma

Department of Engineering Mechanics, Shanghai Jiao Tong University, Shanghai, China

Abstract :

We account both bodies and liquid as a conservative system to explore unsteady motions of two bodies, adopting the Lagrange equation of motion in vector form. Theoretical analysis and numerical results reveal that as two circles initially revolving around each other in an inviscid liquid with a stationary centroid are set to translate, their centroidal trajectory is sinuous and would gradually deviate from the original moving direction. For the rectilinear motion of the bodies in tandem, Galilean relativity is invalid for just two bodies during the hydrodynamic interaction, and their centroid, after initially set in motion, may accelerate and then decelerate, finally approaching a constant speed as far apart from each other.


I6-7 Currents induced by waves in the surf zone and the pollutant transport analysis

Jing-xin Zhang , Hua Liu

School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai, China

Abstract :

As a long-term hydrodynamic phenomenon induced by waves, the currents dominate the pollutant transport in the shoaling and surf zones. The present study was focused on the analysis of the cross-shore currents induced by waves and the pollutant transport. The wave parameters were obtained by solving one phase-averaged wave model including wave shoaling and breaking in the surf zone. The effect of the breaking wave rollers in the surf zone was considered by specifying an additional boundary condition on the mean water level (MWL). We introduced a formula for the vertical varied wave radiation stress into the hydrodynamic model. Furtherly, we simulated and analyzed the pollutant transport in the domain for different pollutant release position.