物理计算的保真与代数动力学算法：Ⅲ辛代数动力学算法与高精度代数动力学算法简介.pdfVIP

物理计算的保真与代数动力学算法： 111．辛代数动力学算法-q高精度代数动力学算法简介 王顺金 张华 四川大学物理学院 理论物理中心成都610064 摘要：本报告包括两项内容：A．介绍辛代数动力学算法设计与计算机实验初步结果：基于 常微分方程的代数动力学的精确解析解，对哈密顿系统，设计出保持局域辛几何结构的辛代 数动力学算法，比较了辛代数动力学算法与常规辛几何算法和龙格．库塔算法的计算机实验 结果。B．简要介绍精度达到机器精度的高精度代数动力学算法的计算机实验的初步结果。 A．辛代数动力学算法设计与计算机实验初步结果 I-常微分方程的代数动力学精确解 n个变量、p个参数的常微分方程可写为 X．=f(口。(f)，X，)，j=1,2，．刀；∥=1,2。P(1a) 口，(，)构成p-维参数空间，口，=巳·自治系统· X，(，)为n-维相空间的曲线，向量场F规定了方程的局域微分结构． 把常微分方程提升为偏微分方程： 型：讶， 三：夕F(x)旦 (1b) 8l 急”。ax? 在收敛半径内，用代数动力学方法得(1b)式的精确解 牙(五，，)一∑≯(xo)Roe反釉‘矗 (2) 月IO，l： 分量形式为 X，(A-o，f)=∑≯(Xo)Xo。=g“㈨‘丘。 (3) ^-O，I： 动力系统时间演化的局域微分性质体现在时间平移算子．生成元三上， 动力系统时间演化的整体积分性质体现在时间演化算子．群元c，(，)上． 基于精确解的近似算法有： (1)朴素代数动力学算法(显式)：以(三△r)：至号手彦(4) 舢y： (2)辛代数动力学算法(保辛、隐式)： (3)高精度代数动力学算法(达到机器精度、隐式)： 对于(1)朴素代数动力学算法，请参考文献[1，2】．下面介绍(2)辛代数动力学算 法设计与计算机实验结果．对于(3)高精度代数动力学算法，只介绍计算机实验结果，以 后有机会再介绍算法设计． 11．低阶辛代数动力学算法j以 与显式代数动力学算法不同，辛代数动力学算法是用时间平移微分算子三去设计一定 阶的辛代数动力学算法，仍然保持代数动力学算法的特点．即保持时间平移算子三的完整性 和对于精确解一定的精度。 2．1设计辛代数动力学算法的原则 设计辛代数动力学算法的基本原则是，以时间平移算子上为基本算子去构造近似的 33 时间演化算子sU^，(△f)，使它在(△，)的N阶精度下逼近精确解，而同时保持辛几何结构． 为此．它必须遵从下述具体原则： (1)L的完整性：sUⅣ(△，)必须是L的函数以保持L的完整性： sU^，(△，)=sUⅣ【￡址】 Oa) (2)精度可控制：sUⅣ(△，)必须是严格的时间演化算子U(△．，)的N阶近似 (5b) sD^，(△f)=U_Ⅳ[Latl+so[(Lat)㈧】 (3)保持局域辛结构(对常数、标准、“平直”辛度规而言，，(石)为局域度规)： Od) M7J(X)M=／(Xo) 对比‘·平直’’保辛： M7JM=‘， M：坚 拟o