湍流公式推导系列——(一) 不可压湍动能方程的推导与含义
不可压流体的动量方程如式(1)所示,对其进行系综平均,平均后的方程为式(2),即RANS方程,对式(1)与式(2)做差,得到脉动运动方程 (3)
∂ui∂t+∂uiuj∂xj=−1ρ∂p∂xi+v∂2ui∂xj2(1)\frac{\partial u_{i}}{\partial t}+\frac{\partial u_{i} u_{j}}{\partial x_{j}}=-\frac{1}{\rho} \frac{\partial p}{\partial x_{i}}+v \frac{\partial^{2} u_{i}}{\partial x_{j}^{2}} (1) ∂t∂ui+∂xj∂uiuj=−ρ1∂xi∂p+v∂xj2∂2ui(1)
∂uˉi∂t+∂uˉiuˉj∂xj=−1ρ∂pˉ∂xi+v∂2uˉi∂xj2−∂ui′uj′‾∂xj(2)\frac{\partial \bar{u}_{i}}{\partial t}+\frac{\partial \bar{u}_{i} \bar{u}_{j}}{\partial x_{j}}=-\frac{1}{\rho} \frac{\partial \bar{p}}{\partial x_{i}}+v \frac{\partial^{2} \bar{u}_{i}}{\partial x_{j}^{2}}-\frac{\partial \overline{u_{i}^{\prime} u_{j}^{\prime}}}{\partial x_{j}} (2) ∂t∂uˉi+∂xj∂uˉiuˉj=−ρ1∂xi∂pˉ+v∂xj2∂2uˉi−∂xj∂ui′uj′(2)
∂ui′∂t+uˉj∂ui′∂xj+uj′∂uˉi∂xj=−1ρ∂p′∂xi+v∂2ui′∂xj2−∂∂xj(ui′uj′−ui′uj′‾)(3)\frac{\partial u_{i}^{\prime}}{\partial t}+\bar{u}_{j} \frac{\partial u_{i}^{\prime}}{\partial x_{j}}+u_{j}^{\prime} \frac{\partial \bar{u}_{i}}{\partial x_{j}}=-\frac{1}{\rho} \frac{\partial p^{\prime}}{\partial x_{i}}+v \frac{\partial^{2} u_{i}^{\prime}}{\partial x_{j}^{2}}-\frac{\partial}{\partial x_{j}}\left(u_{i}^{\prime} u_{j}^{\prime}-\overline{u_{i}^{\prime} u_{j}^{\prime}}\right)(3) ∂t∂ui′+uˉj∂xj∂ui′+uj′∂xj∂uˉi=−ρ1∂xi∂p′+v∂xj2∂2ui′−∂xj∂(ui′uj′−ui′uj′)(3)
ui′u_{i}^{\prime}ui′与式(3)相乘后系综平均,并定义湍动能(TKE)为
k=12uk′uk′‾=12(u′2‾+v′2‾+w′2‾)(4)k=\frac{1}{2} \overline{u_{k}^{\prime} u_{k}^{\prime}}=\frac{1}{2}\left(\overline{u^{\prime 2}}+\overline{v^{\prime 2}}+\overline{w^{\prime 2}}\right)(4) k=21uk′uk′=21(u′2+v′2+w′2)(4)
则湍动能输运方程为
∂k∂t+uˉk∂k∂xk=−ui′uk′‾∂uˉi∂xk−∂∂xk(1ρp′uk′‾+k′uk′‾−v∂k∂xk)−v∂ui′∂xk∂ui′∂xk‾(5)\frac{\partial k}{\partial t}+\bar{u}_{k} \frac{\partial k}{\partial x_{k}}=-\overline{u_{i}^{\prime} u_{k}^{\prime}} \frac{\partial \bar{u}_{i}}{\partial x_{k}}-\frac{\partial}{\partial x_{k}}\left(\frac{1}{\rho} \overline{p^{\prime} u_{k}^{\prime}}+\overline{k^{\prime} u_{k}^{\prime}}-v \frac{\partial k}{\partial x_{k}}\right)-v \overline{\frac{\partial u_{i}^{\prime}}{\partial x_{k}} \frac{\partial u_{i}^{\prime}}{\partial x_{k}}}(5) ∂t∂k+uˉk∂xk∂k=−ui′uk′∂xk∂uˉi−∂xk∂(ρ1p′uk′+k′uk′−v∂xk∂k)−v∂xk∂ui′∂xk∂ui′(5)
其中Pk=−u′iu′k‾∂uˉi∂xk{P_k} = - \overline {{{u'}_i}{{u'}_k}} \frac{{\partial {{\bar u}_i}}}{{\partial {x_k}}}Pk=−u′iu′k∂xk∂uˉi为湍动能生成项,表示了由于剪切作用引起能量从平均动能向脉动动能传递;Dk=∂∂xk(1ρp′u′k‾+k′u′k‾−ν∂k∂xk){D_k} = \frac{\partial }{{\partial {x_k}}}\left( {\frac{1}{\rho }\overline {p'{{u'}_k}} + \overline {k'{{u'}_k}} - \nu \frac{{\partial k}}{{\partial {x_k}}}} \right)Dk=∂xk∂(ρ1p′u′k+k′u′k−ν∂xk∂k)为扩散项,表示了湍流的扩散效应; ε=ν∂u′i∂xk∂u′i∂xk‾\varepsilon = \nu \overline {\frac{{\partial {{u'}_i}}}{{\partial {x_k}}}\frac{{\partial {{u'}_i}}}{{\partial {x_k}}}}ε=ν∂xk∂u′i∂xk∂u′i为耗散项,表示由于粘性而引起的湍动能的耗散率。
湍流公式推导系列——(一) 不可压湍动能方程的推导与含义相关推荐
- 湍动能耗散率ε输运方程的推导
湍动能耗散率ε输运方程的推导 一.脉动运动方程 二.湍动能耗散率输运方程的推导 三.湍动能耗散率输运方程各项解析 四.参考资料 一.脉动运动方程 脉动运动方程的推导可以参考博主的另一篇博文<雷诺 ...
- Black Scholes公式推导及求解 Part 1:BS Equation的推导
Black Scholes公式推导及求解 Black Scholes公式推导及求解 Part 1:BS Equation的推导 Black Scholes公式推导及求解 Part 1:BS Equat ...
- 一二三系列之状压DP——Max Correct Set(一)Neko Rules the Catniverse (Large Version)(二)Make It Ascending(三)
文章目录 一:CF1463F 二:CF1152F2 三:CF1342F 一:CF1463F Max Correct Set 有一个结论:以x+yx+yx+y为周期排列填充一定是不劣于最后的答案的 令p ...
- Flink系列-背压(反压)
目录 了解背压 什么是背压 背压产生的原因 背压导致的影响 定位背压 解决背压 了解背压 什么是背压 在流式处理系统中,如果出现下游消费的速度跟不上上游生产数据的速度,就种现象就叫做背压 (backp ...
- FH7333系列LDO稳压IC系列
详细介绍 资料下载 样品申请 概述 FH73xx-1 是一款采用 CMOS 技术的低压差线性稳压器.最大输出电流为 250mA 且允许的最高输入电压为 21V.具有几个固定的输出电压,范围从 3.0V ...
- 【控制】傅里叶系列(一)傅里叶级数 (Fourier series) 的推导
傅里叶级数 傅立叶级数 (Fourier series) 1.把一个周期函数表示成三角级数: 2.麦克劳林公式中的待定系数法: 3.三角函数的正交性: 4.函数展开成傅里叶级数: 傅立叶级数 (Fou ...
- python中哪个函数能生成集合_神奇的python系列11:函数之生成器,列表推导式
1.生成器 生成器的本质是迭代器. 在python中有三种方式来获取生成器 1.通过生成器函数 2.通过各种推到式来实现生成器 3.通过数据的转换也可以获取生成器 #函数 deffunc():prin ...
- 调制解调系列(1) IQ调制(理论推导+工程实现(FM))
IQ调制:正交调制,IQ调制可以做出所有的调制方式. FM:Ufm=UcmCos(Wct+mfsinWst) :三角函数分解得 FM:Ufm=UcmCos(Wct)Cos(mfsinWst) - Uc ...
- Fluent求解器——湍流参数的设置
一.湍流参数的类型 在FLUENT流场仿真的湍流设置中,需要对边界条件的湍流参数进行设定.虽然流动进出口处的湍流设置对于计算结果常常没有较大影响.但某些情况下(流动进出口处的湍流水平接近于边界层内的湍 ...
- OpenFOAM 中边界条件的设定【转载】
转载自:http://blog.sina.com.cn/s/blog_a0b4201d0102v7jt.html 用习惯了FLUENT的操作界面,再使用OpenFOAM就会觉得非常繁琐.遇到的第一个问 ...
最新文章
- 如何自学python基础-零基础如何自学成为Python高手?我有几点干货
- Linux单网卡多个IP(或者多个网卡多个IP)设置
- 数据结构之数组、链表、栈和队列
- MFC标签页控件的使用
- CoreAnimation编程指南(八)事务 转自:http://www.dreamingwish.com/
- C#绘图工具之Move
- 三伏天里小试牛刀andriod 开发 #华为云·寻找黑马程序员#
- (23)System Verilog旗语解决资源共享需求
- 【转】LoadRunner中事务和集合点的放置顺序问题
- 无pygame写一个python贪吃蛇
- net idautomationhc39m条形码字体生成条形码
- 智能供应链预测的应用
- 一文搞懂HTTPProxy丨含基础、高级路由、服务韧性
- 苹果开发者技术支持电话方式改变
- linux centos7下源码 tar安装mysql5.7.22或mysql5.7.20 图文详解
- OpenBSD身份验证绕过和权限提升漏洞
- android tablayout放图片,Android TabLayout的Indicator如何设置为图片
- Java中的正则表达式 regex
- 写给IT自学者的入门指南
- word快捷键被占用
热门文章
- moment解读常用操作及语句——subtract、add、calendar
- 伍伦贡计算机科学硕士申请,伍伦贡大学电脑科学(网络和信息安全)硕士研究生申请要求及申请材料要求清单...
- vscode eslint beautify 格式化 html
- 江南鹤微信公众号文章采集器,开发完成了!以后再也不发愁采集微信公众号文章了!
- 计算机组装防静电措施,浅谈组装电脑如何防静电与去除机箱静电的方法
- 怎么修改teredo服务器,技术员设置win7系统通过teredo连接IPv6的修复方案
- 解决:unable to find valid certification path to requested target
- NodeJS解析前端请求图片链接,将服务器目录下的图片返回给前端用于页面展示
- c++中rand(),srand()使用
- Axure RP 认识