A theory on the second order wave diffraction by a three dimensional body fixed in a regular sea has been developed in the present paper. By regarding the sinusoidal disturb potential as a stationary solu- tion of an ...A theory on the second order wave diffraction by a three dimensional body fixed in a regular sea has been developed in the present paper. By regarding the sinusoidal disturb potential as a stationary solu- tion of an initial value problem, and using Laplace transformation method and Tauberian theorem, the boundary value problems of stationary solution of the first and second order diffraction potential have been de- rived in this paper. Furthermore, the explicit solution of the second order stationary diffraction potential has been obtained with the method of the double Fourier transformation. It is found that the asymptotic behaviour of the second order stationary solution at far field is dependent on two wave systems, the first is 'free wave', travelling independently of the first order wave system, the other is 'phase locked waves', which accompany the first order waves. At the same time, the radiation conditions of the second order diffraction problems are derived. We also find that one can still pursue a steady state formulation with the inclusion of Rayleigh damping. Finally, as an example, the second order wave forces upon a fixed vertical cir- cular cylinder have been calculated, and the numerical results agree well with the experimental data.展开更多
In this paper, the principle of mirror image is used to transform the problem of wave diffraction from a circular cylinder in front of orthogonal vertical walls into the problem of diffraction of four symmetric incide...In this paper, the principle of mirror image is used to transform the problem of wave diffraction from a circular cylinder in front of orthogonal vertical walls into the problem of diffraction of four symmetric incident waves from four symmetrically arranged circular cylinders, and then the eigenfunction expansion of velocity potential and Grafs addition theorem are used to give the analytical solution to the wave diffraction problem. The relation of the total wave force on cylinder to the distance between the cylinder and orthogonal vertical walls and the incidence angle of wave is also studied by numerical computation.展开更多
This paper presents a study on the motion response of a tension-leg platform(TLP) under first-and second-order wave forces, including the mean-drift force, difference and sum-frequency forces. The second-order wave fo...This paper presents a study on the motion response of a tension-leg platform(TLP) under first-and second-order wave forces, including the mean-drift force, difference and sum-frequency forces. The second-order wave force is calculated using the full-field quadratic transfer function(QTF). The coupled effect of the horizontal motions, such as surge, sway and yaw motions, and the set-down motion are taken into consideration by the nonlinear restoring matrix. The time-domain analysis with 50-yr random sea state is performed. A comparison of the results of different case studies is made to assess the influence of second-order wave force on the motions of the platform. The analysis shows that the second-order wave force has a major impact on motions of the TLP. The second-order difference-frequency wave force has an obvious influence on the low-frequency motions of surge and sway, and also will induce a large set-down motion which is an important part of heave motion. Besides, the second-order sum-frequency force will induce a set of high-frequency motions of roll and pitch. However, little influence of second-order wave force is found on the yaw motion.展开更多
So far many investigations have been made on nonlinear wave diffraction problemfor a large-diameter vertical circular cylinder. However,there are still some problemsworthy to be further discussed. It includes that the...So far many investigations have been made on nonlinear wave diffraction problemfor a large-diameter vertical circular cylinder. However,there are still some problemsworthy to be further discussed. It includes that the second order radiation condition isnot very clear and the inhomogeneous term of second order free surface boundarycondition makes the calculation of second order wave force either not easier to beperformed accurately due to its slowly decaying with radial distance or toocomplicated for practical application. In this paper, the second order radiationcondition is posed of the circumferential Fourier components of second orderpotential, instead of the second order potential. It is found that the circumferenatialFourier cormponents of second order potential have to satisfy Sommerfeld radiationcondition. By means of the mathematical formulae derived in this paper, theinhomogeneous term of second order free surface boundary condition were simplifiedand then an exact expression of second order wave force was obtained, which issimpler in form and easier to be used in practical calculation. The calculation resultsagree well with some experimental data.展开更多
A complete semi-analytical solution is obtained for second-order diffraction of plane bichromatic waves by a fixed truncated circular column.The fluid domain is divided into interior and exterior regions.In the exteri...A complete semi-analytical solution is obtained for second-order diffraction of plane bichromatic waves by a fixed truncated circular column.The fluid domain is divided into interior and exterior regions.In the exterior region,the second-order velocity potential is expressed in terms of‘locked-wave’and‘free-wave’ components,both are solved using Fourier and eigenfunction expansions.The re- sulting‘locked wave’potential is expressed by one-dimensional Green's integrals with oscillating integrands.In order to increase computational efficiency,the far-field part of the integrals are carried out analytically.Solutions in both regions are matched on the interface by the potential and its normal derivative continuity conditions.Based on the present approach,the sum-and difference-frequency potentials are efficiently evaluated and are used to generate the quadratic transfer functions which correlates the incident wave spectrum with second-order forcing spectrum on the column.The sum-frequency QTFs for a TLP column are present,which are compared for some frequency pairs with those from a fully numerical procedure.Satisfactory agreement has been obtained.QTF spectra for a case study TLP column,generated using the semi-analytical solution are presented.Also given are the results for nonlinear wave field around the column.展开更多
By using the theories on Stokes multicolored water waves and taking the two- layer ocean as a basic model of stratified ocean, the paper analyzes the problems related to the effects of the nonlinear water wave on offs...By using the theories on Stokes multicolored water waves and taking the two- layer ocean as a basic model of stratified ocean, the paper analyzes the problems related to the effects of the nonlinear water wave on offshore structures. A mathematical expression is presented to describe second order wave radiation conditions. Using integral principle, the analytical integral solutions are given to evaluate second order scattered wave loads on general vertical circular cylinders in the two-layer ocean, and the special recurrence formulas for infinite integrals over free and stratified surfaces are derived.展开更多
The power performances of a point absorber wave energy converter(WEC)operating in a nonlinear multidirectional random sea are rigorously investigated.The absorbed power of the WEC Power-Take-Off system has been predic...The power performances of a point absorber wave energy converter(WEC)operating in a nonlinear multidirectional random sea are rigorously investigated.The absorbed power of the WEC Power-Take-Off system has been predicted by incorporating a second order random wave model into a nonlinear dynamic filter.This is a new approach,and,as the second order random wave model can be utilized to accurately simulate the nonlinear waves in an irregular sea,avoids the inaccuracies resulting from using a first order linear wave model in the simulation process.The predicted results have been systematically analyzed and compared,and the advantages of using this new approach have been convincingly substantiated.展开更多
This paper studies the random internal wave equations describing the density interface displacements and the velocity potentials of N-layer stratified fluid contained between two rigid walls at the top and bottom. The...This paper studies the random internal wave equations describing the density interface displacements and the velocity potentials of N-layer stratified fluid contained between two rigid walls at the top and bottom. The density interface displacements and the velocity potentials were solved to the second-order by an expansion approach used by Longuet-Higgins (1963) and Dean (1979) in the study of random surface waves and by Song (2004) in the study of second- order random wave solutions for internal waves in a two-layer fluid. The obtained results indicate that the first-order solutions are a linear superposition of many wave components with different amplitudes, wave numbers and frequencies, and that the amplitudes of first-order wave components with the same wave numbers and frequencies between the adjacent density interfaces are modulated by each other. They also show that the second-order solutions consist of two parts: the first one is the first-order solutions, and the second one is the solutions of the second-order asymptotic equations, which describe the second-order nonlinear modification and the second-order wave-wave interactions not only among the wave components on same density interfaces but also among the wave components between the adjacent density interfaces. Both the first-order and second-order solutions depend on the density and depth of each layer. It is also deduced that the results of the present work include those derived by Song (2004) for second-order random wave solutions for internal waves in a two-layer fluid as a particular case.展开更多
In the present research, the study of Song (2004) for random interfacial waves in two-layer fluid is extended to the case of fluids moving at different steady uniform speeds. The equations describing the random displa...In the present research, the study of Song (2004) for random interfacial waves in two-layer fluid is extended to the case of fluids moving at different steady uniform speeds. The equations describing the random displacements of the density interface and the associated velocity potentials in two-layer fluid are solved to the second order, and the wave-wave interactions of the wave components and the interactions between the waves and currents are described. As expected, the extended solutions include those obtained by Song (2004) as one special case where the steady uniform currents of the two fluids are taken as zero, and the solutions reduce to those derived by Sharma and Dean (1979) for random surface waves if the density of the upper fluid and the current of the lower fluid are both taken as zero.展开更多
A previous study (Song. 2004. Geophys Res Lett, 31 (15):L15302) of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into...A previous study (Song. 2004. Geophys Res Lett, 31 (15):L15302) of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song (2004. Geophys Res Lett, 31 (15):L15302). The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves (Thorpe. 1968a.Trans R Soc London, 263A:563~614) and standing interracial waves (Thorpe. 1968b. J Fluid Mech, 32:489-528) for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean (1979.Ocean Engineering Rep 20) if the density of the upper layer is much smaller than that of the lower layer.展开更多
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 commercial 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.展开更多
-In this paper, an analytical solution in the outer region of finite water depth is derived for the second-order diffraction potential, which gives a clear physical meaning of the wave transmission and reflection char...-In this paper, an analytical solution in the outer region of finite water depth is derived for the second-order diffraction potential, which gives a clear physical meaning of the wave transmission and reflection characteristics in the far field. A numerical method-simple Green's function technique-for calculating the second-order diffraction potential in the inner region is also described. Numerical results are provided for the second-order wave forces on a semi-submerged cylinder. It is found that the contribution of second-order diffraction potential to second-order wave forces is important. The effect of water depth and submerged depth on the wave force is also discussed.展开更多
With growing computational power, the first-order wave-maker theory has become well established and is widely used for numerical wave flumes. However, existing numerical models based on the first-order wave-maker theo...With growing computational power, the first-order wave-maker theory has become well established and is widely used for numerical wave flumes. However, existing numerical models based on the first-order wave-maker theory lose accuracy as nonlinear effects become prominent. Because spurious harmonic waves and primary waves have different propagation velocities, waves simulated by using the first-order wave-maker theory have an unstable wave profile. In this paper, a numerical wave flume with a piston-type wave-maker based on the second-order wave-maker theory has been established. Dynamic mesh technique was developed. The boundary treatment for irregular wave simulation was specially dealt with. Comparisons of the free-surface elevations using the first-order and second-order wave-maker theory prove that second-order wave-maker theory can generate stable wave profiles in both the spatial and time domains. Harmonic analysis and spectral analysis were used to prove the superiority of the second-order wave-maker theory from other two aspects. To simulate irregular waves, the numerical flume was improved to solve the problem of the water depth variation due to low-frequency motion of the wave board. In summary, the new numerical flume using the second-order wave-maker theory can guarantee the accuracy of waves by adding an extra motion of the wave board. The boundary treatment method can provide a reference for the improvement of nonlinear numerical flume.展开更多
On the basis of the second order solution of two-dimensional random gravity waves in finite uniform depth,which is rederived by a perturbation expansion method,the analytical expression of the second order spectrum is...On the basis of the second order solution of two-dimensional random gravity waves in finite uniform depth,which is rederived by a perturbation expansion method,the analytical expression of the second order spectrum is strictly deduced, and for infinite depth, the correct form of the kernel function is given for each octant do main.In fact,the present study improves and corrects the generally accepted results obtained by Tick展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
Based on the second order random wave solutions of water wave equations in finite water depth, statistical distributions of the depth integrated local horizontal momentum components are derived by use of the charact...Based on the second order random wave solutions of water wave equations in finite water depth, statistical distributions of the depth integrated local horizontal momentum components are derived by use of the characteristic function expansion method. The parameters involved in the distributions can be all determined by the water depth and the wave number spectrum of ocean waves. As an illustrative example, a fully developed wind generated sea is considered and the parameters are calculated for typical wind speeds and water depths by means of the Donelan and Pierson spectrum. The effects of nonlinearity and water depth on the distributions are also investigated.展开更多
Linear wave theory and Longuet-Higgins and Steward’s (1964) group-induced second-order longweve (GSLW) theory ware used in this study on the grouping effect on wave forces acting on a verticalbreakwater. The calculat...Linear wave theory and Longuet-Higgins and Steward’s (1964) group-induced second-order longweve (GSLW) theory ware used in this study on the grouping effect on wave forces acting on a verticalbreakwater. The calculated variance of total wave pressure on the vertical breakwater was closer tothe measured value if the wave grouping effect was considered.展开更多
A 3-D time-domain seakeeping analysis tool has been newly developed by using a higher-order boundary element method with the Rankine source as the kernel function. An iterative time-marching scheme for updating both k...A 3-D time-domain seakeeping analysis tool has been newly developed by using a higher-order boundary element method with the Rankine source as the kernel function. An iterative time-marching scheme for updating both kinematic and dynamic free-surface boundary conditions is adopted for achieving numerical accuracy and stability. A rectangular computational domain moving with the mean speed of ship is introduced. A damping beach at the outer portion of the truncated free surface is installed for satisfying the radiation condition. After numerical convergence checked, the diffraction unsteady problem of a Wigley hull traveling with a constant forward speed in waves is studied. Extensive results including wave exciting forces, wave patterns and pressure distributions on the hull are presented to validate the efficiency and accuracy of the proposed 3-D time-domain iterative Rankine BEM approach. Computed results are compared to be in good agreement with the corresponding experimental data and other published numerical solutions.展开更多
文摘A theory on the second order wave diffraction by a three dimensional body fixed in a regular sea has been developed in the present paper. By regarding the sinusoidal disturb potential as a stationary solu- tion of an initial value problem, and using Laplace transformation method and Tauberian theorem, the boundary value problems of stationary solution of the first and second order diffraction potential have been de- rived in this paper. Furthermore, the explicit solution of the second order stationary diffraction potential has been obtained with the method of the double Fourier transformation. It is found that the asymptotic behaviour of the second order stationary solution at far field is dependent on two wave systems, the first is 'free wave', travelling independently of the first order wave system, the other is 'phase locked waves', which accompany the first order waves. At the same time, the radiation conditions of the second order diffraction problems are derived. We also find that one can still pursue a steady state formulation with the inclusion of Rayleigh damping. Finally, as an example, the second order wave forces upon a fixed vertical cir- cular cylinder have been calculated, and the numerical results agree well with the experimental data.
基金financially supported by the National Natural Science Foundation of China(Grant No.50025924).
文摘In this paper, the principle of mirror image is used to transform the problem of wave diffraction from a circular cylinder in front of orthogonal vertical walls into the problem of diffraction of four symmetric incident waves from four symmetrically arranged circular cylinders, and then the eigenfunction expansion of velocity potential and Grafs addition theorem are used to give the analytical solution to the wave diffraction problem. The relation of the total wave force on cylinder to the distance between the cylinder and orthogonal vertical walls and the incidence angle of wave is also studied by numerical computation.
基金supported by the National Natural Science Foundation of China(Nos.51239008 and 51279130)
文摘This paper presents a study on the motion response of a tension-leg platform(TLP) under first-and second-order wave forces, including the mean-drift force, difference and sum-frequency forces. The second-order wave force is calculated using the full-field quadratic transfer function(QTF). The coupled effect of the horizontal motions, such as surge, sway and yaw motions, and the set-down motion are taken into consideration by the nonlinear restoring matrix. The time-domain analysis with 50-yr random sea state is performed. A comparison of the results of different case studies is made to assess the influence of second-order wave force on the motions of the platform. The analysis shows that the second-order wave force has a major impact on motions of the TLP. The second-order difference-frequency wave force has an obvious influence on the low-frequency motions of surge and sway, and also will induce a large set-down motion which is an important part of heave motion. Besides, the second-order sum-frequency force will induce a set of high-frequency motions of roll and pitch. However, little influence of second-order wave force is found on the yaw motion.
文摘So far many investigations have been made on nonlinear wave diffraction problemfor a large-diameter vertical circular cylinder. However,there are still some problemsworthy to be further discussed. It includes that the second order radiation condition isnot very clear and the inhomogeneous term of second order free surface boundarycondition makes the calculation of second order wave force either not easier to beperformed accurately due to its slowly decaying with radial distance or toocomplicated for practical application. In this paper, the second order radiationcondition is posed of the circumferential Fourier components of second orderpotential, instead of the second order potential. It is found that the circumferenatialFourier cormponents of second order potential have to satisfy Sommerfeld radiationcondition. By means of the mathematical formulae derived in this paper, theinhomogeneous term of second order free surface boundary condition were simplifiedand then an exact expression of second order wave force was obtained, which issimpler in form and easier to be used in practical calculation. The calculation resultsagree well with some experimental data.
文摘A complete semi-analytical solution is obtained for second-order diffraction of plane bichromatic waves by a fixed truncated circular column.The fluid domain is divided into interior and exterior regions.In the exterior region,the second-order velocity potential is expressed in terms of‘locked-wave’and‘free-wave’ components,both are solved using Fourier and eigenfunction expansions.The re- sulting‘locked wave’potential is expressed by one-dimensional Green's integrals with oscillating integrands.In order to increase computational efficiency,the far-field part of the integrals are carried out analytically.Solutions in both regions are matched on the interface by the potential and its normal derivative continuity conditions.Based on the present approach,the sum-and difference-frequency potentials are efficiently evaluated and are used to generate the quadratic transfer functions which correlates the incident wave spectrum with second-order forcing spectrum on the column.The sum-frequency QTFs for a TLP column are present,which are compared for some frequency pairs with those from a fully numerical procedure.Satisfactory agreement has been obtained.QTF spectra for a case study TLP column,generated using the semi-analytical solution are presented.Also given are the results for nonlinear wave field around the column.
基金National Natural Science Foundation of China (19802023)
文摘By using the theories on Stokes multicolored water waves and taking the two- layer ocean as a basic model of stratified ocean, the paper analyzes the problems related to the effects of the nonlinear water wave on offshore structures. A mathematical expression is presented to describe second order wave radiation conditions. Using integral principle, the analytical integral solutions are given to evaluate second order scattered wave loads on general vertical circular cylinders in the two-layer ocean, and the special recurrence formulas for infinite integrals over free and stratified surfaces are derived.
基金The National Natural Science Foundation of China under contract No.51979165。
文摘The power performances of a point absorber wave energy converter(WEC)operating in a nonlinear multidirectional random sea are rigorously investigated.The absorbed power of the WEC Power-Take-Off system has been predicted by incorporating a second order random wave model into a nonlinear dynamic filter.This is a new approach,and,as the second order random wave model can be utilized to accurately simulate the nonlinear waves in an irregular sea,avoids the inaccuracies resulting from using a first order linear wave model in the simulation process.The predicted results have been systematically analyzed and compared,and the advantages of using this new approach have been convincingly substantiated.
基金Project supported by the National Science Fund for Distinguished Young Scholars (Grant No 40425015), the Cooperative Project of Chinese Academy Sciences and the China National 0ffshore oil Corporation ("Behaviours of internal waves and their roles on the marine structures") and the National Natural Science Foundation of China (Grant No10461005).
文摘This paper studies the random internal wave equations describing the density interface displacements and the velocity potentials of N-layer stratified fluid contained between two rigid walls at the top and bottom. The density interface displacements and the velocity potentials were solved to the second-order by an expansion approach used by Longuet-Higgins (1963) and Dean (1979) in the study of random surface waves and by Song (2004) in the study of second- order random wave solutions for internal waves in a two-layer fluid. The obtained results indicate that the first-order solutions are a linear superposition of many wave components with different amplitudes, wave numbers and frequencies, and that the amplitudes of first-order wave components with the same wave numbers and frequencies between the adjacent density interfaces are modulated by each other. They also show that the second-order solutions consist of two parts: the first one is the first-order solutions, and the second one is the solutions of the second-order asymptotic equations, which describe the second-order nonlinear modification and the second-order wave-wave interactions not only among the wave components on same density interfaces but also among the wave components between the adjacent density interfaces. Both the first-order and second-order solutions depend on the density and depth of each layer. It is also deduced that the results of the present work include those derived by Song (2004) for second-order random wave solutions for internal waves in a two-layer fluid as a particular case.
文摘In the present research, the study of Song (2004) for random interfacial waves in two-layer fluid is extended to the case of fluids moving at different steady uniform speeds. The equations describing the random displacements of the density interface and the associated velocity potentials in two-layer fluid are solved to the second order, and the wave-wave interactions of the wave components and the interactions between the waves and currents are described. As expected, the extended solutions include those obtained by Song (2004) as one special case where the steady uniform currents of the two fluids are taken as zero, and the solutions reduce to those derived by Sharma and Dean (1979) for random surface waves if the density of the upper fluid and the current of the lower fluid are both taken as zero.
基金supported by the National Science Foundation for Distinguished Young Scholars of China under contract No.40425015the Cooperative Project of Chinese Academy Sciences and the China National 0ffshore 0il Corporation("Behaviours of internal waves and their roles on the marine stuctures").
文摘A previous study (Song. 2004. Geophys Res Lett, 31 (15):L15302) of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song (2004. Geophys Res Lett, 31 (15):L15302). The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves (Thorpe. 1968a.Trans R Soc London, 263A:563~614) and standing interracial waves (Thorpe. 1968b. J Fluid Mech, 32:489-528) for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean (1979.Ocean Engineering Rep 20) if the density of the upper layer is much smaller than that of the lower layer.
基金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 commercial 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.
文摘-In this paper, an analytical solution in the outer region of finite water depth is derived for the second-order diffraction potential, which gives a clear physical meaning of the wave transmission and reflection characteristics in the far field. A numerical method-simple Green's function technique-for calculating the second-order diffraction potential in the inner region is also described. Numerical results are provided for the second-order wave forces on a semi-submerged cylinder. It is found that the contribution of second-order diffraction potential to second-order wave forces is important. The effect of water depth and submerged depth on the wave force is also discussed.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51579038,51739010,51490672,51879037)
文摘With growing computational power, the first-order wave-maker theory has become well established and is widely used for numerical wave flumes. However, existing numerical models based on the first-order wave-maker theory lose accuracy as nonlinear effects become prominent. Because spurious harmonic waves and primary waves have different propagation velocities, waves simulated by using the first-order wave-maker theory have an unstable wave profile. In this paper, a numerical wave flume with a piston-type wave-maker based on the second-order wave-maker theory has been established. Dynamic mesh technique was developed. The boundary treatment for irregular wave simulation was specially dealt with. Comparisons of the free-surface elevations using the first-order and second-order wave-maker theory prove that second-order wave-maker theory can generate stable wave profiles in both the spatial and time domains. Harmonic analysis and spectral analysis were used to prove the superiority of the second-order wave-maker theory from other two aspects. To simulate irregular waves, the numerical flume was improved to solve the problem of the water depth variation due to low-frequency motion of the wave board. In summary, the new numerical flume using the second-order wave-maker theory can guarantee the accuracy of waves by adding an extra motion of the wave board. The boundary treatment method can provide a reference for the improvement of nonlinear numerical flume.
文摘On the basis of the second order solution of two-dimensional random gravity waves in finite uniform depth,which is rederived by a perturbation expansion method,the analytical expression of the second order spectrum is strictly deduced, and for infinite depth, the correct form of the kernel function is given for each octant do main.In fact,the present study improves and corrects the generally accepted results obtained by Tick
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
文摘Based on the second order random wave solutions of water wave equations in finite water depth, statistical distributions of the depth integrated local horizontal momentum components are derived by use of the characteristic function expansion method. The parameters involved in the distributions can be all determined by the water depth and the wave number spectrum of ocean waves. As an illustrative example, a fully developed wind generated sea is considered and the parameters are calculated for typical wind speeds and water depths by means of the Donelan and Pierson spectrum. The effects of nonlinearity and water depth on the distributions are also investigated.
文摘Linear wave theory and Longuet-Higgins and Steward’s (1964) group-induced second-order longweve (GSLW) theory ware used in this study on the grouping effect on wave forces acting on a verticalbreakwater. The calculated variance of total wave pressure on the vertical breakwater was closer tothe measured value if the wave grouping effect was considered.
基金supported by the Fundamental Research Developing Association for Shipbuilding and Offshore (REDAS), Japan
文摘A 3-D time-domain seakeeping analysis tool has been newly developed by using a higher-order boundary element method with the Rankine source as the kernel function. An iterative time-marching scheme for updating both kinematic and dynamic free-surface boundary conditions is adopted for achieving numerical accuracy and stability. A rectangular computational domain moving with the mean speed of ship is introduced. A damping beach at the outer portion of the truncated free surface is installed for satisfying the radiation condition. After numerical convergence checked, the diffraction unsteady problem of a Wigley hull traveling with a constant forward speed in waves is studied. Extensive results including wave exciting forces, wave patterns and pressure distributions on the hull are presented to validate the efficiency and accuracy of the proposed 3-D time-domain iterative Rankine BEM approach. Computed results are compared to be in good agreement with the corresponding experimental data and other published numerical solutions.