To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as i...To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as it can form a matrix with Tocplitz and Hankel properties. In this paper, a method is proposed to decompose the finite depth Green function into two terms, which can form matrices with the Toeplitz and a Hankel properties respectively. Then, a pFFT method for finite depth problems is developed. Based on the pFFT method, a numerical code pFFT-HOBEM is developed with the discretization of high order elements. The model is validated, and examinations on the computing efficiency and memory requirement of the new method have also been carried out. It shows that the new method has the same advantages as that for infinite depth.展开更多
Two coaxial vertical cylinders-one is a riding hollow cylinder and the other a solid cylinder of greater radius at some distance above an impermeable horizontal bottom,were considered.This problem of diffraction by th...Two coaxial vertical cylinders-one is a riding hollow cylinder and the other a solid cylinder of greater radius at some distance above an impermeable horizontal bottom,were considered.This problem of diffraction by these two cylinders,which were considered as idealization of a buoy and a circular plate,can be considered as a wave energy device.The wave energy that is created and transferred by this device can be appropriately used in many applications in lieu of conventional energy.Method of separation of variables was used to obtain the analytical expressions for the diffracted potentials in four clearly identified regions.By applying the appropriate matching conditions along the three virtual boundaries between the regions,a system of linear equations was obtained,which was solved for the unknown coefficients.The potentials allowed us to obtain the exciting forces acting on both cylinders.Sets of exciting forces were obtained for different radii of the cylinders and for different gaps between the cylinders.It was observed that changes in radius and the gap had significant effect on the forces.It was found that mostly the exciting forces were significant only at lower frequencies.The exciting forces almost vanished at higher frequencies.The problem was also investigated for the base case of no plate arrangement,i.e.,the case having only the floating cylinder tethered to the sea-bed.Comparison of forces for both arrangements was carried out.In order to take care of the radiation of the cylinders due to surge motion,the corresponding added mass and the damping coefficients for both cylinders were also computed.All the results were depicted graphically and compared with available results.展开更多
In this paper, function characteristics of dispersion of ocean wave in finite depth water were analyzed systematically. The functional form of the fitting function is reasonably proposed, in which the parame- ters are...In this paper, function characteristics of dispersion of ocean wave in finite depth water were analyzed systematically. The functional form of the fitting function is reasonably proposed, in which the parame- ters are optimally determined by the least square method (LSM). For infinitely deep and extremely shallow water, the fitting function fits strictly the dispersion to be fitted. A new technique is presented in application of LSM. An empirical formula with maximum error of less than 0.5% for computing wavelength in finite depth water is presented for practical applications.展开更多
Linear surface gravity waves on Maxwell viscoelastic fluids with finite depth are studied in this paper.A dispersion equation describing the spatial decay of the gravity wave in finite depth is derived.A dimensionless...Linear surface gravity waves on Maxwell viscoelastic fluids with finite depth are studied in this paper.A dispersion equation describing the spatial decay of the gravity wave in finite depth is derived.A dimensionless memory(time)number θ is introduced.The dispersion equation for the pure viscous fluid will be a specific case of the dispersion equation for the viscoelastic fluid as θ=0.The complex dispersion equation is numerically solved to investigate the dispersion relation.The influences of θ and water depth on the dispersion characteristics and wave decay are discussed.It is found that the role of elasticity for the Maxwell fluid is to make the surface gravity wave on the Maxwell fluid behave more like the surface gravity wave on the inviscid fluid.展开更多
An integral equation describing riding waves, i. e. , small-scaleperturbation waves superposed on unperturbed surface waves, in shallow water of finite depth wasstudied via explicit Hamiltonian formulation, and the wa...An integral equation describing riding waves, i. e. , small-scaleperturbation waves superposed on unperturbed surface waves, in shallow water of finite depth wasstudied via explicit Hamiltonian formulation, and the water was regarded as ideal incompressiblefluid of uniform density. The kinetic energy, density of the perturbed fluid motion was formulatedwith Hamiltonian canonical variables, elevation of the free surface and the velocity potential atthe free surface. Then the variables were expanded to the first order at the free surface ofunperturbed waves. An integal equation for velocity potential of perturbed waves on the unperturbedfree surface was derived by conformal mapping and the Fourier transformation. The integral equationcould replace the Hamiltonian canonical equations which are difficult to solve. An explicitexpression of Lagrangian density function could be obtained by solving the integral equation. Themethod used in this paper provides a new path to study the Hamiltonian formulation of riding wavesand wave interaction problems.展开更多
Nonlinear wave-body interactions for a stationary surface-piercing body in water of finite depth with fiat and sloping bottoms are simulated in a two-dimensional numerical wave tank, which is constructed mainly based ...Nonlinear wave-body interactions for a stationary surface-piercing body in water of finite depth with fiat and sloping bottoms are simulated in a two-dimensional numerical wave tank, which is constructed mainly based on the spatially averaged Navier-Stokes equations with the k- ε model for simulating the turbulence. The equations are discretized based on the finite volume method and the scheme of the pressure implicit splitting of operators is employed to solve the Navier-Stokes equations. By using the force time histories, the mean and higher-harmonic force components are calculated. The computational results are shown to be in good agreement with experimental and numerical results of other researchers. Then, the horizontal force, the vertical force and the moment on the surface-piercing body under nonlinear regular waves with flat and sloping bottoms are obtained. The results indicate that the bottom topographies have a significant influence on the wave loads on the surface-piercing body.展开更多
For computation of large amplitude motions of ships fastened to a dock, a fast evaluation scheme is implemented for computation of the time-domain Green function for finite water depth. Based on accurate evaluation of...For computation of large amplitude motions of ships fastened to a dock, a fast evaluation scheme is implemented for computation of the time-domain Green function for finite water depth. Based on accurate evaluation of the Green function directly, a fast approximation method for the Green function is developed by use of Chebyshev polynomials. Examinations are carried out of the accuracy of the Green function and its derivatives from the scheme. It is shown that when an appropriate number of polynomial terms are used, very accurate approximation can be obtained.展开更多
A nonlinear short-crested wave system, consisting of two progressive waves propagating at an oblique angle to each other in a fluid of finite depth, is investigated by means of an analytical approach named the homotop...A nonlinear short-crested wave system, consisting of two progressive waves propagating at an oblique angle to each other in a fluid of finite depth, is investigated by means of an analytical approach named the homotopy analysis method (HAM). Highly convergent series solutions are explicitly derived for the velocity potential and the surface wave elevation. We find that, at every value of water depth, there is little difference between the kinetic energy and the potential energy for nonlinear waves. The nonlinear short-crested waves with a larger angle of incidence always contain the more potential wave energy. With the aid of the HAM, we obtain the dispersion relation for nonlinear short-crested waves. Furthermore, it is shown that the wave elevation tends to be smoothened at the crest and be sharpened at the trough as the water depth increases, and the wave pressure crests and troughs become steeper with increasing incident wave steepness.展开更多
Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Co...Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Compressive re- sidual stress and dent profile are important factors to eval- uate the effectiveness of shot peening process. In this pa- per, the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated. Firstly, dimensionless relations of pro- cessing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were de- duced by dimensional analysis method. Secondly, the in- fluence of each dimensionless parameter on dimensionless variables was investigated by the finite element method. Fur- thermore, related empirical formulas were given for each di- mensionless parameter based on the simulation results. Fi- nally, comparison was made and good agreement was found between the simulation results and the empirical formula, which shows that a useful approach is provided in this pa- per for analyzing the influence of each individual parameter.展开更多
A fully three-dimensional surface gravitycapillary short-crested wave system is studied as two progressive wave-trains of equal amplitude and frequency, which are collinear with uniform currents and doubly-periodic in...A fully three-dimensional surface gravitycapillary short-crested wave system is studied as two progressive wave-trains of equal amplitude and frequency, which are collinear with uniform currents and doubly-periodic in the horizontal plane, are propagating at an angle to each other. The first- and second-order asymptotic analytical solutions of the short-crested wave system are obtained via a perturbation expansion in a small parameter associated with the wave steepness, therefore depicting a series of typical three-dimensional wave patterns involving currents, shallow and deep water, and surface capillary waves, and comparing them with each other.展开更多
Very Large Floating Structures (VLFS) have drawn considerable attention recently due to their potential significance in the exploitation of ocean resources and in the utilization of ocean space. Efficient and accurate...Very Large Floating Structures (VLFS) have drawn considerable attention recently due to their potential significance in the exploitation of ocean resources and in the utilization of ocean space. Efficient and accurate estimation of their hydroelastic responses to waves is very important for the design. Recently, an efficient numerical algorithm was developed by Ertekin and Kim (1999). However, in their analysis, the linear Level I Green-Naghdi (GN) theory is employed to describe fluid dynamics instead of the conventional linear wave (LW) theory of finite water depth. They claimed that this linear level I GN theory provided better predictions of the hydroelastic responses of VLFS than the linear wave theory. In this paper, a detailed derivation is given in the conventional linear wave theory framework with the same quantity as used in the linear level I GN theory framework. This allows a critical comparison between the linear wave theory and the linear level I GN theory. It is found that the linear level I GN theory can be regarded as an approximation to the linear wave theory of finite water depth. The consequences of the differences between these two theories in the predicted hydroelastic responses are studied quantitatively. And it is found that the linear level I GN theory is not superior to the linear wave theory. Finally, various factors affecting the hydroelastic response of VLFS are studied with the implemented algorithm.展开更多
The spectral form of wind waves is investigated based on the ocean wave data observed at three nearshore stations of Taiwan. In this study, the generalized forms of Pierson-Moskowitz spectrum and JONSWAP spectrum are ...The spectral form of wind waves is investigated based on the ocean wave data observed at three nearshore stations of Taiwan. In this study, the generalized forms of Pierson-Moskowitz spectrum and JONSWAP spectrum are used to describe the local wave spectrum by selecting suitable spectral form parameters. It is shown that, at a specific site, the similarity of wave spectral form exists. Thus it is possible to use a representative spectral form for a given nearshore region to describe the wave spectrum at this nearshore. On the other hand, the effects of relative water depth on spectral form are examined. The feasibility of two spectral models in finite water depth is evaluated by using the same field wave data.展开更多
The flexural-gravity wave responses due to a load steadily moving or suddenly accelerated along a rectilinear orbit are analytically studied within the framework of the linear potential theory. A thin viscoelastic pla...The flexural-gravity wave responses due to a load steadily moving or suddenly accelerated along a rectilinear orbit are analytically studied within the framework of the linear potential theory. A thin viscoelastic plate model is used for a very large floating structure. The initially quiescent fluid in the ocean is assumed to be homogenous, incompressible, and inviscid, and the disturbed motion be irrotational. A moving line source on the plate surface is considered as a moving point in the two-dimensional coordinates. Under the assumptions of small-amplitude wave motion and small plate deflection, a linear fluid-plate coupling model is established. The integral solutions for the surface deflections and the wave resistances are analytically obtained by the Fourier transform method. To study the dynamic characteristics of the flexural-gravity wave response, the asymptotic representations of the wave resistances are derived by the residue theorem and the methods of stationary phase. It shows that the steady wave resistance is zero when the speed of moving load is less than the minimal phase speed. The wave resistances due to the accelerate motion consist of two parts, namely the steady and transient wave responses. Eventually the transient wave resistance declines toward zero and the wave resistance approaches the steady component as the time goes to the infinity. Furthermore, the effect of the strain relaxation time for this viscoelastic plate is studied and it exhibits more influence for a high-speed motion.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.51490672 and 51379032)
文摘To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as it can form a matrix with Tocplitz and Hankel properties. In this paper, a method is proposed to decompose the finite depth Green function into two terms, which can form matrices with the Toeplitz and a Hankel properties respectively. Then, a pFFT method for finite depth problems is developed. Based on the pFFT method, a numerical code pFFT-HOBEM is developed with the discretization of high order elements. The model is validated, and examinations on the computing efficiency and memory requirement of the new method have also been carried out. It shows that the new method has the same advantages as that for infinite depth.
文摘Two coaxial vertical cylinders-one is a riding hollow cylinder and the other a solid cylinder of greater radius at some distance above an impermeable horizontal bottom,were considered.This problem of diffraction by these two cylinders,which were considered as idealization of a buoy and a circular plate,can be considered as a wave energy device.The wave energy that is created and transferred by this device can be appropriately used in many applications in lieu of conventional energy.Method of separation of variables was used to obtain the analytical expressions for the diffracted potentials in four clearly identified regions.By applying the appropriate matching conditions along the three virtual boundaries between the regions,a system of linear equations was obtained,which was solved for the unknown coefficients.The potentials allowed us to obtain the exciting forces acting on both cylinders.Sets of exciting forces were obtained for different radii of the cylinders and for different gaps between the cylinders.It was observed that changes in radius and the gap had significant effect on the forces.It was found that mostly the exciting forces were significant only at lower frequencies.The exciting forces almost vanished at higher frequencies.The problem was also investigated for the base case of no plate arrangement,i.e.,the case having only the floating cylinder tethered to the sea-bed.Comparison of forces for both arrangements was carried out.In order to take care of the radiation of the cylinders due to surge motion,the corresponding added mass and the damping coefficients for both cylinders were also computed.All the results were depicted graphically and compared with available results.
基金This study was financially supported by Key Technologies R & D Programme of China (No.2001BA603B-01).
文摘In this paper, function characteristics of dispersion of ocean wave in finite depth water were analyzed systematically. The functional form of the fitting function is reasonably proposed, in which the parame- ters are optimally determined by the least square method (LSM). For infinitely deep and extremely shallow water, the fitting function fits strictly the dispersion to be fitted. A new technique is presented in application of LSM. An empirical formula with maximum error of less than 0.5% for computing wavelength in finite depth water is presented for practical applications.
基金The project supported by the National Natural Science Foundation of China(50279029)
文摘Linear surface gravity waves on Maxwell viscoelastic fluids with finite depth are studied in this paper.A dispersion equation describing the spatial decay of the gravity wave in finite depth is derived.A dimensionless memory(time)number θ is introduced.The dispersion equation for the pure viscous fluid will be a specific case of the dispersion equation for the viscoelastic fluid as θ=0.The complex dispersion equation is numerically solved to investigate the dispersion relation.The influences of θ and water depth on the dispersion characteristics and wave decay are discussed.It is found that the role of elasticity for the Maxwell fluid is to make the surface gravity wave on the Maxwell fluid behave more like the surface gravity wave on the inviscid fluid.
文摘An integral equation describing riding waves, i. e. , small-scaleperturbation waves superposed on unperturbed surface waves, in shallow water of finite depth wasstudied via explicit Hamiltonian formulation, and the water was regarded as ideal incompressiblefluid of uniform density. The kinetic energy, density of the perturbed fluid motion was formulatedwith Hamiltonian canonical variables, elevation of the free surface and the velocity potential atthe free surface. Then the variables were expanded to the first order at the free surface ofunperturbed waves. An integal equation for velocity potential of perturbed waves on the unperturbedfree surface was derived by conformal mapping and the Fourier transformation. The integral equationcould replace the Hamiltonian canonical equations which are difficult to solve. An explicitexpression of Lagrangian density function could be obtained by solving the integral equation. Themethod used in this paper provides a new path to study the Hamiltonian formulation of riding wavesand wave interaction problems.
基金Project supported by the National Natural Science Foundation of China (Grant No. 40776057)the Knowledge Innovation Program of Chinese Academy of Sciences (Grant Nos. KJCX2-YW-L07,KZCX2-YW-212-2)
文摘Nonlinear wave-body interactions for a stationary surface-piercing body in water of finite depth with fiat and sloping bottoms are simulated in a two-dimensional numerical wave tank, which is constructed mainly based on the spatially averaged Navier-Stokes equations with the k- ε model for simulating the turbulence. The equations are discretized based on the finite volume method and the scheme of the pressure implicit splitting of operators is employed to solve the Navier-Stokes equations. By using the force time histories, the mean and higher-harmonic force components are calculated. The computational results are shown to be in good agreement with experimental and numerical results of other researchers. Then, the horizontal force, the vertical force and the moment on the surface-piercing body under nonlinear regular waves with flat and sloping bottoms are obtained. The results indicate that the bottom topographies have a significant influence on the wave loads on the surface-piercing body.
文摘For computation of large amplitude motions of ships fastened to a dock, a fast evaluation scheme is implemented for computation of the time-domain Green function for finite water depth. Based on accurate evaluation of the Green function directly, a fast approximation method for the Green function is developed by use of Chebyshev polynomials. Examinations are carried out of the accuracy of the Green function and its derivatives from the scheme. It is shown that when an appropriate number of polynomial terms are used, very accurate approximation can be obtained.
基金Supported by the National Key Basic Research Development Program of China(973 Program,Grant No.2014CB046203)the National Natural Science Foundation of China(Grant No.11472166)the Natural Science Founda-tion of Shanghai(Grant No.14ZR1416200)
文摘A nonlinear short-crested wave system, consisting of two progressive waves propagating at an oblique angle to each other in a fluid of finite depth, is investigated by means of an analytical approach named the homotopy analysis method (HAM). Highly convergent series solutions are explicitly derived for the velocity potential and the surface wave elevation. We find that, at every value of water depth, there is little difference between the kinetic energy and the potential energy for nonlinear waves. The nonlinear short-crested waves with a larger angle of incidence always contain the more potential wave energy. With the aid of the HAM, we obtain the dispersion relation for nonlinear short-crested waves. Furthermore, it is shown that the wave elevation tends to be smoothened at the crest and be sharpened at the trough as the water depth increases, and the wave pressure crests and troughs become steeper with increasing incident wave steepness.
基金supported by the National Natural Science Foun-dation of China (10972228,11002150,and 91016025)the Basic Research Equipment Project of Chinese Academy of Sciences(YZ200930)
文摘Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Compressive re- sidual stress and dent profile are important factors to eval- uate the effectiveness of shot peening process. In this pa- per, the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated. Firstly, dimensionless relations of pro- cessing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were de- duced by dimensional analysis method. Secondly, the in- fluence of each dimensionless parameter on dimensionless variables was investigated by the finite element method. Fur- thermore, related empirical formulas were given for each di- mensionless parameter based on the simulation results. Fi- nally, comparison was made and good agreement was found between the simulation results and the empirical formula, which shows that a useful approach is provided in this pa- per for analyzing the influence of each individual parameter.
基金The project supported by the Foundation for the Author of National Excellent Doctoral Dissertation of China (200428)the National Natural Science Foundation of China (10272072and 50424913)the Shanghai Natural Science Foundation(05ZR14048)
文摘A fully three-dimensional surface gravitycapillary short-crested wave system is studied as two progressive wave-trains of equal amplitude and frequency, which are collinear with uniform currents and doubly-periodic in the horizontal plane, are propagating at an angle to each other. The first- and second-order asymptotic analytical solutions of the short-crested wave system are obtained via a perturbation expansion in a small parameter associated with the wave steepness, therefore depicting a series of typical three-dimensional wave patterns involving currents, shallow and deep water, and surface capillary waves, and comparing them with each other.
基金by the National Natural Science Foundation of China(50039010)the Science and Technology Development Foundation of Shanghai Municipal Government(00XD14015)
文摘Very Large Floating Structures (VLFS) have drawn considerable attention recently due to their potential significance in the exploitation of ocean resources and in the utilization of ocean space. Efficient and accurate estimation of their hydroelastic responses to waves is very important for the design. Recently, an efficient numerical algorithm was developed by Ertekin and Kim (1999). However, in their analysis, the linear Level I Green-Naghdi (GN) theory is employed to describe fluid dynamics instead of the conventional linear wave (LW) theory of finite water depth. They claimed that this linear level I GN theory provided better predictions of the hydroelastic responses of VLFS than the linear wave theory. In this paper, a detailed derivation is given in the conventional linear wave theory framework with the same quantity as used in the linear level I GN theory framework. This allows a critical comparison between the linear wave theory and the linear level I GN theory. It is found that the linear level I GN theory can be regarded as an approximation to the linear wave theory of finite water depth. The consequences of the differences between these two theories in the predicted hydroelastic responses are studied quantitatively. And it is found that the linear level I GN theory is not superior to the linear wave theory. Finally, various factors affecting the hydroelastic response of VLFS are studied with the implemented algorithm.
文摘The spectral form of wind waves is investigated based on the ocean wave data observed at three nearshore stations of Taiwan. In this study, the generalized forms of Pierson-Moskowitz spectrum and JONSWAP spectrum are used to describe the local wave spectrum by selecting suitable spectral form parameters. It is shown that, at a specific site, the similarity of wave spectral form exists. Thus it is possible to use a representative spectral form for a given nearshore region to describe the wave spectrum at this nearshore. On the other hand, the effects of relative water depth on spectral form are examined. The feasibility of two spectral models in finite water depth is evaluated by using the same field wave data.
基金Project supported by the National Natural Science Foundation of China(Grant No.11472166)the State Key Laboratory of Hydraulic Engineering Simulation and Safety at Tianjin University(Grant No.HESS-1506)the National Key Basic Research Development Program of China(973 Program,Grant No.2014 CB046203)
文摘The flexural-gravity wave responses due to a load steadily moving or suddenly accelerated along a rectilinear orbit are analytically studied within the framework of the linear potential theory. A thin viscoelastic plate model is used for a very large floating structure. The initially quiescent fluid in the ocean is assumed to be homogenous, incompressible, and inviscid, and the disturbed motion be irrotational. A moving line source on the plate surface is considered as a moving point in the two-dimensional coordinates. Under the assumptions of small-amplitude wave motion and small plate deflection, a linear fluid-plate coupling model is established. The integral solutions for the surface deflections and the wave resistances are analytically obtained by the Fourier transform method. To study the dynamic characteristics of the flexural-gravity wave response, the asymptotic representations of the wave resistances are derived by the residue theorem and the methods of stationary phase. It shows that the steady wave resistance is zero when the speed of moving load is less than the minimal phase speed. The wave resistances due to the accelerate motion consist of two parts, namely the steady and transient wave responses. Eventually the transient wave resistance declines toward zero and the wave resistance approaches the steady component as the time goes to the infinity. Furthermore, the effect of the strain relaxation time for this viscoelastic plate is studied and it exhibits more influence for a high-speed motion.