fluent经验
关于FLUENT中Y+的一些讨论
关于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，不过这个值不是绝对的。
Fluent学习总结
FLUENT学习总结1 概述：FLUENT是目前处于世界领先地位的商业CFD软件包之一，最初由FLUENT Inc.公司发行。
2006年2月ANSYS Inc.公司收购FLUENT Inc.公司后成为全球最大的CAE软件公司。
FLUENT是一个用于模拟和分析复杂几何区域内的流体流动与传热现象的专用软件。
FLUENT提供了灵活的网格特性，可以支持多种网格。
用户可以自由选择使用结构化或者非结构化网格来划分复杂的几何区域，例如针对二维问题支持三角形网格或四边形网格；针对三维问题支持四面体、六面体、棱锥、楔形、多面体网格；同时也支持混合网格。
用户可以利用FLUENT提供的网格自适应特性在求解过程中根据所获得的计算结果来优化网格。
FLUENT是使用C语言开发的，支持并行计算，支持UNIX和Windows等多种平台，采用用户/服务器的结构，能够在安装不同操作系统的工作站和服务器之间协同完成同一个任务。
FLUENT通过菜单界面与用户进行交互，用户可以通过多窗口的方式随时观察计算的进程和计算结果。
计算结果可以采用云图、等值线图、矢量图、剖面图、XY散点图、动画等多种方式显示、存贮和打印，也可以将计算结果保存为其他CFD软件、FEM软件或后处理软件所支持的格式。
FLUENT还提供了用户编程接口，用户可以在FLUENT的基础上定制、控制相关的输入输出，并进行二次开发。
1.1 FLUENT软件包的组成针对不同的计算对象，CFD软件都包含有3个主要功能部分：前处理、求解器、后处理。
其中前处理是指完成计算对象的建模、网格生成的程序；求解器是指求解控制方程的程序；后处理是指对计算结果进行显示、输出的程序。
FLUENT软件是基于CFD软件的思想设计的。
FLUENT软件包主要由GAMBIT、Tgrid、Filters、FLUENT几部分组成。
（1）前处理器。
包括GAMBIT、Tgrid和Fliters。
其中GAMBIT是由FLUENT Inc.公司自主开发的专用CFD前置处理器，用于模拟对象的几何建模以及网格生成。
fluent经验之谈(过来人的总结)
fluent经验之谈(过来人的总结).docFluent经验之谈（过来人的总结）引言Fluent作为计算流体动力学（CFD）领域内一款强大的软件工具，被广泛应用于工程设计、科研和教育等多个领域。
它能够帮助工程师和研究人员模拟和分析流动、热传递和化学反应等复杂现象。
本文档将基于个人使用Fluent的经验，提供一些实用的技巧和建议，以帮助新用户更高效地学习和使用Fluent。
Fluent软件概述Fluent的主要功能流动模拟：包括层流、湍流等流动特性的模拟。
热传递分析：涉及导热、对流和辐射等热传递方式。
化学反应模拟：模拟燃烧、化学反应等过程。
Fluent的应用领域航空航天：飞机设计、发动机性能分析等。
汽车工业：汽车空气动力学、冷却系统设计等。
能源领域：风力发电、太阳能热利用等。
环境工程：污染物扩散、室内空气质量等。
Fluent学习路径基础知识流体力学基础：理解流体的基本性质和流动规律。
数值方法：了解有限体积法、有限元法等数值求解方法。
Fluent界面熟悉用户界面：熟悉Fluent的图形用户界面（GUI）。
命令行操作：学习使用Fluent的命令行工具。
实践操作案例练习：通过实际案例练习来加深理解。
参数调整：学习如何调整模型参数以获得更准确的结果。
Fluent建模技巧几何建模精确建模：确保几何模型的准确性，避免简化过度。
边界条件：合理设置边界条件，如入口、出口、壁面等。
网格划分网格质量：生成高质量的网格，避免过度拉伸或扭曲。
网格细化：在关键区域进行网格细化，提高模拟精度。
物理模型选择流动模型：根据流动特性选择合适的流动模型，如k-ε、k-ω等。
湍流模型：选择适合流动特性的湍流模型。
Fluent求解设置求解器配置压力-速度耦合：选择合适的耦合求解器，如SIMPLE、PISO等。
迭代方法：设置适当的迭代方法和收敛标准。
监控和收敛残差监控：监控残差曲线，判断模拟是否收敛。
收敛标准：根据问题特性设置合理的收敛标准。
FLUENT学习经验总结(狠珍贵,学长传授)
1对于刚接触到FLUENT新手来说，面对铺天盖地的学习资料和令人难读的FLUENT help，如何学习才能在最短的时间内入门并掌握基本学习方法呢？答：学习任何一个软件，对于每一个人来说，都存在入门的时期。
认真勤学是必须的，什么是最好的学习方法，我也不能妄加定论，在此，我愿意将我三年前入门FLUENT心得介绍一下，希望能给学习FLUENT的新手一点帮助。
由于当时我需要学习FLUENT来做毕业设计，老师给了我一本书，韩占忠的《FLUENT流体工程仿真计算实例与应用》，当然，学这本书之前必须要有两个条件，第一，具有流体力学的基础，第二，有FLUENT 安装软件可以应用。
然后就照着书上二维的计算例子，一个例子，一个步骤地去学习，然后学习三维，再针对具体你所遇到的项目进行针对性的计算。
不能急于求成，从前处理器GAMBIT，到通过FLUENT进行仿真，再到后处理，如TECPLOT，进行循序渐进的学习，坚持，效果是非常显著的。
如果身边有懂得FLUENT的老师，那么遇到问题向老师请教是最有效的方法，碰到不懂的问题也可以上网或者查找相关书籍来得到答案。
另外我还有本《计算流体动力学分析》王福军的，两者结合起来学习效果更好。
2 CFD计算中涉及到的流体及流动的基本概念和术语：理想流体和粘性流体；牛顿流体和非牛顿流体；可压缩流体和不可压缩流体；层流和湍流；定常流动和非定常流动；亚音速与超音速流动；热传导和扩散等。
A.理想流体（Ideal Fluid）和粘性流体（Viscous Fluid）：流体在静止时虽不能承受切应力，但在运动时，对相邻的两层流体间的相对运动，即相对滑动速度却是有抵抗的，这种抵抗力称为粘性应力。
流体所具备的这种抵抗两层流体相对滑动速度，或普遍说来抵抗变形的性质称为粘性。
粘性的大小依赖于流体的性质，并显著地随温度变化。
实验表明，粘性应力的大小与粘性及相对速度成正比。
当流体的粘性较小（实际上最重要的流体如空气、水等的粘性都是很小的），运动的相对速度也不大时，所产生的粘性应力比起其他类型的力如惯性力可忽略不计。
Fluent后处理原创经验
Fluent后处理原创经验一、前言：由于时间关系，我一直没有学习fluent的相关的后处理软件，有几次虽然想学tocplot，但是总觉得好难上手……本科学习的时候，出于爱好，学了一些网页制作方面的软件photoshop、flash等，上研究生之后学了matlab等，在充分利用我原先学过的软件，并把它们很好地结合之后，我顺利地完成了fluent计算完毕之后的一些简单的后处理任务。
如果你以前的经历和我一样，学过一些图像处理方面的软件，用过matlab，会excell，而且对后处理数据要求不高，仅限于二维的数值处理，也许这篇文章对你会有所帮助。
当然我这个方法比较笨。
大连理工大学:Haiong2004/12/21 二、Fluent自带后处理功能与本方法将要介绍的比较Fluent自带后处理本文将要介绍的方法直接用fluent读出的xy图 Matlab绘出的xy图绘出的图形很粗，如果将这样的多条曲线重叠在一起，可能看不清楚清晰，而且可以充分利用matlab的绘图功能，得到你想要的效果Fluent自带的后处理功能做出的动画效果不清晰而且如果长宽比较大，则不能充分利用视窗好处是比较节省时间能做出较为清晰的图片，动画作出的图片动画可以充分利用视窗，将不必要的剪除方法比较笨，耗时较多,需要撑握多个软件，所以针对的使用者是已经学过在使用此方法时需要用到的软件三、用Matlab绘xy图出的文件的名字（此处为：wantmyplot，后缀默认）和选择你要保存的文件夹。
选择“所有文件”，选中wantmyplot文件，双击，出现下图：观察预览文件框里的数据，可以发现前4行为我们不需要的文本文件，因此在“导入起始行（R）”中输入“5”，点击完成。
出行如下左图：3）打开Matlab.在”Command Window”中输入”x=[]”。
4）切换到Excell窗口，按住ctrl+end，当前单元到最后一行，移动上下左右键使当前单元为A1001，如上图右图黑框。
fluent过来人经验谈之continuity不收敛的问题
continuity不收敛的问题（1）连续性方程不收敛是怎么回事？在计算过程中其它指数都收敛了，就continuity不收敛是怎么回事。
这和fluent程序的求解方法SIMPLE有关。
SIMPLE根据连续方程推导出压力修正方法求解压力。
由于连续方程中流场耦合项被过渡简化，使得压力修正方程不能准确反映流场的变化，从而导致该方程收敛缓慢。
你可以试验SIMPLEC方法，应该会收敛快些。
在计算模拟中，continuity总不收敛，除了加密网格，还有别的办法吗？别的条件都已经收敛了，就差它自己了，还有收敛的标准是什么？是不是到了一定的尺度就能收敛了，比如10－e5具体的数量级就收敛了continuity是质量残差，具体是表示本次计算结果与上次计算结果的差别，如果别的条件收敛了，就差它。
可以点report，打开里面FLUX选项，算出进口与出口的质量流量差，看它是否小于0.5%.如果小于，可以判断它收敛.(2) fluent残差曲线图中continuity是什么含义？是质量守恒方程的反映，也就是连续性的残差。
这个收敛的快并不能说明你的计算就一定正确，还要看动量方程的迭代计算。
表示某次迭代与上一次迭代在所有cells积分的差值，continuty表示连续性方程的残差(3) 正在学习Fluent,模拟圆管内的流动，速度入口，出口outflow 运行后xy的速度很快就到1e-06了，但是continuity老是降不下去，维持在1e-00和1e-03之间，减小松弛因子好像也没什么变化大家有什么建议吗？你查看了流量是否平衡吗？在report->flux里面操作，mass flow rate，把所有进出口都选上，compute一下，看看nut flux是什么水平，如果它的值小于总进口流量的1%，并且其他检测量在继续迭代之后不会发生波动，也可以认为你的解是收敛的。
造成连续方程高残差不收敛的原因主要有以下几点：1.网格质量，主要可能是相邻单元的尺寸大小相差较大，它们的尺寸之比最好控制在1.2以内，不能超过1.4.2.离散格式及压力速度耦合方法，如果是结构网格，建议使用高阶格式，如2阶迎风格式等，如果是非结构网格，除pressure保持standard格式不变外，其他格式改用高阶格式；压力速度耦合关系，如果使用SIMPLE，SIMPLEC，PISO等segerated solver对联系方程收敛没有提高的话，可以尝试使用coupled solver。
fluent-udf-经验
[转]菜鸟学UDF的感觉，希望对UDFers有用光看书，感觉UDF不难。
看例子，有些看个四五遍之后才能差不多看懂。
原来，得靠UDF 帮助。
我主要用的是fluent v6.3自带的html格式的帮助，里面东西很全，当然也包括UDF Manual。
里面自带的search功能相当好，只是要注意用好+或-号（逻辑符号），另外，这个功能似乎有些浏览器支持不太好，不过基本上用IE不太容易出问题。
对于从零开始学习UDF，建议还是先看一下UDF中文帮助，我估计大家知道的都是马世虎翻译的那本吧，感觉挺好。
（没想到马世虎跟我是校友，去年给安世亚太投过一份简历，他给我打过电话，当时一阵兴奋，呵呵。
）对于只涉及到边界条件或物性等的UDF，一般用interpret就可以的，这些我觉得只需要根据例子改一下就是了。
$$ 对于要添加UDS方程的，相对难一点。
我编程用的是三到五个UDS，几十个UDM。
一开始编程时，没有头绪，后来看别人编的，才慢慢发现了一些基本思路。
比如，可以用枚举定义UDS或UDM，这样用起来方便。
enum{NP,RHOH2O_Y_UP_X,RHOH2O_Y_UP_Y,RHOH2O_Y_UP_Z,N_REQUIRED_UDS};//枚举UDS变量名对于UDM，则用N_REQUIRED_UDM代表个数。
然后在INIT与ADJUST函数中，检查变量个数时则比较方便，如：DEFINE_INIT(init_parameter,domain){if (n_uds < N_REQUIRED_UDS)Error(”Not enough user defined scalars!(init)\n”);if (n_udm<N_REQUIRED_UDM)Error(”Not enough user defined memories(init)!\n”);initialise(domain);//代表初始化}DEFINE_ADJUST(adjust_compute,domain){if (n_uds < N_REQUIRED_UDS)Error(”Not enough user defined scalars!(adjust)\n”);if (n_udm<N_REQUIRED_UDM)Error(”Not e nough user defined memories(adjust)!\n”);update_parameter(domain);//代表主函数}初始化时，则可：cell_t c;Thread *t;int i;thread_loop_c(t,d){if(NNULLP(THREAD_STORAGE(t,SV_UDS_I(NP)))&&NNULLP(THREAD_STORAGE( t,SV_UDS_I(NP_R)))) //为各UDS提供存储空间{begin_c_loop(c, t){for (i=0; i<N_REQUIRED_UDS; i++)C_UDSI(c,t,i) = 0.0;}end_c_loop(c, t);}if(NNULLP(THREAD_STORAGE(t,SV_UDM_I))){begin_c_loop(c, t){for (i=0; i<N_REQUIRED_UDM; i++)C_UDMI(c,t,i) = 0.0;}end_c_loop(c, t);}}对于各UDM量，则可：real udm_v;udm_v=0;//用之前对变量进行初始化…//UDM相关运行C_UDMI(c,t,UDM_V)=udm_v;//把值输入给UDM，当然之前要对UDM_V进行定义用UDM有个好处，一是可以在后处理中显示，二是传递变量相当方便，比如在ADJUST 中计算的量用于源项或对流项等，用UDM可以直接调用。
fluent一些笔记
1．提示：Welcome to Fluent 6.3.26Copyright 2006 Fluent IncAll Rights ReservedCannot open dump file "fl_s1119.dmp"Error: Unable to open dump file9Error encountered in critical code sectionHit return to exit.原因：使用优化大师进行系统清理，而优化大师默认是要把lib文件夹下的这三个文件给删除的，这就导致FLUENT启动的时候找不到"fl_s1119.dmp。
处理方法：把lib文件夹做下备份，优化完以后覆盖即可。
（毕竟系统还是要清理的）偶是发现了，备份的时候最好是压缩包，要不然优化大师会把你的备份文件也一起清理掉的。
2．1、输出grid图形2、选择surface---plane,打开plane surface面板3、通过确定三个点来确定平面位置。
单击slect point，出现提示，不点选cancel.在grid 图形的多孔介质区域任意位置右键点选3个点。
4、回到plane surface面板，勾选plane tool，则在grid图形的多孔介质区域出现一个平面。
若出现的平面与我们的预期相差比较大的话，可以单击reset points,可以获得一个特殊位置的平面。
5、打开多孔介质的控制面板，选择porou zone标签，点击update from plane tool按钮，获得方向矢量1，和方向矢量2的原始值，并与左下角的坐标系统比较，确定我们大概的旋转方向。
6、对比grid图形左下角的坐标系统，红线和红色箭头代表的是方向矢量1，绿线和绿色箭头代表的是方向矢量2应该使红线和X正方向平行，绿线和Y正方向平行。
具体的操作应该是：（首先把平面移动到图形外有利于旋转，比较清楚。
平面法线方向的移动是用鼠标右键单击平面阴影部分并拖动，横向移动则需按下shift并进行如上操作。
FLUENT使用基础经验
FLUENT使⽤基础经验Fluent使⽤经验记录基于Fluent全攻略基础教程成功模拟的⼏点建议1、Check和Scale⽹格，⽹格斜度Skewness < 0.98才能求解；2、能量亚松弛因⼦范围选在0.95-1，耦合传热问题时，⾼导热系数下，较⼩的亚松弛因⼦会阻碍收敛；3、对于三⾓形和四⾯体⽹格，采⽤Node-Based平均梯度，⽐默认的Cell-Based平均梯度结果更为精确；4、残差⾄少下降3个量级，对于pressure-based求解器能量⽅程残差要达到10-6，物态⽅程达到10-5；5、收敛的解不⼀定是正确的解，需要⽤⾼阶格式求解；6、⽹格不相关处理，采⽤Adatption；基本概念N-S⽅程：由Navier和Stokes 通过在Euler’s Equation 上加⼊粘性项得到的粘性不可压流体动量⽅程，对于可压流体还需能量⽅程和状态⽅程。
N-S⽅程简化：N-S > RANS > 边界层⽅程> ⽆粘⾮线性⽅程【如⾼Re，粘性⽐重下降，转变为Euler’s Equation，势能⽅程，跨⾳⼩扰动⽅程】> ⽆粘线性⽅程【如Laplace⽅程】。
数值求解本质：将控制⽅程线性化并⽤离散⽅法和格式变为代数⽅程，⽤线性⽅程迭代求解。
何时使⽤湍流模型：流动为湍流时，Re>2320认为是湍流，因为实际中初始流场常为湍流，⽽湍流下临界Re=2320。
问题：⽹格smooth/swap操作和含义？如果不加湍流模型那么求解可以进⾏吗，求解的是什么⽅程？Hydraulic Diameter：⽔⼒直径(Hydraulic diameter)是，⾮圆形截⾯管道等效成圆形截⾯管道的⼀个⼏何尺⼨，⽤于计算雷诺数，判断管道内流体是层流还是湍流状态。
Turbulence Intensity：湍流强度，湍流脉动速度与平均速度之⽐。
湍流强度可以⽤I=0.16Re^(-1/8)来近似计算，如何推导的？Turbulence Kinetic Energy：湍流动能，湍流速度脉动⽅差与质量乘积的1/2。
fluent过来人经验谈之continuity不收敛的问题精品
continuity不收敛的问题（1）连续性方程不收敛是怎么回事？在计算过程中其它指数都收敛了，就continuity不收敛是怎么回事。
这和fluent程序的求解方法SIMPLE有关。
SIMPLE根据连续方程推导出压力修正方法求解压力。
由于连续方程中流场耦合项被过渡简化，使得压力修正方程不能准确反映流场的变化，从而导致该方程收敛缓慢。
你可以试验SIMPLEC方法，应该会收敛快些。
在计算模拟中，continuity总不收敛，除了加密网格，还有别的办法吗？别的条件都已经收敛了，就差它自己了，还有收敛的标准是什么？是不是到了一定的尺度就能收敛了，比如10－e5具体的数量级就收敛了continuity是质量残差，具体是表示本次计算结果与上次计算结果的差别，如果别的条件收敛了，就差它。
可以点report，打开里面FLUX选项，算出进口与出口的质量流量差，看它是否小于0.5%.如果小于，可以判断它收敛.(2) fluent残差曲线图中continuity是什么含义？是质量守恒方程的反映，也就是连续性的残差。
这个收敛的快并不能说明你的计算就一定正确，还要看动量方程的迭代计算。
表示某次迭代与上一次迭代在所有cells积分的差值，continuty表示连续性方程的残差(3) 正在学习Fluent,模拟圆管内的流动，速度入口，出口outflow 运行后xy的速度很快就到1e-06了，但是continuity老是降不下去，维持在1e-00和1e-03之间，减小松弛因子好像也没什么变化大家有什么建议吗？你查看了流量是否平衡吗？在report->flux里面操作，mass flow rate，把所有进出口都选上，compute一下，看看nut flux是什么水平，如果它的值小于总进口流量的1%，并且其他检测量在继续迭代之后不会发生波动，也可以认为你的解是收敛的。
造成连续方程高残差不收敛的原因主要有以下几点：1.网格质量，主要可能是相邻单元的尺寸大小相差较大，它们的尺寸之比最好控制在1.2以内，不能超过1.4.2.离散格式及压力速度耦合方法，如果是结构网格，建议使用高阶格式，如2阶迎风格式等，如果是非结构网格，除pressure保持standard格式不变外，其他格式改用高阶格式；压力速度耦合关系，如果使用SIMPLE，SIMPLEC，PISO等segerated solver对联系方程收敛没有提高的话，可以尝试使用coupled solver。
Fluent网格划分经验
Fluent⽹格划分经验转载⾸先，在⽹格划分之前，你最好从数值仿真的全局出发，⽐如精度要求，计算时间要求，机⼦配置等等，思考⼀下是使⽤结构⽹格，还是⾮结构⽹格，抑或是混合⽹格；因为这关系到接下来的⽹格划分布置和划分策略。
然后，在确定了⽹格类型之后，就是根据模型情况，构思⼀下⽹格拓扑，就是⾃⼰要明确最终想得到什么样的⽹格，⽐如翼型⽹格，是C型，还是O型；⼀个圆⾯是想得到“内⽅外圆”的铜钱币类型的⽹格，还是⼀般的⽹格，等等。
这⼀步有时可能不太清楚，⾃⼰有时都不知道什么样的⽹格拓扑是合适的，那就需要平时多看看这⽅⾯的帖⼦，收集⼀些划分⽐较好的⽹格图⽚，体会体会。
确定了⽹格拓扑之后，对模型进⾏划分⽹格前的准备，⽐如分割啊，对尺度⼩对计算结果影响不⼤的次要⼏何进⾏简化，等等。
接着，划分⽹格。
划分⽹格都是从线⽹格，⾯⽹格，到体⽹格的；线⽹格的划分，也就是⽹格节点的布置，对⽹格的质量影响⽐较⼤，⽐如歪斜，长宽⽐，等等，节点密度在GAMBIT中可以通过很多的⽅法进⾏控制调整，⼤家可以看相关的资料。
⾯⽹格的划分，⾮结构的⽹格咱就不说了，结构⽹格可能有时⽐较⿇烦，这就要求⼤家最好对那⼏种⽹格策略⽐较了解，⽐如Quad-Map划分⽅法所适⽤的模型形状，在划分的时候对顶点类型及⽹格节点数的要求（Quad-Map，适⽤于边数⼤于或等于4的⾯，顶点要求为4个End类型，其他为Side类型，对应边的⽹格节点数必须相等），以此类推，其他的划分⽅法也有这⽅⾯的要求以及适合的形状。
当出现了不能划分的时候，可以根据GAMBIT给的提⽰进⾏修改顶点类型或⽹格节点数来满⾜划分⽅法的要求。
如果实在不能划分，则退⽽求其次，改⽤其他⽅法进⾏划分或者对⾯进⾏分割；等等。
关于体⽹格的划分，与⾯⽹格划分所要注意的东西类似。
另外，根据我个⼈的经验，如果模型⽐较简单规则，⼤家最好尽量使⽤结构⽹格，⽐较容易划分，计算结果也⽐较好，计算时间也相对较短；对于复杂的⼏何，在尽量少的损失精度的前提下，尽量使⽤分块混合⽹格。
fluent过来人经验谈之continuity不收敛的问题
continuity不收敛的问题（1）连续性方程不收敛是怎么回事？在计算过程中其它指数都收敛了，就continuity不收敛是怎么回事。
这和fluent程序的求解方法SIMPLE有关。
SIMPLE根据连续方程推导出压力修正方法求解压力。
由于连续方程中流场耦合项被过渡简化，使得压力修正方程不能准确反映流场的变化，从而导致该方程收敛缓慢。
你可以试验SIMPLEC方法，应该会收敛快些。
在计算模拟中，continuity总不收敛，除了加密网格，还有别的办法吗？别的条件都已经收敛了，就差它自己了，还有收敛的标准是什么？是不是到了一定的尺度就能收敛了，比如10－e5具体的数量级就收敛了continuity是质量残差，具体是表示本次计算结果与上次计算结果的差别，如果别的条件收敛了，就差它。
可以点report，打开里面FLUX选项，算出进口与出口的质量流量差，看它是否小于0.5%.如果小于，可以判断它收敛.(2) fluent残差曲线图中continuity是什么含义？是质量守恒方程的反映，也就是连续性的残差。
这个收敛的快并不能说明你的计算就一定正确，还要看动量方程的迭代计算。
表示某次迭代与上一次迭代在所有cells积分的差值，continuty表示连续性方程的残差(3) 正在学习Fluent,模拟圆管内的流动，速度入口，出口outflow 运行后xy的速度很快就到1e-06了，但是continuity老是降不下去，维持在1e-00和1e-03之间，减小松弛因子好像也没什么变化大家有什么建议吗？你查看了流量是否平衡吗？在report->flux里面操作，mass flow rate，把所有进出口都选上，compute一下，看看nut flux是什么水平，如果它的值小于总进口流量的1%，并且其他检测量在继续迭代之后不会发生波动，也可以认为你的解是收敛的。
造成连续方程高残差不收敛的原因主要有以下几点：1.网格质量，主要可能是相邻单元的尺寸大小相差较大，它们的尺寸之比最好控制在1.2以内，不能超过1.4.2.离散格式及压力速度耦合方法，如果是结构网格，建议使用高阶格式，如2阶迎风格式等，如果是非结构网格，除pressure保持standard格式不变外，其他格式改用高阶格式；压力速度耦合关系，如果使用SIMPLE，SIMPLEC，PISO等segerated solver对联系方程收敛没有提高的话，可以尝试使用coupled solver。
fluent经验
Cooper/hex（hex/wedge）：对指定的源面上的节点模式进行扫掠，从而形成体网格；即把一个体看成是一个
或多个逻辑cylinders，每个cylinders 都包括一个桶状和两个盖（源面）
TGrid/（Tet/Hybrid）：主要采用四面体单元，但是在恰当的地方也用六面体或者锥体，楔体等单元
问题二：在画面网格时已经定义了spacing,那么在画体网格时再定义，会不会有冲突？还有，究竟怎样定义
才是合理的？
第一个问题:在做网格时,所有的错误提示都表明部分网格生成要么有问题没有生成,要么生成的质量很差,要
重新检查一下几何体和网格类型.
Fluent 常见问题回答 摘自清洁能源论坛，Jim King 整理
可采用分块生成网格的办法，或使用边界层网格。
2.二维轴对称问题，想在对称轴附近的狭长区域加密网格。
最简单的就是将模型分块，轴附近的狭长区域是一块，其它是另一块，两块分别分网格。
3. 经常出现在两个面交线上的网格间距不同的现象，也就是两块网格不连续的现象，怎么克服这种情况？
将相邻的线（面）merge 或connect 以后，在公共边及其相对的边上设置相同的节点数就好了。 先将交线
1000000.
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
Fluent 常见问题回答 摘自清洁能源论坛，Jim King 整理
1.这是表示你点的地方是二条等势线的中间吧，不用flood&line 表示，只用line 表现时，就很清楚 了。
FLUENT感悟
1、四面体网格的自动聚合：自动的将四面体网格转化为多面体网格，使网格总数降低，加快收敛。
方法如下：第一个是可以将所有网格都转换；第二个是将指定的转换。
2、Gradient Option 提供了基于控制体中心和控制体节点外，新增“Least Square Cell Based”这是支持新增的支持多面体网格而设计的，它实际上假设变量在网格单元之间线性变化。
3、欧拉多项流模型只能应用于不可压缩流体。
4、由于Fluent软件包最初是在Unix系统下运行的，面向高端用户，随着Windows的产生，逐渐开发应用于Windows操作系统，但是安装前期处理软件Gambit必须在Unix操作系统下才能运行，所以安装之前必须首先在Windows系统下安装一个模拟Unix系统的软件，这就是Exceed。
5、当利用Fluent解决某一问题时，需要详细考虑一下几个问题：A确定计算目标。
预期从CFD模拟计算中获得什么样的结果？获得这些结果需要时间是多久？从模型中得到什么样的精度？B选择计算模型。
所遇到的问题是否可以简化，计算区域如何界定，需要确定什么样的边界条件，各个边界的信息是否充分，什么样的网格拓扑更加有利于解决问题。
C确定物理模型。
确定无粘度还是粘度，层流还是湍流，定常还是非定常，可压缩还是不可压缩，是否需要应用其他模型。
D确定解的程序。
确定是否使用默认解还是参数值，那一次求解格式能够加快收敛，估计得到收敛解需要多长时间，是否需要高阶格式E分析改进计算。
根据收敛结果分析流型是否正确，关键的物理现象是否合理，力、力矩、流量、温度等与实验值比较是否满足实验精度，是否需要做模型方面的改进，计算域是否需要扩大，边界层网格是否需要加密。
6、Re=ρνr/u根据所给参数，计算出雷诺数来确定层流还是湍流。
7、在命令中输入Reset表示可清空此次任务以前的所有操作，从新开始建立模型。
8、Label选项域表示我们为创建的点或者直线起的名字，有利于命名点不至于乱，在视图控制面板中，选择Label小方框选中，就可以将点的名称显示出来。
fluent过来人经验谈之continuity不收敛的问题-推荐下载
continuity不收敛的问题（1）连续性方程不收敛是怎么回事？ 在计算过程中其它指数都收敛了，就continuity不收敛是怎么回事。
这和fluent程序的求解方法SIMPLE有关。
SIMPLE根据连续方程推导出压力修正方法求解压力。
由于连续方程中流场耦合项被过渡简化，使得压力修正方程不能准确反映流场的变化，从而导致该方程收敛缓慢。
你可以试验SIMPLEC方法，应该会收敛快些。
在计算模拟中，continuity总不收敛，除了加密网格，还有别的办法吗？别的条件都已经收敛了，就差它自己了，还有收敛的标准是什么？是不是到了一定的尺度就能收敛了，比如10－e5具体的数量级就收敛了continuity是质量残差，具体是表示本次计算结果与上次计算结果的差别，如果别的条件收敛了，就差它。
可以点report，打开里面FLUX选项，算出进口与出口的质量流量差，看它是否小于0.5%.如果小于，可以判断它收敛.(2) fluent残差曲线图中continuity是什么含义？是质量守恒方程的反映，也就是连续性的残差。
这个收敛的快并不能说明你的计算就一定正确，还要看动量方程的迭代计算。
表示某次迭代与上一次迭代在所有cells积分的差值，continuty表示连续性方程的残差(3) 正在学习Fluent,模拟圆管内的流动，速度入口，出口outflow运行后xy的速度很快就到1e-06了，但是continuity老是降不下去，维持在1e-00和1e-03之间，减小松弛因子好像也没什么变化大家有什么建议吗？你查看了流量是否平衡吗？在report->flux里面操作，mass flow rate，把所有进出口都选上，compute一下，看看nut flux是什么水平，如果它的值小于总进口流量的1%，并且其他检测量在继续迭代之后不会发生波动，也可以认为你的解是收敛的。
造成连续方程高残差不收敛的原因主要有以下几点：1.网格质量，主要可能是相邻单元的尺寸大小相差较大，它们的尺寸之比最好控制在1.2以内，不能超过1.4.2.离散格式及压力速度耦合方法，如果是结构网格，建议使用高阶格式，如2阶迎风格式等，如果是非结构网格，除pressure保持standard格式不变外，其他格式改用高阶格式；压力速度耦合关系，如果使用SIMPLE，SIMPLEC，PISO 等segerated solver对联系方程收敛没有提高的话，可以尝试使用coupled solver。
FLUENT学习经验
FLUENT学习经验1. 什么是结构化网格和非结构化网格1.1结构化网格从严格意义上讲，结构化网格是指网格区域内所有的内部点都具有相同的毗邻单元。
它可以很容易地实现区域的边界拟合，适于流体和表面应力集中等方面的计算。
它的主要优点是：网格生成的速度快。
网格生成的质量好。
数据结构简单。
对曲面或空间的拟合大多数采用参数化或样条插值的方法得到，区域光滑，与实际的模型更容易接近。
它的最典型的缺点是适用的范围比较窄，只适用于形状规则的图形。
尤其随着近几年的计算机和数值方法的快速发展，人们对求解区域的几何形状的复杂性的要求越来越高，在这种情况下，结构化网格生成技术就显得力不从心了。
1.2非结构化网格同结构化网格的定义相对应，非结构化网格是指网格区域内的内部点不具有相同的毗邻单元。
即与网格剖分区域内的不同内点相连的网格数目不同。
从定义上可以看出，结构化网格和非结构化网格有相互重叠的部分，即非结构化网格中可能会包含结构化网格的部分。
2.如果一个几何造型中既有结构化网格，也有非结构化网格，分块完成的，分别生成网格后，也可以直接就调入fluent中计算。
3.在fluent中，对同一个几何造型，如果既可以生成结构化网格，也可生成非结构化网格，当然前者要比后者的生成复杂的多，那么应该选择哪种网格，两者计算结果是否相同，哪个的计算结果更好些呢？一般来说，结构网格的计算结果比非结构网格更容易收敛，也更准确。
但后者容易做。
影响精度主要是网格质量，和你是用那种网格形式关系并不是很大，如果结构话网格的质量很差，结果同样不可靠，相对而言，结构化网格更有利于计算机存储数据和加快计算速度。
结构化网格据说计算速度快一些，但是网格划分需要技巧和耐心。
非结构化网格容易生成，但相对来说速度要差一些。
4.在gambit中，只有map和submap生成的是结构化网格，其余均为非结构化网格。
5 我们经常遇到计算区是对称的问题，如同心圆环内的自然对流，圆柱绕流，我们为了节省计算资源，许多时候都把计算区域趣味一半，但有些问题的真实情况是两步分的流场及物理量的分布并不对称呀，问我们如何判断该不该区一般呢? 对称的问题一般用在流场稳态解..需满足 1.几何图形对秤..2.边界条件对秤.. 也就是物理条件对秤..3.structral网格..所以对秤轴的Flux和properties gradient必须为0...6 按照算例学习了一段时间，有些简单的问题还可以分析对错，但是对于一些头脑里没有概念的问题，是做出了很多图行了，矢量图了，但是如何比较仔细全面的分析其合理性，觉得有些困难答：一般来讲计算应该辅助以高精度的实验作为证明，无法或不容易用实验实现的往往是计算一个经典的或别人算过的例子对比一下。
fluent经验
4.不能用于计算流量分配问题（比如有多个出口的问题）
6.在压力出口中，会要求输入相应的backflow turbulent intensity等值，这些值只有在迭代时产生返流的时候才会使用，
通常设置成一个合理的值。算例14中，设置为intensity 10%，diameter hydraulic按实际模型数值。
器，设定你想要的y+值，它就能给你计算出第一层网格高度，与计算结果的y+很接近。
/APPS/YPlus/
10.消除左手面：/grid/modify-zones/re-fa
11.LES网格尺寸的确定方法：/thread-956892-1-1.html
（2）输入mesh/modify-zones/make-periodic,再根据提示选择相应的面。
5.outflow边界条件不需要给定任何入口的物理条件，但是应用也会有限制，大致为以下四点：
1.只能用于不可压缩流动
2.出口处流动充分发展
7.后处理的时候，显示速度矢量图的时候，箭头的长度可以不按速度的大小给出，而仅由箭头的颜色决定，具体的操作:
Vector options.勾选Fixed Length
8.波尔兹曼数能表达式在第14个例子的最后。
9.不要使用那些书上写的y+与yp的计算公式，那个公式一般只能提供数量级上的参考。推荐大家使用NASA的粘性网格间距计算
（2）求解值不再随迭代发生改变：有时候，残差还在下降，但是某些监视的流动变量不再发生变化即可。
（3）系统的质量、动量、能量达到平衡：利用flux report实现，要求净不平衡量小于0.2%。
学习FLUENT心得笔记
進口排气孔边界条件：用于计算进口排气孔处的损失系数，流动方向，以及周围的温度和压力。
输入：除了一些常见的参数外，主要是一个损失系数（前面的11个和压力边界条件相同）。
对于损失系数，按照公式计算：其中ρ为密度，kL是一个无量纲的经验系数，注意：△p表示流动方向的压力损失，你可以定义为常数或者速的多项式、分段式函数。
定义面板和定义温度相关属性的相同。
进气风扇边界条件用于模型化一个外部的有指定压力升高、流动方向、周围温度和压力的进气风扇。
输入：前11项和压力边界条件的一样。
通过进气风扇的压力上升被认为是流速的函数。
对于逆向流，进气风扇被当作一个带损失系数的出口排气孔。
你可以设置压力上升为常量，或者速度的函数。
压力出口边界条件：需要指定一个静压，这只适用于亚音速流动，对于超音速，这个条件是无用的。
流动的一些特性将由上游推倒得到。
如果在解算过程中流动逆相，需要设置一系列的“逆流”条件。
输入：静压，逆流条件：总温，湍流参数，化学成分质量分数，混合物分数和变迁，过程变量，多相边界条件。
辐射条件，分散相边界条件。
定义静压：注意输入的静压和工作条件面板的工作压力相关，注意关于液体静压的评论。
系统也提供一个关于径向平衡边界条件的选择，选择该项的化，输入的静压只适用于最小半径，其他部分的压力通过下是计算，r为距离回转轴的半径距离，vθ为切向速度。
注意折椅边界条件对于旋转速度是零也适用。
该条件只适用于3D计算和轴对称计算。
定义逆流条件：适用于流体被拖动穿过出口。
定义辐射参数：定义分散相条件：6．9压力far-field边界条件：用于定于无穷远处自由流的压力条件，常被称作特性边界条件，因为这里使用因为这里使用特性信息（黎曼常量）来计算边界的流动变量。
该条件仅适用于利用理想气体公式计算密度的流动，其他的不允许。
该区域必须足够的远。
6．10出口流边界条件：用于模拟结算前流动的速度和方向等都未知的流动，不需要任何设定，但是要能够了解该条件的限制条件。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Fluent 问题集锦问题1:如果体网格做好后,感觉质量不好,然后将体网格删除,在其面上重新作网格,结果发现网格都脱离面,不再附体了,比其先前的网格质量更差了.原因: 删除体网格时，也许连同较低层次的网格都删除了.上面的脱离面可能是需要的体的面.解决方法: 重新生成了面,在重新划分网格问题2:在gambit下做一虚的曲面的网格,结果面上的网格线脱离曲面,由此产生的体网格出现负体积.原因: 估计是曲面扭曲太严重造成的解决方法: 可以试试分区域划分体网格,先将曲面分成几个小面，生成各自的面网格，再划体网格。
问题3:当好网格文件的时候，并检查了网格质量满足要求，但输出*.msh时报错误.原因: 应该不是网格数量和尺寸.可能是在定义边界条件或continuum type时出了问题.解决方法: 先把边界条件删除重新导出看行不行.其二如果有两个几何信息重合在一起, 也可能出现上诉情况,将几何信息合并掉.问题4:当把两个面(其中一个实际是由若干小面组成,将若干小面定义为了group了)拼接在一起,也就是说两者之间有流体通过,两个面各属不同的体,网格导入到fluent时,使用interface时出现网格check的错误,将interface的边界条件删除,就不会发生网格检查的错误.如何将两个面的网格相连.原因: interface后的两个体的交接面,fluent以将其作为内部流体处理（非重叠部分默认为wall，合并后网格会在某些地方发生畸变，导致合并失败.也可能准备合并的两个面几何位臵有误差,应该准确的在同一几何位臵(合并的面大小相等时),在合并之前要合理分块。
解决方法: 为了避免网格发生畸变(可能一个面上的网格跑到另外的面上了),可以一面网格粗,一面网格细,避免; 再者就是通过将一个面的网格直接映射到另一面上的,两个面默认为interior.也可以将网格拼接一起.上述语言有些模糊不清,仅供参考,并希望高手批评指正,^_^/cgi-bin/bbs/topic_show.cgi?id=17102&h=1#136480和/cgi-bin/bbs/topic_show.cgi?id=29648&h=1#2330872、关于入口速度设臵的讨论，见/cgi-bin/bbs/topic_show.cgi?id=22671&h=1#18454233、关于gambit建模出现的问题的讨论，见/cgi-bin/bbs/topic_show.cgi?id=32453&h=1#2476224、关于出口条件的设臵问题，见/cgi-bin/bbs/topic_show.cgi?id=26815&h=1#206126问题：FLUENT的初始化面板中有一项是设臵从哪个地方开始计算(compute from),选择从不同的边界开始计算有很大的区别吗?该怎样根据具体问题选择从哪里计算呢?比如有两个速度入口A和B,还有压力出口等等，是选速度入口还是压力出口?如果选速度入口,有两个,该选哪个呀?有没有什么原则标准之类的东西？答：一般是选取ALL ZONE，即所有区域的平均处理，通常也可选择有代表性的进口(如多个进口时)进行初始化。
对于一般流动问题，初始值的设定并不重要，因为计算容易收敛。
但当几何条件复杂，而且流动速度高变化快(如音速流动)，初始条件要仔细选择。
如果不收敛，还应试验不同的初始条件，甚至逐次改变边界条件最后达到所要求的条件。
问题：要判断自己模拟的结果是否是正确的，似乎解的收敛性要比那些初始条件和边界条件更重要，可以这样理解吗？也就是说，对于一个具体的问题，初始条件和边界条件的设定并不是唯一的，为了使解收敛，需要不断调整初始条件和边界条件直到解收敛为止，是吗？如果解收敛了，是不是就可以基本确定模拟的结果是正确的呢？答：对于一个具体的问题，边界条件的设定当然是唯一的，只不过初始化时可以选择不同的初始条件（指定常流），为了使解的收敛比较好，我一般是逐渐的调节边界条件到额定值（"额定值"是指你题目中要求的入口或出口条件，例如计算一个管内流动，要求入口压力和温度为10MPa和3000K，那么我开始叠代时选择入口压力和温度为1MPa和500K（假设，这看你自己问题了），等流场计算的初具规模、收敛的较好了，再逐渐调高压力和温度，经过好几次调节后最终到达额定值10MPa和3000K，这样比一开始就设为10MPa和3000K收敛的要好些）这样每次叠代可以比较容易收敛，每次调节后不用再初始化即自动调用上次的解为这次的初始解，然后继续叠代。
即使解收敛了，这并不意味着就可以基本确定模拟的结果是正确的，还需要和实验的结果以及理论分析结果进行对比分析。
在fluent中，用courant number来调节计算的稳定性与收敛性。
一般来说，随着courantnumber的从小到大的变化，收敛速度逐渐加快，但是稳定性逐渐降低。
所以具体的问题，在计算的过程中，最好是把courant number从小开始设臵，看看迭代残差的收敛情况，如果收敛速度较慢而且比较稳定的话，可以适当的增加courant number的大小，根据自己具体的问题，找出一个比较合适的courant number，让收敛速度能够足够的快，而且能够保持它的稳定性对于流体力学试验、数值模拟、理论研究的关系，本人以为应当从哲学的方面借以揭示。
因为我们是在蒙着眼睛过河，我们面对的是不清楚的事实，通过观察，测量，得到数据，然后用自己接受的一套规律把它理解，最后依靠这套规律来预测、解释和研究这个事实。
我们都是被动的在接受存在的事物。
这就像哲学上的认识论，人类的知识总是在不断增加，但最终能够完全认识整个宇宙吗？恐怕没有答案。
人们对流体的认识不断加深，但最终能够完全用数学描述流体吗？不可能。
那我们做理论来做什么？为了尽可能的了解，而作数模则是在尽可能了解的基础上做近似，用求解通过观察获得的理论的方法来反演事实。
这肯定是不准确的，但如果方法正确，应该是近似准确的。
就像很简单的两个球，假定他们是绝对刚体，则我们的动量定理可以发挥作用，但刚体并不存在。
所以我们的认识是错误的，但却在一定条件下接近正确。
所以说，做数值模拟的最终结果需要试验来验证，但试验的正确性需要谁来验证？反复的试验？问题：比较两种耦合算法的性能，怎样把二者的RESIDUAL HISTORY 放到一张图上呢?RESIDUAL 中没有WRITE FILE 这一选项呀?答：计算完毕以后，首先把二者的RESIDUAL HISTORY图拷贝下来，然后可以放到一张图上编辑用Tecplot显示Fluent计算的轴对称流场，请问如何把对称轴的下面一般流场显示出来？Fluent保存的流场数据和网格都是只有上半流场的。
答：复制一下上半场（Create Zone），把Y的符号反号一下就可以了(Equation)。
另外最简单的办法，自己写一段小程序，将数据复制一份，之后反号。
问：一对相互平行的斜面，几何结构和尺寸完全相同，用相同的参数控制网格，结果却不相同，相互之间相差将近10个网格。
这样的话，要实现cooper的拓扑功能，只能由一个面映射到对面，同时也会导致体网格严重畸变。
所以想先划分两个平行面的面网格，然后利用cooper实现体网格的生成。
可是两个面的网格结构即使存在细微的差别，也会导致体网格生成失败，不知道大家有没有遇到这种情况答：你说的情形应该是不可能出现的，如果使用了虚体，你仔细检查一下你对应的线是不是正确，我曾经碰到过这种情形，如果两边都是实体，应该是完全一样的。
问：我用的就是实体，一个四分之一圆柱，用一组相互平行的斜面进行Split操作，得到几个相互连接的体，而在所得到的这些切割面上生成面网格，死活都不能一致。
答：你只要严格控制对应边的节点个数，对应面网格自然会相同，也就可以Fluent后处理的匹处理！如果有很多个只有边界条件不同的cas/dat文件要进行处理，可以这样做：file----write----start journal---（命名一个记录文件）; 然后fluent 把你下面的每一步操作记录下来，在你想记录结束的时候，file----write-----stop journal ，对于下一个文件，file----read----journal-----(那个记录文件），你原来的操作在新的cas、dat文件上从来一遍。
计算中的疑问：turbulent viscosity limited to viscosity ratio of 1.000000e+5 in cell **,是说明什么？计算有问题了吗？是指在计算过程中turbulent viscosity超过了设定的阈值.如果实际中的确是超过的话,可以在solve->controls->limits中调整.更详细的说明请看帮助!问：什么是旋流数？答：切向动量和轴向动量比. 强旋流动用混合长，K-e模型都不行，反应不出强旋情况下的各向异型。
用代数应力模型或者雷诺应力模型比较合适。
关于模拟湍流流场适用性的问题具体可以参考周力行的书。
一般来说，混合长，K-e模型（用的最广）算自由射流，剪切流，弱旋，无浮力流等等各向异型不强的情况符合很好。
对强旋流动，也可以采用修正后的K-e模型（有很多修正，对强旋就找针对强旋修正的模型）用于生成流场的边界，已有了这根曲线上足够多点的坐标了。
然后想用CAD做，完了再导入到Gambit里面，但是结果导入后发现曲线被分成了很多段，怎么让它连接成一根光滑的线呢？另外好像CAD里面只能导出region，单单一根线怎么导出的？直接用jou文件导入gambit,格式用vertex create coordinate x y z 用fortran编,把你每个点的坐标放到前面" vertex create coordinate"命令后面,然后文件名字用.jou,最后用gambit运行jou,然后扫描点生成曲线,肯定光滑.关于湍流问题的数值模拟目前工程中常用基于Reynolds平均的方法，用各种封闭方法来获得时均的控制方程，以求得工程上比较受关心的时均流场及湍动信息,但实际上在诸如化工、冶金、能源等诸多过程工业领域内常碰到的混合，扩散等等问题，时均模型在解释这些过程的物理机制时常常显得力不从心.DNS, LES等方法，试图在空域和时域上都能更精确的贴近湍流过程,但以目前大众化的计算能力来看，用来分析工程实际似乎还是有较远的距离近期内这种矛盾如何更好地解决?这个问题，应该说关于湍流模拟的很多领域都存在。
现在对于各领域的湍流模拟，有两条路：一是通过简化的模拟，包括简化的模型、边界条件以及算法等，这样可以以较快的速度较小的代价获得计算结果，仍然保持了数值模拟能获得详细信息的优点，但是这些简化的方法一般都得基于高级的模拟技术或实验，适用性也需要认真考虑。