-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.展开更多
Seismic data processing typically deals with seismic wave reflections and neglects wave diffraction that affect the resolution. As a general rule, wave diffractions are treated as noise in seismic data processing. How...Seismic data processing typically deals with seismic wave reflections and neglects wave diffraction that affect the resolution. As a general rule, wave diffractions are treated as noise in seismic data processing. However, wave diffractions generally originate from geological structures, such as fractures, karst caves, and faults. The wave diffraction energy is much weaker than that of the reflections. Therefore, even if wave diffractions can be traced back to their origin, their energy is masked by that of the reflections. Separating and imaging diffractions and reflections can improve the imaging accuracy of diffractive targets. Based on the geometrical differences between reflections and diffractions on the plane-wave record; that is, reflections are quasi-linear and diffractions are quasi-hyperbolic, we use plane-wave prediction fltering to separate the wave diffractions. First, we estimate the local slope of the seismic event using plane- wave destruction filtering and, then, we predict and extract the wave reflections based on the local slope. Thus, we obtain the diffracted wavefield by directly subtracting the reflected wavefield from the entire wavefield. Finally, we image the diffracted wavefield and obtain high-resolution diffractive target results. 2D SEG salt model data suggest that the plane-wave prediction filtering eliminates the phase reversal in the plane-wave destruction filtering and maintains the original wavefield phase, improving the accuracy of imaging heterogeneous objects.展开更多
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.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green's functions of compressive and shear wave sources are deriv...This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green's functions of compressive and shear wave sources are derived based on Biot's theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by: (1) checking the satisfaction extent of the boundary conditions, (2) comparing the degenerated solutions of a single-phased case with well- known solutions, and (3) examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.展开更多
This paper presents a closed-form solution for diffraction of plane SH waves by a semi-circular cavity in half-space by using wave function expansion method. Accuracy of the solution is checked by the displacement res...This paper presents a closed-form solution for diffraction of plane SH waves by a semi-circular cavity in half-space by using wave function expansion method. Accuracy of the solution is checked by the displacement residual and stress residual along the boundaries. Numerical results show that there are notable differences for response amplitudes between a semi-circular cavity and a whole-circular cavity in a half-space.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional a...This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.展开更多
Diffraction of a two-dimensional (2D) semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called "impro...Diffraction of a two-dimensional (2D) semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called "improved cosine half- range expansion" algorithm exhibits an excellent performance in reducing displacement residual errors at two rim points of concern. The governing equations are developed in a manner that minimizes the residues of the boundary conditions. Detailed derivation and analysis procedures as well as truncation of infinite linear governing equations are presented. The semi-circular cavity model presented in this paper, due to its simple profile, is expected to be used in seismic wave propagation studies as a benchmark for examining the accuracies of various analytical or numerical methods for mixed-boundary wave propagation problems.展开更多
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.展开更多
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.展开更多
This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffract...This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated (either drained or undrained) cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.展开更多
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.展开更多
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.展开更多
A numerical model based on the mild-slope equation of water wave propagation over complicated bathymetry,taking into account the combined effects of refraction,diffraction and dissipation due to wavebreaking is presen...A numerical model based on the mild-slope equation of water wave propagation over complicated bathymetry,taking into account the combined effects of refraction,diffraction and dissipation due to wavebreaking is presented.Wave breaking is simulated by modifying the wave height probability density func-tion and the wave energy dissipation mechanism is parameterized according to that of the hydraulic jumpformulation.Solutions of the wave height,phase function,and the wave direction at every grid point areobtained by finite difference approximation of the governing equations,using Gauss-Seidel Iterative Method(GSIM)row by row.Its computational convenience allows it to be applied to large coast regions tostudy the wave transformation problem.Several case studies have been made and the results compare verywell with the experiment data and other model solutions.The capability and utility of the model forreal coast areas are illustrated by application to a shallow bay of northeast Australia.展开更多
In this paper, we investigate the diffraction tomography for quantitative imaging damages of partly through-thickness holes with various shapes in isotropic plates by using converted and non-converted scattered Lamb w...In this paper, we investigate the diffraction tomography for quantitative imaging damages of partly through-thickness holes with various shapes in isotropic plates by using converted and non-converted scattered Lamb waves generated nu- merically. Finite element simulations are carried out to provide the scattered wave data. The validity of the finite element model is confirmed by the comparison of scattering directivity pattern (SDP) of circle blind hole damage between the finite element simulations and the analytical results. The imaging method is based on a theoretical relation between the one-dimensional (1D) Fourier transform of the scattered projection and two-dimensional (2D) spatial Fourier transform of the scattering object. A quantitative image of the damage is obtained by carrying out the 2D inverse Fourier transform of the scattering object. The proposed approach employs a circle transducer network containing forward and backward projections, which lead to so-called transmission mode (TMDT) and reflection mode diffraction tomography (RMDT), respectively. The reconstructed results of the two projections for a non-converted SO scattered mode are investigated to illuminate the influence of the scattering field data. The results show that Lamb wave diffraction tomography using the combination of TMDT and RMDT improves the imaging effect compared with by using only the TMDT or RMDT. The scattered data of the converted A0 mode are also used to assess the performance of the diffraction tomography method. It is found that the circle and elliptical shaped damages can still be reasonably identified from the reconstructed images while the reconstructed results of other complex shaped damages like crisscross rectangles and racecourse are relatively poor.展开更多
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.展开更多
The noncollinear interaction of guided optical waves with magnetostatic waves under inclined bias magnetic field is theoretically studied in detail. Similar approach can also be applied to the collinear interaction. C...The noncollinear interaction of guided optical waves with magnetostatic waves under inclined bias magnetic field is theoretically studied in detail. Similar approach can also be applied to the collinear interaction. Calculation results indicate that the diffraction efficiency (DE) in magnitude is equal to the mode-conversion efficiency (MCE) under vertical bias magnetic field, but they differ greatly under inclined bias magnetic field. By comparison to the case of vertical magnetization, the DE or the MCE can be greatly increased under inclined magnetic field. The characteristic of the DE curves obtained is basically in agreement with the experimental result.展开更多
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.展开更多
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.展开更多
Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line sourc...Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.展开更多
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.展开更多
文摘-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.
基金funded jointly by the National Natural Science Foundation of China(No.41104069)the National Key Basic Research Program of China(973 Program:2011CB202402)+1 种基金the Shandong University Science and Technology Planning Project(No.J17KA197)the College of Petroleum Engineering in Shengli College China University of Petroleum"Chunhui Project"(No.KY2015003)
文摘Seismic data processing typically deals with seismic wave reflections and neglects wave diffraction that affect the resolution. As a general rule, wave diffractions are treated as noise in seismic data processing. However, wave diffractions generally originate from geological structures, such as fractures, karst caves, and faults. The wave diffraction energy is much weaker than that of the reflections. Therefore, even if wave diffractions can be traced back to their origin, their energy is masked by that of the reflections. Separating and imaging diffractions and reflections can improve the imaging accuracy of diffractive targets. Based on the geometrical differences between reflections and diffractions on the plane-wave record; that is, reflections are quasi-linear and diffractions are quasi-hyperbolic, we use plane-wave prediction fltering to separate the wave diffractions. First, we estimate the local slope of the seismic event using plane- wave destruction filtering and, then, we predict and extract the wave reflections based on the local slope. Thus, we obtain the diffracted wavefield by directly subtracting the reflected wavefield from the entire wavefield. Finally, we image the diffracted wavefield and obtain high-resolution diffractive target results. 2D SEG salt model data suggest that the plane-wave prediction filtering eliminates the phase reversal in the plane-wave destruction filtering and maintains the original wavefield phase, improving the accuracy of imaging heterogeneous objects.
文摘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.
基金Program for New Century Excellent Talents in University Under Grant No. NCET-05-0248the Key Program for Applied Basic Research of Tianjin Municipality Under Grant No. 07JCZDJC10100
文摘This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green's functions of compressive and shear wave sources are derived based on Biot's theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by: (1) checking the satisfaction extent of the boundary conditions, (2) comparing the degenerated solutions of a single-phased case with well- known solutions, and (3) examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.
基金supported by National Natural Science Foundation of China (No. 50978183)Tianjin Natural Science Foundation (No. 07JCZDJC10100)
文摘This paper presents a closed-form solution for diffraction of plane SH waves by a semi-circular cavity in half-space by using wave function expansion method. Accuracy of the solution is checked by the displacement residual and stress residual along the boundaries. Numerical results show that there are notable differences for response amplitudes between a semi-circular cavity and a whole-circular cavity in a half-space.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.
文摘Diffraction of a two-dimensional (2D) semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called "improved cosine half- range expansion" algorithm exhibits an excellent performance in reducing displacement residual errors at two rim points of concern. The governing equations are developed in a manner that minimizes the residues of the boundary conditions. Detailed derivation and analysis procedures as well as truncation of infinite linear governing equations are presented. The semi-circular cavity model presented in this paper, due to its simple profile, is expected to be used in seismic wave propagation studies as a benchmark for examining the accuracies of various analytical or numerical methods for mixed-boundary wave propagation problems.
基金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.
基金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.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated (either drained or undrained) cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.
基金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.
基金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.
文摘A numerical model based on the mild-slope equation of water wave propagation over complicated bathymetry,taking into account the combined effects of refraction,diffraction and dissipation due to wavebreaking is presented.Wave breaking is simulated by modifying the wave height probability density func-tion and the wave energy dissipation mechanism is parameterized according to that of the hydraulic jumpformulation.Solutions of the wave height,phase function,and the wave direction at every grid point areobtained by finite difference approximation of the governing equations,using Gauss-Seidel Iterative Method(GSIM)row by row.Its computational convenience allows it to be applied to large coast regions tostudy the wave transformation problem.Several case studies have been made and the results compare verywell with the experiment data and other model solutions.The capability and utility of the model forreal coast areas are illustrated by application to a shallow bay of northeast Australia.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11474195,11274226,11674214,and 51478258)
文摘In this paper, we investigate the diffraction tomography for quantitative imaging damages of partly through-thickness holes with various shapes in isotropic plates by using converted and non-converted scattered Lamb waves generated nu- merically. Finite element simulations are carried out to provide the scattered wave data. The validity of the finite element model is confirmed by the comparison of scattering directivity pattern (SDP) of circle blind hole damage between the finite element simulations and the analytical results. The imaging method is based on a theoretical relation between the one-dimensional (1D) Fourier transform of the scattered projection and two-dimensional (2D) spatial Fourier transform of the scattering object. A quantitative image of the damage is obtained by carrying out the 2D inverse Fourier transform of the scattering object. The proposed approach employs a circle transducer network containing forward and backward projections, which lead to so-called transmission mode (TMDT) and reflection mode diffraction tomography (RMDT), respectively. The reconstructed results of the two projections for a non-converted SO scattered mode are investigated to illuminate the influence of the scattering field data. The results show that Lamb wave diffraction tomography using the combination of TMDT and RMDT improves the imaging effect compared with by using only the TMDT or RMDT. The scattered data of the converted A0 mode are also used to assess the performance of the diffraction tomography method. It is found that the circle and elliptical shaped damages can still be reasonably identified from the reconstructed images while the reconstructed results of other complex shaped damages like crisscross rectangles and racecourse are relatively poor.
基金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.
文摘The noncollinear interaction of guided optical waves with magnetostatic waves under inclined bias magnetic field is theoretically studied in detail. Similar approach can also be applied to the collinear interaction. Calculation results indicate that the diffraction efficiency (DE) in magnitude is equal to the mode-conversion efficiency (MCE) under vertical bias magnetic field, but they differ greatly under inclined bias magnetic field. By comparison to the case of vertical magnetization, the DE or the MCE can be greatly increased under inclined magnetic field. The characteristic of the DE curves obtained is basically in agreement with the experimental result.
基金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.
基金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.
基金supported by National Natural Science Foundation of China (50978183)
文摘Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.
基金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.