else if(temp>=855.&&temp<905.) {
if(rate>1.0e-4 && rate<1.e5)
trial=12830000./rate*log(pow((rate*4.7063),0.2817)+ pow((1.+pow((rate*4.7063),0.5634)),0.5))*pow(10.,-0.06 *(temp-855.)); else if (rate>=1.e5)
trial=243.654*pow(10.,-0.06*(temp-855.)); else
trial=1.47897e10*pow(10.,-0.06*(temp-855.)); }
else if(temp>=905.) {
if(rate>1.0e-4 && rate<1.e5)
trial=12830./rate*log(pow((rate*4.7063),0.2817)+ pow((1.+pow((rate*4.7063),0.5634)),0.5)); else if (rate>=1.e5) trial=0.24365; else
trial=1.47897e7; }
if(trial<1.e12&&trial>100.) mu_lam=trial; else if(trial<=1.) mu_lam=1.; else
mu_lam=1.e12; return mu_lam; }
自定义粘度与温度和应变率有关，因此无法使用软件自带的粘度定义形式
模型比较复杂，是在pro/E中建的模，然后用igs导入gambit，不过这样就产生了很多碎线和碎面并且在一些面交界的地方还存在尖角。我曾经做成功过把它们统统merge成一个虚面，中间设置了一个可以容忍尖角的参数，也可以划分网格，但把生成的msh文件导入fluent就会出错，这是virtual geometry的原因还是因为尖角的原因？还有，virtual geometry和普通的真实的几何体到底有什么区别？好像最大的区别是virtual geometry不能进行布尔操作，布尔操作（boolean operation）又是什么？使用virtual geometry需要注意哪些问题？
virtual geometry是很头疼的问题。你把它们统统merge成一个虚面 按理说全是虚的也是可以算的。可能是因为尖角的原因，虚实最大差别：是virtual geometry不能进行布尔操作，boolean operation即是并
对于复杂外形的网格生成，不可避免的会用到virtual geometry，virtual
face ,和virtual edge等，
1。作网格的时候，把所有的面全部合成一个虚面的做法不好，特别是对于复杂外形的网格生成，你最好在模型变化剧烈的地方多分几个面，这样会更有效的控制网格能够在模型表面曲率比较大的地方能够生成规则的结构或者非结构网格。
2对于你输入gambit的时候产生很多碎片的问题，你可以适当的把proe里面的模型精度和它的公差降低，因为gambit的建模工具精度本事就不高。 3。布尔运算就是对于面与面，体与体的联合，相减等运算。这个在所有的cad建模过程中是经常见到的问题。
4。对于虚体生成的计算网格，和实体生成的计算网格，在计算的时候没有区别，关键是看你网格生成的质量如何，与实体虚体无关。 我在作复杂模型计算的时候，大部分都是用的虚体，特别是从其他的建模软件里面导进来的复杂模型，基本上不能够生成实体。
至于计算的效果如何，那是你对于fluent的设置问题和网格的质量问题，与模型无关。
可以用gambit里面的check功能检查一下你的网格质量，看看质量怎么样
在fluent中会出现这么几个压力： Static pressure（静压） Dynamic pressure（动压） Total pressure（总压）
这几个压力是空气动力学的概念,它们之间的关系为: Total pressure（总压）= Static pressure（静压z） + Dynamic pressure（动压）
滞止压力等于总压(因为滞止压力就是速度为0时的压力,此时动压为0.) Static pressure（静压）就是你测量的，比如你现在测量空气压力是一个大气压
而在fluent中，又定义了两个压力： Absolute pressure（绝对压力） Relative pressure（参考压力） 还有两个压力：
operating pressure（操作压力） gauge pressure（表压） 它们之间的关系为: Absolute pressure（绝对压力）= operating pressure（操作压力） + gauge pressure（表压）
上面几个压力实际上有些是一一对应的，只是表述上的差别，比如： Static pressure（静压） gauge pressure（表压）
定义操作压力 对于可压缩流动: 把操作压力设为0
把表压看作绝对压力
1)coupled是耦合的意思，指同能量方程一起求解，而segrated是动量方程、压力方程和能量方程分开单独求解，迭代求解。一般能用耦合的时候尽量用，精度高。而分段并行求解一般精度低。
2)coupled solver从rampant发展来的，是与nasa合作开发的，适合于计算高马赫数（跨音和超音），为了克服低速下numeric stiffness，采用了preconditioning方法。couple implict内存耗费大，一般为segregated的1.5倍以上，you need powerful computer!segregated solver从uns发展而来，是基于simple算法，更适合于亚音场计算，收敛速度快，内存少。
3) fluent rampant是密度基时间推进法，故适合求解高速流动，也是这几个商用软件中唯一的密度基算法，但求解低速时一定要用preconditioning。
fluent UNS同其他商用软件如CFX4，5，TASCFLOW，star－cd一样都属于压力基方法，适用于低速不可压流动。
发信人: stipulation (华健哥~闷声大发财), 信区: NumComp 标 题: 建议大家用fluent网格最好用四面体网格
发信站: BBS 水木清华站 (Thu Apr 4 22:27:50 2002)
我最近一直忙于我得论文
是有限体积求解n－s方程的，结构网格
我鼓捣鼓捣程序的核心，特别是计算通量的时候
处理单位元面积，体积的算法的时候，想到fluent的非结构网格 我感觉如果是四面体网格，或者面元是三角形网格 应该计算起来误差小
我不知道事实上是不是这样的阿
简单说，结构化网格几乎处理的是六面体但愿，和四边形面元
计算面元，通过两个相邻边的差积的平均，这样有的网格划分如果质量不是很好 （正交性不好），造成的误差较大 而三角形网格就没有这个问题
因此我推想，如果fluent内处理有限体积的思想也是如此
那么三角形网格的计算误差应该比混和型（存在四边形面元）的要小
不知道是不是酱紫的，我做的用到fluent都在一些工程问题上，需要的精度不是很高 如果那个大侠做一些诸如湍流的模拟，模拟流场细微结构，可以检验一下我说的对不对
--
不一定吧，具我所知，六面体的网格在收敛的性能上明显比四面体网格要好。只要是 其网格质量比较好的话。还有在
后处理时，四面体网格的速度场的方向比较乱。个人认为高质量的六面体网格更加适合
有限体积法的求解。在湍流边界层的处理时，六面体网格的优点更加明显。
但是在对复杂形状的网格划分时，如何生成高质量的六面体网格是比较困难的。四面体
网格的适应性比较强，生成网格比较容易。 乱。
【 在 stipulation (华健哥~闷声大发财) 的大作中提到: 】 : 我最近一直忙于我得论文
: 是有限体积求解n－s方程的，结构网格
: 我鼓捣鼓捣程序的核心，特别是计算通量的时候
: 处理单位元面积，体积的算法的时候，想到fluent的非结构网格 : 我感觉如果是四面体网格，或者面元是三角形网格 : 应该计算起来误差小
: 我不知道事实上是不是这样的阿
: 简单说，结构化网格几乎处理的是六面体但愿，和四边形面元
: 计算面元，通过两个相邻边的差积的平均，这样有的网格划分如果质量不是很好 : （正交性不好），造成的误差较大 : 而三角形网格就没有这个问题 : ................... --
※ 来源:·BBS 水木清华站 smth.org·[FROM: 202.114.22.210] 发信人: stipulation (华健哥~闷声大发财), 信区: NumComp 标 题: Re: 建议大家用fluent网格最好用四面体网格 发信站: BBS 水木清华站 (Fri Apr 5 21:18:57 2002)
是
你说的有道理
可能我一直作比较复杂的外形几何体的模拟
所以我开始说的就默认了很难生成简单的高质量的六面体网格的情况 当然，诸如很规则的边界外形的流场，能做出高质量的六面体网格 肯定是六面体网格为首选
【 在 xiaoruofu (苦啊) 的大作中提到: 】
: 不一定吧，具我所知，六面体的网格在收敛的性能上明显比四面体网格要好。只要是
: 其网格质量比较好的话。还有在
: 后处理时，四面体网格的速度场的方向比较乱。个人认为高质量的六面体网格更加适合
: 有限体积法的求解。在湍流边界层的处理时，六面体网格的优点更加明显。
: 但是在对复杂形状的网格划分时，如何生成高质量的六面体网格是比较困难的。四面体
: 网格的适应性比较强，生成网格比较容易。 : 乱。
一.边界层
边界层是指定与边或者面相邻区域的网格接点的距离，目的是控制网格密度，从而控制感兴趣区域计算模型的有用的信息量。
例如：在一个液体流管中，我们知道靠近壁面处的速度剃度大，而中心处的速度剃度小，为了使得壁面处的网格密而中心处的网格稀疏，我们就在壁面处加一边界层。这样我们就能控制网个密度。
要定义一个边界层，你要定义以下参数： 1）边界曾依附的边或者面 2）指定边界层方向的面或者体 3）第一列网格的高度 4）相邻列之间的比例因子
5）总列数，指定了边界层的深度
同时，你也可以指定一个过度边界层。要指定一个过度边界层，你需要定义以下参数（过度模式，过度的列数）
1)生成边界层 需要定义以下参数：
i)size:包括指定第一列的高度和相邻列的比例因子
ii)internal continuity ：
当在体的某个面上施加边界层时，gambit会把边界层印在与这个面相邻的所有面上，如果在体的两个或者更多的面上施加边界层，那么边界层就有可能重叠，internal continuity 这个参数就决定了边界层如何重叠
当选择internal continuity 时，gambit不会在相邻的面上互相施加边界层 否则就会在相邻的面上互相施加边界层
具体可以看guide的图示：
同时这个参数还影响施加了边界层的体可以采用何种方式划分网格
corner shape ：
gambit 允许你控制conner（即边界层依附的的两条边的连接点处）附近区域网格的形状