The reflection and diffraction of a planar shock wave around a circular cylinder are a typical problem of the complex nonlinear shock wave phenomena in literature.It has long been studied experimentally,analytically a...The reflection and diffraction of a planar shock wave around a circular cylinder are a typical problem of the complex nonlinear shock wave phenomena in literature.It has long been studied experimentally,analytically as well as numerically.Takayama in 1987 obtained clear experimental pictures ofisopycnics in shock tube under the condi- tion that the impinging shock wave propagates as far as 3 diameters away from the cylinder.To know more complete- ly the whole unsteady process,it is desirable to get experimental results in a region which is more than 10 diameters away from the cylinder.This is what has been done in this paper by using the pulsed laser holographic interferometry for several shock Mach numbers of the impinging shock. Results for several moments are shown,giving more know- ledge about the whole unsteady flow field.This is useful for a reliable and complete understanding of the changing force acting on the cylinder,and provides interesting data to check the performance of many recently developed high resolution numerical methods for unsteady shock wave calculation.展开更多
The layout forms of several breakwater structures can be generalized as asymmetrical arrangements in actual engineering.However,the problem of wave diffraction around asymmetrically arranged breakwaters has not been a...The layout forms of several breakwater structures can be generalized as asymmetrical arrangements in actual engineering.However,the problem of wave diffraction around asymmetrically arranged breakwaters has not been adequately investigated.In this study,we propose an analytical method of wave diffraction for regular waves passing through asymmetrically arranged breakwaters,and we use the Nyström method to obtain the analytical solution numerically.We compared the results of this method with those of previous analytical solutions and with numerical results to demonstrate the validity of our approach.We also provided diffraction coefficient diagrams of breakwaters with different layout forms.Moreover,we described the analytical expression for the problem of diffraction through long-wave incident breakwaters and presented an analysis of the relationship between the diffraction coefficients and the widths of breakwater gates.The analytical method presented in this study contributes to the limited literature on the theory of wave diffraction through asymmetrically arranged breakwaters.展开更多
An analytical method is developed to study the sheltering effects on arc-shaped floating perforated breakwaters. In the process of analysis, the tloating breakwater is assumed to be rigid, thin, vertical, and immovabl...An analytical method is developed to study the sheltering effects on arc-shaped floating perforated breakwaters. In the process of analysis, the tloating breakwater is assumed to be rigid, thin, vertical, and immovable and located in water with constant depth. The fluid domain is divided into two regions by imaginary interface. The velocity potential in each region is expanded by eigenfunction in the context of linear theory. By satisfying continuity of pressure and normal velocity across the imaginary fluid interface, a set of linear algebraic equations can be obtained to determine the unknown coefficients for eigenfunction expansions. The accuracy of the present model was verified by a comparison with existing results for the case of arc-shaped floating breakwater. Numerical results, in the form of contour maps of the non-dimensional wave amplitude around the breakwater and diffracted wave amplitude at typical sections, are presented for a range of wave and breakwater parameters. Results show that the sheltering effects on the arc-shaped floating perforated breakwater are closely related to the incident wavelength, the draft and the porosity of the breakwater.展开更多
Diffracted seismic waves may be used to help identify and track geologically heterogeneous bodies or zones.However,the energy of diffracted waves is weaker than that of reflections.Therefore,the extraction of diffract...Diffracted seismic waves may be used to help identify and track geologically heterogeneous bodies or zones.However,the energy of diffracted waves is weaker than that of reflections.Therefore,the extraction of diffracted waves is the basis for the effective utilization of diffracted waves.Based on the difference in travel times between diffracted and reflected waves,we developed a method for separating the diffracted waves via singular value decomposition filters and presented an effective processing flowchart for diffracted wave separation and imaging.The research results show that the horizontally coherent difference between the reflected and diffracted waves can be further improved using normal move-out(NMO) correction.Then,a band-rank or high-rank approximation is used to suppress the reflected waves with better transverse coherence.Following,separation of reflected and diffracted waves is achieved after the filtered data are transformed into the original data domain by inverse NMO.Synthetic and field examples show that our proposed method has the advantages of fewer constraints,fast processing speed and complete extraction of diffracted waves.And the diffracted wave imaging results can effectively improve the identification accuracy of geological heterogeneous bodies or zones.展开更多
New version of SWAN model includes the wave diffraction effect which is the main improvement compared with the previous versions. Experimental data collected in the wave basin of the University of Delaware were used t...New version of SWAN model includes the wave diffraction effect which is the main improvement compared with the previous versions. Experimental data collected in the wave basin of the University of Delaware were used to test its performance. Wave heights were compared in the four cases （with different wave energies and directional spreading spectra）. The results agreed well with the measurements, especially for the broad directional spectra cases. The effect of wave diffraction was analyzed by switching on/off the corresponding tenn. By introducing the diffraction term, the distributions of wave height and wave direction were smoothed, especially obvious for the narrow spectrum cases. Compared with the calculations without diffraction, the model with diffraction effect gave better results.展开更多
The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interf...The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.展开更多
An analytical method is developed to study wave diffraction on arc-shaped and bottom-mounted perforated breakwaters. The breakwater is assumed to be rigid, thin, vertical, immovable and located in water of constant de...An analytical method is developed to study wave diffraction on arc-shaped and bottom-mounted perforated breakwaters. The breakwater is assumed to be rigid, thin, vertical, immovable and located in water of constant depth. The fluid domain is divided into two regions by imaginary interface. The velocity potential in each region is expanded by eigenfunctions. By satisfying the continuity of pressure and normal velocity across the imaginary fluid interface, a set of linear algebraic equations can be obtained to determine the unknown coefficients of eigenfunctions. Numerical results, in the form of contour maps of the relative wave amplitude around the breakwater, are presented for a range of wave and breakwater parameters. Results show that the wave diffraction on the arc-shaped and bottom-mounted perforated breakwater is related to the incident wavelength and the porosity of the breakwater. The porosity of the perforated breakwater may have great effect on the diffracted field.展开更多
A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water ...A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.展开更多
Two different methods for incorporating diffraction effect into wave action balance equation based coastal spectral wave models, WABED and SWAN, are discussed and evaluated with respect to their formulations, numerica...Two different methods for incorporating diffraction effect into wave action balance equation based coastal spectral wave models, WABED and SWAN, are discussed and evaluated with respect to their formulations, numerical implementations, and modeling capabilities. Both models were nm to simulate the wave transformation through a gap between two infinitely long breakwaters and that across an elliptical shoal observed in laboratory studies, with the emphasis laid on the diffraction induced by either obstacles or wave amplitude variations. Calculations of WABED were compared with Sommerfeld＇s analytical solutions, experimental observations and SWAN simulations. It is shown that both methods can predict reasonably wave diffraction for the two eases studied herein, and a fairly better performance is provided by WABED for stronger diffraction ease.展开更多
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.展开更多
Abstract： The scattering of plane SH-waves by topographic features in a layered transversely isotropic （TI） half-space is investigated by using an indirect boundary element method （IBEM）. Firstly, the anti-plane ...Abstract： The scattering of plane SH-waves by topographic features in a layered transversely isotropic （TI） half-space is investigated by using an indirect boundary element method （IBEM）. Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green＇s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green＇s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.展开更多
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.展开更多
The radiation and diffraction problem of a two-dimensional rectangular body with an opening floating on a semi- infinite fluid domain of finite water depth is analysed based on the linearized velocity potential theory...The radiation and diffraction problem of a two-dimensional rectangular body with an opening floating on a semi- infinite fluid domain of finite water depth is analysed based on the linearized velocity potential theory through an analytical solution procedure. The expressions for potentials are obtained by the method of variation separation, in which the unknown coefficients are determined by the boundary condition and matching requirement on the interface. The effects of the position of the hole and the gap between the body and side wall on hydrodynamic characteristics are investigated. Some resonance is observed like piston motion in a moon pool and sloshing in a closed tank because of the existence of restricted fluid domains.展开更多
Compliant offshore structures such as spars,tension leg platforms(TLPs) and semi-submersibles have been dramatically improved in recent years due to their capability for deep water operation. Waves are the most import...Compliant offshore structures such as spars,tension leg platforms(TLPs) and semi-submersibles have been dramatically improved in recent years due to their capability for deep water operation. Waves are the most important environmental phenomenon affecting these offshore structures. Estimation of wave forces is vital in offshore structure design. For large compliant offshore platforms,Morrison's equation is not valid anymore and usually diffraction theory is used. In this research,by using the finite difference method,a detailed analysis of the first-order diffraction of monochromatic waves on a large cylinder as a structural element is performed to solve the radiation and diffraction potentials. The results showed that the developed model is a reliable tool to estimate the wave forces and hydrodynamic coefficients on large structure elements when wave diffraction and radiation are considered.展开更多
A higher-order boundary element method(HOBEM) for simulating the fully nonlinear regular wave propagation and diffraction around a fixed vertical circular cylinder is investigated. The domain decomposition method with...A higher-order boundary element method(HOBEM) for simulating the fully nonlinear regular wave propagation and diffraction around a fixed vertical circular cylinder is investigated. The domain decomposition method with continuity conditions enforced on the interfaces between the adjacent sub-domains is implemented for reducing the computational cost. By adjusting the algorithm of iterative procedure on the interfaces, four types of coupling strategies are established, that is, Dirchlet/Dirchlet-Neumman/Neumman(D/D-N/N), Dirchlet-Neumman(D-N),Neumman-Dirchlet(N-D) and Mixed Dirchlet-Neumman/Neumman-Dirchlet(Mixed D-N/N-D). Numerical simulations indicate that the domain decomposition methods can provide accurate results compared with that of the single domain method. According to the comparisons of computational efficiency, the D/D-N/N coupling strategy is recommended for the wave propagation problem. As for the wave-body interaction problem, the Mixed D-N/N-D coupling strategy can obtain the highest computational efficiency.展开更多
The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field.The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstru...The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field.The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstruction,which in turn is expected to reduce the shock wave strength at the target location.In the present study the interaction of a plane shock front(generated from a shock tube)with various geometric designs such as,1)zig-zag geometric passage,2)staggered cylindrical obstructions and 3)zigzag passage with cylindrical obstructions have been investigated using computational technique.It is seen from the numerical simulation that,among the various designs,the maximum shock attenuation is produced by the zig-zag passage with cylindrical obstructions which is then followed by zig-zag passage and staggered cylindrical obstructions.A comprehensive investigation on the shock wave reflection and diffraction phenomena happening in the proposed complex passages have also been carried out.In the new zig-zag design,the initial shock wave undergoes shock wave reflection and diffraction process which swaps alternatively as the shock front moves from one turn to the other turn.This cyclic shock reflection and diffraction process helps in diffusing the shock wave energy with practically no obstruction to the flow field.It is found that by combining the shock attenuation ability of zig-zag passage(using shock reflection and diffraction)with the shock attenuation ability of cylindrical blocks(by flow obstruction),a drastic attenuation in shock strength can be achieved with moderate level of flow blocking.展开更多
This paper extends a prediction model for multi-directional random wave transformation based on an energy balance equation by Mase with the consideration of wave shoaling, refraction, diffraction, reflection and break...This paper extends a prediction model for multi-directional random wave transformation based on an energy balance equation by Mase with the consideration of wave shoaling, refraction, diffraction, reflection and breaking. This numerical model is improved by 1) introducing Wen's frequency spectrum and Mitsuyasu's directional function, which are more suitable to the coastal area of China; 2) considering energy dissipation caused by bottom friction, which ensures more accurate results for large-scale and shallow water areas; 3) taking into account a non-linear dispersion relation. Predictions using the extended wave model are carried out to study the feasibility of constructing the Ai Hua yacht port in Qingdao, China, with a comparison between two port layouts in design. Wave fields inside the port for different incident wave directions, water levels and return periods are simulated, and then two kinds of parameters are calculated to evaluate the wave conditions for the two layouts. Analyses show that Layout I is better than Layout II. Calculation results also show that the harbor will be calm for different wave directions under the design water level. On the contrary, the wave conditions do not wholly meet the requirements of a yacht port for ship berthing under the extreme water level. For safety consideration, the elevation of the breakwater might need to be properly increased to prevent wave overtopping under such water level. The extended numerical simulation model may provide an effective approach to computing wave heights in a harbor.展开更多
By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem...By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of an arbitrary cavity is obtained. By employing the orthogonal function expansion technique, the dynamic stress problem can be reduced to the solution of an infinite algebraic equation series. As an example, the numerical results for the dynamic stress concentration factor in thick plates with a circular, elliptic cavity are graphically presented. The numerical results are discussed.展开更多
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.展开更多
An exact analytic solution for wave diffraction by wedge or corner with arbitrary angle (rπ) and reflection coefficients (R0 and Rr) is presented in this paper. It is expressed in two forms-series and integral repres...An exact analytic solution for wave diffraction by wedge or corner with arbitrary angle (rπ) and reflection coefficients (R0 and Rr) is presented in this paper. It is expressed in two forms-series and integral representations, corresponding recurrence relation and asymptotic expressions are also derived. The solution is simplified for some special cases of rπ. For Rr= R0,r= 1/N and Rr≠R0,r = 1/2N, the solution can be reduced to linear superpositions of incident and partially reflected waves, hence a nonlinear solution of forth order for this problem can be obtained by using the author's theory of nonlinear interaction among gravity surface waves. The given solution is related to inhomogeneous Robin boundary conditions, which include the Neumann boundary conditions usually accepted in wave diffraction theory.展开更多
基金The project suported partially by National Natural Science Foundation of China
文摘The reflection and diffraction of a planar shock wave around a circular cylinder are a typical problem of the complex nonlinear shock wave phenomena in literature.It has long been studied experimentally,analytically as well as numerically.Takayama in 1987 obtained clear experimental pictures ofisopycnics in shock tube under the condi- tion that the impinging shock wave propagates as far as 3 diameters away from the cylinder.To know more complete- ly the whole unsteady process,it is desirable to get experimental results in a region which is more than 10 diameters away from the cylinder.This is what has been done in this paper by using the pulsed laser holographic interferometry for several shock Mach numbers of the impinging shock. Results for several moments are shown,giving more know- ledge about the whole unsteady flow field.This is useful for a reliable and complete understanding of the changing force acting on the cylinder,and provides interesting data to check the performance of many recently developed high resolution numerical methods for unsteady shock wave calculation.
基金supported by the National Natural Science Foundation of China(Grant No.51679132)the Science and Technology Commission of Shanghai Municipality(Grant No.21ZR1427000)Shanghai Frontiers Science Center of“Full Penetration”Far-Reaching Offshore Ocean Energy and Power.
文摘The layout forms of several breakwater structures can be generalized as asymmetrical arrangements in actual engineering.However,the problem of wave diffraction around asymmetrically arranged breakwaters has not been adequately investigated.In this study,we propose an analytical method of wave diffraction for regular waves passing through asymmetrically arranged breakwaters,and we use the Nyström method to obtain the analytical solution numerically.We compared the results of this method with those of previous analytical solutions and with numerical results to demonstrate the validity of our approach.We also provided diffraction coefficient diagrams of breakwaters with different layout forms.Moreover,we described the analytical expression for the problem of diffraction through long-wave incident breakwaters and presented an analysis of the relationship between the diffraction coefficients and the widths of breakwater gates.The analytical method presented in this study contributes to the limited literature on the theory of wave diffraction through asymmetrically arranged breakwaters.
基金supported by the Natural Science Foundation of Jiangsu Province(Grant No.Bk2006013)the foundation of the State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University(Grant No.0907)the foundation of Engineering Institute of Engineering Corps and PLA University of Science & Technology
文摘An analytical method is developed to study the sheltering effects on arc-shaped floating perforated breakwaters. In the process of analysis, the tloating breakwater is assumed to be rigid, thin, vertical, and immovable and located in water with constant depth. The fluid domain is divided into two regions by imaginary interface. The velocity potential in each region is expanded by eigenfunction in the context of linear theory. By satisfying continuity of pressure and normal velocity across the imaginary fluid interface, a set of linear algebraic equations can be obtained to determine the unknown coefficients for eigenfunction expansions. The accuracy of the present model was verified by a comparison with existing results for the case of arc-shaped floating breakwater. Numerical results, in the form of contour maps of the non-dimensional wave amplitude around the breakwater and diffracted wave amplitude at typical sections, are presented for a range of wave and breakwater parameters. Results show that the sheltering effects on the arc-shaped floating perforated breakwater are closely related to the incident wavelength, the draft and the porosity of the breakwater.
基金supported by the National Natural Science Foundation of China(41874123)Shaanxi Province Natural Science Basic Research Project(2017JZ007)PetroChina Innovation Foundation(2014D-5006-0303)。
文摘Diffracted seismic waves may be used to help identify and track geologically heterogeneous bodies or zones.However,the energy of diffracted waves is weaker than that of reflections.Therefore,the extraction of diffracted waves is the basis for the effective utilization of diffracted waves.Based on the difference in travel times between diffracted and reflected waves,we developed a method for separating the diffracted waves via singular value decomposition filters and presented an effective processing flowchart for diffracted wave separation and imaging.The research results show that the horizontally coherent difference between the reflected and diffracted waves can be further improved using normal move-out(NMO) correction.Then,a band-rank or high-rank approximation is used to suppress the reflected waves with better transverse coherence.Following,separation of reflected and diffracted waves is achieved after the filtered data are transformed into the original data domain by inverse NMO.Synthetic and field examples show that our proposed method has the advantages of fewer constraints,fast processing speed and complete extraction of diffracted waves.And the diffracted wave imaging results can effectively improve the identification accuracy of geological heterogeneous bodies or zones.
基金This study was supported by the National Key Basic Research Project of China (Grant No2002CB412403)the Research Project in Science and Technology Commission of Shanghai Municipality,China (Grant No04DZ12049)
文摘New version of SWAN model includes the wave diffraction effect which is the main improvement compared with the previous versions. Experimental data collected in the wave basin of the University of Delaware were used to test its performance. Wave heights were compared in the four cases （with different wave energies and directional spreading spectra）. The results agreed well with the measurements, especially for the broad directional spectra cases. The effect of wave diffraction was analyzed by switching on/off the corresponding tenn. By introducing the diffraction term, the distributions of wave height and wave direction were smoothed, especially obvious for the narrow spectrum cases. Compared with the calculations without diffraction, the model with diffraction effect gave better results.
基金sponsored by the National Basic Research Program of China(973 Program,Grant No.2014CB046203)the National Natural Science Foundation of China(Grant No.11072140)
文摘The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.
基金This project was supported by the Natural Science Foundation of Jiangsu Province (Grant NoBk2006013)
文摘An analytical method is developed to study wave diffraction on arc-shaped and bottom-mounted perforated breakwaters. The breakwater is assumed to be rigid, thin, vertical, immovable and located in water of constant depth. The fluid domain is divided into two regions by imaginary interface. The velocity potential in each region is expanded by eigenfunctions. By satisfying the continuity of pressure and normal velocity across the imaginary fluid interface, a set of linear algebraic equations can be obtained to determine the unknown coefficients of eigenfunctions. Numerical results, in the form of contour maps of the relative wave amplitude around the breakwater, are presented for a range of wave and breakwater parameters. Results show that the wave diffraction on the arc-shaped and bottom-mounted perforated breakwater is related to the incident wavelength and the porosity of the breakwater. The porosity of the perforated breakwater may have great effect on the diffracted field.
基金The present work was financially supported by the National Natural Science Foundation of China under contract No.50025924.
文摘A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.
基金supported by the Special Fund of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering(Grant No.2009585812)the National Natural Science Foundation of China(Grant No.50979033)the Programfor New Century Excellent Talentsin University of China(Grant No.NCET-07-0255)
文摘Two different methods for incorporating diffraction effect into wave action balance equation based coastal spectral wave models, WABED and SWAN, are discussed and evaluated with respect to their formulations, numerical implementations, and modeling capabilities. Both models were nm to simulate the wave transformation through a gap between two infinitely long breakwaters and that across an elliptical shoal observed in laboratory studies, with the emphasis laid on the diffraction induced by either obstacles or wave amplitude variations. Calculations of WABED were compared with Sommerfeld＇s analytical solutions, experimental observations and SWAN simulations. It is shown that both methods can predict reasonably wave diffraction for the two eases studied herein, and a fairly better performance is provided by WABED for stronger diffraction ease.
文摘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 under Grant Nos.51578373 and 51578372
文摘Abstract： The scattering of plane SH-waves by topographic features in a layered transversely isotropic （TI） half-space is investigated by using an indirect boundary element method （IBEM）. Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green＇s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green＇s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.
文摘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.
基金supported by the Lloyd's Register Educational Trust (The LRET) through the joint centre involving University College London, Shanghai Jiao Tong University and Harbin Engineering University
文摘The radiation and diffraction problem of a two-dimensional rectangular body with an opening floating on a semi- infinite fluid domain of finite water depth is analysed based on the linearized velocity potential theory through an analytical solution procedure. The expressions for potentials are obtained by the method of variation separation, in which the unknown coefficients are determined by the boundary condition and matching requirement on the interface. The effects of the position of the hole and the gap between the body and side wall on hydrodynamic characteristics are investigated. Some resonance is observed like piston motion in a moon pool and sloshing in a closed tank because of the existence of restricted fluid domains.
文摘Compliant offshore structures such as spars,tension leg platforms(TLPs) and semi-submersibles have been dramatically improved in recent years due to their capability for deep water operation. Waves are the most important environmental phenomenon affecting these offshore structures. Estimation of wave forces is vital in offshore structure design. For large compliant offshore platforms,Morrison's equation is not valid anymore and usually diffraction theory is used. In this research,by using the finite difference method,a detailed analysis of the first-order diffraction of monochromatic waves on a large cylinder as a structural element is performed to solve the radiation and diffraction potentials. The results showed that the developed model is a reliable tool to estimate the wave forces and hydrodynamic coefficients on large structure elements when wave diffraction and radiation are considered.
基金supported by the National Natural Science Foundation of China(Grant No.51490673)the Pre-Research Field Fund Project of the Central Military Commission of China(Grant No.61402070201)the Fundamental Research Funds for the Central Universities(Grant No.DUT18LK09)
文摘A higher-order boundary element method(HOBEM) for simulating the fully nonlinear regular wave propagation and diffraction around a fixed vertical circular cylinder is investigated. The domain decomposition method with continuity conditions enforced on the interfaces between the adjacent sub-domains is implemented for reducing the computational cost. By adjusting the algorithm of iterative procedure on the interfaces, four types of coupling strategies are established, that is, Dirchlet/Dirchlet-Neumman/Neumman(D/D-N/N), Dirchlet-Neumman(D-N),Neumman-Dirchlet(N-D) and Mixed Dirchlet-Neumman/Neumman-Dirchlet(Mixed D-N/N-D). Numerical simulations indicate that the domain decomposition methods can provide accurate results compared with that of the single domain method. According to the comparisons of computational efficiency, the D/D-N/N coupling strategy is recommended for the wave propagation problem. As for the wave-body interaction problem, the Mixed D-N/N-D coupling strategy can obtain the highest computational efficiency.
文摘The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field.The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstruction,which in turn is expected to reduce the shock wave strength at the target location.In the present study the interaction of a plane shock front(generated from a shock tube)with various geometric designs such as,1)zig-zag geometric passage,2)staggered cylindrical obstructions and 3)zigzag passage with cylindrical obstructions have been investigated using computational technique.It is seen from the numerical simulation that,among the various designs,the maximum shock attenuation is produced by the zig-zag passage with cylindrical obstructions which is then followed by zig-zag passage and staggered cylindrical obstructions.A comprehensive investigation on the shock wave reflection and diffraction phenomena happening in the proposed complex passages have also been carried out.In the new zig-zag design,the initial shock wave undergoes shock wave reflection and diffraction process which swaps alternatively as the shock front moves from one turn to the other turn.This cyclic shock reflection and diffraction process helps in diffusing the shock wave energy with practically no obstruction to the flow field.It is found that by combining the shock attenuation ability of zig-zag passage(using shock reflection and diffraction)with the shock attenuation ability of cylindrical blocks(by flow obstruction),a drastic attenuation in shock strength can be achieved with moderate level of flow blocking.
基金supported by the National Natural Science Foundation of China (50879085)the Program for New Century Excellent Talents in University(NCET-07-0778)Fundamental Research Funds for the Central Universities (2012QNA4020)
文摘This paper extends a prediction model for multi-directional random wave transformation based on an energy balance equation by Mase with the consideration of wave shoaling, refraction, diffraction, reflection and breaking. This numerical model is improved by 1) introducing Wen's frequency spectrum and Mitsuyasu's directional function, which are more suitable to the coastal area of China; 2) considering energy dissipation caused by bottom friction, which ensures more accurate results for large-scale and shallow water areas; 3) taking into account a non-linear dispersion relation. Predictions using the extended wave model are carried out to study the feasibility of constructing the Ai Hua yacht port in Qingdao, China, with a comparison between two port layouts in design. Wave fields inside the port for different incident wave directions, water levels and return periods are simulated, and then two kinds of parameters are calculated to evaluate the wave conditions for the two layouts. Analyses show that Layout I is better than Layout II. Calculation results also show that the harbor will be calm for different wave directions under the design water level. On the contrary, the wave conditions do not wholly meet the requirements of a yacht port for ship berthing under the extreme water level. For safety consideration, the elevation of the breakwater might need to be properly increased to prevent wave overtopping under such water level. The extended numerical simulation model may provide an effective approach to computing wave heights in a harbor.
文摘By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of an arbitrary cavity is obtained. By employing the orthogonal function expansion technique, the dynamic stress problem can be reduced to the solution of an infinite algebraic equation series. As an example, the numerical results for the dynamic stress concentration factor in thick plates with a circular, elliptic cavity are graphically presented. The numerical results are discussed.
基金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.
文摘An exact analytic solution for wave diffraction by wedge or corner with arbitrary angle (rπ) and reflection coefficients (R0 and Rr) is presented in this paper. It is expressed in two forms-series and integral representations, corresponding recurrence relation and asymptotic expressions are also derived. The solution is simplified for some special cases of rπ. For Rr= R0,r= 1/N and Rr≠R0,r = 1/2N, the solution can be reduced to linear superpositions of incident and partially reflected waves, hence a nonlinear solution of forth order for this problem can be obtained by using the author's theory of nonlinear interaction among gravity surface waves. The given solution is related to inhomogeneous Robin boundary conditions, which include the Neumann boundary conditions usually accepted in wave diffraction theory.
