辐射不透明度计算中含温有界原子结构理论的不确定度量化.pdfVIP

第32卷 第3期 计 算 物 理 Vo1．32．NO．3 2015年 5月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS May，2015 文章编号 ：1001．246X(2015)03-0264．13 辐射不透明度计算中含温有界原子 结构理论的不确定度量化 孟续军， 马智博 ， 王瑞利， 朱希睿 (北京应用物理与计算数学研究所 ，北京 100094) 摘 要 ：以交换势的泛 函表达式为核心 ，分析产生原子结构理论不确定度的主要原因，阐述辐射不透 明度计算 中含 温有界原子结构理论不确定度量化的具体步骤，给出能级、矩阵元的不确定度量化计算公式，以及判定理论模型优 劣的判据．以汞、金、铁元素为例 ，通过数值计算，根据对 比分析，验证该不确定度量化方法的可行性． 关键词 ：原子结构 ；不确定度 ；辐射不透明度 ；自洽场理论 中图分类号 ：0213．2；0562．4 文献标志码 ：A 0 引言 物质辐射不透明度是研究辐射输运问题的关键参数 ，也是辐射流体力学数值模拟的必要输入条件 ，往往 决定着数值模拟的成败．在不能进行系统级综合性试验的条件下，辐射不透明度参数 的不确定度也成为可靠 性认证所必需的输入信息．由于基础实验条件的限制 ，在很大温度密度范围内，辐射不透明度参数只能通过 理论计算来获得 ，且无法得到实验验证．为此 ，国际上发展了多种理论模型．这些模型计算框架相 同，物理建 模时选择的具体公式有所不同．尽管这些具体公式可能近似程度不 同、考虑的物理机制不同也不一定完整 ， 但经过了一定的 “公理”性的检验和部分基础实验的检验 ，被证明具有一定的可靠性． 辐射不透明度理论就是用这些 “可靠”公式建立的．然而，用 “可靠”的公式组成的理论系统却不一定完 全可靠．在同一个环节上 ，公式的选取就能反映出人们对具体问题认知程度．辐射不透明度拥有很多理论模 型，恰好反映出人们在辐射不透明度 问题上的认知度存在差异． 我们要研究的终极问题是辐射不透明度物理建模的不一致性所导致的不确定度 ，这里称之为辐射不透 明度理论的不确定度．它源 自建模时所用的假设条件与计算公式的不 同，属于认知不确定度 ，需要通过共有 的特征参量的分散程度进行描述． 辐射过程 由四个基本物理过程结合热力学统计体系组成，计算非常复杂 ，需要含温有界原子结构模型为 其提供大量最基本的微观参数 一能级 、波函数 、电子 占据数等．辐射不透明度理论的不确定度不仅来 自辐射 过程 自身的不确定性 ，还要受原子结构传递过来的不确定度的影响．就解决辐射不透明度理论的不确定度问 题而言 ，可将其分为原子结构部分的量化和辐射过程的量化两部分． 目前，国内外对辐射不透明度及原子结构参数不确定度问题的研究主要在理论对算 、实验对 比、误差分 析的层次．在 SandiaNationalLaboratory的工作报告…及 国内的一些文献 中，对物理建模和爆轰力学方面的 不确定度量化 。 ，给出了指导原则和技术框架 ，很有借鉴意义．但还没有应用到辐射不透明度理论体系，也 没有发现从剖析该理论体系的角度开展的研究． 对原子结构模型进行不确定度量化 ，不仅能够为辐射流场计算的不确定度传递分析提供输入信息，也能 为复杂系统的数值模拟校准提供合理的模型参数可调范围，具有重要的工程应用价值．所 以，辐射不透明度 理论的不确定度研究就以原子结构的不确定度量化方法研究进行开篇． 收稿 日期 ：2014—06—04；修 回日期 ：2014—09—19 基金项 目：中国工程物理研究院科学技术发展基金(2012B0102010)及国家 自然科学基金 资助项 目 作者简介 ：孟续军 (1965一)，男，博士 ，研究员 ，从事原子物理过程及辐射特性研究，E-mail：mengxujun@sina．corn 通讯作者 ：朱希睿(1978一)，男，硕士，副研究员 ，从事原子微观过程及状态方程研究，E-mail：zxr1978@sina．corn 第 3期 孟续军等：辐射不透明度计算中含温有界原子结构理论的不确定度量化