In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependen...In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependence. Test problems including the single soliton wave motion, interaction of two solitary waves and interaction of three solitary waves will use to validate the proposed method. The three invariants of the motion are evaluated to determine the conservation properties of the generated scheme. Finally, a Maxwellian initial condition pulse is then studied. The L<sub>2</sub> and L<sub>∞</sub> error norms are computed to study the accuracy and the simplicity of the presented method.展开更多
Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harm...Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.展开更多
A base function expressed with Chebyshev polynomials is reached. The relationship between the coefficients of the partial differential equation and the base function is deduced. Using the relationship, one can obtain ...A base function expressed with Chebyshev polynomials is reached. The relationship between the coefficients of the partial differential equation and the base function is deduced. Using the relationship, one can obtain nearly the same results as those calculated by Fast Fourier Transformation (FFT). The pseudo-spectral matrix method is applied in this paper to simulate numerically the incompressible laminar boundary flow on a plate. The simulation proves to be precise and efficient.展开更多
We propose a new scheme for simulation of a high-order nonlinear Schrodinger equation with a trapped term by using the mid-point rule and Fourier pseudospectral method to approximate time and space derivatives, respec...We propose a new scheme for simulation of a high-order nonlinear Schrodinger equation with a trapped term by using the mid-point rule and Fourier pseudospectral method to approximate time and space derivatives, respectively. The method is proved to be both charge- and energy-conserved. Various numerical experiments for the equation in different cases are conducted. From the numerical evidence, we see the present method provides an accurate solution and conserves the discrete charge and energy invariants to machine accuracy which are consistent with the theoretical analysis.展开更多
This paper presents the analysis of open microstrip structures by using diakoptic method of lines (ML) combined with periodic boundary conditions (PBC). The parameters of microstrip patch are obtained from patch curre...This paper presents the analysis of open microstrip structures by using diakoptic method of lines (ML) combined with periodic boundary conditions (PBC). The parameters of microstrip patch are obtained from patch current excited by plane wave. Impedance matrix elements are computed by using fast Fourier transform(FFT), and reduced equation is solved by using diakoptic technique. Consequently, the computing time is reduced significantly. The convergence property of simulating open structure by using PBC and the comparison of the computer time between using PBC and usual absorbing boundary condition (ABC) show the validity of the method proposed in this paper. Finally, the resonant frequency of a microstrip patch is computed. The numerical results obtained are in good agreement with those published.展开更多
针对现阶段随机冰力时域模拟方法单一的问题,将一种经典的时域模拟方法——快速傅里叶逆变换(Inverse Fast Fourier Transform,IFFT)法,引入海上风机单桩基础随机冰载荷研究中。结合随机冰力模型,采用MATLAB软件平台编写冰载荷时域模拟...针对现阶段随机冰力时域模拟方法单一的问题,将一种经典的时域模拟方法——快速傅里叶逆变换(Inverse Fast Fourier Transform,IFFT)法,引入海上风机单桩基础随机冰载荷研究中。结合随机冰力模型,采用MATLAB软件平台编写冰载荷时域模拟程序。基于山东省东营市渤中某海上风电场8.5 MW单桩式风机的结构与尺寸,进行算例分析,考虑入射角度作用,得到随机冰力在单桩基础结构上的实际分布情况,对IFFT法模拟随机冰力的精度和速度进行综合评价。结果表明,IFFT法作为一种计算简单快速、结果精确的时域模拟方法,适用于单桩基础随机冰力的计算,且所编程序可实现输入少量参数,快速高效获得冰力数值与分布的目的,具有实用性和推广性。展开更多
文摘In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependence. Test problems including the single soliton wave motion, interaction of two solitary waves and interaction of three solitary waves will use to validate the proposed method. The three invariants of the motion are evaluated to determine the conservation properties of the generated scheme. Finally, a Maxwellian initial condition pulse is then studied. The L<sub>2</sub> and L<sub>∞</sub> error norms are computed to study the accuracy and the simplicity of the presented method.
基金Project supported by the National Natural Science Foundation of China (No.10172038)
文摘Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.
文摘A base function expressed with Chebyshev polynomials is reached. The relationship between the coefficients of the partial differential equation and the base function is deduced. Using the relationship, one can obtain nearly the same results as those calculated by Fast Fourier Transformation (FFT). The pseudo-spectral matrix method is applied in this paper to simulate numerically the incompressible laminar boundary flow on a plate. The simulation proves to be precise and efficient.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11201169 and 11271195)the Qing Lan Project of Jiangsu Province,China
文摘We propose a new scheme for simulation of a high-order nonlinear Schrodinger equation with a trapped term by using the mid-point rule and Fourier pseudospectral method to approximate time and space derivatives, respectively. The method is proved to be both charge- and energy-conserved. Various numerical experiments for the equation in different cases are conducted. From the numerical evidence, we see the present method provides an accurate solution and conserves the discrete charge and energy invariants to machine accuracy which are consistent with the theoretical analysis.
基金Supported by the National Natural Science Foundation of China
文摘This paper presents the analysis of open microstrip structures by using diakoptic method of lines (ML) combined with periodic boundary conditions (PBC). The parameters of microstrip patch are obtained from patch current excited by plane wave. Impedance matrix elements are computed by using fast Fourier transform(FFT), and reduced equation is solved by using diakoptic technique. Consequently, the computing time is reduced significantly. The convergence property of simulating open structure by using PBC and the comparison of the computer time between using PBC and usual absorbing boundary condition (ABC) show the validity of the method proposed in this paper. Finally, the resonant frequency of a microstrip patch is computed. The numerical results obtained are in good agreement with those published.
文摘针对现阶段随机冰力时域模拟方法单一的问题,将一种经典的时域模拟方法——快速傅里叶逆变换(Inverse Fast Fourier Transform,IFFT)法,引入海上风机单桩基础随机冰载荷研究中。结合随机冰力模型,采用MATLAB软件平台编写冰载荷时域模拟程序。基于山东省东营市渤中某海上风电场8.5 MW单桩式风机的结构与尺寸,进行算例分析,考虑入射角度作用,得到随机冰力在单桩基础结构上的实际分布情况,对IFFT法模拟随机冰力的精度和速度进行综合评价。结果表明,IFFT法作为一种计算简单快速、结果精确的时域模拟方法,适用于单桩基础随机冰力的计算,且所编程序可实现输入少量参数,快速高效获得冰力数值与分布的目的,具有实用性和推广性。
