Wave decomposition phenomenon and spectrum evolution over submerged bars are investigated by a previously developed numerical model. First, the computed free surface displacements of regular waves at various locations...Wave decomposition phenomenon and spectrum evolution over submerged bars are investigated by a previously developed numerical model. First, the computed free surface displacements of regular waves at various locations are compared with the available experimental data to confirm the validity of the numerical model, and satisfactory agreements are obtained. In addition, variations of decomposition characteristics with incident wave parameters and the change of energy spectrum for regular waves are also studied. Then the spectrum evolution of irregular waves over submerged bars, as well as the influence of incident peak wave period and the steepness of the front slope of the bar on spectrum evolution, is investigated. Wave decomposition and spectral shape are found to be significantly influenced by the incident wave conditions. When the upslope of the bar becomes 1：2, the length of the slope becomes shorter and will not benefit the generation of high frequency energy, so spectrum evolution is not significant.展开更多
Ocean boundaries present a significant effect on the vibroacoustic characteristics and sound propagation of an elastic structure in practice.In this study,an efficient finite element/wave superposition method(FE/WSM)f...Ocean boundaries present a significant effect on the vibroacoustic characteristics and sound propagation of an elastic structure in practice.In this study,an efficient finite element/wave superposition method(FE/WSM)for predicting the three-dimen-sional acoustic radiation from an arbitrary-shaped radiator in Pekeris waveguides with a lossy seabed is proposed.The method is based on the FE method(FEM),WSM,and sound propagation models.First,a near-field vibroacoustic model is established by the FEM to obtain vibration information on a radiator surface.Then,the WSM based on the Helmholtz boundary integral is used to pre-dict the far-field acoustic radiation and propagation.Furthermore,the rigorous image source method and complex normal mode are employed to obtain the near-and far-field Green’s function(GF),respectively.The former,which is based on the spherical wave decomposition,is adopted to accurately solve the near-field source strength,and the far-field acoustic radiation is calculated by the latter and perturbation theory.The simulations of both models are compared to theoretical wavenumber integration solutions.Finally,numerical experiments on elastic spherical and cylindrical shells in Pekeris waveguides are presented to validate the accuracy and efficiency of the proposed method.The results show that the FE/WSM is adaptable to complex radiators and ocean-acoustic envi-ronments,and are easy to implement and computationally efficient in calculating the structural vibration,acoustic radiation,and sound propagation of arbitrarily shaped radiators in practical ocean environments.展开更多
Impacts of the Indonesian Throughflow(ITF) on seasonal circulation in the equatorial eastern Indian Ocean are investigated using the ocean-only model LICOM by opening and closing ITF passages. LICOM had daily forcing ...Impacts of the Indonesian Throughflow(ITF) on seasonal circulation in the equatorial eastern Indian Ocean are investigated using the ocean-only model LICOM by opening and closing ITF passages. LICOM had daily forcing from NCEP reanalysis data during 2000–2011. It can reproduce vertical profiles of mean density and buoyancy frequency of World Ocean Atlas 2013 data. The model also simulates well annual oscillation in the central Indian Ocean and semiannual oscillation in the eastern Indian Ocean of sea level anomalies(SLA) using satellite altimeter data, as well as the semiannual oscillation of surface zonal equatorial currents of Ocean Surface Current Analyses Real Time current data in the equatorial Indian Ocean. The wave decomposition method is used to analyze the propagation and reflection of equatorial long waves based on LICOM output. Wave analysis suggests that ITF blockage mainly influences waves generated from the Indian Ocean but not the Pacific Ocean, and eastern boundary reflections play an important role in semiannual oscillations of SLA and zonal current dif ferences in the equatorial Indian Ocean associated with ITF. Reconstructed ITF-caused SLA using wave decomposition coefficient dif ferences between closed and open ITF-passage experiments suggest both Kelvin and Rossby waves from the first baroclinic mode have comparable contributions to the semiannual oscillations of SLA diff erence. However, reconstructed ITFcaused surface zonal currents at the equator suggest that the first meridional-mode Rossby wave has much greater contribution than the first baroclinic mode Kelvin wave. Both reconstructed sea level and zonal currents demonstrate that the first baroclinic mode has a greater contribution than other baroclinic modes.展开更多
To improve localization accuracy, the spherical microphone arrays are used to capture high-order wavefield in- formation. For the far field sound sources, the array signal model is constructed based on plane wave deco...To improve localization accuracy, the spherical microphone arrays are used to capture high-order wavefield in- formation. For the far field sound sources, the array signal model is constructed based on plane wave decomposition. The spatial spectrum function is calculated by minimum variance distortionless response （MVDR） to scan the three-dimensional space. The peak values of the spectrum function correspond to the directions of multiple sound sources. A diagonal loading method is adopted to solve the ill-conditioned cross spectrum matrix of the received signals. The loading level depends on the alleviation of the ill-condition of the matrix and the accuracy of the inverse calculation. Compared with plane wave decomposition method, our proposed localization algorithm can acquire high spatial resolution and better estimation for multiple sound source directions, especially in low signal to noise ratio （SNR）.展开更多
Explicit fomulas for 2-D electroelastic fundamental solutions in general anisotropic piezoelectric media subjected to a line force and a line charge are obtained by using the plane wave decomposition method and a subs...Explicit fomulas for 2-D electroelastic fundamental solutions in general anisotropic piezoelectric media subjected to a line force and a line charge are obtained by using the plane wave decomposition method and a subsequent application of the residue calculus. 'Anisotropic' means that any material symmetry restrictions are not assumed. 'Two dimensional' includes not only in-plane problems but also anti-plane problems and problems in which in-plane and anti-plane deformations couple each other. As a special case, the solutions for transversely isotropic piezoelectric media are given.展开更多
It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation c...It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation coupled with the first and the second order absorbing boundary conditions. The computational scheme is also developed. The approach allows the absorbing boundary conditions to be naturally imposed, which makes it easier for us to construct high order schemes for the absorbing boundary conditions. A thirdorder Lagrange finite element method with mass lumping is applied to obtain the spatial discretization of the wave equation. The resulting scheme is stable and is very efficient since no matrix inversion is needed at each time step. Moreover, we have shown both abstract and explicit conditional stability results for the fully-discrete schemes. The results are helpful for designing computational parameters in computations. Numerical computations are illustrated to show the efficiency and accuracy of our method. In particular, essentially no boundary reflection is seen at the artificial boundaries.展开更多
基金The National Natural Science Foundation of China for Distinguished Young Scholars under contract No. 50025925the Chinese Postdoctoral Science Foundation under contract No. 20040350088
文摘Wave decomposition phenomenon and spectrum evolution over submerged bars are investigated by a previously developed numerical model. First, the computed free surface displacements of regular waves at various locations are compared with the available experimental data to confirm the validity of the numerical model, and satisfactory agreements are obtained. In addition, variations of decomposition characteristics with incident wave parameters and the change of energy spectrum for regular waves are also studied. Then the spectrum evolution of irregular waves over submerged bars, as well as the influence of incident peak wave period and the steepness of the front slope of the bar on spectrum evolution, is investigated. Wave decomposition and spectral shape are found to be significantly influenced by the incident wave conditions. When the upslope of the bar becomes 1：2, the length of the slope becomes shorter and will not benefit the generation of high frequency energy, so spectrum evolution is not significant.
基金financially supported by the National Key Research and Development Plan of China (No. 2016YFC1401203)the National Natural Science Foundation of China (Nos. 42006168 and 11404079)
文摘Ocean boundaries present a significant effect on the vibroacoustic characteristics and sound propagation of an elastic structure in practice.In this study,an efficient finite element/wave superposition method(FE/WSM)for predicting the three-dimen-sional acoustic radiation from an arbitrary-shaped radiator in Pekeris waveguides with a lossy seabed is proposed.The method is based on the FE method(FEM),WSM,and sound propagation models.First,a near-field vibroacoustic model is established by the FEM to obtain vibration information on a radiator surface.Then,the WSM based on the Helmholtz boundary integral is used to pre-dict the far-field acoustic radiation and propagation.Furthermore,the rigorous image source method and complex normal mode are employed to obtain the near-and far-field Green’s function(GF),respectively.The former,which is based on the spherical wave decomposition,is adopted to accurately solve the near-field source strength,and the far-field acoustic radiation is calculated by the latter and perturbation theory.The simulations of both models are compared to theoretical wavenumber integration solutions.Finally,numerical experiments on elastic spherical and cylindrical shells in Pekeris waveguides are presented to validate the accuracy and efficiency of the proposed method.The results show that the FE/WSM is adaptable to complex radiators and ocean-acoustic envi-ronments,and are easy to implement and computationally efficient in calculating the structural vibration,acoustic radiation,and sound propagation of arbitrarily shaped radiators in practical ocean environments.
基金Supported by the National Natural Science Foundation of China(No.41206018)the National Natural Science Foundation of China(Nos.41176019,41421005,U1406401)+4 种基金the Chinese Academy of Sciences(No.XDA11010203)to WANG Jingthe Chinese Academy of Sciences(No.XDA11010301)the National Basic Research Program of China(973 Program)(No.2012CB956001)the Specialized Scientific Research Project for Public Welfare Industries(Meteorology)(No.GYHY201306018)the State Oceanic Administration of China(No.GASI-03-01-01-05)to YUAN Dongliang
文摘Impacts of the Indonesian Throughflow(ITF) on seasonal circulation in the equatorial eastern Indian Ocean are investigated using the ocean-only model LICOM by opening and closing ITF passages. LICOM had daily forcing from NCEP reanalysis data during 2000–2011. It can reproduce vertical profiles of mean density and buoyancy frequency of World Ocean Atlas 2013 data. The model also simulates well annual oscillation in the central Indian Ocean and semiannual oscillation in the eastern Indian Ocean of sea level anomalies(SLA) using satellite altimeter data, as well as the semiannual oscillation of surface zonal equatorial currents of Ocean Surface Current Analyses Real Time current data in the equatorial Indian Ocean. The wave decomposition method is used to analyze the propagation and reflection of equatorial long waves based on LICOM output. Wave analysis suggests that ITF blockage mainly influences waves generated from the Indian Ocean but not the Pacific Ocean, and eastern boundary reflections play an important role in semiannual oscillations of SLA and zonal current dif ferences in the equatorial Indian Ocean associated with ITF. Reconstructed ITF-caused SLA using wave decomposition coefficient dif ferences between closed and open ITF-passage experiments suggest both Kelvin and Rossby waves from the first baroclinic mode have comparable contributions to the semiannual oscillations of SLA diff erence. However, reconstructed ITFcaused surface zonal currents at the equator suggest that the first meridional-mode Rossby wave has much greater contribution than the first baroclinic mode Kelvin wave. Both reconstructed sea level and zonal currents demonstrate that the first baroclinic mode has a greater contribution than other baroclinic modes.
基金Project supported by the National Natural Science Foundation of China (Grant No.61001160)the Doctoral Foundation of Ministry of Education (Grant No.20093108120018)the Shanghai Leading Academic Discipline Project (Grant No.S30108)
文摘To improve localization accuracy, the spherical microphone arrays are used to capture high-order wavefield in- formation. For the far field sound sources, the array signal model is constructed based on plane wave decomposition. The spatial spectrum function is calculated by minimum variance distortionless response （MVDR） to scan the three-dimensional space. The peak values of the spectrum function correspond to the directions of multiple sound sources. A diagonal loading method is adopted to solve the ill-conditioned cross spectrum matrix of the received signals. The loading level depends on the alleviation of the ill-condition of the matrix and the accuracy of the inverse calculation. Compared with plane wave decomposition method, our proposed localization algorithm can acquire high spatial resolution and better estimation for multiple sound source directions, especially in low signal to noise ratio （SNR）.
文摘Explicit fomulas for 2-D electroelastic fundamental solutions in general anisotropic piezoelectric media subjected to a line force and a line charge are obtained by using the plane wave decomposition method and a subsequent application of the residue calculus. 'Anisotropic' means that any material symmetry restrictions are not assumed. 'Two dimensional' includes not only in-plane problems but also anti-plane problems and problems in which in-plane and anti-plane deformations couple each other. As a special case, the solutions for transversely isotropic piezoelectric media are given.
文摘It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation coupled with the first and the second order absorbing boundary conditions. The computational scheme is also developed. The approach allows the absorbing boundary conditions to be naturally imposed, which makes it easier for us to construct high order schemes for the absorbing boundary conditions. A thirdorder Lagrange finite element method with mass lumping is applied to obtain the spatial discretization of the wave equation. The resulting scheme is stable and is very efficient since no matrix inversion is needed at each time step. Moreover, we have shown both abstract and explicit conditional stability results for the fully-discrete schemes. The results are helpful for designing computational parameters in computations. Numerical computations are illustrated to show the efficiency and accuracy of our method. In particular, essentially no boundary reflection is seen at the artificial boundaries.
