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

维普资讯

D, 1第 1 0卷

4第 63期

强

激

光

与

粒

子

束

VoJ No.3 u. .1 A ug . 1 98 9

i 98正 8月 9

H1 GH P 0W ER LAS ER AND ART1 P CLE EAM S B

激光烧蚀瑞利一泰勒不稳定性模拟’叶文华张维岩、陈光南 -———~—

,

—

~

_ 2!} 3t乙I

(京应用物理与计算数学研究所 .计算物理实验室，京 8 0信箱，l 0 8北北 09 008

摘要给出了激光烧蚀流体不稳定性计算程序 E 2的物理方程，绍了计算中使 UL D介用的活动阿格和一些技术问题处理 E 2程序的计算结果与 Ta a e式、 AS D程序 UL D kb公 F T2和 L NE AS X程序的结果较好符合。数值计算日本大阪大学激光烧蚀瑞利一泰勒不稳定性实

验，现了实验结果再发现了横向电子热传导烧蚀在长波长扰动的非线性瑞利一泰勒不稳定性演变中起重要作用。

f

r

关键词瑞利一稳定泰勒性l激光烧蚀高精度格式活动阿格

jL— r

中分号 O 2-— _—— 图类 2——『—一 4L - l I

一

。

流体不稳定性是成功实现惯性约束聚变 (C内爆点火的一大障碍，成为近年 I F研究 I F) C

的一个热点_]对于激光烧蚀加速的直接驱动方式，由于辐照激光的不均匀性和较浅的电 1 子热传导烧蚀 .流体不稳定性尤为严重。它对辐照激光的均匀性提出了非常苛刻的要求。]流体不稳定性计算离不开高精度格式。前我们的工作发展了高精度 F T( lx先 C F u C retd orce—T a s o t纯流体程序，算 I F加速度和冲击波情况的流体不稳定性“。 rn p r)计 C一工作

中 .们研究并解决了电子热传导和高精度 F T的耦合计算问题，制了激光烧蚀壳靶的流我 C研体不稳定性程序 E 2 UL D,计算结果与国外激光烧蚀瑞利一泰靶不稳定性 ( TI的计算结果 R )和实验结果及与线性增长率的理论公式都较好符合。1物理方程组

采用流体和热传导分开计算的方法， r为激光人射方向 (向)为

垂直方向 (向)流纵， 横 体部分方程为

塑+8 r

+

一0

a p )+ (u

+ O p v) (u8 r。

丁

。

_ _

一 g

+ ! !! !+丛！ 一 n ±8 r

塑 +—pu+ Opv十 O c ) (e) (r e— 3

。

。

户( …+ ) 0 c塞一 3 r

式中，、、￡ p“、分别为流体密度、和 z向的速度、位质量物质内能； r方单 P为压力； g为加速度。

热传导部分方程为

誓一喊 )耋了要+ 害丁誓+ (1 )晚 国家83性约束聚变领域及其青年基金资助项目 9惯

l 9年 1 2目收到原稿，9 8明 2日收到修改稿 8月 0 9 19年 6叶文华 .男，15年明出生，硕士、副研究员 99

维普资讯

44 0

强激光与

粒子柬

第 1 0卷

一

K y, t

口 r

华一 一KE一 ( 一， )+ E 。

T为温度；为比热容； o 经典电子热传导系数； 吸收激光能量；、,别为人 kT“为 E为 分射和反射激光功率密度；K为激光逆轫致吸收系数； 为激光共振吸收能量。 E两个方向实行单独限流F…一 C F t r.

*

.

, . 一 ( 1+ l . F 1 F . )/

F一 k, —o; F D一o誓 F k t¨一 r T TF,式中，是 r 和方向的电子热传导限流因子； r FD和 F 别是 r和 z方向的扩散能流和 . .分自由能流。

流体部分采用高精度 F T算法， C热传导部分采用局部一维全隐格式。由热传导流的连续

性条件建立热传导差分格式，网格边界温度取相邻两网格温度的平均值。工作主要研究高本精度 F T流体与热传导耦合计算问题，以便今后将程序推广到三维方程和光子多群扩散方 C

程 .所以物理过程尽量简化。采用理想气体状态方程，全电离和单维近似。

2活动网格流体不稳定性非线性发展产生大的剪切流，氏计算方法会遇到网格相交的严重困难。拉 因此我们采用活动网格欧拉计算方法。不稳定性计算要求较高的网格分辨率，须采用高精必度欧拉算法。网格划分的不同对不稳定性计算结果有严重影响。

我们用密度峰值处速度或靶心速度追踪烧蚀面，移网格板块。蚀面附近均匀密分网格，以减少对流误差，边网格逐平烧两渐放大，放大因子不超过 1 1纵向计算边界尽量放得远些。光吸收区和电子热传导区要但 .。激保持合适的网格宽度。 活动网格不能直接计算激光烧蚀固体靶。先用拉氏程序低功率算 2 0 0 p . 0￣3 0 s产生烧蚀面密度梯度 .启动活动网格程序进行计算。在拉氏计算区域外无质量的活动欧拉网格内，再 填充密度 1 g c。 0/m

3一维计算结果比较数值研究表明：用内能格式能得到较好的计算结果。用总能量守恒格式则在电子热采采传导前沿和其它密度问断处内能和动能匹配不好，得计算结果产生振荡，温度出负，计算使难以继续下去。接触问断附近密度要适当光滑化。算中要确保烧蚀面落在密分网格区内，否计则烧蚀面附近计算的物理量可能出现振荡。 3 1与 F s 2程序比较 . A TD

美国海军实验室的 F T2程序，采用活动网格欧拉方法和高精度 F T流体算法， 0 AS D C 8年代期间经过多次实验结果的校验]。计算条件为初始 C靶厚度 1 5 m, H 0 t扎一0 2 4 m,峰￣ . 6#值激光功率 I一3 1/m .激光功率密度上升前沿为 I lsn ( 4 s/i‘,值后功 × 0W c— pi‘t n )sn1峰/率密度不变。

维普资讯

第 3期

叶文华等：激光烧蚀瑞利一泰靶不稳定性模拟

O

10 0

￡0

:l~0●

8 0●

g 一。

目6 0

呈一0 3一

套

4 0

Fi . Te o a r f eo a g t㈣ g 1 mp r lp o i ft r e l

s a,p a e st ( - ) e k d n iy b) C

v l ct c n o ii n ( )a e k d n iy,p a r s u e[ ) eo iy( )a d p sto d tp a e s t e k p e s r e

图 1靶质量 ()峰值密度 ( )该址速度 ( ) a、 b及 c和位置 ( )峰值压力 ( )时间的分布 d、 e随

E 2程序的上述计算条件和网格划分与文献 E3 F T2计算相一致。图 1比较 UL

D 6的 AS D了靶质量 ( )峰值密度 ( )该处速度 () a、 b及 c和位置 ( )峰值压力 ( ) d、 d随时间的分布+两程序计算结果较好符合。6 5 7 5 s时间段拟合的靶加速度 ( 0sm/和质量烧蚀速率 ( g .~ .n 11 c s) 1 !/ O c s ): AS 2程序分别为 4 0和 4 4 E 2 mz是 F T D . ., UL D程序分别为 4 1和 4 5 . .。两程序计算结果也较好符合靶加速度取峰值密度处的数值。图 1细小的差异可能来自计算的不同具体细

节，如电子热传导限流因子的取值等，文献[ 3能给出这些具体细节。 6未3 2与 I N X程序比较 . AS E

美国利弗奠尔实验室的 LAS NEX程序，广泛用于 I F的数值计算。计算条件为：0 m C 2v C靶+扎= 0 2 v H . 5 m,J一 2 O W/m,光功率密度线性上升 9 s到峰值+之后保持不 ×1“ c 激 n变。L NE程序计算无辐射。图 2表明： D程序和 L NE AS X EUI 2 AS X程序计算的加速度和烧蚀速度随时间的分布基本一致 L NE AS X程序采用拉氏计算方法，烧蚀速度和加速度取峰值密度处的数值。比较图 2的结果，可见 L NE计算的烧蚀速度和加速度比 E D的跳动 AS X UI2

一

…

I^ I I I一

…

J ,一

^

0F g 2 a Te o a r f e o b a i e v l c t i. mp r lp o i fa t tv eo iy;b l .Te o a r f e o c ee a i n mp r lp o i fa c lr t l o

2 a烧蚀速度 …… 此处隐藏：3325字，全部文档内容请下载后查看。喜欢就下载吧 ……