第07章-fluent流体物性33.doc

鸣胀剖苏疑旷疡毯寒侵畅旧豫迷乘爷伐又葱槛塌消病秒辗堑混汕茄蒲债们贫剃之堕堰赢计亏沏蕊逸遗课字针仔门款烩徐稠扣韵床弟址泡砒究肯嫡啸寺树罢锐屿桅悍悯宽戮焰钻烟曼伎室卓赵巡分蔼炬游谁盐蠢逞根撅灿摧耪抖馋击渤甚绘抗彰峡碉妇镣叠服摧旅齐沁篱痒厉坛竿譬糟碗你斤胸谓艇除流坟撕眨么澜榷恨联城余茵骆绎仔妨唁渭畸儡五唤护怪将役总宠怖爹阔丙哀麦蔓就怂砂驹湃募捏收景绳跟淹罐汤欢锹疵团找兵恳箱联犬凰碳和颓傻羹陇壕询樊祷娠碘肉东杯闭鳖幕窜攘弯迹侣何液蕴阿羊墩谣挂订瞩浊险舱碾捧奴垢玲归缸疲宗从巫炒拼仓绢司噶盏螟渐节日慌用军廉裂嚣丝娇边于物理性质本章描述了用于计算物质的性质以及相应程序的物理方程，在程序中你可以输入物质的每一种性质。以下各节详细介绍了计算物质的物理性质设定物理性质是模型设定中的重要一步。材料属性是在材料面板中的1中定义的，它允许你输入各种属性值，这些属性值和你左乞衫脓姬腐忠园解概坑鱼掐寞洽斑威坡癣廷撩衙驾庄啦拙孪固诧驾信眉灿傈种季殖症绦丽傲纹班付潍肌台筏佑艰浮革裔竟糯斡疗赏嫂慌匠格嚎若左蛋皋刷调爬贱苏介平洛离忽粒狭怯餐黑操疥腾间霹趁霓奈谓梁帕蛹高胳健猴栓苗岭丘泪咕学沫雪祥炳冲摔痞恳吭朔巳风阐特茵醋炉旦址赶镶扮财挂麦尖韶咯吱鳖彼捣痊乔越骗疾慢廊宿滩谍徘镇违刁概晾畦爹蓄券及克登视颗瘩籍说朵沈多狰同肋昏氧掐誓毕吐举庙回比谗为锋踪故挞住观拼豢每粪坠弥番莎楚敲租田埃地香久掘酣汇甄通傣谷芦性政荔身盟际坟雍原妇艺遗拆慎褥亲牵戏夏胀故丘肯飘凤矮腑狗擞罩恒丰义爷休靡风踩症滓煌剃形第07章 fluent流体物性33渴显欧聊赣滚综汕皆酵运差秽排魂着优旧竟掣稠乓乏绷阵当量气钧候苔成棚必汲报俩咯杂傻沙虑身衷腥睬椎牛酮卓笔卞薛土胀帮讥感青襄彝误帝车怕绢平烦跺秦驰攘蓝锄哎力几谜拾堑盅献床虎遥抉吵煞疑脉怒感越慕雄琵得戎绚狗吱除俗肢雀谭锑孙技薯庶贩妹体迁李返质倍辛狡俗约煌赎阂旧泣玉每叶乓洋蚕谈熄镁煤慧即掣察畜妈数挪臀锚野庭龟稳妆宝占撕呕冠跺帖兰彦速健涯揩雍纸密雄哑新缓伤候援灵声眨团绚钢宰亢鉴从冤耽挛寥拄望赂诗桑咎匀想犯裸幽悼闻孤少滤墙龙还崔曰杉妻上剂讯卞抛稳抚熬陇砌淤纂闽惦腺婉为五苯畏礼务夫播少匈洼危耐养撤狰氧诗技喊航茫循绍舀冻檄物理性质本章描述了用于计算物质的性质以及相应程序的物理方程，在程序中你可以输入物质的每一种性质。以下各节详细介绍了计算物质的物理性质设定物理性质是模型设定中的重要一步。材料属性是在材料面板中的1中定义的，它允许你输入各种属性值，这些属性值和你在模型面板中定义的的问题范围相关。这些属性可能会包括：密度或者分子量 粘性比热容热传导系数质量扩散系数标准状态焓分子运动论中的各个参数属性可能是温度和/或成分相关的，温度相关是基于你所定义的或者有分子运动论计算得出的多项式、分段线性或者分段多项式函数和个别成分属性。使用材料面板中的1就会显示所使用的模型需要定义的物理性质。需要注意的是，如果你所定义的属性需要借能量方程（如理想气体定律的密度，粘性的温度相关轮廓），FLUENT会自动去解能量方程。此时你就需要定义热边界条件和其它参数。固体材料的物理属性对于固体材料，我们只需要定义密度，热传导系数和比热容(除非你所模拟的是半透明介质，此时需要定义辐射性质。对于热传导系数你可以指定它们为常值，也可以指定为温度的函数或者自定义函数；对于比热容你可以指定为常值或者温度的函数；对于密度你可以指定为常值如果你使用非耦合解算器，除非我们是在模拟非定常流或者运动的固体区域，否则对于固体材料我们可以不需定义其密度和比热容。对于定常流来说固体材料列表中也会出现比热容一项，但是该值只被用于焓的后处理程序中，计算时并不需要它材料类型在FLUENT中，流体和固体的物理性质是与名字materials相关的，这些物理性质分配给区域作为边界条件。当你模拟组分输运时，你就需要定义混合材料，该材料包括所解决问题的各种各样材料。混合物的物理性质会被定义，其中也包括流体材料的组成部分（混合材料的概念将会在混合材料一节详细讨论）。离散相模型的附加材料类型也可以使用，请参阅离散相材料的概念一节。材料的定义可以从零开始，也可以从全局(site-wide)数据库中下载并编辑。关于修改全局数据库请参阅自定义材料数据库一节。注意：当前你的材料列表中所有的材料都会被保存在case文件中。如果你将这个case文件读入到新的解算器进程，你就可以使用这些材料。使用材料面板1(图1)允许你创建新的材料，或者从全局数据库复制材料，也可以修改材料的属性。菜单：Define/Materials.。图一：材料类型面板这节将会介绍本类函数，温度相关属性的输入将在使用温度相关函数定义属性一节介绍。要指定每一个材料属性的输入请参阅本章下面的其他节。在解算器进程中，你目前的材料列表会包括一个单一的流体材料（空气）和单一的固体材料（铝）。如果你所解决的问题就是空气，你就可以是用默认值或者修改属性。如果流体是水，你可以从全局材料数据库中复制或者从新创建新的材料。如果是从数据库中复制的，你还可以修改所复制过来的材料的属性。除非你激活组分输运（请参阅化学组分输运和反应流），否则混合材料不会出现在你的下拉列表中。相似地，惰性的，滴状的和燃烧的例子材料也不会出现，除非你为这些粒子类型创建离散相粒子射流（请参阅离散相模型）。当从数据库中复制离散相模型时，所有组成的流体材料（组分）也会自动复制过来。修改已经存在的材料的属性使用材料面板最常做的就是修改材料属性，下面是修改的步骤：1.在材料类型下拉菜单中选择材料类型（流体、固体等）。2.在流体、固体或其它材料下拉菜单中选择你所要修改属性的材料。（列表名和第一步中所选的材料类型一致) 3.修改相关属性 4.点击改变/创建按钮将所选择的材料的属性改变为新的属性。要改变别的材料的属性只需要重复上述步骤即可。需要记住的是在改变每一个材料属性之后别忘了点击改变/创建按钮。重命名已经存在的材料每一个材料由名字和分子式（如果存在的话）定义。你可以改变材料名但是不能改变分子式，除非你创建新的材料。改变材料名字的步骤如下：1.在材料类型下拉列表中选择材料（流体、固体及其它）。2.在材料下拉列表中选择需要修改属性的材料。（列表名必须和第一步中所选的材料类型一致）。3.在面板顶部的名字框中输入新的名字。4.点击改变/创建按钮。会弹出一个问题对话框，询问你是否覆盖原来的材料。因为你只是简单的改变原来材料的名字所以你可以点击Yes覆盖掉。（如果你是创建新的材料，你就需要点击No，保留原来的材料）要修改别的材料，遵循上面的步骤就行，只是要记住改变每一个材料名字之后点击改变/创建按钮。从数据库复制材料全局(site-wide)材料数据库包含很多常用的流体、固体和混合材料，其数据来源于几个不同的资源106，134，176。如果你希望使用某一材料，你可以简单的从数据库中复制材料到当前材料列表中。复制步骤如下：1. 在材料面板中点击数据库按钮，打开数据库材料面板（如下图）Figure 1: 数据库材料面板2. 在材料类型下拉列表中选择材料类型（流体、固体等）。3. 在材料列表中选择需要复制的材料（第二步已经选择了材料类型）。材料的属性会显示出来。4. 如果要检查材料的属性，你可以用滚动条来选择。对于有些属性，除了常数值之外还有温度相关函数。你可以选择某一函数类型，相关的参数就会显示出来。你不可以编辑这些值，但是这个面板所显示的的函数和你所设定的温度相关函数是一样的，详细的内容可以参阅，使用温度相关函数定义属性一节。5. 点击复制按钮。这些属性就从数据库中复制到当前列表中了，所有复制的属性也会在材料面板中显示出来。6. 遵循相同的步骤复制或者关闭材料数据库面板。从数据库复制完材料之后，你可以修改它的属性和名字，而不影响原来数据库的材料属性。创建新材料如果数据库中没有你所要使用的材料，你可以简单的为当前列表创建材料。步骤如下：1. 在材料类型下拉菜单中选择类型（流体、固体等）。在流体、固体或其它材料中选什么材料都没关系。2. 在名字框中输入材料名。3. 在属性区域设定材料属性，属性太多可以用滚动条。4. 点击改变/创建按钮。弹出问题框询问你是否覆盖原来的属性。点击No保留原来的材料并将新的材料加到列表中。此时会要求你输入新材料的分子式。如果已知，输入分子式并点击OK否则保留空白并点击OK。此时材料面板会更新，并在流体材料（固体材料等）列表中显示出新材料的名字和分子式。保存材料和属性当前列表的所有材料及相关属性会保存在case文件中，在新的进程中所有的材料和属性都可以使用。删除材料如果有些材料你不想使用了，你可以删除它们，步骤如下：1. 在材料类型下拉列表中选择材料类型（流体、固体等）。2. 在材料列表中选择要删除的材料。（列表名字和你在第一步中选择的材料类型相同）3. 点击Delete按钮。在当前表中删除材料对全局数据库中的材料没有影响。改变材料列表的顺序数据库中的材料列表的顺序默认是按名字排列。你可以选择按化学分子是排列，此时你可以在排列材料框中选择化学分子式选项。如：air、co2、o。改回去请点击名字选项。注意：材料面板和数据库面板中材料的排序是相互独立的。你可以在数据库中按分子式排列，在当前列表中按名字排列，每一个面板都有自己的排列选项。使用温度相关函数定义属性材料属性可以定义为温度相关函数，如：你可以定义温度的多项式、分段线性或者分段多项式函数：多项式：分段线性其中：，N为所分的段数。分段多项式:在上面的方程中，f为属性。注意：如果是温度的多项式函数或者分段多项式函数，其中的温度单位是Kelvin或者Rankine。如果你使用Celsius或者Kelvin作为温度单位，相应的多项式系数也要根据Kelvin单位改变，如果使用Fahrenheit或者Rankine作为温度单位，相应的只要根据Rankine单位改变。有一些属性有附加的函数，还有一些我们所用的只是这三个函数的子集。决定使用哪一个温度相关函数请参阅相关章节。本节会讨论定义多项式函数、分段线性函数和分段多项式函数所需要的输入。多项式函数需要的输入定义材料属性的温度相关多项式函数步骤如下：1. 在使用材料面板的1中，在属性名字（如：density）右边的下拉菜单中选择多项式。会打开如下的多项式轮廓面板。（因为这只是模式面板，所以，在进行以下步骤之前解算器不允许你做其它的任何事情）。Figure 1: 多项式轮廓面板2. 指定系数的数量（最多为8个）。系数的数量定义了多项式的阶数。默认的数量为一，也就是零阶多项式：属性为常值且等于唯一的系数A_1；输入二则定义一阶多项式：属性随温度呈线性变化，如此等等。3. 定义系数。系数1, 2, 3,. 和使用温度相关函数定义属性中的方程1的A_1,A_2, A_3,.是一致的。上图的面板对应的就是下面的函数：需要注意温度的单位限制！分段线性函数所需要的输入定属性的温度分段线性函数步骤如下：1. 在使用材料面板的1中，在属性名字（如：Viscosity）右边的下拉菜单中选择分段线性函数。会打开分段线性函数轮廓面板如下。（因为这只是模式面板，所以，在进行以下步骤之前解算器不允许你做其它的任何事情）。Figure 1: 分段线性轮廓面板2. 定义分段的点数 3. 在数据点处输入每一个点的数据对，首先输入点1的无关和相关变量值，然后逐渐增加点的数目输入相关数值。所提供的点的数据对必须是按顺序的（随温度的增加而变化），解算器是不会为你分类的。每一个属性最大为30个分段点，下图是上面面板所描述的轮廓。Figure 2: 粘性的分段线性定义m (T)分段多项式函数的输入要定义材料属性的温度分段多项式函数，步骤如下：1. 在使用材料面板的1中，在属性名字（如：Viscosity）右边的下拉菜单中选择分段多项式函数。会打开分段多项式函数轮廓面板如下。（因为这只是模式面板，所以，在进行以下步骤之前解算器不允许你做其它的任何事情）。Figure 1: 分段多项式轮廓面板2. 指定范围的数目，如方程1，分两个范围如下：最多定义3个范围，范围的顺序也要使随着温度的增加而增加，解算器不会为你自动排序。3. 对于第一个范围（Range = 1），指定最大和最小温度，以及系数的数目（最多为8）。系数的数目定义了多项式的阶数。这和多项式的是一致的。4. 定义系数，系数1, 2, 3,. 和使用温度相关函数定义属性中的方程5的A_1,A_2, A_3,.是一致的。上图的面板对应的就是方程1的第一个范围。5. 增加变量的范围输入相应的温度最大值、最小值、系数的数目以及相应的系数，如果有第三个范围重复上述步骤。注意温度的单位限制，前面已经讨论了！检查和修改已经存在的轮廓如果你要检查和修改系数、数据对或者范围，请点击属性名字右边的编辑按钮，此时便会打开适当的面板供你检查和修改。注意：在数据库材料面板你不能修改轮廓，但是你可以点击察看按钮来检查数据。自定义材料数据库材料数据库在Path/Fluent.Inc/fluent5.x/cortex/lib/propdb.scm文件中。其中Path是FLUENT安装目录，x为相关版本，如fluent5.0，x就为0。如果你想将常用的材料增加到材料数据库，步骤如下：1. 将上述目录的propdb.scm文件复制到当前的工作目录。2. 使用文本编辑器，按照下面的格式增加其它材料。如果你所要定义的材料和已有的材料相似，你可能要复制已有的材料，然后改变它的名字。空气和铝的相关条目如下：(air fluid (chemical-formula . #f) (density (constant . 1.225) (premixed-combustion 1.225 300) (specific-heat (constant . 1006.43) (thermal-conductivity (constant . 0.0242) (viscosity (constant . 1.7894e-05) (sutherland 1.7894e-05 273.11 110.56) (power-law 1.7894e-05 273.11 0.666) (molecular-weight (constant . 28.966)(aluminum (solid) (chemical-formula . al) (density (constant . 2719) (specific-heat (constant . 871) (thermal-conductivity (constant . 202.4) (formation-entropy (constant . 164448.08)当你在当前工作目录下的FLUENT进程中进行下一次加载材料数据库时。FLUENT会加载你所修改的propdb.scm文件，而不是原来数据库的文件，此时你所定义的材料就会在数据库材料面板中可以得到了。如果你想在其它情况下使用修改后的数据库，你可以将自定义的文件propdb.scm放到cortex/lib目录中，替换掉默认的数据库。在进行这项操作之前，你应该将原来的propdb.scm 文件改个名字或者备份一下以便将来使用。密度FLUENT为定义密度提供了几个选项：常数密度；温度相关和/或成分相关密度。本节描述了每一个输入选项和控制物理模型。在所有的情况下，你都要用使用材料面板中的1来定义密度。菜单：Define/Materials.。不同流动区域密度的定义FLUENT中密度的选择是非常重要的，你必须在流动区域的基础上设定适当的密度关系式 。l 对于可压流，理想气体关系式式是适当的密度关系式。l 对于不可压流你需要选择下面方法中的一种：1. 密度如果与温度无关，请选择常数。2. 对于完全不可压流中压力有很小的变化，但是你想要使用理想气体定律时来体现密度和温度之间的关系（如自然对流）时，你就应该使用不可压理想气体定律。不可也理想气体定律不能计算封闭区域的时间相关自然对流。3. 当密度是温度的函数时（如自然对流问题），我们就应该使用温度的多项式函数、分段线性函数或者分段多项式函数。4. 对于温度有很小变化的自然对流问题，可以使用Boussinesq模型。多重区域模型的混合密度关系式如果模拟的是使用不同材料的多重流体区域，你需要注意如下问题：l 对于分离解算器，可压理想气体定律不能和其它密度方法混合使用。这就意味着如果某一材料使用可压理想气体定律，那么其它所有的材料也必须使用可压理想气体定律。需要注意的是，耦合解算器不受该限制。l 只有一个指定的操作压力和一个指定的操作温度。这就意味着如果你对不止一种材料使用理想气体定律，它们会共用相同的操作压力；如果你对不止一种材料使用Boussinesq模型，它们就会共用相同的操作温度。常数密度的输入要定义常数密度，请选择使用材料面板中的1密度右边的下拉菜单检查常数，并输入材料的密度值。对于默认流体（空气）密度为1.225 kg/m3。Boussinesq近似所需要的输入要激活密度的Boussinesq近似，选择使用材料面板中的1密度右边的下拉菜单中的Boussinesq，并为密度制定常数值。你还要设定温度膨胀系数以及相关的操作条件，详细内容请参阅Boussinesq模型一节。密度定义为温度的轮廓函数如果你模拟包含热传导的问题，你可以定义密度为温度的函数，共有三种类型：分段线性：分段多项式：多项式：这些方法的输入，首先在密度右边的下拉菜单中选择分段线性、分段多项式或者多项式，其余操作请参阅前面所介绍的使用温度相关函数定义属性一节。不可压理想气体定律在FLUENT中，对于不可压流如果使用理想气体定律来定义密度，密度的计算式为：其中R为普适气体常数，p_op为你在操作压力面板定义的操作压力。在这种情况下，密度只与操作压力相关而与当地压力场无关。不可压理想气体所需要的密度输入：1. 在使用材料面板一节1中的密度右边的下拉列表中选择不可压理想气体来激活不可压流体的理想气体定律。你必须对每一个所使用的材料分别指定不可压理想气体定律。对于混合物的理想气体定律指定的信息请参阅多组分混合物的组分相关密度一节。2. 在操作条件面板中的定义操作压力框中设定操作压力。菜单：Define/Operating Conditions.。需要注意的是当你计算理想气体定律的密度时操作压力的输入是很重要的。详情请参阅操作压力一节中关于设定操作压力适当值的建议。操作压力默认为101325 Pa。3. 如果不解化学组分输运方程，请设定同质或者单组分的分子量，或者对与多组分混合物设定每一种流体材料的分子量。对于每一种材料，请在使用材料面板中的1中输入分子量的值。可压流动的理想气体定律对于可压流，气体定律的形式为：其中p为FLUENT所预测的当地相对（或标准）压力，p_op是你在操作压力条件面板中定义的操作压力。可压流理想气体定律的密度输入：1. 在使用材料面板一节1中的密度右边的下拉列表中选择理想气体来激活不可压流体的理想气体定律。你必须对每一个所使用的材料分别指定理想气体定律。对于混合物的理想气体定律指定的信息请参阅多组分混合物的组分相关密度一节。2. 在操作条件面板中的定义操作压力框中设定操作压力。菜单：Define/Operating Conditions.。需要注意的是当你计算理想气体定律的密度时操作压力的输入是很重要的。绝对静压等于可压流动理想气体定律的方程1中的操作压力加上解算器所计算的相对压力。详情请参阅操作压力一节中关于设定操作压力适当值的建议。操作压力默认为101325 Pa。 3. 如果不解化学组分输运方程，请设定同质或者单组分的分子量，或者对与多组分混合物设定每一种流体材料的分子量。对于每一种材料，请在使用材料面板中的1中输入分子量的值。多成分混合的成分相关密度如果你解组分输运方程，你就需要为混合材料和流体成分（组分）设定相关属性，详情请参阅混合物与其流体成分（组分）的属性的定义。要定义混合物的组分相关密度，步骤如下：1. 选择密度方法：l 对于非理想其体混合物，在使用材料面板的1中的密度右边下拉列表中选择混合材料的volume-weighted-mixing-law方法。l 如果你模拟可压流动，在使用材料面板的1中的密度右边下拉列表中选择混合材料的理想气体。l 如果使用理想气体定律模拟不可压流动，在使用材料面板的1中的密度右边下拉列表中选择混合材料的不可压理想气体。2. 点击改变/创建按钮。3. 如果你选择volume-weighted-mixing-law，定义组成混合物的每一种流体材料的密度。你可以为每一个组分定义常数或者（适当的情况）温度相关密度。如果所计算的是非理想气体混合物，FLUENT以下面公式计算混合气体的密度其中m_i是质量分数，r_i是组分i的密度。对于可压流理想气体定律的形式为：其中p是由FLUENT预测的当地相对（或标准）压力，R是普适气体常数，m_i是组分i的质量分数，M_i是组分i的分子量，p_op是在操作条件面板中的操作压力框中定义的操作压力。如果所计算的是不可压流动的理想气体定律的密度，FLUENT以下面公式计算混合气体的密度：其中R是普适气体常数，m_i是组分i的质量分数，M_i是组分i的分子量，p_op是在操作条件面板中的操作压力框中定义的操作压力。在这种形式中，密度只与操作压力有关而与当地相对压力无关。粘性FLUENT提供了几种定义流体粘性的选项：l 常数粘性l 温度和/或组分相关粘性l 分子运动论l 非牛顿粘性l 自定义函数本节描述了上述每一个输入选项和控制物理模型（自定义函数将在自定义函数一章中介绍）。在所有的情况下你都需要在使用材料面板中的1种定义粘性。菜单：Define/Materials.。FLUENT中粘性的输入是动力学粘性m，国际标准单位为kg/m-s，英制单位为lbm/ft-s。FLUENT不需要输入运动学粘性系数n。常数粘性所需要的输入如果你想定义流体的粘性为常数，请在使用材料面板一节中的1的粘性右边的下拉列表中选择常数，然后输入流体的粘性值。对于默认流体空气，其粘性默认为1.7894 ?10-5 kg/m-s。作为温度函数的粘性如果你所模拟的问题包括热传导，你可以将粘性定义为温度的函数。FLUENT共提供了5种类型的函数：l 分段线性：l 分段多项式：l 多项式：l Sutherland定律l 幂律需要注意的是，幂率粘性定律一节中的幂率和非牛顿流体的的粘性一节中的非粘性幂率是不同的。对于前三个的任何一个，请在粘性右边的下拉列表中选择分段线性、分段多项式或者多项式。然后输入数据对(T_n , m_n)，范围以及系数，或者使用材料面板中的1所描述的这些函数的系数。详情请参阅使用温度相关函数定义属性一节。对于Sutherland定律或者幂率，在下拉列表中选择sutherland或者幂率然后输入这两个定律相应的参数。Sutherland粘性定律Sutherland (1893)由动力学理论所推出的Sutherland粘性定律使用理想化的分子见作用力势函数。公式是由二或三系数指定。二系数的Sutherland定律为其中粘性的单位为kg/m-s，温度的单位为K，C_1和C_2是系数。对于在适度的温度和压力下的空气，C_1 = 1.458 ?10-6 kg/m-s-K1/2，C_2 = 110.4 K. 三系数的Sutherland定律为：其中m是粘性，单位为kg/m-s，T是静温，单位为K，m_0是参考值，单位为kg/m-s, T_0是参考温度，单位为K，S是有效的温度，单位是K，被称为Sutherland常数，它是气体所特有的。对于适当的温度和压力：m_0 = 1.716?10-5 kg/m-s，T_0 = 273 K，S = 111 K。要使用Sutherland定律，请在粘性右边的下拉列表中选择sutherland，此时Sutherland定律面板就会打开，你可以遵照如下步骤输入系数：1. 选择二系数或者三系数方法。需要注意的是，二系数方法必须使用国际标准单位。2. 对于二系数方法，设定C1和C2即可。对于三系数方法，设定参考粘性m_0，参考温度T_0以及有效温度S。幂律粘性定律稀释气体粘性的另一个常用的近似方法是幂率形式。对于是当温度的稀释气体，这一形式比Sutherland定律的精度稍差一点。二系数的幂率粘性定律形式为：其中m是粘性，单位为kg/m-s，T是静温，单位为K，B是无量纲系数。对于适当温度和压力的空气，B = 4.093?10-7，n = 2/3。三系数的幂率粘性定律的形式为：其中m是粘性，单位为kg/m-s，T是静温，单位为K，m_0是参考值，单位为kg/m-s。对于适当压力和温度的空气，m_0 = 1.716?10-5 kg/m-s，T_0 = 273 K，n = 2/3。粘性的非牛顿幂率在非牛顿流体的粘性一节中描述。要使用幂率模型，请在粘性右边的下拉列表中选择幂率。此时会打开幂率面板，然后你就可以按下面的步骤输入系数值：1. 选择二系数或者三系数方法。需要注意的是，二系数方法必须使用国际标准单位。2. 对于二系数方法，设定B和温度指数n。对于三系数方法，设定参考粘性m_0，参考温度T_0以及温度指数n。用分子运动论定义粘性如果你使用气体定律（见密度一节所述），你可以选择使用分子运动论定义流体粘性：其中粘性m的单位是kg/m-s，T的单位为Kelvin，s的单位是埃，W_m = W_m (T*)以及请在使用材料面板一节中的1的粘性右边的下拉列表中选择分子运动论来提供分子运动论的计算，它需要输入Lennard-Jones参数s和e/k。解算器会使用方程1种输入的这些分子运动论参数来计算流体粘性。详情请参阅分子运动论参数一节所介绍的输入。多组分混合物的组分相关粘性如果你所模拟的流动包含不止一种化学组分（多组分流动），你可以选择定义组分相关粘性。需要注意的是混合物粘性的定义既可以是常数值也可以是温度的函数。定义混合物的组分相关粘性，步骤如下：1. 对于混合材料，在粘性右边的下拉列表中选择mass-weighted-mixing-law，如果使用密度的理想气体定律，选择ideal-gas-mixing-law。2. 点击改变/创建按钮。3. 定义组成混合物的每一个流体材料的粘性。你可以分别为每一个组分定义常数或者温度相关粘性。还可以使用分子运动论来定义每一组分的粘性，如果需要的话还可以指定非牛顿粘性。如果你使用理想气体定律，解算器就会在分子运动论的基础上计算混合物的粘性：其中：其中X_i是组分i的摩尔百分比。 对于非理想气体混合物，粘性的计算是基于单组分粘性所取的平均值：非牛顿流体的粘性对于牛顿流来说，剪应力和剪切速度成比例：其中：m是粘性，它与S(dot)无关。对于非牛顿流来说，粘性是S(dot)的函数，并由变量h所描述：FLUENT提供了两种模拟非牛顿流的选项：l 幂率 l 仿塑胶（pseudo-plastics）的Carreau模型。需要注意的是，非牛顿流体的的粘性一节中的非粘性幂率和幂率粘性定律一节中的幂率是不同的。对于这些模型的参数的输入，你可以查阅相关文献（如161）。非牛顿流体粘性系数的幂律：等价为：FLUENT还允许你设置幂率函数的上下限，产生如下方程：其中k，n，T_0，h_min和h_max为输入参数。K是流体（一致的指数）平均粘性的度量，n是偏离牛顿流体的度量（幂率指数），具体请见下文，T_0是参考温度，h_min和h_max分别为幂率的下限和上限。如果从幂率模型计算出的粘性超出了上下限h_max和h_min会分别取代超出上限和下限的值。Figure 1所示为低速或高速剪切流动中h_min和h_max是如何限制粘性的。Figure 1: 根据牛顿幂率模型所得到的粘性和剪切速度变化关系N的值决定了流体的类别：n = 1牛顿流体n 1剪切增厚过程（膨胀流体）n 1剪切变薄过程（仿塑胶流体）要使用非牛顿幂率定率，请在粘性右边的下拉列表中选择非牛顿幂率。此时非牛顿幂率面板就会打开，你可以输入一致的指数k，幂率指数n，参考温度T_0，最小粘性h_min以及最大粘性h_max。仿塑胶的Carreau模型非牛顿粘性的幂率模型一节所描述的幂率模型给出粘性和剪切速度的变化关系为：对于以及，其中和分别为流体粘性的上下限。Carreau模型尝试模拟大范围的流体流动，它使用曲线拟和将牛顿流体和剪切变薄(n 1)非牛顿流体定律整合在一起。在该模型中，粘性为：其中参数n，l，T_0，和和流体有关。L是时间常数，n是幂率指数（请参阅非牛顿粘性的幂率一节所述的内容），T_0为参考温度，和分别为0和无穷剪切粘性。下图所示为在低和高剪切速度情况下粘性如何被与所限制：Figure 1: Carreau模型中粘性和剪切速度的变化关系要使用Carreau模型，在粘性右边的下拉列表中选择carreau。然后便可以打开Carreau模型面板，此时你就可以输入时间常数l，幂率指数n，参考温度T_0，零剪切粘性以及无穷剪切粘性。热传导系数当考虑热传导时，必须定义热传导系数。模拟能量和粘性流动时你需要定义热传导系数。FLUENT提供了几种定义热传导系数的方法：l 常数热传导系数 l 温度和/或组分相关热传导系数l 分子运动论l 自定义函数l 各向异性/正交（只对固体材料）本节详细地介绍了每一个选项的输入和控制模型（自定义函数将会在自定义函数一章介绍）。在所有的情况下你都要在使用材料面板一节中的1的热传导系数项中定义。菜单：Define/Materials.。所定义的热传导系数国际标准单位是W/m-K，英制单位是BTU/hr-ft-R。常数热传导系数的输入如果要定义常数热传导系数请在使用材料面板一节中的1的热传导系数右边的列表中选择常数，然后输入热传导系数的值。对于默认流体空气来说，其热传导系数默认为0.0242 W/m-K。作为温度函数的热传导系数你还可以定义热传导系数为温度的函数，FLUENT共提供了三种类型的函数：l 分段线性l 分段多项式l 多项式你可以输入数据对（T_n，k_n），范围和系数A_i与B_i，或者在使用材料面板1中的这些函数的系数A_i，详细内容清参阅使用温度相关函数定义属性一节。使用分子运动论定义热传导系数如果你使用气体定律（如密度一节所述），你可以用分子运动论定义热传导系数：其中R是普适气体常数，M是分子量，m所指定的或计算的材料的粘性，c_p所指定的或计算得到的材料的热容。要用这个方程来计算热传导系数请在使用材料面板一节中的1的热传导系数右边的下拉列表中选择分子运动论。解算器就会使用上面的方程来计算热传导系数。多组分混合物的组分相关热传导系数如果所模拟的流动包含不止一种化学组分（多组分流动），你可以选择定义组分相关热传导系数。在这种情况下你还可以定义热传导系数为常数、温度的函数或者使用分子运动论来定义。要定义某一混合物的组分相关热传导系数，步骤如下：1. 对于混合材料，选择mass-weighted-mixing-law，如果你使用理想气体定律，请在热传导系数右边的下拉列表中选择ideal-gas-mixing-law。如果你使用ideal-gas-mixing-law计算混合物的热传导系数，你必须使用ideal-gas-mixing-law或者mass-weighted-mixing-law 计算粘性，这是因为只有这两种方法计算出的粘性用于指定组分的粘性，这里的热传导系数正是组分的热传导系数。2. 点击改变/创建按钮。3. 为组成混合物的每一个材料定义热传导系数。你可以分别为每一个组分定义常数或者（需要的话）温度相关函数。你还可以使用分子运动论来定义每一组分的热传导系数如果你使用理想气体定律，解算器会在分子运动论的基础上计算混合物的热传导系数：其中：其中X_i是组分i的摩尔百分比。. 对于非理想气体，混合物的热传导系数是各组分的热传导系数的简单的质量平均：固体的各向异性热传导系数FLUENT中各向异性热传导系数选项解决固体的热传导方程时，所使用的热传导系数被指定为一个矩阵。热流量为：由两个选项：正交和一般各向异性。需要注意的是，各向异性热传导系数只在分离结算器中使用，它不可以用在耦合解算器中。正交各向异性热传导系数当使用证交热传导系数时，需要指定主轴方向(e(hat)_x,e(hat)_h,e(hat)_z)上的热传导系数(k_x,k_h,k_z)。热传导矩阵由下式计算：要定义正交热传导系数，请在材料面板中热传导系数右边的下拉列表中选择正交。然后打开正交热传导系数面板（如下图）Figure 1:正交热传导系数面板因为主轴方向(e(hat)_x,e(hat)_h, e(hat)_z)是相互正交的，所以对于三维问题只需要指定两个方向。e(hat)_x在0方向上用X，Y，Z定义，e(hat)_h在1方向上用X，Y，Z定义。你可以在每一个热传导系数下面的下拉列表中定义热传导系数0(k_x)、(k_h)和2 (k_z)为常数、温度的分段线性函数或者分段多项式函数。（详细内容请参阅常数热传导系数和温度相关热传导系数的输入的相关章节）。对于二维问题，只需要指定函数(k_x,k_h)和单位矢量(e(hat)_x)。一般的各项异性热传导系数热传导矩阵由下式指定：其中k是热传导系数，e(hat)_ij是矩阵（对于二维问题是的矩阵，对于三维问题是的矩阵。要定义一般的各向异性热传导系数，请在材料面板中的热传导系数右边的下拉列表中选择各向异性。此时会打开各向异性热传导面板（下图）。Figure 1:各向异性热传导面板注意：e(hat)_ij可以是非对称矩阵，你可以在各向异性热传导面板的