Flume experiments play a pivotal role in studying wave propagation,with wave elements typically assumed to remain constant in the perpendicular direction.However,evident cross wave phenomena were observed within flume...Flume experiments play a pivotal role in studying wave propagation,with wave elements typically assumed to remain constant in the perpendicular direction.However,evident cross wave phenomena were observed within flumes under certain conditions.This paper presents new analytical solutions for both primary and cross waves on double shoals in a flume via linear shallow-water equations,which may be used to idealize dynamic experimental configurations of coral reefs.The primary waves on double shoals are described by the associated Legendre functions,whereas the ultimate solutions are derived by considering the incident and reflected waves in front of a bathymetry and the transmitted waves positioned behind it.The effects of the angular frequency and topographic parameters on the primary waves are subsequently analyzed.Cross waves on double shoals constitute a type of topographically trapped wave whose solutions are formulated by combining the first and second types of the associated Legendre functions.The angular frequency is not only determined by the wavenumber but also influenced by the topographic parameters.Numerical experiments are conducted to investigate the generation mechanism of cross waves.The consistency between the numerical results and analytical solutions confirms the validity of the new analytical framework of cross waves on double shoals.展开更多
The paper considers a catastrophic event-the eruption of Hunga Tonga-Hunga Ha’apai volcano on January 15,2022.The process of preparation and eruption of Hunga Tonga volcano generated tsunami waves that were observed ...The paper considers a catastrophic event-the eruption of Hunga Tonga-Hunga Ha’apai volcano on January 15,2022.The process of preparation and eruption of Hunga Tonga volcano generated tsunami waves that were observed throughout the World Ocean.This event was notable for its unprecedented global impact and the early appearance of tsunami waves at distant coastal stations.So,the first waves at tide gauge stations in Chile and Peru were recorded 4 hours earlier than the arrival time of tsunami waves to the tide gauge after the eruption of Tonga volcano.Two mechanisms are possible for the generation of early tsunami waves:acoustic Lamb waves generated by a volcanic explosion and submarine landslides that occurred on the slopes of the volcano during the preparatory phase of the eruption.In this study,numerical simulation of various pre-eruption landslide scenarios on the slope of Hunga Tonga volcano is carried out in an attempt to explain these early tsunami waves.Under computation the elastoplastic model of landslide was taken into account.Wave characteristics of a tsunami on the coast of Chile and Peru generated by a landslide process on a volcanic slope are obtained.A detailed comparison of virtual tide gauge data with observational ones is used to validate this model.The results obtained can be used to improve early warning systems.展开更多
A three-dimensional numerical model of sand wave dynamics,incorporating the interaction of currents and waves at various angles,has been developed using the Regional Ocean Modeling System(ROMS).This model accounts for...A three-dimensional numerical model of sand wave dynamics,incorporating the interaction of currents and waves at various angles,has been developed using the Regional Ocean Modeling System(ROMS).This model accounts for both bedload and suspended load sediment transport under combined waves and current conditions.The investigation examines the influence of several key parameters,including the rotation angle of sand waves relative to the main current,tidal current velocity amplitude,residual current,water depth,wave height,wave period,and wave direction,on sand wave evolution.The growth rate and migration rate of sand waves decrease as their rotation angle increases.For rotation angles smaller than 15°,sand wave evolution can be effectively simulated by a vertical 2D model with an error within 10%.The numerical results demonstrate that variations in tidal current velocity amplitude or residual current affect both vertical growth and horizontal migration of sand waves.As tidal current velocity amplitude and residual current increase,the growth rate initially rises to a maximum before decreasing.The migration rate shows a consistent increase with increasing tidal current amplitude and residual current.Under combined waves and current,both growth and migration rates decrease as water depth increases.With increasing wave height and period,the growth rate and migration rate initially rise to maximum values before declining,while showing a consistent increase with wave height and period.The change rate of sand waves reaches its maximum when wave propagation aligns parallel to tidal currents,and reaches its minimum when wave propagation is perpendicular to the currents.This phenomenon can be explained by the fluctuation of total bed shear stress relative to the angle of interaction between waves and current.展开更多
With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ...With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.展开更多
The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general condit...The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.展开更多
A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational r...A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.展开更多
Numerical simulations based on a high-resolution three-dimensional MIT general circulation model(MITgcm)using realistic topography and tidal forcing are conducted to investigate the generation and propagation of the s...Numerical simulations based on a high-resolution three-dimensional MIT general circulation model(MITgcm)using realistic topography and tidal forcing are conducted to investigate the generation and propagation of the so-called type-a waves(large-amplitude rank-ordered wave packets)and type-b waves(isolated wave packets)in the northern South China Sea.At first,we summarized and analyzed the generation and propagation characteristics of these waves.Then,energy budget at the Luzon Strait is calculated.Energy generation has three local maxima every day,of which the largest one corresponds to the emergence of the type-a wave.Energy flux at the west boundary of the Luzon Strait shows two local maxima each day.The larger one is consistent with the generation of the type-a wave and the smaller one is in correspondence with the generation of the type-b wave.Sensitivity experiments are designed to explore the role of the east and west ridge of the Luzon Strait on the generation and propagation of the type-a and type-b waves.It is found that the east ridge is indispensable on the generation of the type-a wave while the west ridge has little contribution.The west ridge diminishes the type-a waves'amplitude but hardly changes their propagation speed.The type-b waves also come from perturbation signals which originate from the east ridge and are enhanced in amplitude and reduced in propagation speed by the west ridge.展开更多
Prediction of coastal sediment transport is of particularly importance for analyzing coast erosion accurately and solving the corresponding coast protection engineering problems.The present study provided a numerical ...Prediction of coastal sediment transport is of particularly importance for analyzing coast erosion accurately and solving the corresponding coast protection engineering problems.The present study provided a numerical scheme for sediment transport in coastal waves and wave-induced currents.In the scheme,the sand transport model was implemented with wave refraction-diffraction model and near-shore current model.Coastal water wave was simulated by using the parabolic mild-slope equation in which wave refraction,diffraction and breaking effects are considered.Wave-induced current was simulated by using the nonlinear shallow water equations in which wave provides radiation stresses for driving current.Then,sediment transport in waves and wave-induced currents was simulated by using the two-dimensional suspended sediment transport equations for suspended sediment and the bed-load transport equation for bed load.The numerical scheme was validated by experiment results from the Large-scale Sediment Transport Facility at the US Army Corps of Engineer Research and Development Center in Vicksburg.The numerical results showed that the present scheme is an effective tool for modeling coastal sediment transport in waves and near-shore currents.展开更多
This paper introduces a numerical model for studying the evolution of a periodic wave train, shoaling, and breaking in surf zone. The model can solve the Reynolds averaged Navier-Stokes （RANS） equations for a mean f...This paper introduces a numerical model for studying the evolution of a periodic wave train, shoaling, and breaking in surf zone. The model can solve the Reynolds averaged Navier-Stokes （RANS） equations for a mean flow, and the k-ε equations for turbulence kinetic energy k and turbulence dissipation rate ε. To track a free surface, the volume of fluid （VOF） function, satisfying the advection equation was introduced. In the numerical treatment, third-order upwind difference scheme was applied to the convection terms of the RANS equations in order to reduce the effect of numerical viscosity. The shoaling and breaking processes of a periodic wave train on gently sloping beaches were modeled. The computed wave heights of a sloping beach and the distribution of breaking wave pressure on a vertical wall were compared with laboratory data.展开更多
In this paper, the water waves and wave-induced longshore currents in Obak6y coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic ...In this paper, the water waves and wave-induced longshore currents in Obak6y coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic mild-slope equation for coastal water waves and the nonlinear shallow water equation for the wave-induced currents. The wave transformation under the effects of shoaling, refraction, diffraction and breaking is considered, and the wave provides radiation stresses for driving currents in the model. The numerical results for the water wave-induced longshore currents were validated by the measured data to demonstrate the efficiency of the numerical model. Then the water waves and longshore currents induced by the waves from main directions were numerically simulated and analyzed based on the numerical results. The numerical results show that the movement of the longshore currents was different while the wave proDaRated to a coastal zone from different directions.展开更多
A two-dimensional （2D） numerical model is developed for the wave sim- ulation and propagation in a wave flume. The fluid flow is assumed to be viscous and incompressible, and the Navier-Stokes and continuity equatio...A two-dimensional （2D） numerical model is developed for the wave sim- ulation and propagation in a wave flume. The fluid flow is assumed to be viscous and incompressible, and the Navier-Stokes and continuity equations are used as the governing equations. The standard k-e model is used to model the turbulent flow. The Navier- Stokes equations are discretized using the staggered grid finite difference method and solved by the simplified marker and cell （SMAC） method. Waves are generated and propagated using a piston type wave maker. An open boundary condition is used at the end of the numerical flume. Some standard tests, such as the lid-driven cavity, the constant unidirectional velocity field, the shearing flow, and the dam-break on the dry bed, are performed to valid the model. To demonstrate the capability and accuracy of the present method, the results of generated waves are compared with available wave theories. Finally, the clustering technique （CT） is used for the mesh generation, and the best condition is suggested.展开更多
The author’s combined numerical model consisting of a third generation shallow water wave model and a 3 D tide surge model with wave dependent surface wind stress were used to study the influence of waves on tide sur...The author’s combined numerical model consisting of a third generation shallow water wave model and a 3 D tide surge model with wave dependent surface wind stress were used to study the influence of waves on tide surge motion. For the typical weather case, in this study, the magnitude and mechanism of the influence of waves on tide surges in the Bohai Sea were revealed for the first time. The results showed that although consideration of the wave dependent surface wind stresses raise slightly the traditional surface wind stress, due to the accumulated effects, the computed results are improved on the whole. Storm level maximum modulation can reach 0.4 m. The results computed by the combined model agreed well with the measured data.展开更多
The modeling of generation and subsequent propagation of irregular waves in a numerical wave flume is performed by mean of the boundary element method. Random waves are generated by a piston-type wave generator at one...The modeling of generation and subsequent propagation of irregular waves in a numerical wave flume is performed by mean of the boundary element method. Random waves are generated by a piston-type wave generator at one end of the flume with the Mitsuyasu-Bretschneider spectrum used as the target spectrum for the generation. An artificial absorbing beach is placed at the other end of the flume to minimize wave reflection. Surface fluctuations are described by use of the Lagrangian description, and finite difference is adopted for the approximation of time derivative. To monitor the developments of the waves, a number of pseudo wave gauges are installed along the tank. Through comparison of the spectra from those gauges with the target spectrum, satisfactory results can be obtained from the present numerical scheme.展开更多
Among all environmental forces acting on ocean structures and marine vessels, those resulting from wave impacts are likely to yield the highest loads. Being highly nonlinear, transient and complex, a theoretical analy...Among all environmental forces acting on ocean structures and marine vessels, those resulting from wave impacts are likely to yield the highest loads. Being highly nonlinear, transient and complex, a theoretical analysis of their impact would be impossible without numerical simulations. In this paper, a pressure-split two-stage numerical algorithm is proposed based on Volume Of Fluid （VOF） methodology. The algorithm is characterized by introduction of two pressures at each half and full cycle time step, and thus it is a second-order accurate algorithm in time. A simplified second-order Godunov-type solver is used for the continuity equations. The method is applied to simulation of breaking waves in a 2-D water tank, and a qualitative comparison with experimental photo observations is made. Quite consistent results are observed between simulations and experiments. Commercially available software and Boundary Integral Method （BIM） have also been used to simulate the same problem. The results from present code and BIM are in good agreement with respect to breaking location and timing, while the results obtained from the comrnercial software which is only first-order accurate in time has clearly showed a temporal and spatial lag, verifying the need to use a higher order numerical scheme.展开更多
A numerical wave load model based on two-phase(water-air) Reynolds-Averaged Navier Stokes(RANS) type equations is used to evaluate hydrodynamic forces exerted on flatted-bottom seafloor mining tool during its entering...A numerical wave load model based on two-phase(water-air) Reynolds-Averaged Navier Stokes(RANS) type equations is used to evaluate hydrodynamic forces exerted on flatted-bottom seafloor mining tool during its entering ocean waves of deploying process.The discretization of the RANS equations is achieved by a finite volume approach(FV).The volume of fluid method(VOF) is employed to track the complicated free surface.A numerical wave tank is built to generate the ocean waves which are suitable for deploying seafloor mining tool.A typical deploying condition is employed to reflect the process of flatted-bottom body impacting with waves,and the pressure distribution of bottom is also presented.Four different lowering velocities are applied to obtain the time histories of maximum pressure of bottom,and it can be concluded that the pressure coefficient decreases with water velocity increasing,which is similar with ordinary water entry case.The numerical results clearly demonstrate the characteristics of flatted-bottom body entering ocean waves.展开更多
A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course ...A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course in waves. Two experiments are carried out respectively to mea- sure the wave loads and the free motions for a pair of side-by- side arranged ship models advancing with an identical speed in head regular waves. For comparison, each model is also tested alone. Predictions obtained by the present solution are found in favorable agreement with the model tests and are more accurate than the traditional method based on the three dimensional pulsating （3DP） source Green function. Numer- ical resonances and peak shift can be found in the 3DP pre- dictions, which result from the wave energy trapped in the gap between two ships and the extremely inhomogeneous wave load distribution on each hull. However, they can be eliminated by 3DTP, in which the speed affects the free sur- face and most of the wave energy can be escaped from the gap. Both the experiment and the present prediction show that hydrodynamic interaction effects on wave loads and free motions are significant. The present solver may serve as a validated tool to predict wave loads and motions of two ves- sels under replenishment at sea, and may help to evaluate the hydrodynamic interaction effects on the ships safety in replenishment operation.展开更多
A numerical wave model based on the modified four-order nonlinear Schoedinger （NKS） equation in deep water is developed to simulate freak waves. A standard split-step, pseudo-spectral method is used to solve NLS equ...A numerical wave model based on the modified four-order nonlinear Schoedinger （NKS） equation in deep water is developed to simulate freak waves. A standard split-step, pseudo-spectral method is used to solve NLS equation. The validation of the model is firstly verified, and then the simulation of freak waves is perforated by changing sideband condi- tions. Results show that freak waves entirely consistent with the definition in the evolution of wave trains are obtained. The possible occurrence mechanism of freak waves is discussed and the relevant characteristics are also analyzed.展开更多
This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are result...This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are resulted boundary and initial conditions. The method of splitting into physical processes receives system from three equations. Then we define the approximation order and investigate stability conditions of the discrete model. The sweep method was used to calculate the system of equations. This work presents surface gravity wave profiles for different propagation phases.展开更多
To deal with the numerical dispersion problem, by combining the staggeredgrid technology with the compact finite difference scheme, we derive a compact staggered- grid finite difference scheme from the first-order vel...To deal with the numerical dispersion problem, by combining the staggeredgrid technology with the compact finite difference scheme, we derive a compact staggered- grid finite difference scheme from the first-order velocity-stress wave equations for the transversely isotropic media. Comparing the principal truncation error terms of the compact staggered-grid finite difference scheme, the staggered-grid finite difference scheme, and the compact finite difference scheme, we analyze the approximation accuracy of these three schemes using Fourier analysis. Finally, seismic wave numerical simulation in transversely isotropic （VTI） media is performed using the three schemes. The results indicate that the compact staggered-grid finite difference scheme has the smallest truncation error, the highest accuracy, and the weakest numerical dispersion among the three schemes. In summary, the numerical modeling shows the validity of the compact staggered-grid finite difference scheme.展开更多
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 equ...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.展开更多
基金financially supported by the National Natural Science Foundation of China (Grant No. 52071128)the Natural Science Foundation of Jiangsu Basic Research Program (Grant No. BK20220082)
文摘Flume experiments play a pivotal role in studying wave propagation,with wave elements typically assumed to remain constant in the perpendicular direction.However,evident cross wave phenomena were observed within flumes under certain conditions.This paper presents new analytical solutions for both primary and cross waves on double shoals in a flume via linear shallow-water equations,which may be used to idealize dynamic experimental configurations of coral reefs.The primary waves on double shoals are described by the associated Legendre functions,whereas the ultimate solutions are derived by considering the incident and reflected waves in front of a bathymetry and the transmitted waves positioned behind it.The effects of the angular frequency and topographic parameters on the primary waves are subsequently analyzed.Cross waves on double shoals constitute a type of topographically trapped wave whose solutions are formulated by combining the first and second types of the associated Legendre functions.The angular frequency is not only determined by the wavenumber but also influenced by the topographic parameters.Numerical experiments are conducted to investigate the generation mechanism of cross waves.The consistency between the numerical results and analytical solutions confirms the validity of the new analytical framework of cross waves on double shoals.
基金funded by the Ministry of Science and Higher Education of the Russian Federation(agreement No.075-15-2022-1127 dated July 1,2022).
文摘The paper considers a catastrophic event-the eruption of Hunga Tonga-Hunga Ha’apai volcano on January 15,2022.The process of preparation and eruption of Hunga Tonga volcano generated tsunami waves that were observed throughout the World Ocean.This event was notable for its unprecedented global impact and the early appearance of tsunami waves at distant coastal stations.So,the first waves at tide gauge stations in Chile and Peru were recorded 4 hours earlier than the arrival time of tsunami waves to the tide gauge after the eruption of Tonga volcano.Two mechanisms are possible for the generation of early tsunami waves:acoustic Lamb waves generated by a volcanic explosion and submarine landslides that occurred on the slopes of the volcano during the preparatory phase of the eruption.In this study,numerical simulation of various pre-eruption landslide scenarios on the slope of Hunga Tonga volcano is carried out in an attempt to explain these early tsunami waves.Under computation the elastoplastic model of landslide was taken into account.Wave characteristics of a tsunami on the coast of Chile and Peru generated by a landslide process on a volcanic slope are obtained.A detailed comparison of virtual tide gauge data with observational ones is used to validate this model.The results obtained can be used to improve early warning systems.
基金the National Natural Science Foundation of China(Grant Nos.52371289 and 51979192).
文摘A three-dimensional numerical model of sand wave dynamics,incorporating the interaction of currents and waves at various angles,has been developed using the Regional Ocean Modeling System(ROMS).This model accounts for both bedload and suspended load sediment transport under combined waves and current conditions.The investigation examines the influence of several key parameters,including the rotation angle of sand waves relative to the main current,tidal current velocity amplitude,residual current,water depth,wave height,wave period,and wave direction,on sand wave evolution.The growth rate and migration rate of sand waves decrease as their rotation angle increases.For rotation angles smaller than 15°,sand wave evolution can be effectively simulated by a vertical 2D model with an error within 10%.The numerical results demonstrate that variations in tidal current velocity amplitude or residual current affect both vertical growth and horizontal migration of sand waves.As tidal current velocity amplitude and residual current increase,the growth rate initially rises to a maximum before decreasing.The migration rate shows a consistent increase with increasing tidal current amplitude and residual current.Under combined waves and current,both growth and migration rates decrease as water depth increases.With increasing wave height and period,the growth rate and migration rate initially rise to maximum values before declining,while showing a consistent increase with wave height and period.The change rate of sand waves reaches its maximum when wave propagation aligns parallel to tidal currents,and reaches its minimum when wave propagation is perpendicular to the currents.This phenomenon can be explained by the fluctuation of total bed shear stress relative to the angle of interaction between waves and current.
文摘With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.
文摘The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.
基金Trans-Century Training program Fund for the Talent,Ministry of Education of China
文摘A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.
基金The National Key Research and Development Plan of China under contract No.2016YFC1401300the National Natural Science Foundation of China under contract No.41276008the Taishan Scholars Program
文摘Numerical simulations based on a high-resolution three-dimensional MIT general circulation model(MITgcm)using realistic topography and tidal forcing are conducted to investigate the generation and propagation of the so-called type-a waves(large-amplitude rank-ordered wave packets)and type-b waves(isolated wave packets)in the northern South China Sea.At first,we summarized and analyzed the generation and propagation characteristics of these waves.Then,energy budget at the Luzon Strait is calculated.Energy generation has three local maxima every day,of which the largest one corresponds to the emergence of the type-a wave.Energy flux at the west boundary of the Luzon Strait shows two local maxima each day.The larger one is consistent with the generation of the type-a wave and the smaller one is in correspondence with the generation of the type-b wave.Sensitivity experiments are designed to explore the role of the east and west ridge of the Luzon Strait on the generation and propagation of the type-a and type-b waves.It is found that the east ridge is indispensable on the generation of the type-a wave while the west ridge has little contribution.The west ridge diminishes the type-a waves'amplitude but hardly changes their propagation speed.The type-b waves also come from perturbation signals which originate from the east ridge and are enhanced in amplitude and reduced in propagation speed by the west ridge.
基金The National Natural Science Foundation of China under contract Nos 51579036 and 51579030the Fundamental Research Funds for the Central Universities of China under contract No.DUT14YQ10
文摘Prediction of coastal sediment transport is of particularly importance for analyzing coast erosion accurately and solving the corresponding coast protection engineering problems.The present study provided a numerical scheme for sediment transport in coastal waves and wave-induced currents.In the scheme,the sand transport model was implemented with wave refraction-diffraction model and near-shore current model.Coastal water wave was simulated by using the parabolic mild-slope equation in which wave refraction,diffraction and breaking effects are considered.Wave-induced current was simulated by using the nonlinear shallow water equations in which wave provides radiation stresses for driving current.Then,sediment transport in waves and wave-induced currents was simulated by using the two-dimensional suspended sediment transport equations for suspended sediment and the bed-load transport equation for bed load.The numerical scheme was validated by experiment results from the Large-scale Sediment Transport Facility at the US Army Corps of Engineer Research and Development Center in Vicksburg.The numerical results showed that the present scheme is an effective tool for modeling coastal sediment transport in waves and near-shore currents.
基金Supported by the High-Tech Research and Development Program of China (863 Program, No. 2001AA633070 2003AA604040)the National Natural Science Foundation of China (No. 40476015).
文摘This paper introduces a numerical model for studying the evolution of a periodic wave train, shoaling, and breaking in surf zone. The model can solve the Reynolds averaged Navier-Stokes （RANS） equations for a mean flow, and the k-ε equations for turbulence kinetic energy k and turbulence dissipation rate ε. To track a free surface, the volume of fluid （VOF） function, satisfying the advection equation was introduced. In the numerical treatment, third-order upwind difference scheme was applied to the convection terms of the RANS equations in order to reduce the effect of numerical viscosity. The shoaling and breaking processes of a periodic wave train on gently sloping beaches were modeled. The computed wave heights of a sloping beach and the distribution of breaking wave pressure on a vertical wall were compared with laboratory data.
基金The National Basic Research Program of China under contract No.2013CB430403the National Natural Science Foundation of China under contract No.51179025+1 种基金the Open Foundation of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering under contract No.2013491511the Open Foundation of State Key Laboratory of Ocean Engineering under contract No.1305
文摘In this paper, the water waves and wave-induced longshore currents in Obak6y coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic mild-slope equation for coastal water waves and the nonlinear shallow water equation for the wave-induced currents. The wave transformation under the effects of shoaling, refraction, diffraction and breaking is considered, and the wave provides radiation stresses for driving currents in the model. The numerical results for the water wave-induced longshore currents were validated by the measured data to demonstrate the efficiency of the numerical model. Then the water waves and longshore currents induced by the waves from main directions were numerically simulated and analyzed based on the numerical results. The numerical results show that the movement of the longshore currents was different while the wave proDaRated to a coastal zone from different directions.
文摘A two-dimensional （2D） numerical model is developed for the wave sim- ulation and propagation in a wave flume. The fluid flow is assumed to be viscous and incompressible, and the Navier-Stokes and continuity equations are used as the governing equations. The standard k-e model is used to model the turbulent flow. The Navier- Stokes equations are discretized using the staggered grid finite difference method and solved by the simplified marker and cell （SMAC） method. Waves are generated and propagated using a piston type wave maker. An open boundary condition is used at the end of the numerical flume. Some standard tests, such as the lid-driven cavity, the constant unidirectional velocity field, the shearing flow, and the dam-break on the dry bed, are performed to valid the model. To demonstrate the capability and accuracy of the present method, the results of generated waves are compared with available wave theories. Finally, the clustering technique （CT） is used for the mesh generation, and the best condition is suggested.
文摘The author’s combined numerical model consisting of a third generation shallow water wave model and a 3 D tide surge model with wave dependent surface wind stress were used to study the influence of waves on tide surge motion. For the typical weather case, in this study, the magnitude and mechanism of the influence of waves on tide surges in the Bohai Sea were revealed for the first time. The results showed that although consideration of the wave dependent surface wind stresses raise slightly the traditional surface wind stress, due to the accumulated effects, the computed results are improved on the whole. Storm level maximum modulation can reach 0.4 m. The results computed by the combined model agreed well with the measured data.
文摘The modeling of generation and subsequent propagation of irregular waves in a numerical wave flume is performed by mean of the boundary element method. Random waves are generated by a piston-type wave generator at one end of the flume with the Mitsuyasu-Bretschneider spectrum used as the target spectrum for the generation. An artificial absorbing beach is placed at the other end of the flume to minimize wave reflection. Surface fluctuations are described by use of the Lagrangian description, and finite difference is adopted for the approximation of time derivative. To monitor the developments of the waves, a number of pseudo wave gauges are installed along the tank. Through comparison of the spectra from those gauges with the target spectrum, satisfactory results can be obtained from the present numerical scheme.
基金Supported by the National Natural Science Foundation of China under Grant Nos 50679010 and 50579004.
文摘Among all environmental forces acting on ocean structures and marine vessels, those resulting from wave impacts are likely to yield the highest loads. Being highly nonlinear, transient and complex, a theoretical analysis of their impact would be impossible without numerical simulations. In this paper, a pressure-split two-stage numerical algorithm is proposed based on Volume Of Fluid （VOF） methodology. The algorithm is characterized by introduction of two pressures at each half and full cycle time step, and thus it is a second-order accurate algorithm in time. A simplified second-order Godunov-type solver is used for the continuity equations. The method is applied to simulation of breaking waves in a 2-D water tank, and a qualitative comparison with experimental photo observations is made. Quite consistent results are observed between simulations and experiments. Commercially available software and Boundary Integral Method （BIM） have also been used to simulate the same problem. The results from present code and BIM are in good agreement with respect to breaking location and timing, while the results obtained from the comrnercial software which is only first-order accurate in time has clearly showed a temporal and spatial lag, verifying the need to use a higher order numerical scheme.
基金Project(51305463)supported by National Natural Science Foundation of ChinaProject(2012QNZT01601005125)supported by Free Exploration Plan of Central South University,ChinaProject supported by Postdoctoral Foundation of Central South university,China
文摘A numerical wave load model based on two-phase(water-air) Reynolds-Averaged Navier Stokes(RANS) type equations is used to evaluate hydrodynamic forces exerted on flatted-bottom seafloor mining tool during its entering ocean waves of deploying process.The discretization of the RANS equations is achieved by a finite volume approach(FV).The volume of fluid method(VOF) is employed to track the complicated free surface.A numerical wave tank is built to generate the ocean waves which are suitable for deploying seafloor mining tool.A typical deploying condition is employed to reflect the process of flatted-bottom body impacting with waves,and the pressure distribution of bottom is also presented.Four different lowering velocities are applied to obtain the time histories of maximum pressure of bottom,and it can be concluded that the pressure coefficient decreases with water velocity increasing,which is similar with ordinary water entry case.The numerical results clearly demonstrate the characteristics of flatted-bottom body entering ocean waves.
基金supported by the National Natural Science Foundation of China(50879090)the Key Research Program of Hydrodynamics of China(9140A14030712JB11044)
文摘A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course in waves. Two experiments are carried out respectively to mea- sure the wave loads and the free motions for a pair of side-by- side arranged ship models advancing with an identical speed in head regular waves. For comparison, each model is also tested alone. Predictions obtained by the present solution are found in favorable agreement with the model tests and are more accurate than the traditional method based on the three dimensional pulsating （3DP） source Green function. Numer- ical resonances and peak shift can be found in the 3DP pre- dictions, which result from the wave energy trapped in the gap between two ships and the extremely inhomogeneous wave load distribution on each hull. However, they can be eliminated by 3DTP, in which the speed affects the free sur- face and most of the wave energy can be escaped from the gap. Both the experiment and the present prediction show that hydrodynamic interaction effects on wave loads and free motions are significant. The present solver may serve as a validated tool to predict wave loads and motions of two ves- sels under replenishment at sea, and may help to evaluate the hydrodynamic interaction effects on the ships safety in replenishment operation.
基金This paperis part of the research project (104061)supported by the keytechnology program,the Ministry of Educa-tion of China
文摘A numerical wave model based on the modified four-order nonlinear Schoedinger （NKS） equation in deep water is developed to simulate freak waves. A standard split-step, pseudo-spectral method is used to solve NLS equation. The validation of the model is firstly verified, and then the simulation of freak waves is perforated by changing sideband condi- tions. Results show that freak waves entirely consistent with the definition in the evolution of wave trains are obtained. The possible occurrence mechanism of freak waves is discussed and the relevant characteristics are also analyzed.
文摘This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are resulted boundary and initial conditions. The method of splitting into physical processes receives system from three equations. Then we define the approximation order and investigate stability conditions of the discrete model. The sweep method was used to calculate the system of equations. This work presents surface gravity wave profiles for different propagation phases.
基金supported by the National High-Tech Research and Development Program of China(Grant No.2006AA06Z202)the Open Fund of the Key Laboratory of Geophysical Exploration of CNPC(Grant No.GPKL0802)+1 种基金the Graduate Student Innovation Fund of China University of Petroleum(East China)(Grant No.S2008-1)the Program for New Century Excellent Talents in University(Grant No.NCET-07-0845)
文摘To deal with the numerical dispersion problem, by combining the staggeredgrid technology with the compact finite difference scheme, we derive a compact staggered- grid finite difference scheme from the first-order velocity-stress wave equations for the transversely isotropic media. Comparing the principal truncation error terms of the compact staggered-grid finite difference scheme, the staggered-grid finite difference scheme, and the compact finite difference scheme, we analyze the approximation accuracy of these three schemes using Fourier analysis. Finally, seismic wave numerical simulation in transversely isotropic （VTI） media is performed using the three schemes. The results indicate that the compact staggered-grid finite difference scheme has the smallest truncation error, the highest accuracy, and the weakest numerical dispersion among the three schemes. In summary, the numerical modeling shows the validity of the compact staggered-grid finite difference scheme.
基金Supported by the Fund of National Nature Sciences of China
文摘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.
