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.展开更多
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.展开更多
文摘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.
基金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.
