流体仿真中,六面体(Hex)网格与四面体(Tet)网格的争论一直伴随着整个CFD的发展过程,坊间也流传着许许多多关于六面体网格、结构化网格、四面体网格、甚至是Cutcell网格等相关内容的种类繁多的观点。由于CFD的工程性质,很多问题并不需要工程师进行过于深入、过于细致的研究,大家往往还是会把关注的重点放在仿真问题的求解之中,不管六面体还是四面体,只要能解决问题,就都是好网格!

但是值得一提的是,对于绝大多数的流体仿真工作者,都有一个共识,那就是:六面体(结构化)的网格好!

具体来讲,六面体网格(Hex)到底好在哪呢?大概会得到以下几种答案:

● 六面体网格求解效率高

● 六面体网格数量少

● 六面体网格精度高、算的准

● 六面体网格中没有四面体那样的伪扩散

● 只要是六面体网格就一定好用,cutcell也是六面体的一种

● 六面体网格守恒性好

● 四面体属于忽悠人的网格,不能用

● ……

不论对错,关于六面体网格的观点还是非常多、非常明确的;但由于篇幅的关系,我们可能没有办法一一的对上述观点进行分析和讨论;所以,本文仅对上述观点中的第一条进行相对详细的分析和介绍,那就是:在公平的条件下,六面体网格和四面体网格,谁的求解效率高?

问题的分解与明确

● 公平的条件如何界定?

● 求解效率高是通过什么方式体现的呢?

采用案例测试的方式进行,使用相同的计算机和相同的算例设置,仅网格不同(采用replace方式更换),初步估计稳态计算1000个迭代步,通过计算的真实时间来判定求解效率。

● 求解之前的内容(如网格划分、网格输出、求解设置、算例保存、网格替换等),暂不列入判定范围

● 仅比较纯六面体网格和纯四面体网格在 Fluent 中的计算效率

通常的流体仿真工作是必须包含边界层的,因为纯四面体在近壁面附近求解精度相对较低;不过,本次测试为了避免非四面体的层状网格(Prism)对最终的求解产生影响,所以采用了完全一致的网格类型进行对比。


知识要点复习

★ 六面体网格的优势

● 采用六面体网格可以减小网格数量

● 对于流动方向沿着网格的问题,六面体网格可以减小数值误差

通常情况下,六面体网格数数量更少

★ 六面体网格的劣势

● 几何形状复杂时工作量极大(能力 & 耐心)

● 几何曲率较大时,网格质量难以有效控制

需要注意的是:有一些相对复杂的几何模型在生成六面体网格的时候会产生质量很差的网格,而且只能通过网格加密才能有效改善。

随着网格的加密,六面体网格的数量反而超过了自动生成的四面体数量……

当网格方向与流动方向对齐的时候,六面体可以极大的减小数值误差

当网格方向与流动方向不一致的时候,六面体的数值精度与四面体一致


公平的条件

比赛之前,必须要指定规则;同样的道理,六面体与四面体的网格对比也需要一个公平的条件,那么这个公平的条件指的是什么呢?有三个备选选项:

● 网格边 (edge) 的长度一致

● 网格面 ( face ) 的面积一致

● 网格单元的体积 (cell) 一致

★ 边 (edge) 的长度一致

通常情况下,对于边(edge)长度相等的hex与tet ,体网格数量之比约为1:5~1:8

★ 面 ( face ) 的面积一致

当三角形与四边形面积相等时,边长的比例约为1.3~1.5

★ 体积 (cell) 一致

体积一致,也就意味着六面体与四面体的体网格数量相等

为了能够确认这个公平的条件,我们还专门进行了一些问卷的调查,结果如下:

总体来讲,大家的意见还是相对分散的,为了公平起见,我们将选择多数工程师的意见“面 ( face ) 的面积一致”作为相对公平的评判标准。

表面网格(face)面积一致的纯六面体网格与纯四面体网格


求解的结果

对比条件

六面体

四面体

表面网格数量(一致)

105000

105532

Cell

1.25M

3.04M

Face

3.80M

6.22M

Node

1.31M

0.54M

.msh 文件大小

116 MB

131 MB

峰值内存占用

5.85 GB

9.63 GB

1000步计算时间

10637 s

11487 s

当面网格数量一致时,我们可以得到以下结论:

● Hex 与 Tet 体网格数量比约为 1:2.5

● Hex 与 Tet 网格文件大小比例约为 1:1.13

● Hex 与 Tet 求解效率比例约为 1.08:1

● Tet 占据内存更大

因此,在相对“公平”的条件下,六面体网格的求解效率可以认为和四面体一致!!

相信对于大多数流体仿真工程师来讲,这是一个颠覆性的结论;当然,也有很多人会说上面所谓“公平”的前提条件并不完全公平,所以结论并不合理;无论如何,这些都是可以接受的观点。

还是引用作者最常用的一句话,文章的目的并不是要否定别人的观点、也并不希望一定能够说服大家,而仅仅是从某些特定的角度,提供一些信息,希望能够引发大家的思考。

一些简单的结论:

● Fluent 采用的是有限体积法,因此,体网格的数量在很大程度上代表了计算区域的离散精细程度,从这个角度来讲,网格数量少,也不一定是优势。

● 简单的讲,当六面体与四面体体网格数量比例是1:2.3时,两者的计算效率(求解速度)是一致的,但四面体求解时占用内存更高(约1.6倍)。

● 仿真分析要有理有据,提出观点不能人云亦云。

3 设置网格数的大小_流体仿真中,六面体(Hex)网格的求解效率真的比四面体(Tet)高”很多”么?...相关推荐

  1. 流体仿真中,六面体(Hex)网格的求解效率真的比四面体(Tet)高”很多”么?

    作者 | 张杨 流体仿真中,六面体(Hex)网格与四面体(Tet)网格的争论一直伴随着整个CFD的发展过程,坊间也流传着许许多多关于六面体网格.结构化网格.四面体网格.甚至是Cutcell网格等相关内 ...

  2. comsol移动网格_移动网格技术在计算流体动力学数值仿真中的应用

    流体流动现象广泛存在于自然界和各种工程领域中,所有这些过程都要受质量守恒.动量守恒.能量守恒等基本物理定律的支配,即要满足一定的控制微分方程[.计算流体动力学(computational fluid ...

  3. java object大小_在Java中,确定对象object大小的最佳方法是什么? - Break易站

    例如,假设我有一个应用程序,它可以读取包含一堆数据行的CSV文件.我根据数据类型向用户提供行数摘要,但是我想确保不读取太多数据行并导致OutOfMemoryErrors.每一行都会翻译成一个对象.有没 ...

  4. 前端flv.js设置缓冲时间和大小_不可不知的 WEB 前端网站优化 —— 雅虎 34 条军规...

    不可不知的 WEB 前端网站优化 -- 雅虎 34 条军规 不得不说现在依然适用于大部分的网站 当年雅虎推荐了一套优化网站加载速度的34条法则(包括Yslow规则22条),以下是详细说明. 1. Mi ...

  5. css grid随页面大小_前端开发中各种设置CSS间距的优点缺点及实例「实践」

    前言 如果两个或多个元素很接近,那么用户就会认为它们以某种方式属于彼此.当对多个设计元素进行分组时,用户可以根据它们之间的空间大小来决定它们之间的关系.没有间距,用户将很难浏览页面并知道哪些内容相关而 ...

  6. qt的qprinter打印功能怎么设置打印时字体大小_如何打印一份成绩单

    学校考试之后,要看自己班的成绩,要打印自己班的成绩. 需要进行以下操作,设置好后,只打印自己班的那页就可以了. 惯用WPS,本文一些方法可能仅适用WPS.因为EXCEL的打印预览我没有掌握. 一.冻结 ...

  7. 前端flv.js设置缓冲时间和大小_好程序员web前端细解cookie那些事

    好程序员web前端细解cookie那些事,在互联网时代,IT行业飞速发展,带动了web前端开发行业的兴趣.由于行业新兴起时间不久,专业人才缺乏,薪资待遇较高,已成为众多IT学子选择就业的首选,今天就为 ...

  8. python怎么设置画布大小_在Matplotlib中设置绘图画布的大小

    这是我对Matplotlib最大的挫折之一.我经常使用光栅数据,例如我想在其中添加彩色地图.图例和一些标题.matplotlib库中的任何简单示例都会导致分辨率不同,因此会重新采样数据.尤其是在做图像 ...

  9. 3模型大小_分布式训练中数据并行远远不够,「模型并行+数据并行」才是王道...

    选自arXiv 作者:Saptadeep Pal等 机器之心编译参与:魔王.杜伟 数据并行(DP)是应用最广的并行策略,对在多个设备上部署深度学习模型非常有用.但该方法存在缺陷,如随着训练设备数量不断 ...

最新文章

  1. NBT:牛瘤胃微生物组的参考基因组集
  2. 如何建立JSP操作用以提高数据库访问效率
  3. 四边形可以分为几类_学习知识:四边形有几种类型
  4. Java线程的挂起与恢复 wait(), notify()方法介绍
  5. window.open 不显示地址栏_谷歌浏览器Chrome显示「由贵单位管理」怎么解决?非阿里...
  6. c# xml的增删改查操作 xmlDocument 的用法
  7. Difference between a normal BSP application and a WebClient UI component
  8. ZOJ-Crashing Balloon
  9. lucene 多索引目录搜索实现方法
  10. IntelliJ IDEA的安装详解
  11. SegmentFault D-Day 2016 年度总结报告
  12. Lyncee 数字全息显微镜 DHM Digital Holographic Microscopy
  13. pythonopencv算法_OpenCV算法精解:基于Python与C++
  14. cpu多开测试软件,教你用多核CPU多开畅玩大型3D游戏
  15. spring data jpa 使用@Query 不确定参数查询
  16. 运维日志审计是什么意思?用什么工具好?
  17. 用html做自我介绍
  18. Objective C转Swift注意事项(一)合理使用结构体,枚举,extensions
  19. 手机版云盘怎么打开连接服务器,手把手教你简易NAS构建,手机/平板/智能电视随意调取,家庭存储云共享,有了自己的网络云盘后再也不用担心容量不够了!...
  20. rop检查_我科成功实施首例全麻下小儿眼底荧光造影检查!

热门文章

  1. 2.Java 面试题整理(基础篇二)
  2. Spark基础学习笔记08:Scala简介与安装
  3. 【codevs3110】一把鼻涕一把泪的堆排序
  4. english 2012020601
  5. php秒数怎么转换为时间格式,php将秒数转换为时间(年、月、日、小时…)
  6. Intel 64/x86_64/IA-32/x86处理器 - SIMD指令集 - MMX技术(5) - 逻辑指令
  7. DSP编程的一些优化技巧和注意事项
  8. keyshot局部打光_keyshot7耳机渲染打光教程
  9. sysbench 1.0.6 mysql_mysql sysbench 1.0.X
  10. 使用Nsight查找CE3的渲染bug