一种用于冷冻靶温度场预测的定向红外光-热耦合模拟方法
本发明属于冷冻靶系统温度场技术领域，具体涉及一种用于冷冻靶温度场预测的定向红外光-热耦合模拟方法。
背景技术：
惯性约束核聚变是通过激光均匀地辐照在冷冻靶丸的表面，从而达到高温高密度的点火条件从而实现聚变反应。为实现点火要求，避免瑞利-泰勒不稳定性，靶丸内氘-氘燃料冰层厚度均匀性需大于99％，燃料冰层内表面均方根粗糙度需小于1μm。燃料冰层的低模粗糙度主要受靶丸周围的温度场所决定，因此冷冻靶温度场控制的重要性尤为突出。
为了使靶丸内氘-氘燃料冰层质量满足点火要求，需要通过添加外界激励的方式使得燃料冰层更为均匀。现有一种被认为可行的方案是通过向靶丸表面投射定向红外，辐照靶丸表面低温区域，从而改善冰层质量。受冷冻靶装置尺寸限制，靶丸及其周围温度场无法通过传感器测量，因此需要通过数值计算的方法对实验进行预测及指导。
目前冷冻靶系统的结构(文献pointdesigntargets,specifications,andrequirementsforthe2010ignitioncampaignonthenationalignitionfacility)定向红外条件下温度场数值模拟的主要思路为光-热解耦，即红外光场与温度场的计算分开进行(文献computationaldesignofinfraredenhancedlayeringoficfcapsules)，忽略了红外光场和温度场之间的强烈耦合作用，仅仅考虑了红外光场对靶丸温度场的热贡献，忽略了温度场自身的红外效应，使得计算结果失真，无法为实验提供可靠的理论指导。
技术实现要素：
为了克服上述现有技术的缺点，本发明目的在于提高了一种用于冷冻靶温度场预测的定向红外光-热耦合模拟方法，考虑了红外光场和温度场之间的强烈耦合作用，提高了计算结果的准确性。
为了达到上述目的，本发明采取的技术方案为：
一种用于冷冻靶温度场预测的定向红外光-热耦合模拟方法，通过蒙特卡洛方法对冷冻靶系统外置定向红外光场进行求解，将红外热贡献与温度场控制方程相耦合，并引入策略来保证计算的可靠性和收敛速度，包括以下步骤：
1)确定光源粒子的初始位置及能量相关参数；
2)确定光源粒子的下一个碰撞距离l及粒子的新位置：
3)判断光源粒子是否逸出计算域，若是，停止光源粒子追踪；若否，继续进行下一步；
4)判断光源粒子是否与壁面发生碰撞，若是，执行步骤5)；若否，执行步骤6)；
5)根据壁面发射率判断光源粒子是否被壁面吸收，若是，计入热贡献，停止光源粒子追踪；若否，根据壁面漫射系数确定光源粒子反射方向；
6)根据单次反照度判断光源粒子是否被氦气中的杂质吸收，若是，计入热贡献，停止光源粒子追踪；若否，根据相位函数确定光源粒子的散射方向；
7)执行步骤2)，直至光源粒子追踪停止；
8)将计入的光源粒子热贡献带入导热微分方程，并结合动量方程、连续性方程计算出冷冻靶温度场分布。
所述的步骤1)中光源粒子的初始位置和能量相关参数要根据相应的热物理边界来确定。
所述的步骤2)中光源粒子在一次传播过程中的碰撞距离l通过下式进行计算：
其中r1为[0,1]区间内均匀分布的随机数，μa和μs分别为氦气的吸收系数和散射系数。
所述的步骤5)中光源粒子的新位置坐标为：
xn+1＝xn+lun+1
yn+1＝yn+lvn+1
zn+1＝zn+lwn+1
其中，x,y,z表示光源粒子在笛卡尔坐标系下的坐标值；下标n表示上一个状态，下表n+1表示新状态；u,v,w表示粒子前进方向；
光源粒子接触到壁面后，一部分会被壁面所吸收，取[0,1]区间上的均匀分布随机数r4，记壁面吸收率为ε，若r4≤ε，则光源粒子被壁面吸收，每个被吸收的光源粒子对壁面热量的贡献为q；剩余的光源粒子会在壁面处发生反射，反射分为镜面反射和漫反射，记壁面的漫反射系数为df，取[0,1]区间上的均匀分布随机数r5，若r5≤df，光源粒子发生漫反射，反之发生镜面反射；对于漫反射而言，分别取两个[0,1]区间内均匀分布的随机数r6，r7，记壁面的单位法向向量a＝(xn,yn,zn)，随机向量b＝(sin(2πr6)sin(2πr7),sin(2πr6)cos(2πr7),cos(2πr6))，若
则有
un+1＝sin(2πr6)sin(2πr7)
vn+1＝sin(2πr6)cos(2πr7)
wn1＝cos(2πr6)
反之重新取一组r6，r7直至满足上述条件；
对于镜面反射，有
对于灰体模型，壁面发射率等同于壁面吸收率，壁面漫射系数反映了光源粒子与壁面碰撞后发生漫反射的份额占比。
所述的步骤6)中假定杂质粒子为球形粒子，发生碰撞后，光源粒子在垂直传输方向的平面内各个方向散射概率相同，因此在该平面内方位角
光源粒子的前进方向为：
步骤6)中对于纯净的氦气，其吸收系数和散射系数均为0，此时不考虑光源粒子在氦气中的吸收及散射效应。
本发明的有益效果：本发明针对冷冻靶系统中面临的多组分，多结构的实际问题，创造性地提出红外光-热耦合的模拟方法，传统方法在计算光场时无法计算温度场，计算温度场时无法求解考虑光场的热效应，本发明考虑了红外光场和温度场之间的强烈耦合作用，从而提高了计算结果的准确性。
附图说明
图1是本发明实施例的冷冻靶结构示意图。
图2是本发明的流程图。
图3是本发明实施例的计算结果。
具体实施方式
下面对结合附图和实施例对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本实施例采用的点火靶结构，包括圆柱形的黑腔1和置于黑腔1中心的靶丸，黑腔1与靶丸之间为氦气腔2；黑腔1两端分别开设北极激光入射口6和南极激光入射口7；黑腔1的壁面上开有小孔用来插入光纤3；靶丸具有多层结构，最外侧为ch靶壳4，ch靶壳4内壁附着有固体燃料冰层5，固体燃料冰层5内部为空腔。
黑腔1的高度为10mm，内径5.44mm，黑腔1壁厚0.5mm；ch靶壳4外径1.1mm，内径0.91mm，固体燃料冰层5厚度60μm；光纤3内径0.3mm，入射口与黑腔1壁面的夹角为45°，入射口距黑腔1的中平面2.72mm。
参照图2，一种用于冷冻靶温度场预测的定向红外光-热耦合模拟方法，包括以下步骤：
1)根据相应的热物理边界来确定，确定光源粒子的初始位置及能量等相关参数；
2)确定光源粒子的下一个碰撞距离l及光源粒子的新位置；
光源粒子在一次传播过程中的碰撞距离l通过下式进行计算：
其中r1为[0,1]区间内均匀分布的随机数，μa和μs分别为氦气的吸收系数和散射系数；对于无吸收和散射的理想气体而言，光源粒子沿原方向矢量持续传播；
3)判断光源粒子是否通过北极激光入射口6和南极激光入射口7逸出计算域，若是，停止光源粒子追踪；若否，继续进行下一步；
4)判断光源粒子是否与壁面发生碰撞，若是，执行步骤5)；若否，执行步骤6)；
5)根据壁面发射率判断光源粒子是否被壁面吸收，若是，计入热贡献，停止光源粒子追踪；若否，根据壁面漫射系数确定粒子反射方向；
光源粒子的新位置坐标为：
xn+1＝xn+lun+1
yn+1＝yn+lvn+1
zn+1＝zn+lwn+1
其中，x,y,z表示光源粒子在笛卡尔坐标系下的坐标值；下标n表示上一个状态，下表n+1表示新状态；u,v,w表示粒子前进方向；
光源粒子接触到壁面后，一部分会被壁面所吸收，取[0,1]区间上的均匀分布随机数r4，记壁面吸收率为ε，若r4≤ε，则光源粒子被壁面吸收，每个被吸收的光源粒子对壁面热量的贡献为q；剩余的光源粒子会在壁面处发生反射，反射分为镜面反射和漫反射，记壁面的漫反射系数为df，取[0,1]区间上的均匀分布随机数r5，若r5≤df，光源粒子发生漫反射，反之发生镜面反射；对于漫反射而言，分别取两个[0,1]区间内均匀分布的随机数r6，r7，记壁面的单位法向向量a＝(xn,yn,zn)，随机向量b＝(sin(2πr6)sin(2πr7),sin(2πr6)cos(2πr7),cos(2πr6))，若
则有
un+1＝sin(2πr6)sin(2πr7)
vn+1＝sin(2πr6)cos(2πr7)
wn1＝cos(2πr6)
反之重新取一组r6，r7直至满足上述条件；
对于镜面反射，有
对于灰体模型，壁面发射率等同于壁面吸收率，壁面漫射系数反映了光源粒子与壁面碰撞后发生漫反射的份额占比；
6)根据单次反照度判断光源粒子是否被氦气中的杂质吸收，若是，计入热贡献，停止光源粒子追踪；若否，根据相位函数确定光源粒子的散射方向；
假定杂质粒子为球形粒子，发生碰撞后，光源粒子在垂直传输方向的平面内各个方向散射概率相同，因此在该平面内方位角
光源粒子的前进方向为：
步骤6)中对于纯净的氦气，其吸收系数和散射系数均为0，此时不考虑光源粒子在氦气中的吸收及散射效应；
7)执行步骤2)，直至光源粒子追踪停止；
8)将计入的光源粒子热贡献带入导热微分方程，并结合动量方程、连续性方程计算出冷冻靶温度场分布。
本实施例最终求得的靶丸外表面温度分布云图如图3所示，在靶丸正对光斑处温度会出现明显的抬升。
- 还没有人留言评论。精彩留言会获得点赞！