激光烧蚀瑞利——泰勒不稳定性模拟.pdfVIP

维普资讯 D ，14第6 第 10卷 3期 强 激 光 与 粒 子 束 VoJ．1u．No．3 i998正 8月 H1GH P0W ER LASER AND PART1CLE BEAM S Aug ．1998 激光烧蚀瑞利一泰勒不稳定性模拟’ 叶文华 张维岩、陈光南 - — — — ～ — ， — ～ _}乙32t! (北京应用物理与计算数学研究所 ．计算物理实验室，北京8009信箱，l00088 I 摘 要 给 出了激光烧蚀流体不稳定性计算程序 EUL2D 的物理方程 ，介绍了计算 中使 用 的活动阿格和一些技术 问题处理 EUL2D程序 的计算结果与 Takabe公式、FAST2D程序 和LASNEX程序 的结果较好符合。数值计算 日本大阪大学激光烧蚀瑞利一泰勒不稳定性实 验，再现了实验结果 发现了横向电子热传导烧蚀在长波长扰动的非线性瑞利 一泰勒不稳定 性演变 中起重要作用 。 f r 关键词 瑞利一泰勒 稳定性l激光烧蚀 高精度格式 活动阿格 j—L r 中图分类号O242L-———_-l『———一 一 。 I 流体不稳定性是成功实现惯性约束聚变 (ICF)内爆点火 的一大障碍，成为近年 ICF研究 的一个热点_1]对于激光烧蚀加速的直接驱动方式，由于辐照激光的不均匀性和较浅的电 子热传导烧蚀．流体不稳定性尤为严重。它对辐照激光的均匀性提出了非常苛刻的要求 ]。 流体不稳定性计算离不开高精度格式。先前我们 的工作发展 了高精度 FCT(Flux Corrected—Transport)纯流体程序 ，计算 ICF加速度和冲击波情况的流体不稳定性 “一。工作 中．我们研究并解决了电子热传导和高精度FCT的耦合计算问题，研制了激光烧蚀壳靶的流 体不稳定性程序 EUL2D，计算结果与国外激光烧蚀瑞利一泰靶不稳定性 (RTI)的计算结果 和实验结果及与线性增长率的理论公式都较好符合。 1 物 理方程 组 采用流体和热传导分开计算的方法，r为激光人射方 向(纵 向)， 为垂直方向(横向)流 体部分方程为 塑 + + 一0 8r a(pu)+ + O(puv) 丁 。 8r 。 一 g __ + !!!!+丛 !± 一n 8r 塑 +—O(p—eu)+O(pev) c3r 十 户( 。 。 …c3r+塞)一0 式中，p、“、、￡分别为流体密度 、r和z方 向的速度 、单位质量物质内能；P为压力 ；g为加速度。 热传导部分方程为 誓一害喊丁誓)+耋(了1要)+晚 ·国家893惯性约束聚变领域及其青年基金资助项 目 l998年1月20目收到原稿 ，1998年明 26日收到修改稿 叶文华．男 ，1959年明 出生 ，硕士 、副研究员