激光自聚焦和成丝模拟中的并行计算方法.pdfVIP

第 25 卷 第 2 期 计 　　算 　　物 　　理 Vol. 25 ,No. 2 　 2008 年 3 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Mar. , 2008 ( ) [文章编号] 　1001246X 2008 激光自聚焦和成丝模拟中的并行计算方法 陈　军 , 　莫则尧 , 　郑春阳 , 　李 　斌 (北京应用物理与计算数学研究所 , 北京 　100088) [摘 　要] 　已编制完成三维非线性流体力学与激光传播耦合并行模拟程序 , 用以研究高强度激光在稀疏等离子 体中的自聚焦和成丝不稳定性形成机制. 介绍该三维并行程序中所使用的数值方法和并行算法. [ 关键词] 　并行计算 ; 激光成丝不稳定性 ; 自聚焦 ; 区域分解 [ 中图分类号] 　O5523 [文献标识码] 　A 0 　引言 在激光等离子体相互作用研究中,激光成丝是一类严重影响激光能量吸收的不稳定性过程. 激光在等离 子体中传播时 , 由于激光本身的空间不均匀性 ,如果某处等离子体密度低 ,将导致激光在此处的折射率增加 , 激光将向此聚集. 在该区域 ,它通过有质动力排开等离子体 ,进一步降低该区域的等离子体密度 ,从而使激光 进一步聚集. 另一方面 ,如果激光束本身有局部的小尺度热斑 ,它将通过增加碰撞吸收使得当地等离子体膨 胀 ,产生密度凹坑 ,引起激光能量向此聚集. 这个正反馈过程最终导致激光被等离子体分割成丝状结构 ,此 时 ,称为激光成丝现象. 极端情形下 ,整束激光将排开所有等离子体而自聚焦到一个非常小的局部区域 ,即激 光自聚焦. 这种情况下 ,激光将无法沿预定方向传播[1 ,2 ] . 激光自聚焦和成丝不稳定性对激光等离子体相互作 用影响较大 , 如局部激光能量的聚集将引入或加强其它的参量不稳定性 ,影响激光能量在等离子体中的沉 积 ,从而影响对称压缩和内爆稳定性 , 因此需要对它们加以严格控制. 在此情形下 , 需要研制相关的数值模 拟工具来研究激光自聚焦和成丝机制. 另外 ,激光成丝模拟涉及到的计算量大 ,原理性研究也需要上亿个网 格 ,因此必须通过大规模并行来实现模拟. 相关研究方面 ,美国LLNL 实验室的 PF3D 程序[3 ,4 ] 未考虑塌缩现象. 另一程序 ALPS[5 ] , 引入块结构网格 自适应方法 ,但未考虑非局域电子热传导. 国内对激光成丝问题的探索停留在一维和二维问题的理论分析 上[6 ] ,还未见三维大规模并行数值模拟研究该问题. 我们已编制完成三维非线性流体力学与激光传播耦合的并行数值模拟程序 LAP3D ,用以研究高强度激 光在稀疏等离子体区域传播中的激光自聚焦和成丝形成机制. 下面介绍其数值方法和并行算法. 1 　数学物理方程 在三维笛卡尔坐标系上 ,假设等离子体满足准中性条件 ,且激光沿 z 向传播 ,并忽略电子惯性 ,所求解的 方程组如下 : ρ Δ ρ ( ) + ·u = 0 , 1 t S Δ Δ Δ Δ ( ) ( ) ρ ψ ( )