[发明专利]一种凸边界的内腔真空辐射仿真的高效并行计算方法

摘要

本发明属于内腔真空辐射仿真方法技术领域，具体是涉及一种基于有限元法的内腔真空辐射仿真的高效并行计算方法。本发明技术方案通过一种基于有限元法的内腔真空辐射仿真的高效并行算法，通过分区并行策略与遍历优化，实现了内腔真空辐射问题的高效仿真。与传统的循环遍历法相比，本发明基于分区并行计算的思想，设计了一种新的针对内腔真空辐射的高效并行计算方法，极大的提高了内腔真空辐射问题的仿真效率。因此，本发明具有较高的实际应用价值和很高的工程应用价值。

主权项

1.一种凸边界的内腔真空辐射仿真的高效并行计算方法，其特征在于：包括以下步骤：步骤1、采集具有内腔问题的固体结构的实际外形和物理尺寸，建立三维实体单元的有限元网格；采集本步骤建立的有限元网格的温度边界条件、热流边界条件、对流换热边界条件的边界面网格；步骤2、采集固体结构的热传导系数；采集温度边界条件、热流边界条件、对流换热系数、流体温度值；其中，温度边界条件为边界温度值，热流边界条件为边界热流值；步骤3、采集内腔辐射边界条件的边界面网格；采集内腔边界的发射率，记为ε；步骤4、设置瞬态热传递计算的时间步长；设置计算总时间；步骤5、基于步骤1建立的三维实体单元的有限元网格，对该网格进行均匀分区；将分区数记为F；步骤6、基于步骤3采集的内腔辐射边界条件与步骤5建立的有限元分区网格，采集每个网格分区中对应的内腔边界条件；将每个分区中的内腔边界面记为{U1,U2,…,UF}；采集每个分区内腔边界面网格的面积A；步骤7、基于步骤5的分区数量F，开启F个并行计算进程，每一个网格分区对应一个计算进程；计算进程号与网格分区号对应；基于步骤1建立的三维实体单元的有限元网格、步骤1采集的各个边界条件的边界面网格、步骤2采集的固体热传导系数与各个边界条件的数值、步骤4设置的瞬态热传递计算的时间步长与步骤5采集的分区网格，进行分区并行的不含内腔辐射值的有限元计算；步骤8、基于步骤6采集的每个分区的内腔边界面{U1,U2,…,UF}，采集步骤7中计算得到的每个分区的内腔边界面上温度值集合，记作{T1,T2,…,TF}；其中，T_i为第i个内腔边界面上的每个边界面网格的集合，i＝1、2…F；i＝1、2…F，ni为第i个网格分区的内腔边界面上的面单元数量；步骤9、基于步骤8得到的每个内腔边界单元上的温度值i＝1、2…F；计算每个边界单元通过辐射损失的能量定义其为辐射能量值，i＝1、2…F,j＝1、2…ni，ni为第i个网格分区的内腔边界面上的面单元数量；其中，ε为步骤3采集的内腔边界发射率，A为步骤6中采集的每个分区内腔边界面网格的面积，σ为Stefan‑Boltzmann常量；步骤10、将步骤9计算得到的每个边界单元的辐射能量值，存储到该单元对应的计算进程的内存空间中；步骤11、基于步骤6采集的每个分区中的内腔边界面，在其对应的计算进程上的所有内腔边界面上的网格单元进行遍历；步骤12、基于步骤9计算的每个边界单元的辐射能量值，在步骤11的内腔边界上的单元遍历过程中，每个边界单元遍历时接收其它所有边界单元所计算得到的辐射能量值；在该步骤中：1)若发送辐射能量值的单元与正在遍历的单元位于同一块分区网格中，那么说明此时，辐射能量值的传递是在同一个计算进程下完成的；此时直接通过内存的计算实现辐射能量值的传递；2)若发送辐射能量值的单元与正在遍历的单元不在同一块分区网格中，那么说明此时辐射能量值的传递需要从一个计算进程传递到另一个计算进程；此时通过跨进程的信息发送和信息接收来实现辐射能量值的传递；跨进程的信息发送与信息接收通过MPI并行编程实现；步骤13、完成步骤12后，每个边界单元接收完所有其它边界单元的辐射能量值，得到其它单元因为辐射传递到本单元而增加的能量；将步骤9计算得到的本边界单元的辐射能量值，加上其它单元因为辐射传递到本单元而增加的能量，得到内腔辐射后本单元获取的能量；步骤14、将步骤13得到的边界单元获取的能量作为能量输入，作为下一个时间步中步骤7计算的边界条件；步骤15、结束步骤11开始的边界单元遍历；步骤16、更新计算时刻；若第n个时间步时的计算时刻为t，则n+1时间步的计算时刻为t+Δt；Δt表示第n个时间步时和第n+1个时间步时之间的时间间隔；步骤17、重复步骤7～16，直到当前计算时刻超过步骤4设置的计算总时间，计算结束。