真空球对称引力场中测地线方程的Maple求法.pdfVIP

第29卷第6期 泰 山 学 院 学报 V01．29 NO．6 2007年11月 JOURNAL OF TAISI／AN UNIVERSITY NOV． 聊 真 空 球 对 称 引 d-J场 中 测 地 线 方 程 的 Maple求 法 皂 链 奄 (泰山学院物理与电子科学系，山东 泰安 271021) [摘 要] 用Maple软件求出了真空球对称引力场中粒子的运动学方程，文章中给出的Maple程序可以计 算各种形式的引力场中的Christoffel记号、Riemann曲率张量、Ricci张量和曲率标量以及爱因斯坦张量，借助这 些物理量可以研究引力场中粒子的行为、引力效应和空间的几何性质等． [关键词] Maple软件；张量分析；真空球对称引力场；测地线方程 [中图分类号] O41 [文献标识码] A [文章编号] 1672—259o(2oo7)o6一oo44—04 0 引言 自从爱因斯坦于1915年发表广义相对论以来，张量分析一直在理论物理领域起着无可代替的重要 作用．随着现代工程技术的发展，工程技术人员使用的数学工具逐渐从单变量发展到多变量，从线性理 论发展到非线性理论，而张量方法的广泛使用，就是这种趋势的一个重要方面．应用张量记号，不仅使工 程理论的表达大大简化，而且因为摆脱了坐标系的局限性，使理论的物理意义变得更为清晰．但是，因为 张量运算的规则比较复杂，尤其对于运算量大的问题，若采用手算的方法不仅费时费力而且容易出错， 因此有必要研究现代化的计算软件在张量分析中的应用．本文以广义相对论中的一个典型计算问题为 例说明科学计算软件Maple在张量分析中的应用．Maple是由加拿大Waterloo Maple Inc在1984年开始推 出的一款工程计算软件，目前最高版本为MaplelO．与Mathematica一样，Maple属于计算机代数系统，因 此二者均适合大量的解析推导运算．但是相比较而言Maple的核心部分所占的内存比Mathematica少，其 所有的数学库和附带的包都是用Maple语言写的，源程序可打开来查看和改写，而Mathematica的大部分 都是用c语言编写的，用户无法打开源程序进行检查和修改，所以Maple比Mathematica更胜一筹． 在广义相对论中涉及大量的张量运算，其中一个最典型的计算过程为：根据被研究对象的物理条 件给出某种形式的度规张量(引力场的引力势)，从该度规张量出发，计算Christoffel记号，把Christoffel 记号代入短程线方程，得到测地线方程(即运动微分方程)，最后利用测地线方程方程(组)研究该引力场 中粒子的运动规律．在上述计算过程中，所需要的数学技巧只是偏微分运算(张量的协变微商)，应该说 没有什么大的障碍．麻烦在于计算量太大，度规和Christoffel记号都是4阶2维以上的张量，如果需要求 爱因斯坦场方程还会碰到Riemann曲率张量、Ricci张量和曲率标量，其中Riemann曲率张量是4阶4维 张量(共有256个分量)．如此大量的符号计算若采用手工推导是相当费时，而且不能保证一次计算的准 确性，往往会因为一个小符号的疏忽而前功尽弃．但是若采用用一个短小的Maple程序，如此复杂的符 号运算便可在瞬间完成． [收稿日期]2O07一o7一l8 [基金项目]泰山学院科研资助立项项目(Y95—2—12) [作者简介]宫衍香(1977一)，男，山东泰安人，泰山学院物理与电子科学系讲师，硕士 第6期 宫衍香：真空球对称弓I力场中测地线方程的Maple求法 45 1 计算步骤分析 真空球对称引力场的协变度规张量 为n gn：[1—2kM 1一 (1) 2 g ：r