之前一直没弄懂为什么很多流化床模拟里面都不去提及的问题。这还是个比较复杂的问题,暂且搁置。但是我们起码要知道:加和不加湍流模型,区别到底体现在哪?这个影响大不大,量级如何?

从工程上来说,我们只要知道这个影响的量级,就可以大致判定可不可以忽略该影响。假如同时存在其他主导规律的时候(比如颗粒对气体的曳力),该影响的量级远远小于曳力的量级,那么就可以放心地忽略该影响。

所以湍流的影响是多大呢?如何估算他的影响的量级呢?

首先,我们必须从CFD计算的NS方程上找源头。也就是说,要知道加不加湍流,在NS方程上的区别

这里就直接摘抄教科书里的内容了。
所采用的教科书为Versteeg编写的 An Introduction to Computational Fluid Dynamics: The Finite Volume Method 第二版
参考章节为3.5节,从第61页开始。

首先写出NS方程。

第一个是连续方程,后三个是分解成xyz三个方向的动量方程。

我们知道,流体在湍流状态下是有无规则的脉动运动的。但是经过时间平均后,会消除这些脉动。如图所示为流体速度随时间变化。

可以将其分解为脉动速度(u‘)和时均速度(U)

上式中,也包含压力的分解。压力也是同理的。

于是将上面分解后的式子代入到NS方程中去。

然后再对他们取时间平均。取平均的方法就是先对时间积分,然后再除以时间。

对于连续方程,时间平均后直接把脉动速度平均没了。于是得到

对于动量方程,先以x方向动量方程为例,对每一项都采取时间平均

这里要注意几点

  1. 首先,单独对脉动速度平均,直接就平均没了。 u′‾=0\overline{u'}= 0 u′=0
  2. 其次,对平均速度平均,相当于没平均。 U‾=U\overline{U}= U U=U
  3. 然后,对平均速度和脉动速度平均,可以把U看作常数提出来 Uu′‾=Uu′‾\overline{Uu'}=U \overline{u'} Uu′=Uu′
  4. 最后,两个脉动速度的乘积的平均,是没法有任何简化的。 u′u′‾=u′u′‾u′v′‾=u′v′‾\overline{u'u'}= \overline{u'u'}\\ \overline{u'v'}= \overline{u'v'} u′u′=u′u′u′v′=u′v′

另外,div代表散度,即 div(u)=∂ux+∂vy+∂wzdiv(\bold u)=\frac{\partial u}{x}+\frac{\partial v}{y}+\frac{\partial w}{z} div(u)=x∂u​+y∂v​+z∂w​ grad表示梯度,即 div(T)=(∂Tx,∂Ty,∂Tz)div(\bold T)=(\frac{\partial T}{x}, \frac{\partial T}{y},\frac{\partial T}{z}) div(T)=(x∂T​,y∂T​,z∂T​)

于是带入到NS方程中,得到

那么(III)就是额外增加的项,称为雷诺应力。
所以
湍流对NS方程的影响就是附加了一个雷诺应力项!

把它看作一个附加在动量方程的源项,相当于附加了一个额外的作用力。因此通常把它放到右侧和压力以及粘性力并列:

最后那项就是多出来的雷诺应力
这里又把密度提出来了。是因为通常来说密度也是不恒定的,也是一个变量。

所以说,有时候采用密度加权的时间平均,又被称之为favre平均。这个平均用上方小波浪线代替小横线。最后写出来,考虑了湍流之后的NS方程为:

【流体力学】加和不加湍流模型在NS方程上的体现相关推荐

  1. 湍流参数计算c语言,常用的湍流模型

    原文作者:Shawn Wasserman,ENGINEERING.com网站仿真编辑. 1 为正确的湍流模型使用正确的CFD模型 一个虚构的故事.海森堡说"当我遇到上帝的时候,我会问他两个问 ...

  2. @RequestParam加与不加的区别

    最简单的两种写法,加或不加@RequestParam注解 1 @RequestMapping("/list") 2 public String test(int userId) { ...

  3. href 里面 链接前面加/与不加的区别?(绝对路径与相对路径)

    在写href链接时,有绝对路径与相对路径,href 里面 链接前面加/与不加的区别? href="/cp/images/lis.jpg" 相对路径 cp前面/会获取当前路径,组合成 ...

  4. 【深度学习】Keras加载权重更新模型训练的教程(MobileNet)

    [深度学习]Keras加载权重更新模型训练的教程(MobileNet) 文章目录 1 重新训练 2 keras常用模块的简单介绍 3 使用预训练模型提取特征(口罩检测) 4 总结 1 重新训练 重新建 ...

  5. 深度优先搜索 和问题 简单函数递归 “加 还是不加”

    Foj 1106 题目:Sum of Factorials 代码: #include<cstdio> #include<iostream> #include<stdlib ...

  6. 超图桌面版加载obj 3D模型 - 2

    在 https://blog.csdn.net/bcbobo21cn/article/details/109041525 里,加载obj格式模型没有出来效果: 下面来看一下其他方法:当前用的版本是10 ...

  7. 超图桌面版加载obj 3D模型

    根据目前看到的资料帮助,使用超图桌面版加载obj 3D模型的步骤如下: 新建场景,图层管理器---普通图层右键,新建KML----开启KML可编辑,右键添加模型. 下面来操作一下: 下一个obj模型: ...

  8. opengl加载显示3D模型UC类型文件

    opengl加载显示3D模型UC类型文件 前置条件 注意事项 项目展示 项目完整C++源代码 前置条件 opengl三方库freeglut,显示uc格式的三维模型文件, 注意事项 源代码 model_ ...

  9. opengl加载显示3D模型3d类型文件

    opengl加载显示3d模型UC类型文件 前置条件 注意事项 项目展示 项目完整C++源代码 前置条件 opengl三方库freeglut,显示3d格式的三维模型文件, 注意事项 源代码 model_ ...

  10. opengl加载显示3DS模型3DS类型文件

    opengl加载显示3DS模型UC类型文件 前置条件 注意事项 项目展示 项目完整C++源代码 前置条件 opengl三方库freeglut,显示3DS格式的三维模型文件 注意事项 源代码 model ...

最新文章

  1. JIT Code Generation代码生成
  2. 由浅入深理解java集合(二)——集合 Set
  3. python哨兵循环_Python通用循环的构造方法实例分析
  4. 广域网设备:远程路由器/交换机/服务器—Vecloud
  5. Python基础教程【读书笔记】 - 2016/7/5
  6. PyQt5 笔记4 -- 多线程
  7. Sping5——响应式编程
  8. [转]C#操作XML方法详解
  9. 离线使用iPhone SDK文档的方法
  10. '800a0005' 图片上传出现写入文件失败的错误 -- 修改pload_5xsoft.inc
  11. java 内存调优_JVM内存模型以及性能调优
  12. 通过缓存数据库结果提高PHP性能
  13. java 执行oracle命令_利用oracle存储过程执行操作系统命令
  14. redis 学习笔记(6)-cluster集群搭建
  15. Oracle P6培训系列:17建立基线
  16. 网管员不得不关注的九大热门网络技术
  17. css内联样式!important
  18. android密度计算器,密度计算器
  19. 服务器2012系统登录密码忘记6,Server2012忘记管理员密码的处理方法
  20. Smith预估控制实验 simulink模型

热门文章

  1. java 汉字区位码表_汉字编码解析
  2. 汉字区位码转换为“汉字ASCII码“
  3. 花了一天的时间给粉丝做了一个小米官网(高仿)
  4. 云计算平台是什么意思?可以划分为哪三类?
  5. 达摩院的地球云计算平台AI Earth使用体验
  6. VB 连接加密码的ACCESS数据库*.mdb,*.accdb
  7. directx9 截屏
  8. GPS之家论坛最新地图下载精选 汇集论坛精华 不断更新中(2013 3 30)
  9. 使用Java完成Socket文件传输
  10. Optimus双显卡笔记本上用MediaCoder转换iPhone/iPod 4视频