matlab速度梯度,[转载]关于FLUENT中Y+的一些讨论
一、关于 fluent计算时壁面函数法和网格的关系,还有一个小问题
1:各位用 fluent的同仁和高手们,我想要比较好的使用 fluent软件,最重要的就是要学好理 论,在这里我想请教各位一个问题,在使用标准 k-eplison和一些其他的封闭模型时,对于近壁区的流动要使用壁面函数法求解。那么在划分网格时,是不是一定要把把第一个内节点布置在湍流充分发展的区域内呢?我们如果自动生成网格时,如果说第一个节点在壁面的粘性底层内,是不是对计算有一定的影响呢?还有一个问题就是在
gambit中设置的
wall
壁面,怎么到 fluent设置为内部表面 interior,好像在边界条件设置时没有这个边界呀。
2:为什么要用壁面函数??就是因为,k-epsilon模型中,k的 boundary
condition已知,在壁面上为零,而
epsilon的 boundary
condition在壁面上为一未知的非零量,如此如何来解
两方程模型???所以,我们就需要壁面函数来确定至少第一内节点上的值,当然也包括壁面上的值。实际上就是把
epsilon方程的 boundary
condition放到了流体内部。至于壁面函数的应用范围,要看它是如何获得的,简单说,他们都是由于,靠近壁面,雷诺应力在粘性底层内基本消失,所以,navier-stokes变为可解,而求得。所以,凡是应用壁面函数求得的节点,都应设置在粘性底层(y+=5-8)或者至少为线性底层(y+>30?具体数值忘记了),当然你放得越低,精度越高,但是网格越小。我在
matlab内自己写的
code,在 y+=5-8内放
10层,fluent应该可以更高。放在 fully developed
region是完全错误的。
4:二楼的兄弟,谢谢!我的意思是壁面函数法和 k-epsilon混合使用,是不是它只计算壁面到第一个节点线之间的区域?如果是这样的话,划分网格是不是要计算这个距离呢?Y+这个值是我们控制,还是
fluent在求解时自动计算呢?y+的临界值好像是
11.63,不过这个值 不是绝对的。
为什么要使用壁面函数呢?首先,在 CFD中应用湍流模型并不一定需要使用壁面函数,在粘性支层中可以对 N-S方程直接求解。在粘性支层中,速度梯度很大,vorticity不为零,所以要直接求解,就必须在粘性支层中布置较多节点,一般要 10层以上,这就是一般的低 Re数湍流模型。当然这样将占用较多的计算资源。而在边界层中,是存在解析解的,如果
在粘性支层内不求解三维 N-S方程,而用一维数学模型代替,将大大降低计算资源的使用,这就是壁面函数。一般高
Re数湍流模型都使用壁面函数。第一层网格节点布置在粘性支层之外。那么你如何判断你的边界层网格节点布置是否合适呢?这就要检查你的
y+,y+就是第一层网格质心到壁面的无量纲距离,与速度、粘度、剪应力等等都有关系。对于
y+的值,
各个学者推荐的范围是不一样的,但一般在 30-60之内肯定是没有问题的。也有推荐 10-110甚至 200的。y+的值合理,意味着你的第一层边界网格布置比较合理,如果
y+不合理,就要调整你的边界层网格。
二、关于 Y+的几点笔记
y+普遍存在于湍流问题中,Y+是由 solver解出來的結果,网格划分时,底层网格一般布置到对数分布律成立的范围内,即 11.5~30<=y+<=200~400。在计算开始时,y+并不知道,这些值需要在计算过程中加以调整。数值计算实践表明,y+对传热特性的影响比较大,往往存在一个合适的取值范围,在该范围内数值计算结果与实验数据的符合较好。算每个模型都要先大概算一下,然后得到
y+,然后再算第一层高度,重新画网格,貌似像是一个迭代的过程。
估算可以采取以下的方法:
1,专门计算器
最好是使用 NASA Viscous Grid Space
Calculator计算,网址为:
http://geolab.larc.nasa.gov/APPS/YPlus/
2,
3,一些经验之谈
Laminar
sublayer (y+ < 5)
Buffer
region (5 < y+ < 30)
Turbulent
region (y+ > 30)
增强的壁面函数(Enhanced Wall
Treatment)对于
y
plus大于 30的有很好的作用,对于 5到 30的也能求解,但对于 laminar底层,对求解非常敏感,特别是在高雷诺数的情况下。所以对于高雷诺数的湍流流动,可以先选择增强的壁面函数设置,然后计算完以后,通过
计算的 y plus对网格进行修改。
湍流边界层分为内区和外区,而内区有分成:
粘性子层,粘性作用主导
过渡层,湍流作用和粘性作用相当
对数律层,湍流作用主导
外区:惯性力主导
距离壁面最近节点定义了 y+的值,在他们之间应用壁面函数,一般应该使该节点位于粘性子层之外,即
y+>=11.06好像,同时在湍流边界层内应该布置一定数量的网格节点,而
y+的上限由雷诺数判定,我觉得和湍流边界层的厚度有关,因为从经验上,对数律层只能到总的边界层厚度的
20%。
问题是:具体在生成网格时,如何来控制网格满足一定的 y+要求呢?
对大空间的钝体绕流计算,合理的 y+应该在什么范围内?而据我所知,y+是结果算出来之后查看的,如果不是很关心剪切力,湍流边界层需要准确计算吗?一刚开始,壁面网格生成凭经验确定。计算后,根据流场查看
y+值是否满足事前的假定,不满足则需要修改网格,再计算、查看y+,如此继续下去。因而我感觉事先是无法准确确定网格
y+是与速度相关的量,不仅仅取决于第一层网格的尺寸,因此必须进行实际运算才能知道。不过如果做的活多了,就会知道第一层网格该怎么取了!
主要是计算后检验,需要反复几次。
一般有 y+在 30至 100左右,太密了,壁函数就作用不大了;过大,壁面处精度不够。对于
粗糙壁面,yp一般为 4倍的粗糙度。30—100是壁面函数的敏感区,大于 100时必须要采用壁面函数修正,但是精度有较大影响,
在小于 30时,可以不用壁面函数,目前一般的求解,尽量不要用壁面函数,为了求解精度的需要,取值小于
30,就可以了。
可以使用公式进行估计
近壁 y+估计值(经验公式) 4 2 2
近壁 y+实际值(计算结果) 3~8 1~6 1~5
如果使用标准壁面函数,真如上面所说,过大和过小,都不好,一般为
30左右最好!
如果使用低雷诺数模型,及 enhance wall
function,y+需要是 1的量界!
至于划分网格的过程,需要多次的反复
ANSYS 14
FLUENT用户教程中二维翼型的例子
The values of
Y+ are dependent
on the resolution of the mesh and the Reynolds
number of the flow, and are
defined only in wall-adjacent cells. The value of
Y+in
the wall-adjacent cells dictates how wall shear stress is
calculated.When you use the Spalart-Allmaras model, you should
check that Y+ of the
wall-adjacent cells is either very small
(on the order of
Y+), or
approximately 30 or greater. Otherwise, you
shouldmodify
your mesh.
The equation
for Y+ is
Where Y is the distance from the wall to the
cell center, μ
is the molecular
viscosity,ρis the density of the air,
and τWis the wall shear
stress.
Figure 5.11 (p. 254) indicates that, except for a few small regions
(notably at the shockand the
trailing edge), Y+ >30and for much of these regions it does not drop
significantlybelow
30. Therefore, you can conclude that the near-wall mesh resolution
isacceptable.
Figure
5.11 XY Plot of y+ Distribution
二维管流教程
matlab速度梯度,[转载]关于FLUENT中Y+的一些讨论相关推荐
- 转载 Fluent中的压力分类
大气压 ( atmospheric pressure )指的大气对浸在它里面的物体产生的压强.一标准大气压(1atm)=760毫米汞柱(mmHg)=101325Pa. 表压(gauge pressur ...
- Matlab中N是什么意思,MATLAB中y=FFT(X,N)中的N是什么意思
matlab傅里叶变换中fft(x,n),x,n分别是什么含义? fft(x,n)是一维快速傅里叶变换,x相当于信号,n是变换点数.离散傅里叶变换DFT的快速算法就是FFT. matlab中FFT函数 ...
- 连加在matlab中的表述,matlab中y(1)+.y(i)的连加程式怎么写啊?
matlab中y(1)+.y(i)的连加程式怎么写啊?以下文字资料是由(历史新知网www.lishixinzhi.com)小编为大家搜集整理后发布的内容,让我们赶快一起来看一下吧! matlab中y( ...
- matlab中 y =ft(x)的意思,matlab中y=fft(x)语句的意思
matlab傅里叶变换中fft(x,n),x,n分别是什么含义? fft(x,n)是一维快速傅里叶变换,x相当于信号,n是变换点数.离散傅里叶变换DFT的快速算法就是FFT. matlab中FFT函数 ...
- Fluent中的网格自适应技术
Fluent中的网格自适应技术 Fluent中的网格自适应技术可以允许我们根据数据计算结果来修改网格梳密布置或网格走向.在fluent中非结构化网格特征的两个重要优点是: 1.跟结构化网格相比减少创建 ...
- fortran调用matlab画图,[转载]Matlab与Fortran的交互--基本概念
如果需要进行Matlab环绕Fortran语言的交互,必须通过Matlab提供的应用编程接口(API)来实现与外部的接口,从而可以在Matlab环境中实现调用Fortran程序.输入或输出数据以及和其 ...
- Fluent Ribbon项目出现“命名空间“clr-namespace:Fluent;assembly=Fluent”中不存在“RibbonWindow”名称”的解决方法...
之前在学习@aganqin的Fluent Ribbon项目(http://www.cnblogs.com/aganqin/p/3269384.html).但是一直都有引用了Fluent.dll后仍旧出 ...
- fluent周期边界_在Fluent中采用TUI设置周期性边界的方法
在Fluent中采用TUI设置周期性边界的方法 By Zjz (1) 首先,指定计算域的Rotational Axis(Cell Zone Conditions中,Frame Motion) (2) ...
- matlab读Excel表格数据画图,matlab读Excel表格数据画图-matlab如何从excel表格中读取数据?...
如何将excel表格中大量数据导入matlab中并作图 1.打开matlab,点击主页下面的数据导入,你可以导入excel数据,这里导入你自己命名的huitushuju文件. 2.单击"打开 ...
最新文章
- Python的subprocess子进程和管道进行交互
- close和shutdown的区别
- 在c语言中scanf什么时候用,scanf什么时候用??c语言?
- 从零开始学视觉Transformer(4):ViT模型全流程解析
- Oracle关联查询-数据类型不一致问题 ORA-01722: 无效数字
- android e图标一会有一会没有反应,android解决菜单Menu添加Icon无效问题,让菜单同时显示图标和文字...
- java json 修改字段_我们如何使用Java中的Jackson来更改JSON中的字段名称?
- 网络校时(非NTP)
- .net 中使用rabbitmq
- html显示时间,并且实时刷新(详细代码)
- 2021华为软件精英挑战赛
- oracle 安装包 下载、plsql 64位 安装包下载 、 plsql注册码
- SkeyePlayer RTSP/RTMP低延迟播放器源码解析系列之效率优化方案
- vivado实现cameralink接口
- 用php写出一个网站的模板,25 个 PHP 的 Web 开发程序网站模板
- lineout是什么接口
- t检验自由度的意义_统计学常用概念:T检验、F检验、卡方检验、P值、自由度...
- 华东理工大学计算机考研资料汇总
- UE4编译错误:Unable to delete XXX
- 漏洞复现-Wordpress 4.6 PwnScriptum RCE命令执行
热门文章
- 恩施机器人编程_真慧学机器人编程:挖掘孩子的自身潜力,开发孩子智力
- IT小常识---局域网---互联网---计算机
- Premature end of file 错误解决
- Keynote幻灯片制作三
- [20140724]kitty常用快捷键.txt
- 掺镝的钒磷酸钇(Y(P,V)O4:Dy3+)稀土发光材料|镝掺杂钒磷酸钇(Y(P,V)O4:Dy^3+)纳米晶-齐岳生物wyf
- iceberg的调研报告
- 优傲机器人 官网CB3 在线学习 模块3-设置工具
- elementUI分页组件实现表格数据序号累加功能
- 王者荣耀米莱狄契约魔法皮肤特效一览