CN101964024A - 一种快速确定固相颗粒所在气相非结构网格的方法 - Google Patents
一种快速确定固相颗粒所在气相非结构网格的方法 Download PDFInfo
- Publication number
- CN101964024A CN101964024A CN 201010512932 CN201010512932A CN101964024A CN 101964024 A CN101964024 A CN 101964024A CN 201010512932 CN201010512932 CN 201010512932 CN 201010512932 A CN201010512932 A CN 201010512932A CN 101964024 A CN101964024 A CN 101964024A
- Authority
- CN
- China
- Prior art keywords
- structured grid
- grid
- structured
- particle
- unstructured
- 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种快速确定固相颗粒所在气相非结构网格的方法。该方法包括:(1)在气相计算所划分的非结构网格中,另建立一套结构网格对整个计算区域加以覆盖;其中,结构网格尺度大于非结构网格尺度;(2)扫描计算区域内的所有非结构网格,将每个结构网格所覆盖的非结构网格的序号进行统计注册,建立数组以备后继计算使用;(3)确定需要搜寻所在非结构网格的离散颗粒所在的结构网格,求解离散颗粒与该结构网格内所有非结构网格节点的距离,其中最小距离者即为该颗粒所处的非结构网格。本发明的方法,大大提高数值模拟的准确性,能显著减少颗粒所处非结构网格的搜索次数,从而大大提高了气固两相数值计算的速度和精度。
Description
技术领域
本发明属于气固两相流动计算机数值模拟领域,具体涉及一种快速确定固相颗粒所在气相非结构网格的方法。
背景技术
气固两相流动广泛存在于能源、化工、环保等领域,随着计算机速度的大幅提高,气固两相流动的计算机数值模拟技术也随之快速发展,并在生产实际中发挥出越来越重要的作用。
在气固两相流动计算机数值模拟中,一种行之有效的措施是将气相场和固相颗粒场分开处理,分别采用欧拉方法和拉格朗日方法进行数值模拟,这样可以充分提高计算效率和增加数值模拟的准确性。在计算离散颗粒所受的气体曳力以及颗粒对气相流场反作用的求解过程中,必须获得颗粒所处当地气相网格的信息,如当地气体流速等。对于复杂外形几何体的非结构网格,如何准确、快速判断颗粒所处的当前网格,是一个难题。目前大多采用对多个非结构网格进行数据平均的方法加以解决,该过程将不可避免地牺牲了大量气相场信息,降低了计算精度。若能找到一种快速确定固相颗粒所在气相非结构网格的方法,则可使数值模拟的效率大大提高。
发明内容
发明目的:针对现有技术中存在的不足,本发明的目的是提供一种快速确定固相颗粒所在气相非结构网格的方法,能显著减少颗粒所处非结构网格的搜索次数,从而大大提高了气固两相数值计算的速度和精度。
技术方案:为了实现上述发明目的,本发明采用的技术方案如下:
一种快速确定固相颗粒所在气相非结构网格的方法,包括以下步骤:
(1)在气相计算所划分的非结构网格中,另建立一套结构网格对整个计算区域加以覆盖;其中,结构网格尺度大于非结构网格尺度;
(2)扫描计算区域内的所有非结构网格,将每个结构网格所覆盖的非结构网格的序号进行统计注册,建立数组以备后继计算使用;
(3)确定需要搜寻所在非结构网格的离散颗粒所在的结构网格,求解离散颗粒与该结构网格内所有非结构网格节点的距离,其中最小距离者即为该颗粒所处的非结构网格。
气相流场的数值模拟通常采用二维或三维进行,当为二维情形进行时,步骤(3)中,当离散颗粒位于结构网格边沿附近时,求解离散颗粒与非结构网格节点距离的范围扩大到该结构网格相邻的8个结构网格。
当为三维情形进行时,步骤(3)中,当离散颗粒位于结构网格边沿附近时,求解离散颗粒与非结构网格节点距离的范围扩大到该结构网格相邻的26个结构网格。
有益效果:本发明的快速确定固相颗粒所在气相非结构网格的方法,结合邻域搜索法以及计算颗粒-网格节点最小距离法,对非结构网格无需进行平均处理,而是通过高效分割搜索直接确定颗粒所处的非结构网格,大大提高数值模拟的准确性,能显著减少颗粒所处非结构网格的搜索次数,从而大大提高了气固两相数值计算的速度和精度。
附图说明
图1是气相计算所划分的非结构网格和及其所覆盖的结构网格示意图。
图2是数值计算中搜索离散颗粒所处非结构网格示意图。
具体实施方式
下面结合具体实施例对本发明做进一步的解释。
实施例1
气相流场的数值模拟通常采用二维或三维进行,现以二维情形进行叙述,三维情形可依此类推:
(1)对于气相计算所划分的非结构网格,首先另建立一套结构网格对整个计算区域加以覆盖,如图1所示,椭圆区域为气相计算所划分的非结构网格求解区域,四边形为所建立的结构网格覆盖区域;(2)扫描计算区域内的所有非结构网格,将每个结构网格所覆盖的非结构网格的序号进行统计注册,建立数组以备后继计算使用;(3)对于需要搜寻所在非结构网格的离散颗粒,首先确定其所在的附加结构网格,然后求解颗粒与该结构网格内所有非结构网格节点的距离,其中最小距离者即为该颗粒所处的非结构网格。(4)考虑到当颗粒位于结构网格边沿附近时,其所处的非结构网格有可能在结构网格的相邻网格内产生,如图2所示,离散颗粒b所处1号非结构网格的节点位于该颗粒所处结构网格A的相邻结构网格B中。将上述(3)中所需求解颗粒与非结构网格节点距离的范围扩大到该结构网格相邻的8个结构网格(对于三维计算,扩大到相邻的26个结构网格区域)。
本发明的方法利用所建立的结构网格,快速确定颗粒所处的空间位置,并使得对非结构网格的搜索范围大大缩小,以2,000,000个非结构网格计算空间为例,如果采用的结构网格数为2,000个,平均而言,每个结构网格覆盖1,000个非结构网格,则计算颗粒到非结构网格节点距离的次数可从计算2,000,0000次减少到9,000次(三维情况下减少到27,000次),从而大大提高了计算效率。在实施上述方法过程中,要求结构网格尺度大于非结构网格尺度。总体而言,结构网格覆盖的非结构网格数目越少,计算颗粒到非结构网格节点距离的次数越少,所花费的计算时间也越短。
Claims (3)
1.一种快速确定固相颗粒所在气相非结构网格的方法,其特征在于,包括以下步骤:
(1)在气相计算所划分的非结构网格中,另建立一套结构网格对整个计算区域加以覆盖;其中,结构网格尺度大于非结构网格尺度;
(2)扫描计算区域内的所有非结构网格,将每个结构网格所覆盖的非结构网格的序号进行统计注册,建立数组以备后继计算使用;
(3)确定需要搜寻所在非结构网格的离散颗粒所在的结构网格,求解离散颗粒与该结构网格内所有非结构网格节点的距离,其中最小距离者即为该颗粒所处的非结构网格。
2.根据权利要求1所述的快速确定固相颗粒所在气相非结构网格的方法,其特征在于:步骤(3)中,当离散颗粒位于结构网格边沿附近时,对于二维计算,求解离散颗粒与非结构网格节点距离的范围扩大到该结构网格相邻的8个结构网格。
3.根据权利要求1所述的快速确定固相颗粒所在气相非结构网格的方法,其特征在于:步骤(3)中,当离散颗粒位于结构网格边沿附近时,对于三维计算,求解离散颗粒与非结构网格节点距离的范围扩大到该结构网格相邻的26个结构网格。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105129326A CN101964024B (zh) | 2010-10-20 | 2010-10-20 | 一种快速确定固相颗粒所在气相非结构网格的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105129326A CN101964024B (zh) | 2010-10-20 | 2010-10-20 | 一种快速确定固相颗粒所在气相非结构网格的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101964024A true CN101964024A (zh) | 2011-02-02 |
CN101964024B CN101964024B (zh) | 2012-05-09 |
Family
ID=43516895
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105129326A Expired - Fee Related CN101964024B (zh) | 2010-10-20 | 2010-10-20 | 一种快速确定固相颗粒所在气相非结构网格的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101964024B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036126A (zh) * | 2014-06-04 | 2014-09-10 | 广东电网公司电力科学研究院 | 肺泡收缩和扩展过程中颗粒物运动的数值模拟方法和系统 |
CN112784401A (zh) * | 2020-12-31 | 2021-05-11 | 南京玻璃纤维研究设计院有限公司 | 一种纤维预制体cvi致密化过程的仿真方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101532906A (zh) * | 2009-04-27 | 2009-09-16 | 东南大学 | 风力发电机叶片的流体动力学和结构力学分析方法 |
US20100112490A1 (en) * | 2006-07-20 | 2010-05-06 | Zhijian Lu | Optical Diffusers, Photomasks and their Methods of Fabrication |
-
2010
- 2010-10-20 CN CN2010105129326A patent/CN101964024B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100112490A1 (en) * | 2006-07-20 | 2010-05-06 | Zhijian Lu | Optical Diffusers, Photomasks and their Methods of Fabrication |
CN101532906A (zh) * | 2009-04-27 | 2009-09-16 | 东南大学 | 风力发电机叶片的流体动力学和结构力学分析方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036126A (zh) * | 2014-06-04 | 2014-09-10 | 广东电网公司电力科学研究院 | 肺泡收缩和扩展过程中颗粒物运动的数值模拟方法和系统 |
CN112784401A (zh) * | 2020-12-31 | 2021-05-11 | 南京玻璃纤维研究设计院有限公司 | 一种纤维预制体cvi致密化过程的仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101964024B (zh) | 2012-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107392875A (zh) | 一种基于k近邻域划分的点云数据去噪方法 | |
KR102044354B1 (ko) | 데이터베이스 구축을 위한 자기장 맵을 생성하는 방법 및 장치 | |
CN108834077B (zh) | 基于用户移动特性的跟踪区划分方法、装置及电子设备 | |
CN103220777A (zh) | 一种移动设备定位系统 | |
CN105973246A (zh) | 一种地磁地图的绘制方法、装置及机器人 | |
CN107703480A (zh) | 基于机器学习的混合核函数室内定位方法 | |
CN104869541A (zh) | 一种室内定位追踪方法 | |
CN104038901A (zh) | 一种减少指纹数据采集工作量的室内定位方法 | |
CN104036095A (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
CN105101090B (zh) | 一种环境监测无线传感网的节点定位方法 | |
CN102890703A (zh) | 一种网络异质多维标度方法 | |
CN105404898A (zh) | 一种松散型点云数据分割方法和设备 | |
JP5930551B2 (ja) | 滞留点抽出方法、滞留点抽出装置及び滞留点抽出プログラム | |
CN105334497A (zh) | 一种三维定位方法和装置 | |
CN103778191A (zh) | 一种顾及空间邻近关系的矢量等高线数据划分方法 | |
CN101964024B (zh) | 一种快速确定固相颗粒所在气相非结构网格的方法 | |
Huang et al. | The analysis of anchor placement for self-localization algorithm in wireless sensor networks | |
CN107506572B (zh) | 获取目标点的高度的方法和装置 | |
CN104778260B (zh) | 一种动态雷达环境知识库建模方法 | |
CN107613502B (zh) | 一种传感器网络非规则区域节点定位方法及其装置 | |
CA2894863A1 (en) | Indoor localization using crowdsourced data | |
CN106251364A (zh) | 目标跟踪方法及装置 | |
Li et al. | An indoor positioning method based on RSSI probability distribution | |
CN102156806A (zh) | 一种快速获取固相颗粒所处非结构气相网格信息的方法 | |
CN113626731A (zh) | 基于人工智能的位置预测的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120509 Termination date: 20141020 |
|
EXPY | Termination of patent right or utility model |