公开(公告)号：CN113761762A

公开(公告)日：2021-12-07

申请号：CN202110886391.1

申请日：2021-08-03

Applicant: 西北核技术研究所

Inventor： 韩峰 ,  王建国 ,  刘钰 ,  付梅艳 ,  闫凯 ,  姚成宝

IPC: G06F30/23 ,  G06T17/20 ,  G06F17/12 ,  G06F111/10

Abstract: 为克服平衡通量构造方法较为复杂、某些方法仅适用于二维问题的缺陷，本发明提出了一种用于有限元数值模拟后验误差估计的平衡通量构造方法，包括步骤：1)在集合Th中每个网格单元上定义矢量基函数；Th表示对计算区域Ω进行剖分所得到的网格单元的集合；2)在每个网格单元上计算所述矢量基函数的系数；3)在相邻网格单元公共边/面的两侧更新矢量基函数的系数；4)基于更新后的矢量基函数的系数和步骤1)定义的矢量基函数，计算平衡通量。本发明仅需要在Th中每个单元上求解一个三阶线性方程组和二阶线性方程组(d＝2)，或仅需要在Th中每个单元上求解一个四阶线性方程组和三阶线性方程组(d＝3)。

公开(公告)号：CN113761762B

公开(公告)日：2023-10-20

申请号：CN202110886391.1

申请日：2021-08-03

Applicant: 西北核技术研究所

Inventor： 韩峰 ,  王建国 ,  刘钰 ,  付梅艳 ,  闫凯 ,  姚成宝

IPC: G06F30/23 ,  G06T17/20 ,  G06F17/12 ,  G06F111/10

Abstract: 为克服平衡通量构造方法较为复杂、某些方法仅适用于二维问题的缺陷，本发明提出了一种用于有限元数值模拟后验误差估计的平衡通量构造方法，包括步骤：1)在集合Th中每个网格单元上定义矢量基函数；Th表示对计算区域Ω进行剖分所得到的网格单元的集合；2)在每个网格单元上计算所述矢量基函数的系数；3)在相邻网格单元公共边/面的两侧更新矢量基函数的系数；4)基于更新后的矢量基函数的系数和步骤1)定义的矢量基函数，计算平衡通量。本发明仅需要在Th中每个单元上求解一个三阶线性方程组和二阶线性方程组(d＝2)，或仅需要在Th中每个单元上求解一个四阶线性方程组和三阶线性方程组(d＝3)。