一种面向冲击系统的二维九速离散玻尔兹曼方法及装置

技术特征：
1.一种面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，包括以下步骤，s1:根据冲击系统的实际工况，初始化密度、速度和温度的物理场分布，以获得离散玻尔兹曼模型的离散化的平衡态速度分布函数，并将离散化的平衡态速度分布函数的值赋给离散化的速度分布函数；s2:计算所述离散玻尔兹曼模型中的碰撞项和对流项；s3:将所述碰撞项和所述对流项的数值结果代入所述离散玻尔兹曼模型中，并更新所述速度分布函数；s4:根据所述速度分布函数，计算密度、速度、温度和压强，并确定动理学矩关系；s5:基于计算的密度、速度、温度和所述动理学矩关系,计算所述平衡态速度分布函数的值；s6:判断时间迭代步数是否满足预设的条件；s7:若所述时间迭代步数满足条件，则输出预先设定的物理量的值。2.根据权利要求1所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述判断时间迭代步数是否满足预设的条件后，还包括以下步骤，若所述时间迭代步数不满足条件，则返回至步骤s2进行迭代循环。3.根据权利要求1所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述离散玻尔兹曼模型的建立步骤包括，对玻尔兹曼方程的速度空间进行离散化，得到包含有限个速度的离散玻尔兹曼模型；基于所述离散玻尔兹曼模型的动理学矩关系，确定宏观物理量的值；在连续性极限条件下，结合所述动理学矩关系，确定所述离散玻尔兹曼模型的流体力学方程组和非平衡物理量。4.根据权利要求3所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述离散玻尔兹曼模型的控制方程包括，其中，t为时间，τ为松弛时间，为hamilton算子，f
i
为离散化的速度分布函数，为离散化的平衡态速度分布函数，v
i
为离散速度，i＝0,1,2,
…
,8。5.根据权利要求4所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述离散玻尔兹曼模型的控制方程中v
i
的表达式为，其中，v
a
为常数。6.根据权利要求4所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述宏观物理量包括流体密度、流速和能量，所述流体密度ρ、所述流速u
α
和所述能量e的数学表达式分别为，ρ＝∑f
i
＝∑f
ieq
其中，v
i
为离散速度的大小，v
iα
为离散速度在α方向的分量，u
α
为宏观流速在α方向的分量，η
i
描述系统额外自由度所对应的内能，并且α＝x,y。7.根据权利要求6所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述离散玻尔兹曼模型的流体力学方程组包括，述离散玻尔兹曼模型的流体力学方程组包括，述离散玻尔兹曼模型的流体力学方程组包括，其中，压强p＝ρt，温度u＝|u|为宏观流速的大小，比热比d＝2为几何空间维度数，i表示额外自由度，并且α＝x,y；β＝x,y。8.根据权利要求6所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述离散玻尔兹曼模型的非平衡物理量包括非组织动量通量和非组织能量通量，所述非组织动量通量和所述非组织能量通量的表达式分别为，动量通量和所述非组织能量通量的表达式分别为，其中，对应非组织动量通量，对应非组织能量通量，为离散速度相对于流速u的本动速度，为本动速度在α方向的分量，为本动速度在β方向的分量。9.根据权利要求1所述的面向冲击系统的二维九速离散玻尔兹曼方法，其特征在于，所述基于计算的密度、速度、温度和所述动理学矩关系,计算所述平衡态速度分布函数的值的步骤包括，将所述平衡态速度分布函数的所述动理学矩关系统一为矩阵形式；通过矩阵求逆的方法求解所述平衡态速度分布函数。10.一种面向冲击系统的二维九速离散玻尔兹曼装置，其特征在于，包括，初始化模块，用于根据冲击系统的实际工况，初始化密度、速度和温度的物理场分布，以获得离散玻尔兹曼模型的离散化的平衡态速度分布函数，并将离散化的平衡态速度分布函数的值赋给离散化的速度分布函数；碰撞项模块，用于计算所述离散玻尔兹曼模型中的碰撞项；对流项模块，用于计算所述离散玻尔兹曼模型中的对流项；更新模块，用于将所述碰撞项和所述对流项的数值结果代入所述离散玻尔兹曼模型
中，并更新所述速度分布函数；动理学矩关系模块，用于计算更新的所述速度分布函数的密度、速度、温度和压强，并确定动理学矩关系；平衡态速度分布函数模块，用于基于计算的密度、速度、温度和所述动理学矩关系,计算所述平衡态速度分布函数的值；判断模块，用于判断时间迭代步数是否满足预设的条件；第一执行模块，用于若所述时间迭代步数满足条件，则输出预先设定的物理量的值。11.一种计算机设备，其特征在于，包括存储器、处理器及存储在存储器上的计算机程序，所述处理器执行所述计算机程序以实现权利要求1至9任意一项所述方法的步骤。12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1至9任意一项所述方法的步骤。

技术总结
本申请涉及一种面向冲击系统的二维九速离散玻尔兹曼方法及装置，其方法包括根据冲击系统的实际工况，初始化密度、速度和温度的物理场分布；计算离散玻尔兹曼模型中的碰撞项和对流项；将碰撞项和对流项的数值结果代入离散玻尔兹曼模型中，并更新速度分布函数；根据所述速度分布函数，计算密度、速度、温度和压强，并确定动理学矩关系；基于计算的密度、速度、温度和所述动理学矩关系,计算平衡态速度分布函数的值；判断时间迭代步数是否满足预设的条件；若时间迭代步数满足条件，则输出预先设定的物理量的值。本申请利于全面研究具有跨尺度特点的复杂冲击系统的数值模拟与仿真,具有更广泛的适用范围和更简介的操作流程。广泛的适用范围和更简介的操作流程。广泛的适用范围和更简介的操作流程。


技术研发人员：林传栋
受保护的技术使用者：中山大学
技术研发日：2022.07.06
技术公布日：2022/11/15