In this paper,we propose a numerical method to estimate the unknown order of a Riemann-Liouville fractional derivative for a fractional Stokes' first problem for a heated generalized second grade fluid.The implicit n...In this paper,we propose a numerical method to estimate the unknown order of a Riemann-Liouville fractional derivative for a fractional Stokes' first problem for a heated generalized second grade fluid.The implicit numerical method is employed to solve the direct problem.For the inverse problem,we first obtain the fractional sensitivity equation by means of the digamma function,and then we propose an efficient numerical method,that is,the Levenberg-Marquardt algorithm based on a fractional derivative,to estimate the unknown order of a Riemann-Liouville fractional derivative.In order to demonstrate the effectiveness of the proposed numerical method,two cases in which the measurement values contain random measurement error or not are considered.The computational results demonstrate that the proposed numerical method could efficiently obtain the optimal estimation of the unknown order of a RiemannLiouville fractional derivative for a fractional Stokes' first problem for a heated generalized second grade fluid.展开更多
A method to compute the numerical derivative of eigenvalues of parameterized crystal field Hamiltonian matrix is given, based on the numerical derivatives the general iteration methods such as Levenberg-Marquardt, New...A method to compute the numerical derivative of eigenvalues of parameterized crystal field Hamiltonian matrix is given, based on the numerical derivatives the general iteration methods such as Levenberg-Marquardt, Newton method, and so on, can be used to solve crystal field parameters by fitting to experimental energy levels. With the numerical eigenvalue derivative, a detailed iteration algorithm to compute crystal field parameters by fitting experimental energy levels has also been described. This method is used to compute the crystal parameters of Yb^3+ in Sc2O3 crystal, which is prepared by a co-precipitation method and whose structure was refined by Rietveld method. By fitting on the parameters of a simple overlap model of crystal field, the results show that the new method can fit the crystal field energy splitting with fast convergence and good stability.展开更多
基金supported by the National Natural Science Foundation of China(Grants 11472161,11102102,and 91130017)the Independent Innovation Foundation of Shandong University(Grant 2013ZRYQ002)the Natural Science Foundation of Shandong Province(Grant ZR2014AQ015)
文摘In this paper,we propose a numerical method to estimate the unknown order of a Riemann-Liouville fractional derivative for a fractional Stokes' first problem for a heated generalized second grade fluid.The implicit numerical method is employed to solve the direct problem.For the inverse problem,we first obtain the fractional sensitivity equation by means of the digamma function,and then we propose an efficient numerical method,that is,the Levenberg-Marquardt algorithm based on a fractional derivative,to estimate the unknown order of a Riemann-Liouville fractional derivative.In order to demonstrate the effectiveness of the proposed numerical method,two cases in which the measurement values contain random measurement error or not are considered.The computational results demonstrate that the proposed numerical method could efficiently obtain the optimal estimation of the unknown order of a RiemannLiouville fractional derivative for a fractional Stokes' first problem for a heated generalized second grade fluid.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.50772112 and 50872135)the Natural Science Foundation of Anhui Province of China(Grant No.08040106820)the Knowledge Innovation Program of the Chinese Academy of Sciences(Grant No.YYYJ-1002)
文摘A method to compute the numerical derivative of eigenvalues of parameterized crystal field Hamiltonian matrix is given, based on the numerical derivatives the general iteration methods such as Levenberg-Marquardt, Newton method, and so on, can be used to solve crystal field parameters by fitting to experimental energy levels. With the numerical eigenvalue derivative, a detailed iteration algorithm to compute crystal field parameters by fitting experimental energy levels has also been described. This method is used to compute the crystal parameters of Yb^3+ in Sc2O3 crystal, which is prepared by a co-precipitation method and whose structure was refined by Rietveld method. By fitting on the parameters of a simple overlap model of crystal field, the results show that the new method can fit the crystal field energy splitting with fast convergence and good stability.