In this paper, the solitary wave deformation along a gentle slope and the impact pressure, on the wall are investigated experimentally and the results are compared with numerical results obtained based on the volume o...In this paper, the solitary wave deformation along a gentle slope and the impact pressure, on the wall are investigated experimentally and the results are compared with numerical results obtained based on the volume of fluid (VOF) method. The topography used in the experiment consists of three segments. The left segment is a 1:4 slope, the middle segment 1 :SO slope and the right segment a horizontal bed. Both the wave heights and breaking points obtained from numerical simulation and experiments are in good agreement. Numerical results give reasonable pressure distributions of breaking waves on the wall.展开更多
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 wate...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.展开更多
An internal solitary wave of elevation in a two-layer density stratified system of an incompressible, viscous and homogeneous fluid was studied. The run-up of a wave of elevation encountering different slopes was inve...An internal solitary wave of elevation in a two-layer density stratified system of an incompressible, viscous and homogeneous fluid was studied. The run-up of a wave of elevation encountering different slopes was investigated numerically based on solving the continuity, Navier-Stokes and convective-diffusion equations within the Boussinesq approximation. The commercial software COMSOL Multiphysics was used to conduct the numerical simulations. For gradual shoals, a bolus formed that transported dense fluid up the shoal. The bolus disappeared when it reached its maximum height on the slope due to the draining of the dense fluid. Various shoal angles were simulated to detect the critical angle above which a bolus does not form. An angle of 30 or less resulted in the formation of a bolus. In addition, the simulations demonstrated that the size of the bolus induced by shallower slopes was larger and that the vertical height traveled by the bolus was insensitive to the slope of the shoal.展开更多
Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-...Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-SJTU is used to solve the Reynolds-Averaged Navier–Stokes(RANS) equations with the SST k ?? turbulence model. The PISO algorithm is utilized for the pressure-velocity coupling. The air-water interface is captured via Volume of Fluid(VOF) technique. The present numerical model is validated by simulating the solitary wave run-up and reflected against a vertical wall, and solitary wave run-up on a vertical circular cylinder. Comparisons between numerical results and available experimental data show satisfactory agreement. Furthermore, simulations are carried out to study the solitary wave run-up on the cylinder with different incident wave height H and different cylinder radius a. The relationships of the wave run-up height with the incident wave height H, cylinder radius a are analyzed. The evolutions of the scattering free surface and vortex shedding are also presented to give a better understanding of the process of nonlinear wave-cylinder interaction.展开更多
The concept of candidate particle set is introduced in the MPS gridless numerical method to generate neighboring particle set ma- trix, which can reduce the CPU time to 1/11 of that before introduction. The Bi-CGSTAB ...The concept of candidate particle set is introduced in the MPS gridless numerical method to generate neighboring particle set ma- trix, which can reduce the CPU time to 1/11 of that before introduction. The Bi-CGSTAB (bi-conjugate gradient stabilized) algorithm is applied to solving the Poisson pressure equation, by which the solving speed is significantly accelerated. The process of solitary waves propagating over a numerical flume and interacting with a vertical wall is simulated. The simulated results of water surface elevation are in good agreement with the analytical solution as well as the measured data. The predicted maximum values of the run-up of solitary waves with various relative incident wave heights agree well with the measured results.展开更多
The interaction between structure and wave is a typical phenomenon in naval architecture and ocean engineering.In this paper,numerical simulation is carried out to study the interaction between a two-dimensional subme...The interaction between structure and wave is a typical phenomenon in naval architecture and ocean engineering.In this paper,numerical simulation is carried out to study the interaction between a two-dimensional submerged,fixed,horizontal rigid plate and solitary wave with our in-house meshless particle CFD solver MLParticle-SJTU.First,the in-house CFD solver is verified by experimental results conducted at the State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology.During the verification,the plate is submerged under water and the solitary wave with a given amplitude is generated by a piston-type wave maker.Free surface elevation of the wave and the pressure impacting on the plate is recorded and compared with experimental data respectively.The predicted pressure and surface elevation agree well with the experimental results.Then in order to further investigate factors affecting wave-structure interaction,wave height,submerged depth and plate length are analyzed.展开更多
Using a polynomial expansion method, the general exact solitary wave solution and singular one areconstructed for the non-linear KS equation. This approach is obviously applicable to a large variety of nonlinear evolu...Using a polynomial expansion method, the general exact solitary wave solution and singular one areconstructed for the non-linear KS equation. This approach is obviously applicable to a large variety of nonlinear evolution equation.展开更多
Studies of the breaking criteria for solitary waves on a slope are presented in this paper. The boundary element method is used to model the processes of shoaling and breaking of solitary waves on various slopes. Empi...Studies of the breaking criteria for solitary waves on a slope are presented in this paper. The boundary element method is used to model the processes of shoaling and breaking of solitary waves on various slopes. Empirical formulae that can be used to characterize the breaking of solitary waves are presented. These include the breaking index, the wave height, the water depth, and the maximum particle velocity at the point of breaking. Comparisons with the results of other researches are given.展开更多
In this paper, exact and numerical solutions are calculated for discrete complex Ginzburg-Landau equation with initial condition by considering the modified Adomian decomposition method (mADM), which is an efficient...In this paper, exact and numerical solutions are calculated for discrete complex Ginzburg-Landau equation with initial condition by considering the modified Adomian decomposition method (mADM), which is an efficient method and does not need linearization, weak nonlinearity assumptions or perturbation theory. The numerical solutions are also compared with their corresponding analytical solutions. It is shown that a very good approximation is achieved with the analytical solutions. Finally, the modulational instability is investigated and the corresponding condition is given.展开更多
The propagation, shoaling and breaking of solitary waves on mild slopes are simulated by boundary element method. In this paper, the criterion of breaking solitary waves on mild slopes is discussed. The criterion is t...The propagation, shoaling and breaking of solitary waves on mild slopes are simulated by boundary element method. In this paper, the criterion of breaking solitary waves on mild slopes is discussed. The criterion is that the ratio of horizontal velocity of water particles on the wave crest to wave celerity equals one. However, the case that the ratio of horizontal velocity of water particles on the wave crest to wave celerity is below one but the front face of wave profile becomes vertical is also considered as a breaking criterion. According to the above criteria, the breaking index for slopes 1:10 to 1:25 is studied. The result is compared to other researchers'. The deformation of solitary waves on slopes is discussed and the distribution of fluid velocities at breaking is shown.展开更多
In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation (RLW) is integrated fully by using an exponentia...In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation (RLW) is integrated fully by using an exponential B-spline Galerkin method in space together with Crank-Nicolson method in time. Three numerical examples related to propagation of sin- gle solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.展开更多
As a new method, the Level Set method had been developed to compute the interface of two-phase flow. The basic mathematical theory and the detailed method to solve the free surface hydrodynamic problem had been invest...As a new method, the Level Set method had been developed to compute the interface of two-phase flow. The basic mathematical theory and the detailed method to solve the free surface hydrodynamic problem had been investigated. By using the Level Set method, the transformation of a solitary wave over a front step was simulated. The results were in good agreement with laboratory experiments.展开更多
In this research work,we constructed the optical soliton solutions of nonlinear complex Kundu-Eckhaus(KE)equation with the help of modified mathematical method.We obtained the solutions in the form of dark solitons,br...In this research work,we constructed the optical soliton solutions of nonlinear complex Kundu-Eckhaus(KE)equation with the help of modified mathematical method.We obtained the solutions in the form of dark solitons,bright solitons and combined dark-bright solitons,travelling wave and periodic wave solutions with general coefficients.In our knowledge earlier reported results of the KE equation with specific coefficients.These obtained solutions are more useful in the development of optical fibers,dynamics of solitons,dynamics of adiabatic parameters,dynamics of fluid,problems of biomedical,industrial phenomena and many other branches.All calculations show that this technique is more powerful,effective,straightforward,and fruitfulness to study analytically other higher-order nonlinear complex PDEs involves in mathematical physics,quantum physics,Geo physics,fluid mechanics,hydrodynamics,mathematical biology,field of engineering and many other physical sciences.展开更多
A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid (VOF) method. An active 'absorbing wave-maker based on linear wave theory is set on the left boun...A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid (VOF) method. An active 'absorbing wave-maker based on linear wave theory is set on the left boundary of the wave flume. The progressive waves and the absorbing waves are generated simultaneously at the active wave generating-absorbing boundary. The absorbing waves are generated to eliminate the waves coming back to the generating boundary due to reflection from the outflow boundary and the structures. SIRW method proposed by Frigaard and Brorsen (1995) is used to separate the incident waves and reflected waves. The digital filters are designed based on the surface elevation signals of the two wave gauges. The corrected velocity of the wave-maker paddle is the output from the digital filter in real time. The numerical results of regular and irregular waves by the active absorbing-generating boundary are compared with the numerical results by the ordinary generating boundary to verify the performance of the active absorbing-generator boundary. The differences between the initial incident waves and the estimated incident waves are analyzed.展开更多
Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to e...Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.展开更多
Base on the Lagrangian description and finite differencing of the time step, the generation, propagation and deformation of solitary waves are simulated numerically by means of boundary element method. An algorithm to...Base on the Lagrangian description and finite differencing of the time step, the generation, propagation and deformation of solitary waves are simulated numerically by means of boundary element method. An algorithm to generate waves with any prescribed form is also implanted in the scheme, but in this research, solitary waves are studied. The numerical model is first verified by studying the case of a solitary wave impinging against a vertical wall. Time histories of evolution of a soliton running up on a sloping beach, as well as over a submerged obstacle are then presented. Applications to scattering problems due to solitary waves progressing over a slope onto a shelf are also shown.展开更多
Based on the time-domain higher-order boundary element method(HOBEM), a two-dimensional numerical wave flume is developed to investigate solitary wave interaction with two rectangular boxes with a narrow gap.In the nu...Based on the time-domain higher-order boundary element method(HOBEM), a two-dimensional numerical wave flume is developed to investigate solitary wave interaction with two rectangular boxes with a narrow gap.In the numerical model, the fully nonlinear boundary conditions are satisfied on the free surface, the mixed Eulerian-Lagrangian method is adopted to track the transient water surface and the fourth-order Runga-Kutta method is used to predict the velocity potential and wave elevation on the free surface. The acceleration potential technique is used to compute the transient wave forces along the wetted object surface. A piston-type wavemaker is used to generate solitary waves. The proposed model is validated by comparing the simulated wave run-up and the wave loads with the published experimental and numerical results of the reflection of a solitary wave from a vertical wall. Then, numerical experiments are performed to study the effects of the narrow gap and the size of each box on the wave run-ups at the two sides of the two-box system and in the narrow gap between two boxes, and the wave loads on the two boxes. The interaction between double solitary waves with a time interval between them with a two-box system is also investigated.展开更多
A nonhydrostatic numerical model was developed and numerical experiments performed on the interaction of an internal solitary wave (ISW) with a sill, for a two-layer fluid with a diffusive interface. Based on the bl...A nonhydrostatic numerical model was developed and numerical experiments performed on the interaction of an internal solitary wave (ISW) with a sill, for a two-layer fluid with a diffusive interface. Based on the blocking parameter (Br), the flow was classified into three cases: (1) when bottom topography has little influence on the propagation and spatial structure of the ISW (Br〈0.5), (2) where the ISW is distorted significantly by the blocking effect of the topography (though no wave breaking occurs, (0.5〈Br〈0.7), and (3) where the ISW is broken as it encounters and passes over the bottom topography (0.7〈Br). The numerical results obtained here are consistent with those obtained in laboratory experiments. The breaking process of the incident ISW when Br=0.7 was completely reproduced. Dissipation rate was linearly related to the blocking parameter when B,〈0.7, and the maximum dissipation rate could reach about 34% as Br raised to about 1.0. After that, instead of breaking, more reflection happened. Similarly, breaking induced mixing was also most effective during Br around 1.0, and can be up to 0.16.展开更多
This paper contributes to the diffraction of a solitary Wave by a cylinder, governed by the Boussinesq eq. The Spectrum Method is used to transform 2-D eqs. into a set of 1-D eqs., which is solved by F. D. M. The exam...This paper contributes to the diffraction of a solitary Wave by a cylinder, governed by the Boussinesq eq. The Spectrum Method is used to transform 2-D eqs. into a set of 1-D eqs., which is solved by F. D. M. The example given in the paper shows that this method is both accurate and cost ef- fective.展开更多
Developing serpent-type wave generators to generate solitary waves in a 3D-basin was investigated in this study. Based on the Lagrangian description with time-marching procedures and finite differences of the time der...Developing serpent-type wave generators to generate solitary waves in a 3D-basin was investigated in this study. Based on the Lagrangian description with time-marching procedures and finite differences of the time derivative, a 3D multiple directional wave basin with multidirectional piston wave generators was developed to simulate ocean waves by using BEM with quadrilateral elements, and to simulate wave-caused problems with fully nonlinear water surface conditions. The simulations of perpendicular solitary waves were conducted in the first instance to verify this scheme. Furthermore, the comparison of the waveform variations confirms that the estimation of 3D solitary waves is a feasible scheme.展开更多
文摘In this paper, the solitary wave deformation along a gentle slope and the impact pressure, on the wall are investigated experimentally and the results are compared with numerical results obtained based on the volume of fluid (VOF) method. The topography used in the experiment consists of three segments. The left segment is a 1:4 slope, the middle segment 1 :SO slope and the right segment a horizontal bed. Both the wave heights and breaking points obtained from numerical simulation and experiments are in good agreement. Numerical results give reasonable pressure distributions of breaking waves on the wall.
基金The National Natural Science Foundation of China under contract No.51279023the Public Science and Technology Research Funds Projects of Ocean under contract No.201205023+1 种基金the Special Funds for Postdoctoral Innovative Projects of Liaoning Province of China under contract No.2011921018the Special Funds for Talent Projects of Dalian Ocean University under contract No.SYYJ2011004
文摘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.
文摘An internal solitary wave of elevation in a two-layer density stratified system of an incompressible, viscous and homogeneous fluid was studied. The run-up of a wave of elevation encountering different slopes was investigated numerically based on solving the continuity, Navier-Stokes and convective-diffusion equations within the Boussinesq approximation. The commercial software COMSOL Multiphysics was used to conduct the numerical simulations. For gradual shoals, a bolus formed that transported dense fluid up the shoal. The bolus disappeared when it reached its maximum height on the slope due to the draining of the dense fluid. Various shoal angles were simulated to detect the critical angle above which a bolus does not form. An angle of 30 or less resulted in the formation of a bolus. In addition, the simulations demonstrated that the size of the bolus induced by shallower slopes was larger and that the vertical height traveled by the bolus was insensitive to the slope of the shoal.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51379125,51411130131,and 11432009)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(Grant No.2013022)the National Basic Research Program of China(973 Program,Grant No.2013CB036103)
文摘Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-SJTU is used to solve the Reynolds-Averaged Navier–Stokes(RANS) equations with the SST k ?? turbulence model. The PISO algorithm is utilized for the pressure-velocity coupling. The air-water interface is captured via Volume of Fluid(VOF) technique. The present numerical model is validated by simulating the solitary wave run-up and reflected against a vertical wall, and solitary wave run-up on a vertical circular cylinder. Comparisons between numerical results and available experimental data show satisfactory agreement. Furthermore, simulations are carried out to study the solitary wave run-up on the cylinder with different incident wave height H and different cylinder radius a. The relationships of the wave run-up height with the incident wave height H, cylinder radius a are analyzed. The evolutions of the scattering free surface and vortex shedding are also presented to give a better understanding of the process of nonlinear wave-cylinder interaction.
基金This work was supported by the Natural Science Foundation of China under No.50479047.
文摘The concept of candidate particle set is introduced in the MPS gridless numerical method to generate neighboring particle set ma- trix, which can reduce the CPU time to 1/11 of that before introduction. The Bi-CGSTAB (bi-conjugate gradient stabilized) algorithm is applied to solving the Poisson pressure equation, by which the solving speed is significantly accelerated. The process of solitary waves propagating over a numerical flume and interacting with a vertical wall is simulated. The simulated results of water surface elevation are in good agreement with the analytical solution as well as the measured data. The predicted maximum values of the run-up of solitary waves with various relative incident wave heights agree well with the measured results.
基金the National Natural Science Foundation of China(Grant Nos.51909160 and 51879159)the National Key Research and Development Program of China(Grant Nos.2019YFB1704200 and 2019YFC0312400)+2 种基金Chang Jiang Scholars Program(Grant No.T2014099)Shanghai Excellent Academic Leaders Program(Grant No.17XD1402300)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(Grant No.2016-23/09).
文摘The interaction between structure and wave is a typical phenomenon in naval architecture and ocean engineering.In this paper,numerical simulation is carried out to study the interaction between a two-dimensional submerged,fixed,horizontal rigid plate and solitary wave with our in-house meshless particle CFD solver MLParticle-SJTU.First,the in-house CFD solver is verified by experimental results conducted at the State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology.During the verification,the plate is submerged under water and the solitary wave with a given amplitude is generated by a piston-type wave maker.Free surface elevation of the wave and the pressure impacting on the plate is recorded and compared with experimental data respectively.The predicted pressure and surface elevation agree well with the experimental results.Then in order to further investigate factors affecting wave-structure interaction,wave height,submerged depth and plate length are analyzed.
文摘Using a polynomial expansion method, the general exact solitary wave solution and singular one areconstructed for the non-linear KS equation. This approach is obviously applicable to a large variety of nonlinear evolution equation.
基金The work was financially supported by the Taiwan Science Council(Project No.NSC-89-2611-E-019-064)
文摘Studies of the breaking criteria for solitary waves on a slope are presented in this paper. The boundary element method is used to model the processes of shoaling and breaking of solitary waves on various slopes. Empirical formulae that can be used to characterize the breaking of solitary waves are presented. These include the breaking index, the wave height, the water depth, and the maximum particle velocity at the point of breaking. Comparisons with the results of other researches are given.
基金supported by National Natural Science Foundation of China under Grant No. 10672147
文摘In this paper, exact and numerical solutions are calculated for discrete complex Ginzburg-Landau equation with initial condition by considering the modified Adomian decomposition method (mADM), which is an efficient method and does not need linearization, weak nonlinearity assumptions or perturbation theory. The numerical solutions are also compared with their corresponding analytical solutions. It is shown that a very good approximation is achieved with the analytical solutions. Finally, the modulational instability is investigated and the corresponding condition is given.
文摘The propagation, shoaling and breaking of solitary waves on mild slopes are simulated by boundary element method. In this paper, the criterion of breaking solitary waves on mild slopes is discussed. The criterion is that the ratio of horizontal velocity of water particles on the wave crest to wave celerity equals one. However, the case that the ratio of horizontal velocity of water particles on the wave crest to wave celerity is below one but the front face of wave profile becomes vertical is also considered as a breaking criterion. According to the above criteria, the breaking index for slopes 1:10 to 1:25 is studied. The result is compared to other researchers'. The deformation of solitary waves on slopes is discussed and the distribution of fluid velocities at breaking is shown.
基金supported by the Scientific and Technological Research Council of Turkey(Grant No.113F394)
文摘In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation (RLW) is integrated fully by using an exponential B-spline Galerkin method in space together with Crank-Nicolson method in time. Three numerical examples related to propagation of sin- gle solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.
文摘As a new method, the Level Set method had been developed to compute the interface of two-phase flow. The basic mathematical theory and the detailed method to solve the free surface hydrodynamic problem had been investigated. By using the Level Set method, the transformation of a solitary wave over a front step was simulated. The results were in good agreement with laboratory experiments.
文摘In this research work,we constructed the optical soliton solutions of nonlinear complex Kundu-Eckhaus(KE)equation with the help of modified mathematical method.We obtained the solutions in the form of dark solitons,bright solitons and combined dark-bright solitons,travelling wave and periodic wave solutions with general coefficients.In our knowledge earlier reported results of the KE equation with specific coefficients.These obtained solutions are more useful in the development of optical fibers,dynamics of solitons,dynamics of adiabatic parameters,dynamics of fluid,problems of biomedical,industrial phenomena and many other branches.All calculations show that this technique is more powerful,effective,straightforward,and fruitfulness to study analytically other higher-order nonlinear complex PDEs involves in mathematical physics,quantum physics,Geo physics,fluid mechanics,hydrodynamics,mathematical biology,field of engineering and many other physical sciences.
基金supported by the‘New Century Outstanding Talent’Scheme of the Ministry of Education of China(Grant No.NCET-07-0135)
文摘A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid (VOF) method. An active 'absorbing wave-maker based on linear wave theory is set on the left boundary of the wave flume. The progressive waves and the absorbing waves are generated simultaneously at the active wave generating-absorbing boundary. The absorbing waves are generated to eliminate the waves coming back to the generating boundary due to reflection from the outflow boundary and the structures. SIRW method proposed by Frigaard and Brorsen (1995) is used to separate the incident waves and reflected waves. The digital filters are designed based on the surface elevation signals of the two wave gauges. The corrected velocity of the wave-maker paddle is the output from the digital filter in real time. The numerical results of regular and irregular waves by the active absorbing-generating boundary are compared with the numerical results by the ordinary generating boundary to verify the performance of the active absorbing-generator boundary. The differences between the initial incident waves and the estimated incident waves are analyzed.
基金Supported by the National Natural Science Foundation of China(51379125,51490675,11432009,51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of the Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for Doctoral Students
文摘Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.
文摘Base on the Lagrangian description and finite differencing of the time step, the generation, propagation and deformation of solitary waves are simulated numerically by means of boundary element method. An algorithm to generate waves with any prescribed form is also implanted in the scheme, but in this research, solitary waves are studied. The numerical model is first verified by studying the case of a solitary wave impinging against a vertical wall. Time histories of evolution of a soliton running up on a sloping beach, as well as over a submerged obstacle are then presented. Applications to scattering problems due to solitary waves progressing over a slope onto a shelf are also shown.
基金The National Natural Science Foundation of China under contract Nos 51679036,51222902 and 51490672the Program for New Century Excellent Talents in University under contract No.NCET-13-0076
文摘Based on the time-domain higher-order boundary element method(HOBEM), a two-dimensional numerical wave flume is developed to investigate solitary wave interaction with two rectangular boxes with a narrow gap.In the numerical model, the fully nonlinear boundary conditions are satisfied on the free surface, the mixed Eulerian-Lagrangian method is adopted to track the transient water surface and the fourth-order Runga-Kutta method is used to predict the velocity potential and wave elevation on the free surface. The acceleration potential technique is used to compute the transient wave forces along the wetted object surface. A piston-type wavemaker is used to generate solitary waves. The proposed model is validated by comparing the simulated wave run-up and the wave loads with the published experimental and numerical results of the reflection of a solitary wave from a vertical wall. Then, numerical experiments are performed to study the effects of the narrow gap and the size of each box on the wave run-ups at the two sides of the two-box system and in the narrow gap between two boxes, and the wave loads on the two boxes. The interaction between double solitary waves with a time interval between them with a two-box system is also investigated.
基金The National Natural Science Foundation of China under contract Nos 41528601 and 41376029the Youth Innovation Promotion Association of Chinese Academy of Sciences under contract No.Y4KY07103Lthe Strategic Priority Research Program of the Chinese Academy of Sciences under contract No.XDA11020101
文摘A nonhydrostatic numerical model was developed and numerical experiments performed on the interaction of an internal solitary wave (ISW) with a sill, for a two-layer fluid with a diffusive interface. Based on the blocking parameter (Br), the flow was classified into three cases: (1) when bottom topography has little influence on the propagation and spatial structure of the ISW (Br〈0.5), (2) where the ISW is distorted significantly by the blocking effect of the topography (though no wave breaking occurs, (0.5〈Br〈0.7), and (3) where the ISW is broken as it encounters and passes over the bottom topography (0.7〈Br). The numerical results obtained here are consistent with those obtained in laboratory experiments. The breaking process of the incident ISW when Br=0.7 was completely reproduced. Dissipation rate was linearly related to the blocking parameter when B,〈0.7, and the maximum dissipation rate could reach about 34% as Br raised to about 1.0. After that, instead of breaking, more reflection happened. Similarly, breaking induced mixing was also most effective during Br around 1.0, and can be up to 0.16.
文摘This paper contributes to the diffraction of a solitary Wave by a cylinder, governed by the Boussinesq eq. The Spectrum Method is used to transform 2-D eqs. into a set of 1-D eqs., which is solved by F. D. M. The example given in the paper shows that this method is both accurate and cost ef- fective.
基金supported by the Science Council under the Project Nos.NSC-95-2221-E-019-075-MY3(CRC)andNSC-97-2221-E-236-011-(RSS)
文摘Developing serpent-type wave generators to generate solitary waves in a 3D-basin was investigated in this study. Based on the Lagrangian description with time-marching procedures and finite differences of the time derivative, a 3D multiple directional wave basin with multidirectional piston wave generators was developed to simulate ocean waves by using BEM with quadrilateral elements, and to simulate wave-caused problems with fully nonlinear water surface conditions. The simulations of perpendicular solitary waves were conducted in the first instance to verify this scheme. Furthermore, the comparison of the waveform variations confirms that the estimation of 3D solitary waves is a feasible scheme.