一种快速获取固相颗粒所处非结构气相网格信息的方法技术

本发明专利技术公开了一种快速获取固相颗粒所处非结构气相网格信息的方法。该方法包括：（1）建立大型和小型两套结构网格对非结构气相网格加以覆盖；（2）扫描气相网格，将各个大型结构网格所覆盖的气相网格进行统计注册；（3）当计算小型结构网格与其所处大型结构网格内的气相网格之间距离取得最小值时，即为该小型结构网格所处的气相网格，并将它们之间的关系建立数组；（4）?后续计算中，根据颗粒位置确定其所处的小型结构网格，由步骤（3）中建立的数组来获取该颗粒所处非结构气相网格的信息。本发明专利技术的方法，在保证气固两相数值模拟准确性的同时，除去了传统方法中因搜索颗粒所处气相网格而带来的巨大计算量，从而大大提高了数值计算效率。

【技术实现步骤摘要】

本专利技术属于气固两相流动计算机数值模拟领域，具体涉及一种快速获取气固两相数值模拟中固相颗粒所处非结构气相网格信息的方法。
技术介绍
气固两相流动广泛存在于能源、化工、环保等领域，随着计算机速度的大幅提高， 气固两相流动的计算机数值模拟技术也随之快速发展，并在生产实际中发挥出越来越重要的作用。在气固两相流动计算机数值模拟中，一种行之有效的措施是将气相场和固相颗粒场分开处理，分别采用欧拉方法和拉格朗日方法进行数值模拟，这样可以充分提高计算效率和增加数值模拟的准确性。在计算离散颗粒所受的气体曳力以及颗粒对气相流场反作用的求解过程中，必须获得颗粒所处当地气相网格的信息，如当地气体流速等。对于复杂外形几何体的非结构网格，如何准确、快速判断颗粒所处的当前网格，是一个难题。如果采用对多个非结构气相网格进行数据平均的方法加以解决，则该过程将不可避免地牺牲大量气相场信息，从而降低了计算精度；然而采取在每个时间步长中都求解颗粒到计算区域内所有非结构气相网格中心点之间距离，来获得颗粒所处的当前非结构气相网格的方法，又会带来极其巨大的计算量，从而大大降低了计算速度。若能找到一种快速确定固相颗粒所在非结构气相网格的方法，则可使气固两相数值模拟的效率大大提高。
技术实现思路
专利技术目的针对现有技术中存在的不足，本专利技术的目的是提供一种快速获取气固两相数值模拟中固相颗粒所处非结构气相网格信息的方法，该方法除去了传统方法中因搜索颗粒所处非结构气相网格而带来的巨大计算量，从而大大提高了气固两相数值计算的速度和精度。技术方案为了实现上述专利技术目的，本专利技术采用快速获取固相颗粒所处非结构气相网格信息的方法的技术方案如下步骤1)在气相计算所划分的非结构网格中，分别建立大型、小型两套结构网格对整个计算区域加以覆盖；其中，大型结构网格尺度大于非结构气相网格尺度，小型结构网格尺度小于非结构气相网格尺度；步骤2)扫描计算区域内的所有非结构气相网格，将每个大型结构网格所覆盖的非结构气相网格的序号进行统计注册，建立数组以备后继计算使用；步骤3)确定每个小型结构网格所在的大型结构网格，根据步骤2)中建立的数组，求解小型结构网格的中心点与其所处的大型结构网格内所有非结构气相网格中心点的距离，其中最小距离者即为该小型结构网格所处的非结构气相网格，之后对每个小型结构网格所属的非结构气相网格序号进行统计注册，建立数组以备后继计算使用；步骤4)在每个时间步长的计算中，根据颗粒位置确定其所处的小型结构网格，由步骤33)中建立的数组从而得到该颗粒所处的非结构气相网格信息。在步骤1)中，大型结构网格尺度要大于非结构气相网格尺度，其尺度越大，则建立的大型结构网格个数越少，单个大型结构网格所含的非结构气相网格个数越多，因而步骤 3)中求解小型结构网格的中心点与其所处的大型结构网格内所有非结构气相网格中心点距离的次数就越多，在数值计算中，大型结构网格尺度应根据以上所述原则和计算机软、硬件条件来设定。在步骤1)中，小型结构网格尺度要小于非结构气相网格尺度，其尺度越小，则建立的小型结构网格个数越多，单个非结构气相网格所含的小型结构网格的个数越多，因而步骤4)中获取颗粒所处的非结构气相网格信息越精确，步骤3)中求解小型结构网格的中心点与其所处的大型结构网格内所有非结构气相网格中心点距离的次数就越多，在数值计算中，小型结构网格尺度根据以上所述原则和计算机软、硬件条件来设定。在步骤3)中，当小型结构网格位于大型结构网格边沿附近时，对于二维计算，求解小型结构网格与非结构气相网格中心点距离的范围扩大到该小型结构网格所处的大型结构网格相邻的8个大型结构网格。在步骤3)中，当小型结构网格位于大型结构网格边沿附近时，对于三维计算，求解小型结构网格与非结构气相网格中心点距离的范围扩大到该小型结构网格所处的大型结构网格相邻的沈个大型结构网格。有益效果本专利技术的快速确定固相颗粒所在气相非结构网格的方法，对非结构网格无需进行平均处理，而是首先通过分割计算区域、邻域搜索以及计算小型结构网格-非结构气相网格中心点最小距离法确定小型结构网格所处的非结构气相网格，之后在迭代计算中通过固相颗粒所处小型结构网格来直接获取该颗粒所处的非结构气相网格信息，因而增强了数值模拟的准确性；除此之外，本专利技术无需求解颗粒到计算区域内非结构气相网格中心点之间距离，因此除去了因搜索颗粒所处非结构气相网格而带来的计算量，从而又大大提高了计算的速度。附图说明图1是气相计算所划分的非结构网格和及其所覆盖的大型结构网格示意图。图2是气相计算所划分的非结构网格和及其所覆盖的小型结构网格示意图。图3是迭代计算前搜索各个小型结构网格所处的非结构气相网格示意图。图4是数值计算中由离散颗粒所处的小型结构网格来确定该颗粒所处非结构气相网格示意图。具体实施例方式本专利技术的快速获取气固两相数值模拟中固相颗粒所处非结构气相网格信息的方法，包括以下步骤(1)在气相计算所划分的非结构网格中，分别建立大型、小型两套结构网格对整个计算区域加以覆盖。其中，大型结构网格尺度大于非结构气相网格尺度，其尺度越大，则建立的大型结构网格个数越少，单个大型结构网格所含的非结构气相网格个数越多，因而下文步骤(3)中求解小型结构网格的中心点与其所处的大型结构网格内所有非结构气相网格中心点距离的次数就越多。另外，小型结构网格尺度要小于非结构气相网格尺度，其尺度越小， 则建立的小型结构网格个数越多，单个非结构气相网格所含的小型结构网格的个数越多， 因而下文步骤(4)获取颗粒所处的非结构气相网格信息越精确，但是下文步骤(3)中求解小型结构网格的中心点与其所处的大型结构网格内所有非结构气相网格中心点距离的次数就越多；(2)扫描计算区域内的所有非结构气相网格，将每个大型结构网格所覆盖的非结构气相网格的序号进行统计注册，建立数组以备后继计算使用；(3)确定每个小型结构网格所在的大型结构网格，根据步骤(2)中建立的数组，求解小型结构网格的中心点与其所处的大型结构网格内所有非结构气相网格中心点的距离，其中最小距离者即为该小型结构网格所处的非结构气相网格，对每个小型结构网格所属的非结构气相网格序号进行统计注册，建立数组以备后继计算使用。当小型结构网格位于大型结构网格边沿附近时，对于二维计算，求解小型结构网格与非结构气相网格中心点距离的范围扩大到该小型结构网格所处的大型结构网格相邻的8个大型结构网格。而对于三维计算，则应扩大到该小型结构网格所处的大型结构网格相邻的26个大型结构网格。(4)在每个时间步长的计算中，根据颗粒位置确定其所处的小型结构网格，由步骤(3)中建立的数组从而得到该颗粒所处的非结构气相网格信息。下面结合具体实施例对本专利技术做进一步的解释。实施例1气相流场的数值模拟通常采用二维或三维进行，现以二维情形进行叙述，三维情形可依此类推(1)对于气相计算所划分的非结构网格，首先分别建立大型、小型两套结构网格对整个计算区域加以覆盖。其中，大型结构网格尺度大于非结构气相网格尺度，小型结构网格尺度小于非结构气相网格尺度，如图1、2所示。椭圆区域为气相计算所划分的非结构网格求解区域，大、小两种正方形为所建立的两套结构网格覆盖区域；(2)扫描计算区域内的所有非结构网格，将每个大型结构网格所覆盖的非结构气相网格本文档来自技高网...

【技术保护点】
１．一种快速获取固相颗粒所处非结构气相网格信息的方法，其特征在于该方法包括以下步骤：步骤１）在气相计算所划分的非结构网格中，分别建立大型、小型两套结构网格对整个计算区域加以覆盖；其中，大型结构网格尺度大于非结构气相网格尺度，小型结构网格尺度小于非结构气相网格尺度；步骤２）扫描计算区域内的所有非结构气相网格，将每个大型结构网格所覆盖的非结构气相网格的序号进行统计注册，建立数组以备后继计算使用；步骤３）确定每个小型结构网格所在的大型结构网格，根据步骤２）中建立的数组，求解小型结构网格的中心点与其所处的大型结构网格内所有非结构气相网格中心点的距离，其中最小距离者即为该小型结构网格所处的非结构气相网格，之后对每个小型结构网格所属的非结构气相网格序号进行统计注册，建立数组以备后继计算使用；步骤４）　在每个时间步长的计算中，根据颗粒位置确定其所处的小型结构网格，由步骤３）中建立的数组从而得到该颗粒所处的非结构气相网格信息。

【技术特征摘要】

【专利技术属性】
技术研发人员：袁竹林，朱立平，
申请(专利权)人：东南大学，
类型：发明
国别省市：84