CN112992294B - 多孔介质lbm计算网格生成方法 - Google Patents

多孔介质lbm计算网格生成方法 Download PDF

Info

Publication number
CN112992294B
CN112992294B CN202110416484.8A CN202110416484A CN112992294B CN 112992294 B CN112992294 B CN 112992294B CN 202110416484 A CN202110416484 A CN 202110416484A CN 112992294 B CN112992294 B CN 112992294B
Authority
CN
China
Prior art keywords
lbm
space
porous medium
grid
geometric
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110416484.8A
Other languages
English (en)
Other versions
CN112992294A (zh
Inventor
肖光明
杜雁霞
刘磊
张超
魏东
杨肖峰
向静
邱芷薇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202110416484.8A priority Critical patent/CN112992294B/zh
Publication of CN112992294A publication Critical patent/CN112992294A/zh
Application granted granted Critical
Publication of CN112992294B publication Critical patent/CN112992294B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了多孔介质LBM计算网格生成方法，包括步骤：S1，获取多孔介质几何数模文件；S2，判断几何空间域的范围；S3，计算LBM物理求解域范围及各方向的网格节点数；S4，计算多孔介质几何数模文件的各固体骨架表面单元在LBM计算空间的区域范围；S5，判断固体骨架表面单元是否存在交点并标记；S6，进行流/固求解域与边界的判断；S7，根据交点排列的奇偶性对流/固区域进行区分；S8对LBM计算空间内的流/固区域与边界进行标识，进而完成LBM计算网格的生成等；本发明使得计算网格的生成能同时兼顾微细结构重建的真实性与网格空间分辨率的可调整性，能够提升LBM求解的灵活性与可靠性等。

Description

多孔介质LBM计算网格生成方法
技术领域
本发明涉及多孔介质及计算材料设计领域，更为具体的，涉及多孔介质LBM计算网格生成方法。
背景技术
多孔介质是由固体骨架和空隙中流体等多相物质体系共同占据空间的一种特殊组合体，而含多孔介质特征的工程材料，如高温陶瓷、泡沫金属等，通常具有体积密度小、相对质量轻、比强度高等优点，被广泛应用于对结构重量要求较为苛刻的航空航天领域。在具体的工程使用过程中，一般将多孔介质等价为均匀分布的虚拟连续体，并通过实验测量的方式对材料等效性能进行表征。然而，由于多孔介质在微细观尺度上具有大小不一、分布随机的真实孔隙结构，使得不同批次、加工工艺以及不同服役环境下的材料性能均存在不确定性，极大地制约了该类材料的研发与应用进程。
随着计算机性能的快速增长，计算材料设计（Computational Material Design,CMD）逐渐成为加速新材料研发的重要技术。该技术主要借助从微观到宏观不同尺度的数值模拟方法，研究分析材料多个尺度的物化现象与各类参数的影响规律，从而实现对材料各种力学与热学性质的准确预测，并为材料的优化设计提供指导方向。在众多计算材料数值模拟方法中，格子玻尔兹曼方法（Lattice Boltzmann Method, LBM）作为一种基于非平衡统计物理学Boltzmann方程的计算方法，对于复杂流体系统中的一些多尺度行为具有物理描述上的自适应性，加上边界条件处理简单以及程序易于并行求解等优点，近年来在多孔介质研究领域得到了广泛关注。
LBM计算网格的生成质量是影响多孔介质材料性能预测准确性的重要因素。当前，公开文献提出的多孔介质LBM计算网格生成方式可以概括为规则建模、随机生成与X-CT扫描等3种主要手段。其中，规则建模方式主要通过提取材料孔隙结构的重要几何特征，采用体心、面心、四面体以及八面体等规则模型进行简单近似，由于规则模型的描述通常可以解析表达，使得对应LBM计算网格的生成较为简单，但无法反映材料微观结构的随机分布特征；随机生成方式是采用高斯随机场对材料微细结构进行重建，该方法可实现LBM计算网格的同步生成，但与真实微细结构仍存在较大差异；X-CT扫描方式通过微CT成像技术获得多孔介质材料的三维空间点云数据，并经过一定的滤波处理后直接作为LBM计算网格进行求解，能真实反映多孔介质孔隙尺度的微细结构特征，但其空间分辨率会直接受到设备硬件的限制，无法根据实际计算状态进行网格调整，从而影响LBM求解的灵活性与可靠性。
发明内容
本发明的目的在于克服现有技术的不足，提供多孔介质LBM计算网格生成方法，使得计算网格的生成能同时兼顾微细结构重建的真实性与网格空间分辨率的可调整性，能够根据实际计算状态进行网格调整，提升LBM求解的灵活性与可靠性等。
本发明的目的是通过以下方案实现的：
多孔介质LBM计算网格生成方法，包括步骤：
S1，获取多孔介质几何数模文件；
S2，根据多孔介质几何数模文件的所有表面单元各顶点的坐标来判断几何空间域的范围；
S3，根据LBM计算的网格规模和几何空间域的范围，计算LBM网格空间分辨率、物理求解域范围及各方向的网格节点数；
S4，根据多孔介质几何数模文件和LBM网格空间分辨率设置，计算多孔介质几何数模文件的各表面单元在LBM计算空间的区域范围；
S5，根据各表面单元在LBM计算空间j-k平面的自适应区域范围，然后将该区域范围内的方向射线还原为真实几何空间的方向射线，利用射线相交来判断与射线对应表面单元是否存在交点并标记；
S6，根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判断；
S7，根据交点排列的奇偶性对流/固区域进行区分；
S8，根据各表面单元在LBM计算空间k-i平面与i-j平面的自适应区域范围，分别再从j方向和k方向对LBM计算空间内的流/固区域与边界进行标识，进而完成LBM计算网格的生成。
进一步地，在步骤S1中，通过X-CT扫描设备获得多孔介质几何数模文件；所述多孔 介质几何数模文件具有所有表面单元各顶点的空间坐标信息
Figure 906641DEST_PATH_IMAGE001
；其中，
Figure 679425DEST_PATH_IMAGE002
为表面单元编号，
Figure 889302DEST_PATH_IMAGE003
为表面单元顶点编号，下标
Figure 939298DEST_PATH_IMAGE004
和
Figure 71202DEST_PATH_IMAGE005
分 别表示第
Figure 941069DEST_PATH_IMAGE006
个表面单元及该表面单元对应的第
Figure 816621DEST_PATH_IMAGE007
个顶点。
进一步地，在步骤S2中，通过寻找多孔介质几何数模文件的所有表面单元各顶点 的x,y,z坐标的最大、最小值
Figure 455544DEST_PATH_IMAGE008
来判断几何空间域的范围；所 述几何空间域的范围指的是x, y, z各方向的最长距离
Figure 23929DEST_PATH_IMAGE009
，采用如下公式计算：
Figure 794755DEST_PATH_IMAGE011
Figure 819343DEST_PATH_IMAGE012
Figure 168416DEST_PATH_IMAGE013
进一步地，在步骤S3中，根据LBM计算的网格规模
Figure 137509DEST_PATH_IMAGE014
以及步骤S2中的几何空间域 的范围
Figure 230230DEST_PATH_IMAGE015
，计算LBM物理求解域范围
Figure 233958DEST_PATH_IMAGE016
及各方向的网格节点数
Figure 751003DEST_PATH_IMAGE017
，在LBM网格空间分辨率
Figure 207392DEST_PATH_IMAGE018
满足
Figure 103804DEST_PATH_IMAGE019
的基础上，几何空间域的范围包含于LBM物理求解域范围，即
Figure 837404DEST_PATH_IMAGE020
并使LBM计算网格总节点数接近于设置的网格规模，即
Figure 387334DEST_PATH_IMAGE021
此时，对应LBM计算空间的区域范围为
Figure 471965DEST_PATH_IMAGE022
、
Figure 31122DEST_PATH_IMAGE023
及
Figure 619230DEST_PATH_IMAGE024
；其中，
Figure 746586DEST_PATH_IMAGE025
、
Figure 911988DEST_PATH_IMAGE026
Figure 884623DEST_PATH_IMAGE027
分别为LBM物理求解域x，y，z方向的总长度，
Figure 717450DEST_PATH_IMAGE028
、
Figure 15707DEST_PATH_IMAGE029
、
Figure 668405DEST_PATH_IMAGE030
分别为LBM计算 空间i，j，k方向的网格节点数。
进一步地，在步骤S4中，根据多孔介质几何数模文件及LBM网格空间分辨率设置， 依次计算各表面单元在LBM计算空间的区域范围
Figure 707381DEST_PATH_IMAGE031
、
Figure 270081DEST_PATH_IMAGE032
及
Figure 598294DEST_PATH_IMAGE033
；其中，
Figure 879233DEST_PATH_IMAGE034
Figure 318305DEST_PATH_IMAGE035
Figure 1090DEST_PATH_IMAGE036
同时，使得表面单元的各顶点包含于对应的LBM物理求解域，即
Figure 375571DEST_PATH_IMAGE037
Figure 268441DEST_PATH_IMAGE038
Figure 386569DEST_PATH_IMAGE039
从而实现复杂三维几何模型在LBM计算空间的区域自适应；其中，
Figure 923861DEST_PATH_IMAGE040
为表 面单元编号,
Figure 593877DEST_PATH_IMAGE041
为表面单元顶点编号，下标
Figure 849409DEST_PATH_IMAGE004
和
Figure 33878DEST_PATH_IMAGE005
分别表示第
Figure 550310DEST_PATH_IMAGE006
个表面单元及 该表面单元对应的第
Figure 532172DEST_PATH_IMAGE007
个顶点，
Figure 9421DEST_PATH_IMAGE042
、
Figure 859565DEST_PATH_IMAGE043
分别为表面单元
Figure 371449DEST_PATH_IMAGE004
对应i方向最小及最 大网格节点编号，
Figure 258634DEST_PATH_IMAGE044
、
Figure 613392DEST_PATH_IMAGE045
分别为表面单元
Figure 611435DEST_PATH_IMAGE004
对应j方向最小及最大网格节点编 号，
Figure 977825DEST_PATH_IMAGE046
、
Figure 691703DEST_PATH_IMAGE047
分别为表面单元
Figure 143544DEST_PATH_IMAGE004
对应k方向最小及最大网格节点编号，
Figure 207928DEST_PATH_IMAGE048
、
Figure 694404DEST_PATH_IMAGE049
、
Figure 48025DEST_PATH_IMAGE050
分别表示表面单元
Figure 721583DEST_PATH_IMAGE004
对应顶点
Figure 451642DEST_PATH_IMAGE005
的x，y，z坐标。
进一步地，在步骤S5中，根据各表面单元在LBM计算空间j-k平面的自适应区域范围
Figure 792624DEST_PATH_IMAGE051
将该区域范围内的方向射线还原为真实几何空间的方向射线，对应的y、z坐标分别为
Figure 317146DEST_PATH_IMAGE052
并利用射线相交法判断与对应表面单元
Figure 478000DEST_PATH_IMAGE004
是否存在交点；若存在交点，则计算相 关交点的x坐标位置
Figure 887116DEST_PATH_IMAGE053
，并对LBM计算空间
Figure 941660DEST_PATH_IMAGE054
位置的i方向射线进行交点标记，即交 点总数
Figure 778029DEST_PATH_IMAGE055
；若不存在交点，
Figure 816392DEST_PATH_IMAGE056
，交点坐标
Figure 498040DEST_PATH_IMAGE057
；按照上述方式依次完成所有表面单元适应区域范围内i方向射线 的交点标记；其中，
Figure 672670DEST_PATH_IMAGE058
、
Figure 700448DEST_PATH_IMAGE059
分别表示在j-k平面选取i方向射线与表面单元
Figure 835894DEST_PATH_IMAGE004
进行交点判 断的j编号与k编号，
Figure 445867DEST_PATH_IMAGE060
、
Figure 615948DEST_PATH_IMAGE061
分别表示在j-k平面选取i方向射线j、k编号对应的实际y， z坐标。
进一步地，在步骤S6中，根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判断，包括如下步骤：
首先，将交点坐标
Figure 918754DEST_PATH_IMAGE062
(
Figure 275917DEST_PATH_IMAGE063
)按从小到大的顺序进行重新排 序，即
Figure 830526DEST_PATH_IMAGE064
其次，将交点坐标转换为LBM计算空间的i编号，即
Figure 714168DEST_PATH_IMAGE065
并指定对应LBM空间点
Figure 63241DEST_PATH_IMAGE066
为流/固边界。
进一步地，在步骤S1中，获得多孔介质几何数模文件，保存为STL或OBJ文件格式。
本发明的有益效果是：
本发明实施例使得计算网格的生成能同时兼顾微细结构重建的真实性与网格空间分辨率的可调整性，能够根据实际计算状态进行网格调整，提升LBM求解的灵活性与可靠性等。具体的，包括如下：
1. 本发明实施例不仅能真实反映多孔介质孔隙尺度的微细结构特征，而且空间分辨率能根据LBM网格规模的设置进行自适应调整。
2. 本发明实施例通过复杂三维几何模型表面单元在LBM计算空间的区域自适应，有效减少了各方向射线与表面单元的交点计算量，优化了LBM流/固求解域与边界的识别算法，使得多孔介质LBM计算网格的生成效率得到显著提升。在本发明的实施例中，如图2~图4所示，经算例测试，对于含89354个表面单元的多孔介质模型，百万量级LBM计算网格的生成时间约3秒，十亿量级LBM计算网格的生成时间也仅为425秒，可见显著提升了生成效率。
3. 本发明实施例通过多个方向的共同识别，可以有效降低了LBM计算区域及边界的误判概率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中自适应空间分辨率的多孔介质LBM计算网格生成流程图；
图2为本发明实施例中基于X-CT扫描的多孔介质几何数模（表面单元数：89354个）；
图3为本发明实施例中低空间分辨率的多孔介质LBM网格生成（网格规模：100×100×100；生成时间：约3秒）；
图4为本发明实施例中高空间分辨率的多孔介质LBM网格生成（网格规模：1000×1000×1000；生成时间：约425秒）；
图5为本发明方法步骤流程图。
具体实施方式
本说明书中所有实施例公开的所有特征，或隐含公开的所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以以任何方式组合和/或扩展、替换。
如图1~5所示，其中，图5为多孔介质LBM计算网格生成方法的步骤流程图，包括步骤：
S1，获取多孔介质几何数模文件；
S2，根据多孔介质几何数模文件的所有表面单元各顶点的坐标来判断几何空间域的范围；
S3，根据LBM计算的网格规模和几何空间域的范围，计算LBM网格空间分辨率、物理求解域范围及各方向的网格节点数；
S4，根据多孔介质几何数模文件和LBM网格空间分辨率设置，计算多孔介质几何数模文件的各表面单元在LBM计算空间的区域范围；
S5，根据各表面单元在LBM计算空间j-k平面的自适应区域范围，然后将该区域范围内的方向射线还原为真实几何空间的方向射线，利用射线相交来判断与射线对应表面单元是否存在交点并标记；
S6，根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判断；
S7，根据交点排列的奇偶性对流/固区域进行区分；
S8，根据各表面单元在LBM计算空间k-i平面与i-j平面的自适应区域范围，分别再从j方向和k方向对LBM计算空间内的流/固区域与边界进行标识，进而完成LBM计算网格的生成。
进一步地，在步骤S1中，通过X-CT扫描设备获得多孔介质几何数模文件；所述多孔 介质几何数模文件具有所有表面单元各顶点的空间坐标信息
Figure 297913DEST_PATH_IMAGE001
；其中，
Figure 125055DEST_PATH_IMAGE002
为表面单元编号，
Figure 269729DEST_PATH_IMAGE003
为表面单元顶点编号，下标
Figure 914337DEST_PATH_IMAGE004
和
Figure 508741DEST_PATH_IMAGE005
分 别表示第
Figure 264208DEST_PATH_IMAGE006
个表面单元及该表面单元对应的第
Figure 997809DEST_PATH_IMAGE007
个顶点。
进一步地，在步骤S2中，通过寻找多孔介质几何数模文件的所有表面单元各顶点 的x,y,z坐标的最大、最小值
Figure 688684DEST_PATH_IMAGE008
来判断几何空间域的范围；所 述几何空间域的范围指的是x, y, z各方向的最长距离
Figure 897949DEST_PATH_IMAGE009
，采用如下公式计算：
Figure 186159DEST_PATH_IMAGE011
Figure 906990DEST_PATH_IMAGE012
Figure 478917DEST_PATH_IMAGE013
进一步地，在步骤S3中，根据LBM计算的网格规模
Figure 576186DEST_PATH_IMAGE014
以及步骤S2中的几何空间域 的范围
Figure 549958DEST_PATH_IMAGE015
，计算LBM物理求解域范围
Figure 317057DEST_PATH_IMAGE016
及各方向的网格节点数
Figure 235334DEST_PATH_IMAGE017
，在LBM网格空间分辨率
Figure 274310DEST_PATH_IMAGE018
满足
Figure 837009DEST_PATH_IMAGE019
的基础上，几何空间域的范围包含于LBM物理求解域范围，即
Figure 165223DEST_PATH_IMAGE020
并使LBM计算网格总节点数接近于设置的网格规模，即
Figure 180583DEST_PATH_IMAGE021
此时，对应LBM计算空间的区域范围为
Figure 150813DEST_PATH_IMAGE022
、
Figure 568019DEST_PATH_IMAGE023
及
Figure 67134DEST_PATH_IMAGE024
；其中，
Figure 569790DEST_PATH_IMAGE025
、
Figure 953498DEST_PATH_IMAGE026
Figure 349844DEST_PATH_IMAGE027
分别为LBM物理求解域x，y，z方向的总长度，
Figure 895226DEST_PATH_IMAGE028
、
Figure 275392DEST_PATH_IMAGE029
、
Figure 931633DEST_PATH_IMAGE030
分别为LBM计算 空间i，j，k方向的网格节点数。
进一步地，在步骤S4中，根据多孔介质几何数模文件及LBM网格空间分辨率设置， 依次计算各表面单元在LBM计算空间的区域范围
Figure 448065DEST_PATH_IMAGE031
、
Figure 426997DEST_PATH_IMAGE032
及
Figure 904246DEST_PATH_IMAGE033
；其中，
Figure 488811DEST_PATH_IMAGE034
Figure 735116DEST_PATH_IMAGE035
Figure 12514DEST_PATH_IMAGE036
同时，使得表面单元的各顶点包含于对应的LBM物理求解域，即
Figure 711479DEST_PATH_IMAGE037
Figure 99735DEST_PATH_IMAGE038
Figure 466126DEST_PATH_IMAGE039
从而实现复杂三维几何模型在LBM计算空间的区域自适应；其中，
Figure 524212DEST_PATH_IMAGE040
为表 面单元编号,
Figure 100686DEST_PATH_IMAGE041
为表面单元顶点编号，下标
Figure 902420DEST_PATH_IMAGE004
和
Figure 247951DEST_PATH_IMAGE005
分别表示第
Figure 742517DEST_PATH_IMAGE006
个表面单元及 该表面单元对应的第
Figure 540709DEST_PATH_IMAGE007
个顶点，
Figure 143204DEST_PATH_IMAGE042
、
Figure 218608DEST_PATH_IMAGE043
分别为表面单元
Figure 8709DEST_PATH_IMAGE004
对应i方向最小及最 大网格节点编号，
Figure 169563DEST_PATH_IMAGE044
、
Figure 437733DEST_PATH_IMAGE045
分别为表面单元
Figure 633223DEST_PATH_IMAGE004
对应j方向最小及最大网格节点编 号，
Figure 594225DEST_PATH_IMAGE046
、
Figure 976796DEST_PATH_IMAGE047
分别为表面单元
Figure 189603DEST_PATH_IMAGE004
对应k方向最小及最大网格节点编号，
Figure 364232DEST_PATH_IMAGE048
、
Figure 105923DEST_PATH_IMAGE049
、
Figure 366003DEST_PATH_IMAGE050
分别表示表面单元
Figure 116922DEST_PATH_IMAGE004
对应顶点
Figure 146058DEST_PATH_IMAGE005
的x，y，z坐标。
进一步地，在步骤S5中，根据各表面单元在LBM计算空间j-k平面的自适应区域范围
Figure 55720DEST_PATH_IMAGE051
将该区域范围内的方向射线还原为真实几何空间的方向射线，对应的y、z坐标分别为
Figure 944042DEST_PATH_IMAGE052
并利用射线相交法判断与对应表面单元
Figure 357706DEST_PATH_IMAGE004
是否存在交点；若存在交点，则计算相 关交点的x坐标位置
Figure 116714DEST_PATH_IMAGE053
，并对LBM计算空间
Figure 590421DEST_PATH_IMAGE054
位置的i方向射线进行交点标记，即交 点总数
Figure 700459DEST_PATH_IMAGE055
；若不存在交点，
Figure 652235DEST_PATH_IMAGE056
，交点坐标
Figure 531329DEST_PATH_IMAGE057
；按照上述方式依次完成所有表面单元适应区域范围内i方向射线 的交点标记；其中，
Figure 51303DEST_PATH_IMAGE058
、
Figure 507692DEST_PATH_IMAGE059
分别表示在j-k平面选取i方向射线与表面单元
Figure 138525DEST_PATH_IMAGE004
进行交点判 断的j编号与k编号，
Figure 403284DEST_PATH_IMAGE060
、
Figure 218794DEST_PATH_IMAGE061
分别表示在j-k平面选取i方向射线j、k编号对应的实际y， z坐标。
进一步地，在步骤S6中，根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判断，包括如下步骤：
首先，将交点坐标
Figure 300494DEST_PATH_IMAGE062
(
Figure 594073DEST_PATH_IMAGE063
)按从小到大的顺序进行重新排 序，即
Figure 447759DEST_PATH_IMAGE064
其次，将交点坐标转换为LBM计算空间的i编号，即
Figure 575115DEST_PATH_IMAGE065
并指定对应LBM空间点
Figure 147042DEST_PATH_IMAGE066
为流/固边界。
进一步地，在步骤S1中，获得多孔介质几何数模文件，保存为STL或OBJ文件格式。
近年来，采用LBM方法对多孔介质材料微细结构内部的流动与传热过程进行精细化数值模拟，已经成为多孔介质传热传质与计算材料设计研究领域的共同热点之一。影响多孔介质性能预测准确性的主要因素除了LBM模型与方法自身的计算精度外，对应计算网格的生成质量也会产生重要作用。目前，多孔介质计算网格的生成主要涉及材料微细结构的建模技术与LBM求解域及边界的识别技术。其中，多孔介质材料微细结构的重建方法主要包括规则建模、随机生成与X-CT扫描等。
对于规则建模与随机生成重建方法，LBM求解域及边界的判别算法较为简单，并且网格规模可以根据计算状态的不同进行任意调整，但这两类方法都属于近似建模方法，无法准确反映材料微细结构的几何特征。相较而言，X-CT扫描方法能真实地还原材料的孔隙结构特征，而且通过专业的CT图像处理软件，可以得到较为精细的三维CAD数模。然而，多孔介质材料的三维数模往往包含数十万甚至上百万三角形或多边形几何单元，采用传统的射线法进行求解域与边界的判断，计算量将随着LBM网格规模的扩大剧烈增长。在传统的研究中，一般将CT扫描获得的三维空间点云数据经过适当的滤波处理后直接用于求解，这就导致LBM计算网格的空间分辨率受到了设备硬件的限制，无法根据实际计算状态进行网格调整，使得LBM求解的灵活性与可靠性受到影响。
目前尚无直接通过X-CT扫描得到的几何数模进行不同空间分辨率多孔介质LBM计算网格生成的专利及论文。
针对上述存在的问题，本发明实施例的目的是基于多孔介质复杂三维几何模型的自适应空间分辨率技术，寻找一种LBM求解域及边界的高效判别方法，使得计算网格的生成能同时兼顾微细结构重建的真实性与网格空间分辨率的可调整性。
在本发明的其他实施例中，包括如下步骤：
步骤1，通过X-CT扫描设备获得多孔介质几何数模文件；所述多孔介质几何数模文 件包含表面单元的空间坐标信息
Figure 244311DEST_PATH_IMAGE067
；其中，
Figure 952504DEST_PATH_IMAGE068
为表面单 元编号,
Figure 250761DEST_PATH_IMAGE069
为表面单元顶点编号，下标
Figure 169039DEST_PATH_IMAGE004
和
Figure 945365DEST_PATH_IMAGE005
分别表示第
Figure 632698DEST_PATH_IMAGE006
个表面单元及 该表面对应的第
Figure 101857DEST_PATH_IMAGE007
个顶点。
步骤2，根据多孔介质几何数模文件，判断几何空间域的范围；所述几何空间域的 范围包括x, y, z各方向的最长距离
Figure 114287DEST_PATH_IMAGE070
。
步骤3，根据LBM计算的网格规模
Figure 818938DEST_PATH_IMAGE071
以及步骤S2中的几何空间域的范围
Figure 501723DEST_PATH_IMAGE070
，自适应计算LBM物理求解域范围
Figure 838DEST_PATH_IMAGE072
及各方向的网格节 点数
Figure 503495DEST_PATH_IMAGE073
，在LBM网格空间分辨率
Figure 480678DEST_PATH_IMAGE018
严格满足
Figure 17970DEST_PATH_IMAGE074
的基础上，一方面要求几何空间域的范围包含于LBM物理求解域范围，即
Figure 828931DEST_PATH_IMAGE075
另一方面使得LBM计算网格总节点数接近于设置的网格规模，即
Figure 943517DEST_PATH_IMAGE076
此时，对应LBM计算空间的区域范围为
Figure 865337DEST_PATH_IMAGE077
、
Figure 381769DEST_PATH_IMAGE078
及
Figure 363631DEST_PATH_IMAGE079
。
步骤4，根据所述多孔介质几何数模文件及LBM网格空间分辨率设置，依次计算各 表面单元在LBM计算空间的区域范围
Figure 840880DEST_PATH_IMAGE080
、
Figure 425445DEST_PATH_IMAGE081
及
Figure 957837DEST_PATH_IMAGE082
；其中：
Figure 969655DEST_PATH_IMAGE083
Figure 934200DEST_PATH_IMAGE084
Figure 197823DEST_PATH_IMAGE085
同时，使得表面单元的各顶点包含于对应的LBM物理求解域，即
Figure 423268DEST_PATH_IMAGE086
Figure 481353DEST_PATH_IMAGE087
Figure 57828DEST_PATH_IMAGE088
从而实现复杂三维几何模型在LBM计算空间的区域自适应；其中，
Figure 859562DEST_PATH_IMAGE068
为表面单元编号,
Figure 205093DEST_PATH_IMAGE069
为表面单元顶点编号，下标
Figure 699659DEST_PATH_IMAGE004
和
Figure 497851DEST_PATH_IMAGE005
分别表示第
Figure 103276DEST_PATH_IMAGE006
个表 面单元及该表面单元对应的第
Figure 178679DEST_PATH_IMAGE007
个顶点；
步骤5，根据各表面单元在LBM计算空间j-k平面的自适应区域范围
Figure 968781DEST_PATH_IMAGE089
Figure 126705DEST_PATH_IMAGE090
将该区域范围内的方向射线还原为真实几何空间的方向射线，对应的y、z坐标分别为
Figure 660455DEST_PATH_IMAGE091
Figure 855944DEST_PATH_IMAGE092
并利用射线相交法判断与对应表面单元
Figure 426733DEST_PATH_IMAGE004
是否存在交点；若存在交点，则计算相 关交点的x坐标位置
Figure 199517DEST_PATH_IMAGE093
，并对LBM计算空间
Figure 412324DEST_PATH_IMAGE094
位置的i方向射线进行交点标记，即交 点总数
Figure 321374DEST_PATH_IMAGE095
；若不存在交点，
Figure 328644DEST_PATH_IMAGE096
，交点坐标
Figure 323145DEST_PATH_IMAGE097
；按照上述方式依次完成所有表面单元适应区域范围内i方向射线 的交点标记。
步骤6，根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判 断；首先，将交点坐标
Figure 339643DEST_PATH_IMAGE098
(
Figure 244145DEST_PATH_IMAGE099
)按从小到大的顺序进行重新 排序，即
Figure 281371DEST_PATH_IMAGE100
其次，将交点坐标转换为LBM计算空间的i编号，即
Figure 169693DEST_PATH_IMAGE101
并指定对应LBM空间点
Figure 583356DEST_PATH_IMAGE102
为流/固边界；
步骤7，根据交点排列的奇偶性对流/固区域进行区分；例如指定
Figure 605015DEST_PATH_IMAGE103
区间为流体区域，则
Figure 688508DEST_PATH_IMAGE104
区间为固体区域；其 中，
Figure 657601DEST_PATH_IMAGE105
；
步骤8，根据各表面单元在LBM计算空间k-i平面与i-j平面的自适应区域范围，类似于步骤5至步骤7，分别再从j方向和k方向对LBM计算空间内的流/固区域与边界进行标识，进而完成LBM计算网格的生成。
在步骤3中，LBM计算网格的空间分辨率可以通过指定网格规模来确定，也可直接指定各方向的网格节点数；在步骤6中，交点坐标既可以按从小到大的顺序进行排序，也可以按从大到小的顺序进行排序；步骤4至步骤8中，i，j，k方向交点的判断顺序没有严格要求，可以任意变换，同时为了进一步减少计算量，也可以只取两个方向。
在本发明的实施例中，直接根据X-CT扫描得到的几何数模以及LBM计算网格规模进行不同空间分辨率网格生成。
在本发明的实施例中，通过指定LBM计算总的网格规模，并根据几何数模的空间域范围对LBM物理求解域范围及各方向的网格节点数进行自适应计算。
在本发明的实施例中，在计算各表面单元在LBM计算空间各方向区域范围的基础上，根据某一特定方向的自适应区域范围进行射线与表面单元的交点判断与求解。
在本发明的实施例中，分别从x，y，z三个方向对计算区域及边界进行多次识别。
本发明功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，在一台计算机设备（可以是个人计算机，服务器，或者网络设备等）以及相应的软件中执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、或者光盘等各种可以存储程序代码的介质，进行测试或者实际的数据在程序实现中存在于只读存储器（Random Access Memory，RAM）、随机存取存储器（Random Access Memory，RAM）等。

Claims (8)

1.多孔介质LBM计算网格生成方法，其特征在于，包括步骤：
S1，获取多孔介质几何数模文件；
S2，根据多孔介质几何数模文件的所有表面单元各顶点的坐标来判断几何空间域的范围；
S3，根据LBM计算的网格规模和几何空间域的范围，计算LBM网格空间分辨率、物理求解域范围及各方向的网格节点数；
S4，根据多孔介质几何数模文件和LBM网格空间分辨率设置，计算多孔介质几何数模文件的各表面单元在LBM计算空间的区域范围；
S5，根据各表面单元在LBM计算空间j-k平面的自适应区域范围，然后将该区域范围内的方向射线还原为真实几何空间的方向射线，利用射线相交来判断与射线对应表面单元是否存在交点并标记；
S6，根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判断；
S7，根据交点排列的奇偶性对流/固区域进行区分；
S8，根据各表面单元在LBM计算空间k-i平面与i-j平面的自适应区域范围，分别再从j方向和k方向对LBM计算空间内的流/固区域与边界进行标识，进而完成LBM计算网格的生成。
2.根据权利要求1所述的多孔介质LBM计算网格生成方法，其特征在于，在步骤S1中，通 过X-CT扫描设备获得多孔介质几何数模文件；所述多孔介质几何数模文件具有所有表面单 元各顶点的空间坐标信息
Figure 304199DEST_PATH_IMAGE001
；其中，
Figure 384151DEST_PATH_IMAGE002
为表面单元编号，
Figure 3351DEST_PATH_IMAGE003
为表面单元顶点编号，下标
Figure 552275DEST_PATH_IMAGE004
和
Figure 282334DEST_PATH_IMAGE005
分别表示第
Figure 951213DEST_PATH_IMAGE006
个表面单元及该表面单 元对应的第
Figure 741314DEST_PATH_IMAGE007
个顶点。
3.根据权利要求1或2任一所述的多孔介质LBM计算网格生成方法，其特征在于，在步骤 S2中，通过寻找多孔介质几何数模文件的所有表面单元各顶点的x,y,z坐标的最大、最小值
Figure 777534DEST_PATH_IMAGE008
来判断几何空间域的范围；所述几何空间域的范围指的是 x, y, z各方向的最长距离
Figure 45705DEST_PATH_IMAGE009
，采用如下公式计算：
Figure 326830DEST_PATH_IMAGE011
Figure 319188DEST_PATH_IMAGE012
Figure 391049DEST_PATH_IMAGE013
。
4.根据权利要求3所述的多孔介质LBM计算网格生成方法，其特征在于，在步骤S3中，根 据LBM计算的网格规模
Figure 565679DEST_PATH_IMAGE014
以及步骤S2中的几何空间域的范围
Figure 432004DEST_PATH_IMAGE015
，计算LBM物理求解 域范围
Figure 439886DEST_PATH_IMAGE016
及各方向的网格节点数
Figure 49859DEST_PATH_IMAGE017
，在LBM网格空间分辨率
Figure 78995DEST_PATH_IMAGE018
满足
Figure 116221DEST_PATH_IMAGE019
的基础上，几何空间域的范围包含于LBM物理求解域范围，即
Figure 614330DEST_PATH_IMAGE020
并使LBM计算网格总节点数接近于设置的网格规模，即
Figure 293573DEST_PATH_IMAGE021
此时，对应LBM计算空间的区域范围为
Figure 911636DEST_PATH_IMAGE022
、
Figure 650922DEST_PATH_IMAGE023
及
Figure 105168DEST_PATH_IMAGE024
；其中，
Figure 322523DEST_PATH_IMAGE025
、
Figure 60672DEST_PATH_IMAGE026
Figure 439700DEST_PATH_IMAGE027
分别为LBM物理求解域x，y，z方向的总长度，
Figure 646822DEST_PATH_IMAGE028
、
Figure 402288DEST_PATH_IMAGE029
、
Figure 260523DEST_PATH_IMAGE030
分别为LBM计算空间i， j，k方向的网格节点数。
5.根据权利要求4所述的多孔介质LBM计算网格生成方法，其特征在于，
在步骤S4中，根据多孔介质几何数模文件及LBM网格空间分辨率设置，依次计算各表面 单元在LBM计算空间的区域范围
Figure 341612DEST_PATH_IMAGE031
、
Figure 767520DEST_PATH_IMAGE032
及
Figure 61098DEST_PATH_IMAGE033
；其中，
Figure 773839DEST_PATH_IMAGE034
Figure 760250DEST_PATH_IMAGE035
Figure 207543DEST_PATH_IMAGE036
同时，使得表面单元的各顶点包含于对应的LBM物理求解域，即
Figure 304812DEST_PATH_IMAGE037
Figure 606480DEST_PATH_IMAGE038
Figure 763792DEST_PATH_IMAGE039
从而实现复杂三维几何模型在LBM计算空间的区域自适应；其中，
Figure 432802DEST_PATH_IMAGE040
为表面单 元编号,
Figure 68183DEST_PATH_IMAGE041
为表面单元顶点编号，下标
Figure 755516DEST_PATH_IMAGE004
和
Figure 83729DEST_PATH_IMAGE005
分别表示第
Figure 974456DEST_PATH_IMAGE006
个表面单元及该表 面单元对应的第
Figure 679107DEST_PATH_IMAGE007
个顶点，
Figure 220946DEST_PATH_IMAGE042
、
Figure 720061DEST_PATH_IMAGE043
分别为表面单元
Figure 95154DEST_PATH_IMAGE004
对应i方向最小及最大网 格节点编号，
Figure 72337DEST_PATH_IMAGE044
、
Figure 734263DEST_PATH_IMAGE045
分别为表面单元
Figure 404279DEST_PATH_IMAGE004
对应j方向最小及最大网格节点编号，
Figure 269598DEST_PATH_IMAGE046
、
Figure 50472DEST_PATH_IMAGE047
分别为表面单元
Figure 363641DEST_PATH_IMAGE004
对应k方向最小及最大网格节点编号，
Figure 955291DEST_PATH_IMAGE048
、
Figure 557174DEST_PATH_IMAGE049
、
Figure 345001DEST_PATH_IMAGE050
分别表示表面单元
Figure 715939DEST_PATH_IMAGE004
对应顶点
Figure 727758DEST_PATH_IMAGE005
的x，y，z坐标。
6.根据权利要求5所述的多孔介质LBM计算网格生成方法，其特征在于，
在步骤S5中，根据各表面单元在LBM计算空间j-k平面的自适应区域范围
Figure 302090DEST_PATH_IMAGE051
将该区域范围内的方向射线还原为真实几何空间的方向射线，对应的y、z坐标分别为
Figure 690346DEST_PATH_IMAGE052
并利用射线相交法判断与对应表面单元
Figure 915791DEST_PATH_IMAGE004
是否存在交点；若存在交点，则计算相关交 点的x坐标位置
Figure 98510DEST_PATH_IMAGE053
，并对LBM计算空间
Figure 422788DEST_PATH_IMAGE054
位置的i方向射线进行交点标记，即交点总 数
Figure 349156DEST_PATH_IMAGE055
；若不存在交点，
Figure 429107DEST_PATH_IMAGE056
，交点坐标
Figure 48307DEST_PATH_IMAGE057
；按照上述方式依次完成所有表面单元适应区域范围内i方向射线 的交点标记；其中，
Figure 393969DEST_PATH_IMAGE058
、
Figure 858449DEST_PATH_IMAGE059
分别表示在j-k平面选取i方向射线与表面单元
Figure 58486DEST_PATH_IMAGE004
进行交点判 断的j编号与k编号，
Figure 599320DEST_PATH_IMAGE060
、
Figure 619228DEST_PATH_IMAGE061
分别表示在j-k平面选取i方向射线j、k编号对应的实际y， z坐标。
7.根据权利要求6所述的多孔介质LBM计算网格生成方法，其特征在于，
在步骤S6中，根据LBM计算空间方向射线的交点标记结果进行流/固求解域与边界的判断，包括如下步骤：
首先，将交点坐标
Figure 152978DEST_PATH_IMAGE062
(
Figure 207521DEST_PATH_IMAGE063
)按从小到大的顺序进行重新排序，即
Figure 653677DEST_PATH_IMAGE064
其次，将交点坐标转换为LBM计算空间的i编号，即
Figure DEST_PATH_IMAGE065
并指定对应LBM空间点
Figure 957620DEST_PATH_IMAGE066
为流/固边界。
8.根据权利要求1所述的多孔介质LBM计算网格生成方法，其特征在于，在步骤S1中，获得多孔介质几何数模文件，保存为STL或OBJ文件格式。
CN202110416484.8A 2021-04-19 2021-04-19 多孔介质lbm计算网格生成方法 Active CN112992294B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110416484.8A CN112992294B (zh) 2021-04-19 2021-04-19 多孔介质lbm计算网格生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110416484.8A CN112992294B (zh) 2021-04-19 2021-04-19 多孔介质lbm计算网格生成方法

Publications (2)

Publication Number Publication Date
CN112992294A CN112992294A (zh) 2021-06-18
CN112992294B true CN112992294B (zh) 2021-08-10

Family

ID=76341040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110416484.8A Active CN112992294B (zh) 2021-04-19 2021-04-19 多孔介质lbm计算网格生成方法

Country Status (1)

Country Link
CN (1) CN112992294B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115455794B (zh) * 2022-08-10 2024-03-29 武汉理工大学 基于连通孔隙划分计算区域的lbm并行优化方法、装置及存储介质
CN116502373B (zh) * 2023-06-30 2023-10-20 陕西空天信息技术有限公司 流道计算网格生成方法及装置
CN116894282B (zh) * 2023-09-07 2023-11-24 中国空气动力研究与发展中心计算空气动力研究所 空间点集与多连通网格区域拓扑关系的识别方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104268401A (zh) * 2014-09-25 2015-01-07 中国科学技术大学 一种粘性指进现象模拟方法及装置
CN104933225A (zh) * 2015-05-25 2015-09-23 中国科学院过程工程研究所 实现计算流体力学大规模实时模拟的方法
CN107966472A (zh) * 2017-12-05 2018-04-27 中国空气动力研究与发展中心计算空气动力研究所 一种高温接触热阻的无损快速测量方法
US10648292B2 (en) * 2017-03-01 2020-05-12 International Business Machines Corporation Cognitive enhanced oil recovery advisor system based on digital rock simulator
US10762252B2 (en) * 2013-07-31 2020-09-01 Dassault Systemes Simulia Corp. Temperature coupling algorithm for hybrid thermal lattice boltzmann method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103337097A (zh) * 2013-07-02 2013-10-02 上海大学 适用于格子Boltzmann方法的多重笛卡尔网格生成方法
CN103970715A (zh) * 2014-04-21 2014-08-06 中国石油大学(北京) 一种计算网格生成方法及装置
CN108615259B (zh) * 2018-04-27 2022-06-24 苏州数算软云科技有限公司 建立血液流动计算模型用以评价冠状动脉血流情况的方法
CN111489447B (zh) * 2020-04-14 2022-04-29 西北工业大学 一种适用于格子Boltzmann方法的直角网格自适应建模方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10762252B2 (en) * 2013-07-31 2020-09-01 Dassault Systemes Simulia Corp. Temperature coupling algorithm for hybrid thermal lattice boltzmann method
CN104268401A (zh) * 2014-09-25 2015-01-07 中国科学技术大学 一种粘性指进现象模拟方法及装置
CN104933225A (zh) * 2015-05-25 2015-09-23 中国科学院过程工程研究所 实现计算流体力学大规模实时模拟的方法
US10648292B2 (en) * 2017-03-01 2020-05-12 International Business Machines Corporation Cognitive enhanced oil recovery advisor system based on digital rock simulator
CN107966472A (zh) * 2017-12-05 2018-04-27 中国空气动力研究与发展中心计算空气动力研究所 一种高温接触热阻的无损快速测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
岩石孔隙结构表征与流体输运可视化研究进展;宋睿;《西南石油大学学报》;20181126;第86-100页 *

Also Published As

Publication number Publication date
CN112992294A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
CN112992294B (zh) 多孔介质lbm计算网格生成方法
Landry et al. Nonparametric inference of the neutron star equation of state from gravitational wave observations
CN112382352B (zh) 基于机器学习的金属有机骨架材料结构特征快速评估方法
WO2022007398A1 (zh) 一种基于gpu加速技术的非结构网格洪水过程模拟系统
Hao et al. Effective visualization of temporal ensembles
TWI767233B (zh) 用於使用基於機器學習的編碼在電子電路佈局中的圖案分類之方法、非暫時性電腦可讀媒體、及系統
CN108540988B (zh) 一种场景划分方法及装置
Liu et al. Memory-efficient modeling and slicing of large-scale adaptive lattice structures
CN115457492A (zh) 目标检测方法、装置、计算机设备及存储介质
Zhu et al. An improved tolerance analysis method based on skin model shapes of planar parts
CN116258840B (zh) 层级细节表示树的生成方法、装置、设备及存储介质
Pieritz et al. 3D tomography analysis of the inner structure of pebbles and pebble beds
CN117333637A (zh) 三维场景的建模及渲染方法、装置及设备
CN115062551B (zh) 一种基于时序神经网络的湿物理过程参数化方法
CN111751403A (zh) 一种热障涂层数值重构模型测试方法和装置
CN111508075B (zh) 一种闭孔泡沫铝的三维真实有限元模型建模方法
CN114139482A (zh) 一种基于深度度量学习的eda电路失效分析方法
Heidig et al. Influence of resolution of rasterized geometries on porosity and specific surface area exemplified for model geometries of porous media
CN108764514B (zh) 一种基于并行运算的光伏发电功率预测方法
Al-Mestarehi et al. Creating a Complete Model of the Wooden Pattern from Laser Scanner Point Clouds Using Alpha Shapes.
Campbell et al. Mesh-free simulation of aerosol evolution
CN116977592B (zh) 基于卷绕数的三维结构化重建方法、装置和计算机介质
Gao et al. Urban wind field prediction based on sparse sensors and physics‐informed graph‐assisted auto‐encoder
CN116882255B (zh) 一种基于傅里叶级数随机生成多孔介质模型的方法及系统
CN116310226B (zh) 三维对象层级模型生成方法、装置、设备和存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant