二维非平衡辐射流体力学程序研制与应用研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
X光激光(XRL)以其独特的优点使得它在物理学、生物学、化学、材料学以及惯性约束聚变(ICF)等领域有着巨大的应用前景。经过20多年的发展,XRL输出已经能够达到饱和,现在各实验室纷纷开展应用演示研究。在XRL的深入研究过程中,我们发现1D研究越来越不能适应XRL理论研究的需要;XRL诊断ICF等离子体研究也期待一个2D程序。只有通过2D研究,我们才能更清楚地认识XRL复杂的场图结构;才能更好地完成XRL单靶、多靶的优化设计,以得到更强的XRL输出;才能在XRL密度诊断中更好完成待测等离子体的预先设计;才能更好的解释密度诊断实验结果,达到校验理论模型的目的。为适应这种需要,在博士论文期间,我把原有1D程序JB19发展成了非平衡辐射流体力学二维程序XRL2D,并用它来进行了系列物理问题的研究。
     XRL2D程序采用单流体模型,电子离子温度脱离,电子离子能量输运采用限流扩散近似,辐射输运采用多群限流扩散近似,用平均原子模型计算原子过程。程序采用分裂格式,把扩散过程与电离、加热等局域过程分开求解,2D扩散计算采用Kershaw的九点差分格式。文章对物理过程和计算方法进行了详细介绍。
     (1)计算方法研究。为克服Lagrange方法的网格易变形问题,借鉴三角形抗变形能力强的特点,本研究创新性地提出了改进IGA格式,有效地抑制了网格变形。但是该方法对物理上真实存在的流体卷曲、混合等也显得无能为力,因此程序采用了ALE方法。为考察ALE方法网格重分、物理量重映导致的数值扩散对计算精度的影响,进行了一些计算比较,结果显示,虽然存在一定程度的抹平现象,但是网格重分带来的数值耗散并不大,对物理图象并没有显著的影响。文中还给出了XRL2D同LASNEX程序的对比计算,以考察程序的可靠性;论文还探讨了用双向一维方法替代九点差分格式的可能性,因为前者能够大大降低计算量。
     (2)物理问题研究。首先对神光Ⅱ上的准稳态Ni-like Ag XRL实验进行了模拟,计算结果很好的解释了XRL实验场图出现弯月形分布的原因。其次,在XRL等离子体密度诊断方面,对近期神光Ⅱ上的三个实验进行了模拟对比研究;用XRL M-Z干涉仪对点聚焦CH等离子体电子密度的诊断;用XRL对倍频光烧蚀CH薄片的背景照相;用XRL对射流的阴影成像。模拟结果同实验结果有些定量符合,有些定性符合,也有存在明显差异的地方。论文还给出了在神光Ⅱ上近期将进行的XRL诊断应用实验的理论设计,用激光驱动半开腔Au靶,诊断腔口喷出等离子体的电子密度,模拟显示了XRL2D程序的模拟能力。
     本论文的完成得到了国家高技术863计划项目资助课题(803-804-7-1,803-804-1-10-1)
For its unique advantages,x-ray laser has significant prospects of applications in physics,biology,chemistry,materials science and ICF fields.After 20 years' development,the saturated output of XRL has been obtained,and now all laboratories all over the world have beginning the research of its application demonstration.In the study of XRL,we more and more find that the 1D study can't catch up with the step of XRL theoretical deep study,and the study of diagnosing of ICF plasma by XRL also desires a 2D code.Only through 2D study,we can understand complicated filed image structures of XRL,and carefully optimize XRL experiments of single- and multi-targets to get higher intensity of XRL,we can more professionally design plasma to be probed in XRL density diagnosis,and explain the experimental results of electron density diagnosis,to fulfill the goal of testing theoretical models.In order to adapt those requirements,in the period of my PhD,I have developed the 1D code JB19 to a 2D non-equilibrium radiation hydrodynamic code,named XRL2D,and have completed some physical studies using this code.
     Single hydrodynamic model is used in XRL2D code,electron and ion have different temperatures.Flux-limit heat-conduction approximation and flux-limit multi-group diffusion approximation are adopted,and average atomic model is used in code.Split scheme is used,the diffusion processes is departed from local processes such as ionization and heating.The 9-piont differential scheme suggested by Kershaw is used in the calculation of 2D diffusions.The physical processes and calculating methods are detailed in this thesis.
     (1)The studies of calculating methods.In order to overcome the trouble of mesh distorting in Lagrange calculation,taking the advantage of strong anti-distorting for triangle meshes,a modified IGA scheme is proposed innovatively in thesis,and the attempts show that the mesh distorting are restrained effectively.But this method can also do nothing while real physical convolution and mixture of fluid existing,so the ALE method is adopted in code.In order to survey the numerical diffusion induced by mesh rezoning and variables remap in ALE method,some calculations had been done,and comparisons show that the smoothing does exist,but the numerical dissipation is not serious,it do no harm to physical picture.Furthermore,thesis includes the comparisons of simulation results between XRL2D and LASNEX code for same models to survey the reliability of code.And the possibility of displacing 9-point differential scheme by bidirectional 1D method is also discussed in thesis,for the latter can greatly decrease the calculation amount.
     (2)The studies of physical problem.Firstly,the experiment of quasi-steady-state Ni-like Ag XRL on ShenGuang-Ⅱfacility is simulated,the calculated results can well explain why the experimental field image of XRL show a crescent-shaped structure.Secondly,three recent experiments of electron density diagnosis by XRL on ShenGuang-Ⅱalso are simulated,those experiments are;electron density diagnosis of spot focused CH plasma by XRL M-Z interferometer,measurement of CH film driven by 2w laser using a XRL backlighter and measurement of fluid jet by XRL backlighter.Some simulation results agree with experimental results quantitatively, some results agree experiment only qualitatively,and some disagreements also exist. A theoretical design of electron density diagnosis application experiment by XRL is also given out in thesis,in which the plasma jetted out of Au half-hohlraum driven by laser will be measured.This simulation shows the calculating ability of XRL2D code.
     This thesis has been supported by the National High Technology 863 project of China(Grant No.863-804-7-1 and 803-804-1-10-1).
引文
[1]彭惠民,王世绩,X射线激光,北京:国防工业出版社,1997
    [2]Jr.C.Robinson,Advance Made on High-Energy Laser,AWST,114(8),25-27,1981
    [3]S.Suckewer,C.H.Skinner,H.Milchberg,et al.Phys.Rev.Lett.,55,1753-1756,1985
    [4]D.L.Matthews,P.L.Hagelstein,M.D.Rosen,et al.phys.Rev.Lett.,54,110-113,1985
    [5]S.J.Wang,Y.Gu,G.L.Zhou,et al,Chin.Phys.Lett.,8(12),618-621,1991
    [6]A.Carillon,H.Z.Chen,P.Dhez,et al,Phys.Rev.Lett.,68(19),2917-2920,1992
    [7]J.A.Koch,B.J.MacGowan,L.B.DaSilva,et al,Phys.Rev.Lett.,68(22),3291-3294,1992
    [8]L.B.DaSilva,B.J.MacGowan,S.Mrowka,et al,Opt.Lett.,18(14),1174-1176,1993
    [9]J.J.Rocca,J.L.A.Chilla,D.P.Clark,et al,Phys.Rev.Lett.,77(8),1476-1479,1996
    [10]J.Zhang,P.J.Warwick,E.Wolfrum,et al,Phys.Rev.A,54(6),R4653-4656,1996
    [11]M.P.Kalachnikov,P.V.Nickles,M.Schnurer,et al,Phys.Rev.A,57(6),R4778-4783,1998
    [12]J.Zhang,A.G.Macphee,J.Nilsen,et al,Phys.Rev.Lett.,78(20),3856-3859,1997
    [13]R.Tommasini,F.Loewenthal,J.E.Palmer,et al,J.Opt.Soc.Am.B,16(10),1664-1667,1999
    [14]J.Y.Lin,G.J.Tallents,J.Zhang,et al,Opt.Commun.,158,55-60,1998
    [15]J.Zhang,A.G.Macphee,J.Lin,et al,Science,276,1097-1100,1997
    [16]T.N.Lee,E.A.Mclean and R.C.Elton,Phys.Rev.Lett.,59(11),1185-1188,1987
    [17]J.Nilsen and J.C.Moreno,Phys.Rev.Lett,74(17),3376-3379,1995
    [18]H.Daido,S.Ninomiya,T.Imani,et al,Opt.Lett.,21(13),958-960,1996
    [19]S.J.Wang,Y.Gu,G.L.Zhou,et al,J.Opt.Soc.Am.B,9(9),360-368,1992
    [20]王世绩,顾援,周关林等,强激光与粒子束,4(2),165-174,1992
    [21]王世绩,顾援,周关林等,强激光与粒子束,5(4),557-563,1993
    [22]S.J.Wang,Y.Gu,C.S.Mao,et al,Chin.J.Lasers B,2,481-484,1993
    [23]S.J.Wang,Y.Gu,G.L.Zhou,et al,Chin.J.Lasers B,3,507-510,1994
    [24]R.Kodama,D.Neely,Y.Kato,et al,Phys.Rev.Lett.,73,3215-3218,1994
    [25]张国平,盛家田,彭惠民等,强激光与粒子束,2(3),289-310,1990
    [26]张国平,盛家田,杨明伦等,强激光与粒子束,4(4),521-526,1992
    [27]张国平,张覃鑫,吴建周,强激光与粒子束,10,352-357,1998
    [28]张国平,张覃鑫,郑无敌,强激光与粒子束,16(1),35-39,2004
    [29]张国平,张覃鑫,郑无敌,强激光与粒子束,16(2),171-176,2004
    [30]王琛,顾援,王世绩等,强激光与粒子束,14(5),381-384,2002
    [31]王琛,王伟,吴江等,物理学报,53,3752-3755,2004
    [32]P.V.Nickles,V.N.Shlyaptsev,M.Kalachnikov,et al,Phys.Rev.Lett.,78,2748,1997
    [33]J.Dunn,A.L.Osterheld,R.Shepherd,et al,Phys.Rev.Lett.,80,2825,1998
    [34]J.Dunn,Y.Li,A.L.Osterheld,et al,Phys.Rev.Lett.,84,4834,2000
    [35]H.Fiedorowicz,A.Bartnik,R.Jarocki,et al,in "X-ray laser 2002",AIP Conf.Proc.,641,21-25,2002
    [36]K.A.Janulewicz,A.Lucianetti,G.Priebe,et al,in "X-ray laser 2002",AIP Conf.Proc.,641,26-30,2002
    [37]Y.Kato,H.Daido,K.Nagashima,et al,in "X-ray laser 2002",AIP Conf.Proc.,641,31-39,2002
    [38]T.Ozaki,H.Kuroda,R.A.Ganeev,et al,in "X-ray laser 2002",AIP Conf.Proa.,641,65-68,2002
    [39]D.Ros,A.Klisnick,D.Joyeux,et al,in "X-ray laser 2002",AIP Conf.Proa.,641,69-76,2002
    [40]S.Sebban,L.M.Upcraft,P.Balcou,et al,J.Opt.Soc.Am.B.,20,195,2003
    [41]Y.Nagata,K.Midorikawa,S.kubodera,et al,Phys.Rev.Lett.,71,3774,1993
    [42]C.Reinhardt,M.Born,A.Egbert,et al,J.Phys.IV(France),11,201,2001
    [43]H.Fiedorowicz,A.Bartnik,M.Szczurek,et al,J.Phys.IV(France),11,197,2001
    [44]R.Keenan,J.Dunn,S.Moon,et al.,Proa SPIE,5197,213,2003
    [45]V.N.Shlyaptsev,J.Dunn,S.Moon,et al.Proc.SPIE,5197,221,2003
    [46]G.P.Zhang,T.X.Zhang and W.D.Zheng,in "X-ray laser 2004",AIP Conf.Proc.,186,373-376,2004
    [47]D.L.Matthews,in "X-ray laser 1996",AIP Conf.Proc.,151,32-39,1996
    [48]R.A.London,X-ray Laser,chapter 7,Sam Diego,Academic Press,1990
    [49]L.B.Da Silva,T.W.Barbee,Jr.R.Cauble,et al.,Phy.Rev.Lett.,74,3991-3994,1995
    [50]A.S.Wan,T.W.Barbee,R.Cauble,et al.,Phy.Rev.E.,55,6293-6296,1997
    [51]王琛,王伟,孙今人等,物理学报,54,0202-0205,2002
    [52]W.D.Zheng,G.P.Zhang,C.Wang,et.al.,in "X-ray laser 2004",AIP conf.proc.,186,575-578,2004
    [53]J.Filevich,K.Kanizay,M.C.Marconi,et al.,Opt.Lett.,25,356-358,1999
    [54]J.Filevich,J.J.Rocca,E.Jankowska,et al.,Phy.Rev.E,67,056409-1-056409-6,2003
    [55]J.J.Rocca,E.C.Hammarsten,E.Jankowska,et al.,Phys.plasmas,10,2031-2038,2003
    [56]R.F.Smith,J.Dunn,J.Nilsen,et al.,Phy.Rev.Lett.,89,065004-1-4,2002
    [57]R.F.Smith,J.Dunn,J.Nilsen,et al.,J.Opt.Soc.Am.B,20,254-259,2003
    [58]J.Filevich,J.J.Rocca,M.C.Marconi,et al.,Phy.Rev.Lett.,94,035005,2005
    [59]D.H.Kalantar,M.H.Key,L.B.DaSilva,et al.,Rhy.Rev.Lett.,76,3574-3577,1995
    [60]E.Wolfrum,J.S.Wark,J.Zhang,et al.,X-ray lasers Conf.,Kyoto,Japan,159,657,1998
    [61]J.Svatos,D.Joyeux,D.phalippou,et al.,Opt.Lett.,18,1367-1369,1993
    [62]Ph.Zeitoun,M.Boussoukaya,A.Zeitoun-Fakiris,et al.,in "X-ray laser 1998",AIP Conf.Proc.,159,669,1998
    [63]Ph.Zeitoun,F.Albert,P.Jaegle,et al,Nucl.Instr.Meth.in Phys.Res.A.,416,189-191,1998
    [64]C.H.Moreno,M.C.Marconi,J.L.A.Chilla,et al,Phys.Rev.E,60,911-917,1999
    [65]J.J.Rocca,C.H.Moreno,M.C.Marconi,et al..Opt.Lett.,24,420-422,1998
    [66]王琛,顾援,傅思祖等,物理学报,51,847-851,2002
    [67]L.B.Da Silva,T.W.Barbee,Jr.R.Cauble,et al.,Phy.Rev.Lett.,74,3991-3994,1995
    [68]C.D.Decker,R.A.London,J.A.Harte,et al.,Phys.Rev.E,57,5952-5957,1998
    [69]L.B.DaSilva,J.E.Trebes,R.Balhorn,et al,Science,258,269-271,1992
    [70]D.T.Attwood,E.Anderson,G.Denbeaux,et al,in "X-ray lasers 2002" AIP Conf.Proc.,641,461-468,2002
    [71]H.A.Baldis,J.Dunn,M.E.Foord,et al.,Rev.Sci.Instr.,73(12),4223-4228,2002
    [72]D.Riley,R.Keenan,S.J.Topping,et al.,IEEE Transactions On Plasma Science,31(5),1016-1022,2003
    [73]J.Nilsen,J.C.Moreno,L.B.DaSilva,et al.,Phy.Rev.A,55,827-830,1996
    [74]J.Nilsen,J.Zhang,A.G.MacPhee,et al.,Phy.Rev.A,56,3161-3165,1997
    [75]Ch.Siegel,M.Braud,J.E.Balmer,et al.,Optics Communications,210,305 - 312,2002
    [76]王琛,论文题名,X光激光实验研究[博士毕业论文],保存地,中国工程物理研究院研 究生部,(博士毕业论文),42-43页,2004
    [77]M.D.Rosen,D.W.Phillion,V.C.Rupet,et al.,Phys.Fluids,22(10)2020-2031,1979
    [78]W.C.Mead,E.M.Campbell,K.G.Estabrook,et al.,Phys.Fluids,26(8),2316-2331,1983
    [79]D.C.Slater,G.E.Busch,G.Charatis,et al.,Phy.Rev.Lett,46(18),1199-1203,1981
    [80]R.P.Drake,R.L.Kauffman and B.F.Lasinsk,Phys.Fluids B,3(12),3477-3484,1991
    [81]W.C.Mead and E.K.Stover,Phys.Rev.A,38(10),5275~5288,1988
    [82]W.C.Mead and E.M.Campbell,Phys.Rev.Lett,47(18),1289-1292,1981
    [83]S.R.Goldman and W.C.mead,Nuclear Fusion,26(6),813-826,1986
    [84]K.Eidmann,F.Amiranoff,R.Fedosejevs,et al.,Phys.Rev.A,30(5),2568-2589,1984
    [85]C.G.Labaune,E.Fabre,C.E.Max,et al.,Phys.Rev.Lett,48(15),1018-1021,1981
    [86]B.K.F.Young,R.E.Stewart and C.J.Cerjan,Phy.Rev.Lett,61(25),2851-2854,1988
    [87]J.A.Tarvin,W.B.Fechner,J.T.Larsen,et al.,Phys.Rev.Lett,51(15),1355-1358,1983
    [88]B.Yaakobi,J.Delettrez,L.M.Goldman,et al.,Phys.Fluids,27(2),516-526,1983
    [89]P.D.Goldstone,S.R.Goldman,W.C.Mead,et al.,Phys.Rev.Lett,59(1),56-59,1987
    [90]T.Boehly,K.A.Tanaka,T.Mochizuki,et al.,J.Appl.Phys,60(11),3840-3844,1986
    [91]D.Ress,L.J.Suter,E.F.Babl,et al.,Phys.Fluids.B,2(10),2448-2455,1990
    [92]E.F.Gabl,B.H.Farlor,G.E.Busch,et al.,Phys.Fluids.B,2(10),2437-2447,1990
    [93]W.B.Fechner,C.L.Shepard,G.E.Busch,et al.,Phys.Fluids,27(6),1552-1560,1984
    [94]R.Kodama,K.Okada,N.Ikeda,et al.,J.Appl.Phys,59(9),3050-3052,1986
    [95]T.Mochizuki,T.Yabe,K.Okada,et al.,Phys.Rev.A,33(1),525-537,1985
    [96]L.Spitzer,完全电离气体的物理学(中译本),科学出版社,1959
    [97]K.A.Brueckner and S.Jorna,Reviews of Modern Physics,46,235,1974
    [98]常铁强,张均,张家泰等,激光等离子体相互作用与激光聚变,湖南科学技术出版社,1991
    [99]李世昌,高温辐射物理与量子理论,国防工业出版社,1992
    [100]段庆生,常铁强,张维岩等,计算物理,19,57-61,2002
    [101]朱少平和古培俊,强激光与粒子束,11,687-691,1999
    [102]王建国,邹宇,常铁强,物理学报,46(11),2146-2153,1997
    [103]张家泰,激光等离子体相互作用,河南科学技术出版社,1999
    [104]汤文辉,张若棋,物态方程理论及计算概论,国防科技大学出版社,1999
    [105]戴维逊,中子迁移理论,科学出版社,1961
    [106]李德元,徐国荣,水鸿寿等,二维非定常流体力学数值方法,科学出版社,1987
    [107]D.S.Kershaw,J.Comp.Phys.,39,375-395,1981
    [108]段庆生和袁国兴,计算物理,14(2),148-154,1997
    [109]D.J.Benson,Journal of Computational Physics,100,143-162,1992
    [110]张国平,张覃鑫和郑无敌,强激光和粒子束,16,1375,2004
    [111]张国平,张覃鑫和郑无敌,强激光和粒子束,17,1817,2005
    [112]张国平,张覃鑫和郑无敌,强激光和粒子束,18,210,2006
    [113]江少恩,刘忠礼,唐道源等,光学精密工程,8,181-184,2000

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700