在间接驱动惯性约束聚变中,靶丸X光辐照的均匀性是关键问题之一。要细致分析黑腔靶中X光的辐射传输过程,数值模拟是重要的研究手段。由于黑腔靶内部光性薄辐射扩散近似不成立,而在腔壁(Au)内部存在光性厚区域,采用输运数值模拟不仅要求输运方程的数值离散格式具有离散扩散极限[1],而且输运方程的源迭代算法很难收敛[2],将大大提高计算复杂性和难度。一个自然的选择是采用输运和扩散模拟各自的优势进行多尺度建模,在光性薄区域求解多群辐射输运方程,在光学厚区域对多群辐射扩散方程,以取得满意的计算精度和计算效率。
本文基于JASMIN 框架的 “联邦计算”等支撑技术[3],在ICF腔靶耦合数值模拟二维总体程序——LARED集成程序[4]中实现了对辐射传输过程的多群输运与多群扩散多尺度耦合数值模拟能力。64群输运与扩散耦合数值模拟神光II实验标准黑腔完整物理过程,得到计算结果与全输运数值模拟结果吻合很好,较全区扩散计算更能反应黑腔辐射场的不均匀性。给定相匹配计算条件,耦合计算与全区扩散计算的CPU时间开销接近,约为全区输运计算CPU时间的2/3.
致谢:本工作受国家自然科学基金委项目91130002资助
参考文献:
[1]M. L. Adams, P. F. Nowak, Asymptotic analysis of a computational method for time- and frequency- dependent radiative transfer[J], Journal of Computational Physics, 1998, 146, 366-403.
[2] J. E. Morel, T.-Y. Brian Yang and J. S. Warsa, Linear Multifrequency-grey Acceleration Recast for
Preconditioned Krylov iterations, J. Comp. Phys. 227,244-263(2007)
[3] 莫则尧、张爱清主编,JASMIN框架1.8版用户指南,北京应用物理与计算数学研究所,2010.2.
[4] Pei W B. The construction of simulation algorithms for laser fusion[J]. Commun Comput Phys,2007,2:255-270.
通讯地址:北京应用物理与计算数学研究所,北京,100094
TEL : 13681009417
内容所属范围:数值计算及其应用