The vibration suppression of the finite plate with square steel beams is studied using traveling wave method. The finite plate with square beams is modeled as the coupling systems between the plate flexural motion and...The vibration suppression of the finite plate with square steel beams is studied using traveling wave method. The finite plate with square beams is modeled as the coupling systems between the plate flexural motion and the flexural and torsional motions for the square beams. The vibration response at any position of the coupling structure can be obtained by wave method. Numerical results show that comparing to finite element method (FEM), not only the low frequency but also the medium-high frequency vibration response of the finite plate with square beam can be effectively calculated by wave method. The suppression effect can be increased as the square beam is located at one-third of the length of plate or increasing the height of the beam. The study provides reference for arranged square beams applying to vibration suppression of ship and train structures.展开更多
Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the...Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m under- ground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.展开更多
The plane wave numerical technique is recast from Ampere’s and Faraday’s laws for materials that are characterized with a bianisotropic form of the constitutive relations. The populating expressions are provided for...The plane wave numerical technique is recast from Ampere’s and Faraday’s laws for materials that are characterized with a bianisotropic form of the constitutive relations. The populating expressions are provided for the eigenvalue matrix system that can be directly solved for the angular frequencies and field profiles when bianisotropy is included. To demonstrate the computation process and expected state diagrams and field profiles, numerical computation examples are provided for a bianisotropic Bragg Array with central defect. It is shown that the location of the magnetoelectric tensor elements has a significant effect on the eigenstates of an equivalent isotropic (anisotropic) structure. One form of the magnetoelectric tensor (diagonal elements only) leads to the observation of merging states and the formation of exceptional points. The numerical approach presented can be implemented as an add-on to the familiar plane wave numerical technique.展开更多
The wave method is introduced to vibration analysis of the fluid-conveying carbon nanotube. The constitutive relation of carbon nanotube on micro-scale is founded using the non- local elastic theory. The governing equ...The wave method is introduced to vibration analysis of the fluid-conveying carbon nanotube. The constitutive relation of carbon nanotube on micro-scale is founded using the non- local elastic theory. The governing equation on micro-scale is obtained. And the first five orders of the natural frequency of the carbon nanotube conveying fluid with various speeds are calculated through the wave method. Besides, the critical flow velocity when the carbon nanotube loses stability is obtained. Meanwhile, a contrast is made between the result obtained through tile wave method and that in previous researches.展开更多
A low-frequency multi-mode ultrasonic Lamb wave method suitable for character- izing the thickness, the density and the elastic constants of the ultra-thin transversely isotropic laminate composite is presented. The &...A low-frequency multi-mode ultrasonic Lamb wave method suitable for character- izing the thickness, the density and the elastic constants of the ultra-thin transversely isotropic laminate composite is presented. The 'ultra-thin' here means that the thickness of the plate is much less than the wavelength of the ultrasonic wave so that the echoes from the front and back faces of the plate can't be separated in the time domain. The dispersion equations for the low frequency ultrasonic Lamb waves with the propagation directions parallel and vertical to the fiber direction are derived. In conjunction with the least square algorithm method, the secant algorithm is used to estimate the parameters of the ultra-thin fiber-reinforced composite layer. The evaluation errors and the sensitivity of the method to different paramters of the thin composite are analyzed. The technique has been used to characterize the ultra-thin grass fiber reinforced PES composite with thickness down to ten percents of the ultrasonic wavelength. It is observed that the agreement between the nominal and the estimation values is reasonably good.展开更多
The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6×6 Luttinger Kohn model.The effect of the number and period of pla...The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6×6 Luttinger Kohn model.The effect of the number and period of plane waves used for expansion on the stability of energy eigenvalues is examined.For practical calculation,it should choose the period large sufficiently to ensure the envelope functions vanish at the boundary and the number of plane waves large enough to ensure the energy eigenvalues keep unchanged within a prescribed range.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident...The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.展开更多
By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at t...By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at the eigenfrequencies of the corresponding interior problems is proposed. By adding appropriate damp to the virtual source system of the wave superposition method, the unique solutions for all wave numbers can be ensured. Based on this thought, a novel method-wave superposition method with complex radius vector is constructed. Not only is the computational time of this method approximately equal to that of the standard wave superposition method, but also the accuracy is much higher compared with other correlative methods. Finally, by taking the pulsating sphere and oscillating sphere as examples, the results of calculation show that the present method can effectively overcome the non-uniqueness problem.展开更多
In order to monitor the basic mechanical properties and interior damage of concrete structures,the piezoelectric actuator/sensor based wave propagation method was investigated experimentally in the laboratory using a ...In order to monitor the basic mechanical properties and interior damage of concrete structures,the piezoelectric actuator/sensor based wave propagation method was investigated experimentally in the laboratory using a specifically designed test setup.The energy attenuation of stress waves was measured by the relative index between the output voltage of sensors and the excitation voltage at the actuator.Based on the experimental results of concrete cube and cylinder specimens,the effect of excitation frequencies,excitation amplitude,wave propagation paths and the curing age on the output signals of sensors are evaluated.The results show that the relative voltage attenuation coefficient RVAC is an effective indicator for measuring the attenuation of stress waves through the interior of concrete.展开更多
Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orien...Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orientations, etc. of the three array arrangements, i. e., linear array, T-type array and pentagon array, on the resolution of the directional spreading of waves, are investigated experimentally. This study can be used as a reference in the experimental study and the field measurement of directional waves.展开更多
Wave energy resource is a very important ocean renewable energy. A reliable assessment of wave energy resources must be performed before they can be exploited. Compared with wave model, altimeter can provide more accu...Wave energy resource is a very important ocean renewable energy. A reliable assessment of wave energy resources must be performed before they can be exploited. Compared with wave model, altimeter can provide more accurate in situ observations for ocean wave which can be as a novel method for wave energy assessment.The advantage of altimeter data is to provide accurate significant wave height observations for wave. In order to develop characteristic and advantage of altimeter data and apply altimeter data to wave energy assessment, in this study, we established an assessing method for wave energy in local sea area which is dedicated to altimeter data.This method includes three parts including data selection and processing, establishment of evaluation indexes system and criterion of regional division. Then a case study of Northwest Pacific was performed to discuss specific application for this method. The results show that assessing method in this paper can assess reserves and temporal and spatial distribution effectively and provide scientific references for the siting of wave power plants and the design of wave energy convertors.展开更多
A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid (VOF) method. An active 'absorbing wave-maker based on linear wave theory is set on the left boun...A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid (VOF) method. An active 'absorbing wave-maker based on linear wave theory is set on the left boundary of the wave flume. The progressive waves and the absorbing waves are generated simultaneously at the active wave generating-absorbing boundary. The absorbing waves are generated to eliminate the waves coming back to the generating boundary due to reflection from the outflow boundary and the structures. SIRW method proposed by Frigaard and Brorsen (1995) is used to separate the incident waves and reflected waves. The digital filters are designed based on the surface elevation signals of the two wave gauges. The corrected velocity of the wave-maker paddle is the output from the digital filter in real time. The numerical results of regular and irregular waves by the active absorbing-generating boundary are compared with the numerical results by the ordinary generating boundary to verify the performance of the active absorbing-generator boundary. The differences between the initial incident waves and the estimated incident waves are analyzed.展开更多
The (G'/G, 1/G)-expansion method for finding exact travelling wave solutions of nonlinear evolution equations, which can be thought of as an extension of the (G'/G)-expansion method proposed recently, is present...The (G'/G, 1/G)-expansion method for finding exact travelling wave solutions of nonlinear evolution equations, which can be thought of as an extension of the (G'/G)-expansion method proposed recently, is presented. By using this method abundant travelling wave so- lutions with arbitrary parameters of the Zakharov equations are successfully obtained. When the parameters are replaced by special values, the well-known solitary wave solutions of the equations are rediscovered from the travelling waves.展开更多
Nozzle damping is one of the most important factors in the suppression of combustion instability in solid rocket motors.For an engineering solid rocket motor that experiences combustion instability at the end of burni...Nozzle damping is one of the most important factors in the suppression of combustion instability in solid rocket motors.For an engineering solid rocket motor that experiences combustion instability at the end of burning,a wave attenuation method is proposed to assess the nozzle damping characteristics numerically.In this method,a periodic pressure oscillation signal which frequency equals to the first acoustic mode is superimposed on a steady flow at the head end of the chamber.When the pressure oscillation is turned off,the decay rate of the pressure can be used to determine the nozzle attenuation constant.The damping characteristics of three other nozzle geometries are numerically studied with this method under the same operating condition.The results show that the convex nozzle provides more damping than the conical nozzle which in turn provides more damping than the concave nozzle.All the three nozzles have better damping effect than that of basic nozzle geometry.At last,the phase difference in the chamber is analyzed,and the numerical pressure distribution satisfies well with theoretical distribution.展开更多
The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, w...The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.展开更多
In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for th...In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for the corresponding problem is obtained.展开更多
A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedur...A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.展开更多
We have successfully ported an arbitrary highorder discontinuous Galerkin method for solving the threedimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units...We have successfully ported an arbitrary highorder discontinuous Galerkin method for solving the threedimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units (GPUs) using the Compute Unified Device Architecture (CUDA) of NVIDIA and Message Passing Interface (MPI) and obtained a speedup factor of about 28.3 for the single-precision version of our codes and a speedup factor of about 14.9 for the double-precision version. The GPU used in the comparisons is NVIDIA Tesla C2070 Fermi, and the CPU used is Intel Xeon W5660. To effectively overlap inter-process communication with computation, we separate the elements on each subdomain into inner and outer elements and complete the computation on outer elements and fill the MPI buffer first. While the MPI messages travel across the network, the GPU performs computation on inner elements, and all other calculations that do not use information of outer elements from neighboring subdomains. A significant portion of the speedup also comes from a customized matrix-matrix multiplication kernel, which is used extensively throughout our program. Preliminary performance analysis on our parallel GPU codes shows favorable strong and weak scalabilities.展开更多
基金National Natural Science Foundation of China ( No. 10972065) Natural Science Foundation of Heilongjiang Province of China( No. ZD200905)
文摘The vibration suppression of the finite plate with square steel beams is studied using traveling wave method. The finite plate with square beams is modeled as the coupling systems between the plate flexural motion and the flexural and torsional motions for the square beams. The vibration response at any position of the coupling structure can be obtained by wave method. Numerical results show that comparing to finite element method (FEM), not only the low frequency but also the medium-high frequency vibration response of the finite plate with square beam can be effectively calculated by wave method. The suppression effect can be increased as the square beam is located at one-third of the length of plate or increasing the height of the beam. The study provides reference for arranged square beams applying to vibration suppression of ship and train structures.
基金supported by National Natural Science Foundation of China (No. 41174085)Chinese Academy of Sciences (KZZD-EW-TZ-19)China Geological Survey (12120113101400)
文摘Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m under- ground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.
文摘The plane wave numerical technique is recast from Ampere’s and Faraday’s laws for materials that are characterized with a bianisotropic form of the constitutive relations. The populating expressions are provided for the eigenvalue matrix system that can be directly solved for the angular frequencies and field profiles when bianisotropy is included. To demonstrate the computation process and expected state diagrams and field profiles, numerical computation examples are provided for a bianisotropic Bragg Array with central defect. It is shown that the location of the magnetoelectric tensor elements has a significant effect on the eigenstates of an equivalent isotropic (anisotropic) structure. One form of the magnetoelectric tensor (diagonal elements only) leads to the observation of merging states and the formation of exceptional points. The numerical approach presented can be implemented as an add-on to the familiar plane wave numerical technique.
基金the support of a grant from Aeronautical Science Foundation of China(2010ZA53013 and 2011ZA53014)the open funds of Key Laboratory of Advanced Design and Intelligent Computing(Dalian University),Ministry of Education(ADIC2010007)Northwestern Polytechnical University Basic Research Fund(JC201114 andJC20110255)
文摘The wave method is introduced to vibration analysis of the fluid-conveying carbon nanotube. The constitutive relation of carbon nanotube on micro-scale is founded using the non- local elastic theory. The governing equation on micro-scale is obtained. And the first five orders of the natural frequency of the carbon nanotube conveying fluid with various speeds are calculated through the wave method. Besides, the critical flow velocity when the carbon nanotube loses stability is obtained. Meanwhile, a contrast is made between the result obtained through tile wave method and that in previous researches.
基金the National Natural Science Foundation of China (No. 69631020) and theOffice of Naval Research of America (00014-93-1-0340).
文摘A low-frequency multi-mode ultrasonic Lamb wave method suitable for character- izing the thickness, the density and the elastic constants of the ultra-thin transversely isotropic laminate composite is presented. The 'ultra-thin' here means that the thickness of the plate is much less than the wavelength of the ultrasonic wave so that the echoes from the front and back faces of the plate can't be separated in the time domain. The dispersion equations for the low frequency ultrasonic Lamb waves with the propagation directions parallel and vertical to the fiber direction are derived. In conjunction with the least square algorithm method, the secant algorithm is used to estimate the parameters of the ultra-thin fiber-reinforced composite layer. The evaluation errors and the sensitivity of the method to different paramters of the thin composite are analyzed. The technique has been used to characterize the ultra-thin grass fiber reinforced PES composite with thickness down to ten percents of the ultrasonic wavelength. It is observed that the agreement between the nominal and the estimation values is reasonably good.
文摘The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6×6 Luttinger Kohn model.The effect of the number and period of plane waves used for expansion on the stability of energy eigenvalues is examined.For practical calculation,it should choose the period large sufficiently to ensure the envelope functions vanish at the boundary and the number of plane waves large enough to ensure the energy eigenvalues keep unchanged within a prescribed range.
基金Project supported by the National Basic Research Program of China (973Project) (No.2002CB412709) and the National Natural Science Foundation of China (Nos.50278012,10272027,19832010)
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
基金National Natural Science Foundation of China under Grant No.51478247National Key Research and Development Program of China under Grant No.2016YFC1402800
文摘The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.
基金Project supported by the National Natural Science Foundation of China (No. 10172038).
文摘By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at the eigenfrequencies of the corresponding interior problems is proposed. By adding appropriate damp to the virtual source system of the wave superposition method, the unique solutions for all wave numbers can be ensured. Based on this thought, a novel method-wave superposition method with complex radius vector is constructed. Not only is the computational time of this method approximately equal to that of the standard wave superposition method, but also the accuracy is much higher compared with other correlative methods. Finally, by taking the pulsating sphere and oscillating sphere as examples, the results of calculation show that the present method can effectively overcome the non-uniqueness problem.
基金Funded by the National Natural Science Foundation of China (No.50708065)the National High-tech R&D Program(863 Program )(No.2007-AA-11-Z-113)the Key Projects in the Science and Technology Pillar Program of Tianjin(No.11ZCKFSF00300)
文摘In order to monitor the basic mechanical properties and interior damage of concrete structures,the piezoelectric actuator/sensor based wave propagation method was investigated experimentally in the laboratory using a specifically designed test setup.The energy attenuation of stress waves was measured by the relative index between the output voltage of sensors and the excitation voltage at the actuator.Based on the experimental results of concrete cube and cylinder specimens,the effect of excitation frequencies,excitation amplitude,wave propagation paths and the curing age on the output signals of sensors are evaluated.The results show that the relative voltage attenuation coefficient RVAC is an effective indicator for measuring the attenuation of stress waves through the interior of concrete.
文摘Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orientations, etc. of the three array arrangements, i. e., linear array, T-type array and pentagon array, on the resolution of the directional spreading of waves, are investigated experimentally. This study can be used as a reference in the experimental study and the field measurement of directional waves.
基金The Dragon III Project of ESA-MOST Dragon Cooperation under contract No.10412the Ocean Renewable Energy Special Fund Project of State Oceanic Administration under contract No.GHME2011ZC07the National Natural Science Foundation of China(NSFC)under contract No.41176157
文摘Wave energy resource is a very important ocean renewable energy. A reliable assessment of wave energy resources must be performed before they can be exploited. Compared with wave model, altimeter can provide more accurate in situ observations for ocean wave which can be as a novel method for wave energy assessment.The advantage of altimeter data is to provide accurate significant wave height observations for wave. In order to develop characteristic and advantage of altimeter data and apply altimeter data to wave energy assessment, in this study, we established an assessing method for wave energy in local sea area which is dedicated to altimeter data.This method includes three parts including data selection and processing, establishment of evaluation indexes system and criterion of regional division. Then a case study of Northwest Pacific was performed to discuss specific application for this method. The results show that assessing method in this paper can assess reserves and temporal and spatial distribution effectively and provide scientific references for the siting of wave power plants and the design of wave energy convertors.
基金supported by the‘New Century Outstanding Talent’Scheme of the Ministry of Education of China(Grant No.NCET-07-0135)
文摘A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid (VOF) method. An active 'absorbing wave-maker based on linear wave theory is set on the left boundary of the wave flume. The progressive waves and the absorbing waves are generated simultaneously at the active wave generating-absorbing boundary. The absorbing waves are generated to eliminate the waves coming back to the generating boundary due to reflection from the outflow boundary and the structures. SIRW method proposed by Frigaard and Brorsen (1995) is used to separate the incident waves and reflected waves. The digital filters are designed based on the surface elevation signals of the two wave gauges. The corrected velocity of the wave-maker paddle is the output from the digital filter in real time. The numerical results of regular and irregular waves by the active absorbing-generating boundary are compared with the numerical results by the ordinary generating boundary to verify the performance of the active absorbing-generator boundary. The differences between the initial incident waves and the estimated incident waves are analyzed.
基金Supported by the International Cooperation and Exchanges Foundation of Henan Province (084300510060)the Youth Science Foundation of Henan University of Science and Technology of China (2008QN026)
文摘The (G'/G, 1/G)-expansion method for finding exact travelling wave solutions of nonlinear evolution equations, which can be thought of as an extension of the (G'/G)-expansion method proposed recently, is presented. By using this method abundant travelling wave so- lutions with arbitrary parameters of the Zakharov equations are successfully obtained. When the parameters are replaced by special values, the well-known solitary wave solutions of the equations are rediscovered from the travelling waves.
文摘Nozzle damping is one of the most important factors in the suppression of combustion instability in solid rocket motors.For an engineering solid rocket motor that experiences combustion instability at the end of burning,a wave attenuation method is proposed to assess the nozzle damping characteristics numerically.In this method,a periodic pressure oscillation signal which frequency equals to the first acoustic mode is superimposed on a steady flow at the head end of the chamber.When the pressure oscillation is turned off,the decay rate of the pressure can be used to determine the nozzle attenuation constant.The damping characteristics of three other nozzle geometries are numerically studied with this method under the same operating condition.The results show that the convex nozzle provides more damping than the conical nozzle which in turn provides more damping than the concave nozzle.All the three nozzles have better damping effect than that of basic nozzle geometry.At last,the phase difference in the chamber is analyzed,and the numerical pressure distribution satisfies well with theoretical distribution.
基金supported by the National Natural Science Foundation of China (11171208)Shanghai Leading Academic Discipline Project (S30106)
文摘The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.
基金Supported by the National Natural Sciences Foundation of China(40676016 and 10471039)the National Key Project for Basic Research(2003CB415101-03 and 2004CB418304)+2 种基金the Key Project of the Chinese Academy of Sciences(KZCX3-SW-221)in part by E-Institutes of Shanghai Municipal Education Commission(N.E03004)the Natural Science Foundation of Zeijiang,China(Y606268).
文摘In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for the corresponding problem is obtained.
基金The project was financially supported by the National Natural Science Foundation of China under the Grant No. 19732004 the National Science Fund for Distinguished Young Scholars under the Grant No. 50029002
文摘A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.
基金supported by the School of Energy Resources at the University of WyomingThe GPU hardware used in this study was purchased using the NSF Grant EAR-0930040
文摘We have successfully ported an arbitrary highorder discontinuous Galerkin method for solving the threedimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units (GPUs) using the Compute Unified Device Architecture (CUDA) of NVIDIA and Message Passing Interface (MPI) and obtained a speedup factor of about 28.3 for the single-precision version of our codes and a speedup factor of about 14.9 for the double-precision version. The GPU used in the comparisons is NVIDIA Tesla C2070 Fermi, and the CPU used is Intel Xeon W5660. To effectively overlap inter-process communication with computation, we separate the elements on each subdomain into inner and outer elements and complete the computation on outer elements and fill the MPI buffer first. While the MPI messages travel across the network, the GPU performs computation on inner elements, and all other calculations that do not use information of outer elements from neighboring subdomains. A significant portion of the speedup also comes from a customized matrix-matrix multiplication kernel, which is used extensively throughout our program. Preliminary performance analysis on our parallel GPU codes shows favorable strong and weak scalabilities.