申请/专利权人:中国空气动力研究与发展中心计算空气动力研究所
申请日:2022-01-07
公开(公告)日:2023-03-21
公开(公告)号:CN114357907B
主分类号:G06F30/28
分类号:G06F30/28;G06F30/25;G06F111/10;G06F113/08;G06F119/14
优先权:
专利状态码:有效-授权
法律状态:2023.03.21#授权;2022.05.03#实质审查的生效;2022.04.15#公开
摘要:本发明适用于流体力学技术领域,提供了一种适用于拉格朗日型粒子类数值模拟的并行方法,包括获取计算区域内的粒子,将计算区域划分为多个子计算区域,每个子计算区域内的粒子数量相同;每个子计算区域对应的线程获取当前子计算区域中的粒子信息、计算参数信息,以及线程属性信息;设置每个子计算区域的联结区;获取当前时刻与各子计算区域相邻子计算区域的联结区内粒子信息,并对下一时刻各子计算区域中的粒子信息进行更新;当下一时刻的各子计算区域粒子越过当前各子计算区域的边界时,将越界的粒子发送至对应的子计算区域,该方法提高了粒子类数值模拟算法的计算效率和内存供给量。
主权项:1.一种适用于拉格朗日型粒子类数值模拟的并行方法,其特征在于,包括如下步骤:步骤S100:获取计算区域内的粒子,将计算区域划分为多个子计算区域,每个子计算区域内的粒子数量相同;每个子计算区域对应的线程获取当前子计算区域中的粒子信息、计算参数信息,以及线程属性信息;步骤S200:设置每个子计算区域的联结区;步骤S300:获取当前时刻与各子计算区域相邻子计算区域的联结区内粒子信息,并对下一时刻各子计算区域中的粒子信息进行更新;步骤S400:当下一时刻的流体力学领域各子计算区域粒子越过当前各子计算区域的边界时,将更新后下一时刻各子计算区域中越界的粒子发送至对应的子计算区域中;步骤S200中:每个子计算区域的联结区的宽度B不小于阈值宽度R;宽度阈值R的计算方法如下:R=k×rn,m其中,k为预设常数,rn,m表示两个相邻子计算区域的联结区的宽度值,n表示第n个子计算区域,m表示第m个子计算区域,第n个子计算区域与第m个子计算区域相邻,n≠m,n=0~N-1,m=0~N-1,N表示子计算区域的总数;rn,m的计算方法如下: 其中,maxhin表示第n个子计算区域中所有粒子中的最大光滑长度,hin表示第n个子计算区域中的第i个粒子的光滑长度,i表示第n个子计算区域中的第i个粒子;maxhjm表示第m个子计算区域中所有粒子中的最大光滑长度,hjm表示第m个子计算区域中的第j个粒子的光滑长度,j表示第m个子计算区域中的第j个粒子。
全文数据:
权利要求:
百度查询: 中国空气动力研究与发展中心计算空气动力研究所 一种适用于拉格朗日型粒子类数值模拟的并行方法
免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。