A perturbation method is introduced in the context of dynamical system for solving the nonlinear Korteweg-de Vries (KdV) equation. Best efficiency is obtained for few perturbative corrections. It is shown that, the qu...A perturbation method is introduced in the context of dynamical system for solving the nonlinear Korteweg-de Vries (KdV) equation. Best efficiency is obtained for few perturbative corrections. It is shown that, the question of convergence of this approach is completely guaranteed here, because a limited number of term included in the series can describe a sufficient exact solution. Comparisons with the solutions of the quintic spline, and finite difference are presented.展开更多
A B-spline with the symplectic algorithm method for the solution of time-dependent Schrodinger equations （TDSEs） is introduced. The spatial part of the wavefunction is expanded by B-spline and the time evolution is ...A B-spline with the symplectic algorithm method for the solution of time-dependent Schrodinger equations （TDSEs） is introduced. The spatial part of the wavefunction is expanded by B-spline and the time evolution is given in a symplectic scheme. This method allows us to obtain a highly accurate and stable solution of TDSEs. The effectiveness and efficiency of this method is demonstrated by the high-order harmonic spectra of one-dimensional atoms in comparison with other references.展开更多
This paper present an implementation of＂modified cubic B-spline differential quadrature method （MCB-DQM）＂ proposed by Arora ＆ Singh （Applied Mathematics and Computation Vol. 224（1） （2013） 161-177） for numer...This paper present an implementation of＂modified cubic B-spline differential quadrature method （MCB-DQM）＂ proposed by Arora ＆ Singh （Applied Mathematics and Computation Vol. 224（1） （2013） 161-177） for numerical computation of Fokker-Planck equations. The modified cubic B-splines are used as set of basis functions in the differential quadrature to compute the weighting coefficients for the spatial derivatives, which reduces Fokker-Planck equation into system of first-order ordinary differential equations （ODEs）, in time. The well known SSP-RK43 scheme is then applied to solve the resulting system of ODEs. The efficiency of proposed method has been confirmed by three examples having their exact solutions. This shows that MCB-DQM results are capable of achieving high accuracy. Advantage of the scheme is that it can be applied very smoothly to solve the linear or nonlinear physical problems, and a very less storage space is required which causes less accumulation of numerical errors.展开更多
This paper is to discuss an approach which combines B-spline patches and transfinite interpolation to establish a linear algebraic system for solving partial differential equations and modify the WEB-spline method dev...This paper is to discuss an approach which combines B-spline patches and transfinite interpolation to establish a linear algebraic system for solving partial differential equations and modify the WEB-spline method developed by Klaus Hollig to derive this new idea. First of all, the authors replace the R-function method with transfinite interpolation to build a function which vanishes on boundaries. Secondly, the authors simulate the partial differential equation by directly applying differential opera- tors to basis functions, which is similar to the RBF method rather than Hollig＇s method. These new strategies then make the constructing of bases and the linear system much more straightforward. And as the interpolation is brought in, the design of schemes for solving practical PDEs can be more flexi- ble. This new method is easy to carry out and suitable for simulations in the fields such as graphics to achieve rapid rendering. Especially when the specified much faster than WEB-spline method. precision is not very high, this method performs展开更多
Instead of most existing postprocessing schemes, a new preprocessing approach, called multi- neighboring grids （MNG）, is proposed for solving PDE eigen-problems on an existing grid G（A）. The linear or multi-linear...Instead of most existing postprocessing schemes, a new preprocessing approach, called multi- neighboring grids （MNG）, is proposed for solving PDE eigen-problems on an existing grid G（A）. The linear or multi-linear element, based on box-splines, are taken as the first stage Khuh -λh/1Mh/1Uh. In this paper, the j-th stage neighboring-grid scheme is defined as Khuh λh/j Mh/j Uh = λh/j Mh/j Uh , where gh ：- Mh/j-1 Kh/1 and Mhuh is to be found as a better mass distribution over the j-th stage neighboring-grid G（/k）, and Kh/1 can be seen as an expansion of Kh on the j-th neighboring-grid with respect to the （j - 1）-th mass distribution Mh_l. It is shown that for an ODE model eigen-problem, the j-th stage scheme with 2j-th order B-spline basis can reach 2j-th order accuracy and even （2j ＋ 2）-th order accuracy by perturbing the mass matrix. The argument can be extended to high dimensions with separable variable cases. For Laplace eigen-problems with some 2-D and 3-D structured uniform grids, some 2j-th order schemes are presented for j ≤ 3.展开更多
基金Project (No. D0701/01/05) supported by Ministry of the Educationand Scientific Research (M.E.S.R), Algeria
文摘A perturbation method is introduced in the context of dynamical system for solving the nonlinear Korteweg-de Vries (KdV) equation. Best efficiency is obtained for few perturbative corrections. It is shown that, the question of convergence of this approach is completely guaranteed here, because a limited number of term included in the series can describe a sufficient exact solution. Comparisons with the solutions of the quintic spline, and finite difference are presented.
基金Supported by the National Natural Science Foundation of China under Grant No 10374119, and the 0ne-Hundred-Talents Project of Chinese Academy of Science. ACKN0WLEDGMENTS： We gratefully acknowledge Professors Ding Peizhu and Liu Xueshen for their hospitality and help with the symplectic al- gorithm.
文摘A B-spline with the symplectic algorithm method for the solution of time-dependent Schrodinger equations （TDSEs） is introduced. The spatial part of the wavefunction is expanded by B-spline and the time evolution is given in a symplectic scheme. This method allows us to obtain a highly accurate and stable solution of TDSEs. The effectiveness and efficiency of this method is demonstrated by the high-order harmonic spectra of one-dimensional atoms in comparison with other references.
文摘This paper present an implementation of＂modified cubic B-spline differential quadrature method （MCB-DQM）＂ proposed by Arora ＆ Singh （Applied Mathematics and Computation Vol. 224（1） （2013） 161-177） for numerical computation of Fokker-Planck equations. The modified cubic B-splines are used as set of basis functions in the differential quadrature to compute the weighting coefficients for the spatial derivatives, which reduces Fokker-Planck equation into system of first-order ordinary differential equations （ODEs）, in time. The well known SSP-RK43 scheme is then applied to solve the resulting system of ODEs. The efficiency of proposed method has been confirmed by three examples having their exact solutions. This shows that MCB-DQM results are capable of achieving high accuracy. Advantage of the scheme is that it can be applied very smoothly to solve the linear or nonlinear physical problems, and a very less storage space is required which causes less accumulation of numerical errors.
基金partially supported by NKBRSF under Grant No.2011CB302404NSFC under Grant Nos. 10871195,10925105,60821002,and 50875027
文摘This paper is to discuss an approach which combines B-spline patches and transfinite interpolation to establish a linear algebraic system for solving partial differential equations and modify the WEB-spline method developed by Klaus Hollig to derive this new idea. First of all, the authors replace the R-function method with transfinite interpolation to build a function which vanishes on boundaries. Secondly, the authors simulate the partial differential equation by directly applying differential opera- tors to basis functions, which is similar to the RBF method rather than Hollig＇s method. These new strategies then make the constructing of bases and the linear system much more straightforward. And as the interpolation is brought in, the design of schemes for solving practical PDEs can be more flexi- ble. This new method is easy to carry out and suitable for simulations in the fields such as graphics to achieve rapid rendering. Especially when the specified much faster than WEB-spline method. precision is not very high, this method performs
基金supported by National Natural Science Foundation of China(Grant Nos.6097008961170075 and 91230109)
文摘Instead of most existing postprocessing schemes, a new preprocessing approach, called multi- neighboring grids （MNG）, is proposed for solving PDE eigen-problems on an existing grid G（A）. The linear or multi-linear element, based on box-splines, are taken as the first stage Khuh -λh/1Mh/1Uh. In this paper, the j-th stage neighboring-grid scheme is defined as Khuh λh/j Mh/j Uh = λh/j Mh/j Uh , where gh ：- Mh/j-1 Kh/1 and Mhuh is to be found as a better mass distribution over the j-th stage neighboring-grid G（/k）, and Kh/1 can be seen as an expansion of Kh on the j-th neighboring-grid with respect to the （j - 1）-th mass distribution Mh_l. It is shown that for an ODE model eigen-problem, the j-th stage scheme with 2j-th order B-spline basis can reach 2j-th order accuracy and even （2j ＋ 2）-th order accuracy by perturbing the mass matrix. The argument can be extended to high dimensions with separable variable cases. For Laplace eigen-problems with some 2-D and 3-D structured uniform grids, some 2j-th order schemes are presented for j ≤ 3.
