CFD分析中如何选择有效网格系统

《CFD分析中如何选择有效网格系统》由会员分享，可在线阅读，更多相关《CFD分析中如何选择有效网格系统（16页珍藏版）》请在金锄头文库上搜索。

1、Mechanical analysis W h i t e P a P e r w w w m e n t o r c o m 白皮书 CFD 分析中如何选择有效网格系统 海基科技专业代理同步CFD模拟软件 FloEFD 白皮书 CFD 分析中如何选择有效网格系统 2 16 摘要 同步 CFD 是 CFD 中一类新的工具 它帮助结构工程师在三维结构 CAD 模型中仿真现今产品的流体流动和传热 情况 对于三维仿真和分析而言 最重要的步骤就是网格和创建有效的网格系统 本文讨论了为什么矩形自适应 网格是先进技术 以及如何有效为新设计选择网格 从而极大降低精确分析所需的时间 提高产品设计效率 目录 1

2、 网格的需要和选择 3 1 1 为何首先需要一个网格系统 3 1 2 网格系统如何进行选择 3 2 影响网格系统选择的因素和需要考虑的事项 5 2 1 网格形状对于网格质量的影响 5 2 2 非矩形几何体的描述 6 2 2 1 Patankar 和其同事 参考 2 6 2 2 2 Spalding 和其同事 参考 3 7 2 2 3 NASA Ames 的工作 参考 4 和参考 5 7 2 2 4 剑桥大学的工作 参考 6 8 2 3 网格排列 结构化或非结构化 8 2 4 网格生成 10 3 FloEFD 的矩形自适应网格技术 11 3 1 初始网格 11 3 2 求解自适应网格 13 4

3、参考 14 附录 1 非正交网格中的扩散通量和压力梯度 15 海基科技专业代理同步CFD模拟软件 FloEFD 白皮书 CFD 分析中如何选择有效网格系统 3 16 1 网格的需要和选择 1 1 为何首先需要一个网格系统 在进行任何 CFD 分析之前 考虑所需的网格系统是非常有必要的 所有的 CFD 分析都是建立在控制流体动力学现象的微分方程之上 这些微分方程有 Navier Stokes 方 程 能量守恒方程等 众所周知 这些微分方程是无法获得解析解的 除非进行大量的简化 因此 只有采用 离散化 才能进行求解 通过在整个分析区域上覆盖一个虚拟的网格系统的方式 将所考虑的区域划分成许多小的体积

4、或单元格 对小体积内和小体积之间所考虑特性的变量 速度 压力和温度等 进行假设 因此可以推导得出这些微分控制方程的近似形式 也就是所谓的有限体积法 只要这个体积足够小 这一体积内的控制方程就足够有效 从而在整个区域内的控制方程也足够有效 最后通过迭代的方式求解这些代数方程 从而获得相应的结果 很明显 网格划分是最终获得控制微分方程合理精确解的一种方法 所选择的网格大小和细密程度对求解的精确度有很大影响 网格系统类型的选择 网格的形状和排列可以是任意的 只要定义的网格能方便可靠的获取精确结果 这一 网格就是良好的网格 然而 这一 只要 字眼是非常重要的限定 经验表明 对于任何实际应用 为 CFD

5、 计算选择网格系统时 必须考 虑以下影响因数 定义问题和以后做相应修改所需的时间 易于获得良好 精确结果 解的强壮性和可靠性 计算速度和存储 这就是为什么 CFD 计算网格系统的选择是一项重要的工作 1 2 网格系统如何进行选择 在用于 CFD 分析的网格系统选择时有两个非常重要的方面 1 网格的形状 主要的选择有 笛卡儿 立方体网格 并且网格面与笛卡儿坐标系中的 X Y Z 轴相平行 六面体 六面体网格 是笛卡儿网格的某种扭曲 可以是 笛卡儿网格拓补 也就是类似笛卡儿网格 但 是网格被扭曲 或者 适体网格 通过扭曲笛卡儿网格 使其很好的与物体的表面贴合 四面体 四个面的网格 例如三棱锥形网格

6、 海基科技专业代理同步CFD模拟软件 FloEFD 白皮书 CFD 分析中如何选择有效网格系统 4 16 2 网格的排列 主要的选择有 结构化网格 网格中节点排列有序 邻点间的关系明确 非结构化网格 节点位置无法用一个固定的法则予以有序的命名 部分非结构化网格 partially unstructured 在某一区域内结构化网格与其它结构化网格以某种方式 结合的网格 并非所有的网格形状与网格排列都具有现实意义 最为常用的网格如下 笛卡儿 无论是结构化还是部分结构化都被广泛的应用到 CFD 的诸多领域 六面体网格 结构化和部分结构化 经常用于 适体 常用于 空气动力学 方面的应用 燃气轮机叶 片

7、 机翼 流线型物体 这主要是因为可以将网格很好的贴合在物体表面 完全非结构化六面体和四面体网格 最初被用于有限元 而不是有限体积法 的 CFD 分析 现在被广泛 的用于有限体积法 通常是棱柱或棱锥形式 这些网格如后一页所示 下面利用笛卡儿网格对正交网格进行 进一步的说明 严格来说 许多对笛卡 儿网格所作的注解也可以应用于 正 交 网格 那就是网格线与正交坐标轴 方向对齐 其中坐标轴互相成 90 度 角 在实际使用中 笛卡儿网格最常用 见的正交网格 基于圆柱坐标系的正交网格也比较常 见 但是使用并不普遍 此外 笛卡儿 网格比其它非正交有更多的优势 我们 会在以后的章节中做进一步的讨论 这一白皮书

8、中考虑了诸多可以选择的网 格形状和排列 但主要集中在第一和第 三两种网格 也就是笛卡儿网格和完全 非结构化 六面体和四面体 网格 第 二种网格 结构化四面体 适体网格 是一种介于以上两者之间的方法 仅仅 适用于空气动力学的应用 海基科技专业代理同步CFD模拟软件 FloEFD 白皮书 CFD 分析中如何选择有效网格系统 5 16 2 影响网格系统选择的因数和需要考虑的事项 2 1 网格形状对于网格质量的影响 为什么 笛卡儿网格形状成为许多应用场合的首选 可以方便的在笛卡儿参考系中对控制方程进行推导和明确的表达 求解的速度分量几乎总是和笛卡儿参考系坐标方向对齐 笛卡尔网格比非正交网格具有更高的网

9、格质量 与笛卡儿网格差异 也就是更大角度的扭曲 越大的非正交网 格 其网格质量方面的 降低 也越明显 网格质量是进行 CFD 分析时 选择网格系统所着重考虑的方面 网格形状 特别是正交性网格的扭曲 对于有限 体积法微分方程推导假设和求解结果方法有很大的影响 附录 1 对最核心的问题进行了详细的阐述 其中考虑了有限体积方程中的两个典型项 它们描述了穿过网格面 的扩散通量以及作为某一方向上速度源的压力梯度 对高度非正交网格中的这两项推导进行了推导 最需要注意的一点是 非正交网格会比笛卡儿网格多产生一个 二次 项 附录 1 考虑了两维的情况 对于每一种考虑的方法仅仅出现了两项 然而 在完全三维的情况

10、中 对于 非正交网格的推导会比笛卡儿网格产生几倍的 二次 项 这些 二次 项的出现会产生很多后果 更多的计算时间 二次 项的计算需要耗费更多的计算时间 由于需要很多项将非正交网格描述成类 似笛卡儿网格 所以可能需要几倍的时间 并且由于计算在求解的过程中迭代进行 所以对时间的影响很 大 更多的存储空间 这可能是最主要的影响 通常情况下都要对关于每一个非正交网格主要几何参数进行 存储 而不是连续的进行计算 这就是为什么非结构化的六面体或四面体网格比笛卡儿网格需要更多的计 算存储空间 实际上在大型复杂计算的过程中 这已成为这种方法 非结构化网格 使用的限制 降低精度和减少迭代求解的强壮性 为了计算这

11、些 二次 项引入了辅助的 cross linkages 也就是说 不是仅仅两个位置的温度被用于热流的计算 远处其它位置的温度也会被用于热流的计算 这会有两个后 果 引入额外的错误 这就意味着 在所有其它条件相同的情况下 高度非正交网格要比正交网格的计算 精度低 换而言之 要实现相同的数值计算精度 非正交网格比正交网格需要更细密 第二个影响方面是有限体积方程系统的收敛稳定性 由于在迭代计算过程中几乎无法直接处理 二次 项 所以使它们具有很大的主导性 从而使迭代求解的可靠性变差 可能会出现不可靠的收敛或发散 这些非正交网格的缺点会随着网格扭曲 非正交性 的增大而变得更明显 所以其结果严重的依赖于实

12、际的应用 问题 至此 非正交网格的不利影响已经被阐述 并且很好被了解 这就是为什么 CFD 的使用者尽可能的要采用笛卡儿网格系统 或其它的正交网格系统 非正交网格系统的用户被要求去阻止差质量网格的产生 通常需要对自动生成的网格进行手动的 调整 这成为整个 CFD 分析过程中最为耗时的工作 海基科技专业代理同步CFD模拟软件 FloEFD 白皮书 CFD 分析中如何选择有效网格系统 6 16 2 2 非矩形几何体的描述 如果笛卡儿网格的优点是那么明显 那么 CFD 的使用者为何还要使用非正交网格 这主要是由于复杂系统的需要 特别是那些非矩形的固体边界 正是由于这个原因 非正交网格系统在机翼等物理

13、外形的贴合方面具有很大的优势 它可以使网格面与物理边界 很好的贴合 然而 在过去十年出现了一些不错的新方法 其中就采用笛卡儿网格 并且非矩形固体形状可以以任意形式穿过 网格 在网格中出现的固体采用合适的 cut cell 技术进行描述 这种方法的优点 可以确保良好的网格质量 具体的方面先前已经阐述 可以避免在自动生成网格之后 再进行手动调整 对于耦合热交换问题 包扩固体区域内存在流动的导热和流体的热交换 常出现在电子散热领域 由于需 要进行耦合求解 很自然网格系统会覆盖流体和固体区域 对于复杂几何外形地问题有不少相关经验 以下引用了四个相关地例子 2 2 1 Patankar 和其同事 下图显

14、示了 Patankar 和其同事所做地工作 上图证明了使用具有流体 固体网格描述的笛卡儿网格所获得结果的精度 将通过圆柱体 Re 26 的绕流流动 与实验流动结果进行了比较 并且与具有相同网格密度的非正交适体网格所得结果进行了比较 两种网格的计算结果均与实验结果相吻合 通过一些其它的 简单 测试 可以得到相同的结论 采用 Cut cell 技术的笛卡儿网格可以与复杂的非正交适体网格获得一样好的计算结果 海基科技专业代理同步CFD模拟软件 FloEFD 白皮书 CFD 分析中如何选择有效网格系统 7 16 2 2 2 Spalding 和其同事 参考 3 下图显示了 Spalding 和其同事所

15、做地工作 参考 3 上图表明置于风洞中的汽车周围是湍流 流动 在这个例子中 笛卡儿网格被嵌 套 也称之为部分非结构化网格 后续 章节会讨论 所得到的重要结论是 通过笛卡儿网格 所获得汽车表面压力变化的结果与实验 测量值相一致 尽管这里没有显示其它 网格系统的结果 但是笛卡儿网格的结 果足以和其它更复杂的适体网格结果相 媲美 2 2 3 NASA Ames 的工作 参考 4 和参考 5 下图参考了 NASA Ames 的工作 主要是飞机和飞行器外部的空气动力学 上图展示了部分非结构化笛卡儿网格 octree structured 在军事直升机空 气动力学方面的应用 参考4 这一网 格系统也被用于

16、NASA Ames 机翼 整 个飞机机身和航空器周围的流动计算 NASA Ames 还利用嵌套的笛卡儿网格 参考 overset structured grids 对机身周围和后部的流动进行计算 参 考 5 采用这类基于笛卡儿网格的技术可以方 便的 相对而言 生成网格 并且与非正 交的网格系统相比在数值计算方面更具 优势 海基科技专业代理同步CFD模拟软件 FloEFD 白皮书 CFD 分析中如何选择有效网格系统 8 16 2 2 4 剑桥大学的工作 参考 6 Dawes 教授的论文回顾了叶轮机械的 CFD 仿真 着重关注了具有特殊几何形状的应用问题 叶轮机械流动是 仿真模拟方面的一个很大挑战 早期的叶轮机械仿真采用结构化的六面体网格 但这限制了叶轮方面 CFD 进入 到 适体 网格的发展 现今诸多通用型 CFD 软件采用 适体 网格 Dawes 教授认为这阻碍了 CFD 软件的使 用 并且由于生成网格时间的原因限制了 CFD 在设计方面的使用 网格生成的转变势在必行 应该对复杂几何 模型采用转变的方式 而不是直接进行处理 Dawes 教授介绍了在计算机图形方面的最新进展 Level