×

重温托马斯·费尔米方程:通过坐标变换加速有理切比雪夫级数。 (英语) Zbl 1507.65260号

小结:我们用(u(0)=1和(u(infty)=0)重新讨论了中性原子Thomas-Ferm问题的光谱解,以说明在存在复杂情况时求解微分方程的一些主题,并对我们早期的处理进行了改进。我们所说的“复杂性”是指问题的特征,这些特征要么破坏了标准切比雪夫级数的指数精度,要么使经典切比雪夫法不适用。Thomas-Fermi问题有四个复杂问题:(i)一个半无限域(在[0,infty]\中)(ii)原点处的(u(r)\中的平方根奇异性(iii)分数幂非线性和(iv)渐近衰减为\(r\rightarrow\infty\),其中包括\(r)的负幂和分数指数。我们早期的处理方法确定了原点处的斜率到小数点后25位,但需要不少于600个基函数来近似一个处处单调的单变量解,并且所有早期的技巧都无法在谱序列的截断中恢复指数收敛率,但只有负幂的高阶收敛性。这里,使用坐标(z\equiv\sqrt{r})中和平方根奇异性,我们表明通过求解原始未知(u(r))而不是之前使用的修改未知(v(r)=\sqrt},可以显著提高精度和收敛速度。在不进一步改变坐标的情况下,有理切比雪夫基(T L_n(z;L))产生原点斜率的十二位小数精度,(u_r(0))具有70个基函数,二十四位小数具有截断(n=100)。
通过适当改变坐标,可以恢复真正的指数精度,其中,(z=mathcal{G}(z)是指数的某些种类。然而,Thomas-Fermi函数的各种Chebyshev级数和Fourier级数具有“复数渐近性”,即对于某个正常数(n_I),(a_n\sim a^{mathit{intermediate}}(n))具有“复式渐近性”。“远非症状”作为(n\rightarrow\infty)通常没有实际意义。对于这个问题,理论上最适合于巨型(n)的TL-with-sinh方法受到了数值病态条件的困扰,只有当目标精度至少为40位小数时,它的渐近优越性才得以实现。
这对于工程来说是荒谬的,但对于基准测试来说也许是有用的。到小数点后六十位,\(u_r(0)=-1.5880710226113753127186845094239501094527466216748256167677 \)。将Coulson-March渐近级数中的常数改进为小数点后的19位(F=13.2709738480269351535)。

MSC公司:

65号35 偏微分方程边值问题的谱、配置及相关方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Allgower,E.L。;Georg,K.,《数值延拓方法导论》,《应用数学经典》,第45卷,(2003),费城工业和应用数学学会·Zbl 1036.65047号
[2] Baker,E.B.,《费米-托马斯统计模型在正离子电位分布计算中的应用》,《物理学》。修订版,36,630-647,(1930)·JFM 56.1327.06号机组
[3] Bayatbabolghani,F。;Parand,K.,《使用Hermite函数求解Thomas-Fermi方程》,《国际数学杂志》。计算。科学。工程,8123-126,(2014)
[4] 本德,C.M。;Milton,K.A。;Pinsky,S.S。;Simmons,L.M.,非线性问题的一种新的微扰方法,数学杂志。物理。,30, 1447-1455, (1989) ·Zbl 0684.34008号
[5] Bethe,H.,Thomas-Fermi核理论,物理学。修订版,167879-907(1968)
[6] Boyd,J.P.,《二维bratu方程的分析和数值研究》,J.Sci。计算。,1, 183-206, (1985) ·Zbl 0649.65057号
[7] Boyd,J.P.,半无限区间上的正交有理函数,J.Compute。物理。,70, 63-88, (1987) ·2013年6月14日Zbl
[8] Boyd,J.P.,在无限区间上使用有理基函数的谱方法,J.Compute。物理。,69, 112-142, (1987) ·Zbl 0615.65090号
[9] Boyd,J.P.,Chebyshev域截断在求解无限区间问题时不如Fourier域截断,J.Sci。计算。,3, 109-120, (1988) ·Zbl 0666.65099号
[10] Boyd,J.P.,具有对数端点奇点的函数的渐近Chebyshev系数,应用。数学。计算。,29,49-67,(1989年)·Zbl 0668.65017号
[11] Boyd,J.P.,应用于无限区间伪谱计算的weideman-Clot-sinh-mapping的无限级整函数的Fourier系数的收敛速度,J.Compute。物理。,110, 360-372, (1994) ·兹伯利0806.65146
[12] Boyd,J.P.,《弱非局部孤立波和超全阶渐近:广义孤子和超共鸣微扰理论》,《数学及其应用》,第442卷,(1998),Kluwer Amsterdam,608页·Zbl 0905.76001号
[13] Boyd,J.P.,Chebyshev和Fourier光谱方法,(2001年),纽约多佛·米诺拉,665页·Zbl 0994.65128号
[14] Boyd,J.P.,删除残差、QR系数牛顿迭代和其他方法,用于孤立波、椭圆余弦波和冲击波非线性特征问题的形式上超定确定离散化,J.Compute。物理。,179, 216-237, (2002) ·Zbl 1006.65116号
[15] Boyd,J.P.,解决Chandrasekhar H函数对数端点奇异性的光谱精确求积,J.Quant。光谱学。辐射。传输。,94, 467-475, (2005)
[16] Boyd,J.P.,《为什么牛顿的方法很难处理行波:小分母,KAM理论,阿诺德的线性傅里叶问题,非唯一性,约束和不稳定故障,数学》。计算。模拟。,74, 72-81, (2007) ·Zbl 1107.65351号
[17] Boyd,J.P.,Thomas-Fermi函数的有理切比雪夫级数:端点奇异性和谱方法,J.Compute。申请。数学。,231, 90-101, (2013) ·Zbl 1260.65065号
[18] Boyd,J.P.,《求解超越方程:切比雪夫多项式代理和其他数值寻根器、扰动级数和预言》,(2014),费城SIAM出版社,460页·Zbl 1311.65047号
[19] Boyd,J.P.,《追踪非线性常微分和积分方程的多解分支:切比雪夫谱方法和傅里叶谱方法以及递增谱同伦DISH》,J.Sci。计算。,1115-1143年,(2016年)·Zbl 1371.65070号
[20] Boyd,J.P.,孤立波和生物的强非线性微扰理论,Evol。等于。《控制理论》(2018),出版
[21] 博伊德,J.P。;Gheorghiu,C.I.,《所有根谱方法:约束、浮点算法和根排除》,应用。数学。莱特。,67, 28-32, (2017) ·Zbl 1358.65049号
[22] 博伊德,J.P。;Petschek,R.,Chebyshev、Legendre和Jacobi多项式之间的关系:Chebyshef多项式的一般优越性和三个重要例外,J.Sci。计算。,59, 1-27, (2014) ·Zbl 1301.41003号
[23] 博伊德,J.P。;Rangan,C。;Buchsbaum,P.H.,《应用于氢原子的半无限区间上的伪谱方法:映射傅里叶方法与拉盖尔级数和有理切比雪夫展开式的比较》,J.Compute。物理。,188, 56-74, (2003) ·Zbl 1028.65086号
[24] Cloot,A.,用于无界域计算的等分布映射和谱方法,Appl。数学。莱特。,4, 23-27, (1991) ·Zbl 0724.65091号
[25] Cloot,A。;Weideman,J.A.C.,无界域上谱计算的自适应算法,J.Comput。物理。,102, 398-406, (1992) ·Zbl 0760.65095号
[26] 科尔森,C.A。;新罕布什尔州March,《使用托马斯·费尔米方法的原子动量》,Proc。R.物理。社会学硕士,63367-374,(1950)·Zbl 0035.42302号
[27] 丹尼斯,J.E。;Schnabel,R.B.,非线性方程和无约束优化的数值方法,应用数学经典,第16卷,(1987),费城SIAM
[28] Deufhard,P.,《非线性问题的牛顿方法:仿射不变性和自适应算法》,(2004),施普林格-弗拉格-海德堡·Zbl 1056.65051号
[29] Fermi,E.,Eine统计方法,Z.Phys。,48, 73-79, (1928) ·JFM 54.0984.02号
[30] Fernandez,F.M.,Thomas-Fermi方程的有理逼近,应用。数学。计算。,217, 6433-6436, (2011) ·Zbl 1210.81133号
[31] 费曼,R.P。;北卡罗来纳州大都会。;Teller,E.,基于广义费米-托马斯理论的元素状态方程,物理学。修订版,751561-1573,(1949)·Zbl 0036.43007号
[32] Fornberg,B.,《伪谱方法实用指南》,(1996),剑桥大学出版社,纽约·Zbl 0844.65084号
[33] 基里斯曼,A。;哈希姆,I。;M.T.卡贾尼。;Maleki,M.,关于求解无限区间上Thomas-Fermi方程的有理第二类Chebyshev伪谱方法,J.Compute。申请。数学。,257, 79-85, (2014) ·Zbl 1294.65080号
[34] 小林,S。;松山,T。;Nagi,S。;Umeda,K.,普通T-F函数初始斜率的精确值,J.Phys。Soc.Jpn.公司。,10, 759-762, (1955)
[35] Lund,J。;Bowers,K.L.,求积和微分方程的Sinc方法,(1992),费城工业和应用数学学会,304页·Zbl 0753.65081号
[36] 曼德尔茨威格,V.B。;Tabakin,F.,《物理中非线性问题的拟线性化方法及其在非线性模型中的应用》,计算。物理学。社区。,141, 268-281, (2001) ·Zbl 0991.65065号
[37] Mason,J.C.,普通Thomas-Fermi函数及其导数的有理逼近,Proc。物理学。《社会学杂志》,84,357-360,(1964)
[38] 麦克劳德,A.J.,托马斯·费米方程的切比雪夫级数解,计算机。物理学。社区。,67, 389-391, (1992) ·Zbl 0875.65075号
[39] Mori,M.,双指数变换的发现及其发展,Publ。京都大学RIMS,41,897-935,(2005)·Zbl 1098.41031号
[40] 莫里,M。;Sugihara,M.,《数值分析中的双指数变换》,J.Compute。申请。数学。,127, 287-296, (2001) ·Zbl 0971.65015号
[41] Nurmuhammad,A。;穆罕默德,M。;莫里,M。;Sugihara,M.,四阶常微分方程边值问题sinc-配置方法中的双指数变换,J.Compute。申请。数学。,182, 32-50, (2001) ·Zbl 1073.65064号
[42] (Olver,F.W.J.;Lozier,D.W.;Boisvert,R.F.;Clark,C.W.,NIST数学函数手册,(2010),剑桥大学出版社,纽约)·Zbl 1198.00002号
[43] Parand,K。;Dehghan,M。;Pirkhedri,A.,求解Thomas-Fermi方程的正弦配置法,J.Compute。申请。数学。,237, 244-252, (2013) ·Zbl 1253.65111号
[44] Parand,K。;Delkhosh,M.,使用有理Chebyshev函数的分数阶精确求解Thomas-Fermi方程,J.Compute。申请。数学。,317, 624-642, (2017) ·Zbl 1362.34039号
[45] Parand,K。;Delkhosh,M.,求解非线性奇异Thomas-Fermi微分方程的新数值解,Bull。贝尔格。数学。Simon Stevin律师,24岁,457-466岁,(2017)·Zbl 1377.65096号
[46] Parand,K。;Ghaderi,A。;Yousefi,H。;Delkhosh,M.,基于有理贝塞尔函数分数阶的非线性Thomas-Fermi方程的新求解方法,Electron。J.差异。Equ.、。,(2016) ·Zbl 1357.34035号
[47] Parand,K。;Mazaheri,P。;Yousefi,H。;Delkhosh,M.,解非线性奇异Thomas-Fermi方程的有理Jacobi函数的分数阶,《欧洲物理学》。J.Plus,132,(2017)
[48] Parand,K。;Shahini,M.,求解Thomas-Fermi方程的有理切比雪夫伪谱方法,物理学。莱特。A、 373210-213(2009)·Zbl 1227.49050号
[49] Parand,K。;Yousefi,H。;Delkhosh,M.,获得托马斯·费尔米方程精确解的一种新的数值技术,《欧洲物理学》。J.Plus,131(2016)
[50] Schwinger,J.,Thomas-Fermi模型:领先修正,Phys。修订版A,22,1827-1832,(1980年)
[51] Schwinger,J.,Thomas-Fermi模型:第二修正,Phys。版本A,242353-2361,(1981)
[52] Seydel,R.,《实用分岔和稳定性分析:从平衡到混沌》(1994),Springer-Verlag-Heidelberg,第428页·Zbl 0806.34028号
[53] 斯莱特,J.C。;Krutter,H.M.,金属的Thomas-Fermi方法,物理学。修订版,47,559-568,(1935)·Zbl 0011.28408号
[54] Sommerfeld,A.,差分渐近积分-gleichung des Thomas-fermischen atmos,Z.Phys。,78, 283-308, (1932) ·JFM 58.1353.03标准
[55] Stenger,F.,基于Whittaker基数或sinc函数的数值方法,SIAM Rev.,23,165-224,(1981)·Zbl 0461.65007号
[56] Stenger,F.,Sinc方法,(1993),纽约施普林格-弗拉格出版社,500页·Zbl 0803.65141号
[57] Stenger,F.,《正弦数值方法概述》,J.Compute。申请。数学。,121, 379-420, (2000) ·Zbl 0964.65010号
[58] Stenger,F.,《sinc数值方法手册》(2010),查普曼和霍尔/CRC博卡拉顿,佛罗里达州,482页
[59] 托马斯,L.H.,《原子场的计算》,数学。程序。外倾角。菲洛斯。《社会学杂志》,第23期,第542-548页,(1927年)·JFM 53.0868.04号
[60] Trefethen,L.N.,Matlab中的光谱方法,(2000),费城工业和应用数学学会·兹比尔0953.68643
[61] Trefethen,L.N.,近似理论和近似实践,(2013),费城工业和应用数学学会(SIAM)·Zbl 1264.41001号
[62] Trefethen,L.N。;Weideman,J.A.C.,指数收敛梯形法则,SIAM Rev.,56,385-458,(2014)·Zbl 1307.65031号
[63] Tuan,P.D。;Elliott,D.,某些函数类的级数展开系数,数学。计算。,26, 213-232, (1972) ·Zbl 0256.42012号
[64] 魏德曼,J.A.C。;Cloot,A.,无限线上演化方程的谱方法和映射,计算。方法应用。机械。工程师,80,467-481,(1990)·Zbl 0732.65095号
[65] 魏德曼,J.A.C。;Trefethen,L.N.,Fejér和Clenshaw-Curtis正交中的扭结现象,数值。数学。,107, 707-727, (2007) ·Zbl 1142.41010号
[66] 新泽西州扎伊采夫。;Matyushkin,I.V。;沙蒙诺夫,D.V.,中心对称原子托马斯·费尔米方程的数值解,俄罗斯微电子。,3, 303-309, (2004)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。