多孔介质边界条件.docx
- 文档编号:24175119
- 上传时间:2023-05-25
- 格式:DOCX
- 页数:22
- 大小:156.83KB
多孔介质边界条件.docx
《多孔介质边界条件.docx》由会员分享,可在线阅读,更多相关《多孔介质边界条件.docx(22页珍藏版)》请在冰豆网上搜索。
多孔介质边界条件
7.19多孔介质边界条件
多孔介质模型适用的范围非常广泛,包括填充床,过滤纸,多孔板,流量分配器,还有管群,管束系统。
当使用这个模型的时候,多孔介质将运用于网格区域,流场中的压降将由输入的条件有关,见Section7.19.2.同样也可以计算热传导,基于介质和流场热量守恒的假设,见Section7.19.3.
通过一个薄膜后的已知速度/压力降低特性可以简化为一维多孔介质模型,简称为“多孔跳跃”。
多孔跳跃模型被运用于一个面区域而不是网格区域,而且也可
以代替完全多孔介质模型在任何可能的时候,因为它更加稳定而且能够很好地收
敛。
见Section7.22.
7.19.1多孔介质模型的限制和假设
多孔介质模型就是在定义为多孔介质的区域结合了一个根据经验假设为主的流
动阻力。
本质上,多孔介质模型仅仅是在动量方程上叠加了一个动量源项。
这种
情况下,以下模型方面的假设和限制就可以很容易得到:
•因为没有表示多孔介质区域的实际存在的体,所以fluent默认是计算基
于连续性方程的虚假速度。
做为一个做精确的选项,你可以适用fluent
中的真是速度,见section7.19.7。
•多孔介质对湍流流场的影响,是近似的,见7.19.4。
•当在移动坐标系中使用多孔介质模型的时候,fluent既有相对坐标系也可以使用绝对坐标系,当激活相对速度阻力方程。
这将得到更精确的源项。
相关信息见section7.19.5和7.19.6。
•当需要定义比热容‘’的时候,'’必须是常数。
7.19.2多孔介质模型动量方程
多孔介质模型的动量方程是在标准动量方程的后面加上动量方程源项。
源项包含
两个部分:
粘性损失项(达西公式项,方程7.19-1右边第一项),和惯性损失
项(方程7.19-1右边第二项)
式中,si是i(x,y,z)动量方程的源项,U是速度大小,D和C是矩阵。
动量源项对多孔介质区域的压力梯度有影响,生成一个与速度大小(速度平方)成正比的压降。
对于各向同性多孔介质简单情况下:
(7.19-2)
式中门是渗透性系数,(是惯性阻力系数,也就是将D,C矩阵简化为对角矩阵,
对角上的系数分别为,其它元素都是0.
同样fluent也可以将源项设定为速度的幕函数型:
多孔介质中的达西定律
通过多孔介质的层流,典型的压降是与速度大学成正比,常数C2可以认为是0
忽略对流加速度和扩散,多孔介质的动量方程源项就可以化简为达西定律:
-(7.19-4)
a
坐标轴三个不同方向的压降fluent计算如下:
1/砌是方程7.19-1中矩阵D的项,勺是£,V,andZ方向速度,
住琢,and埶是-,"',and-方向的多孔介质厚度.
多孔介质材料惯性损失
在咼速流动时,方程7.19-1中常数C2是对多孔介质中关于惯性损失的修正。
这个常数被认为是流动方向单位长度的损失系数,压降定义为动水头的函数。
如果模拟多孔板或者管束系统,有时候可以忽略渗透项,只使用惯性损失项,就得到如下的多孔介质压降方程:
“三-字4(訓闸)(7.19-6)
如下x,y,z方向的压力损失项:
~31
LC%加陶例亦|
thicknessyouhavedefinedinyourmodel.
7.19.3多孔介质能量方程
Fluent在多孔介质区域求解能量输运方程,并且修改了传导通量和瞬态项。
在多孔介质区域,传导通量使用一个有效的传导率,瞬态项包含了多孔介质中的固体区域的热惯性。
伽厨+(1-勺诅)+V何丹£沖))7切可-(工抵Ji)+(诃+於(7.19-8)
式中
,-=流体总能量尽=固体区域总能量
=多孔介质孔隙率
.=多孔介质的有效热传导率
=流体焓源项
多孔介质的有效热传导率多孔介质的有效热传导率:
是根据区域中流体热传导和固体热传导根据体积平均得到:
兔匸'池[n吶址(7.19-9)
式中
=多孔介质的孔隙率
k},&)流体项热传导率(包含湍流影响)
=固体区域热传导率
流体,固体区域热传导都可以通过udf自定义。
非各向同性热传导同样也可以通过udf进行定义。
在这种情况下,流体各向同性
的影响‘‘加到固体各向异性矩阵的对角元素上。
7.19.4多孔介质模型中的湍流模型
默认情况下,fluent会求解标准守恒湍流方程。
这时,固体区域对湍流的生成和耗散率没有影响。
当多孔介质的渗透性很好而且其几何尺度相比于湍流涡的尺度很小的时候是可行的。
其它一些算例,可能还需要抑制多孔介质区域的湍流影响。
如果使用某种湍流模型,除了大涡模型,可以通过设定湍流粘性率为0来消除湍流的影响。
当使用这种功能的时候,fluent将会把入口湍流输运通过多孔介质,而他们对流体混合和动量的影响将会被忽略。
另外多孔介质区域的湍流生成率也被设定为0.选定fluidpanel面板的laminarzone选项,这个模型将会被激活。
激活这个模型将暗示着■是0而且多孔介质区域湍流生成率也为0。
不选这个选项(默认)将会按照计算主流区域的计算方法来计算多孔介质的湍流。
Laminar
zone选项的详细介绍见section7.17.1。
7.19.5多孔介质模型的非定常项
,式
非定常多孔介质计算,多孔介质对时间微分项的影响包含了所有的标量输运方程
和连续性方程。
当考虑孔隙率影响的时候,时间微分方程项变成
孔隙率的影响是自动加入的,默认状况时孔隙率是设定为1.
Dij
inEquation
7.19.6多孔介质模型的用户输入
1.
定义多孔介质区域。
2.
定义多孔介质速度方程(可选)。
3.
选择流过多孔介质区域的流体材料。
4.
激活多孔介质区域的化学反应,如果有合适的反应机理,
选择化学反应机
理。
5.
默认这个选项是激活的,而且可以考虑移动的多孔介质。
见
section
7.19.6。
当使用多孔介质模型的时候,需要另外输入的问题部分如下。
选项如下:
6.
设定粘性阻力系数(
7.19-1,or
"in
Equation7.19-2),还有定义他们应用的方向矢量.或者定义幕函数模型的系数.
7.设定多孔介质区域的孔隙率。
8.为多孔介质选择材料,这种情况只可能发生在热传导模型中。
而且材料比热容只能设定为常数。
9.(optional)设定固体多孔介质部分体积热生成率,或者其它源项如动量,
质量。
10.(optional).设定流体区域的固定值。
11.如果合适,抑制多孔介质区域湍流度。
12.如果有必要,设定旋转轴,或者区域运动等。
设定阻力系数或者渗透系数方法如下。
如果选择幕函数近似来定义多孔介质动量源项,你需要输入C0,C1两个系数而不是阻力系数与流动方向。
7.1.4)的
可以在BoundaryConditionspanel(asdescribedinSectionFluidpanel(Figure7.19.1)设定多孔介质的所有参数,.
Fluid
ZaneName
*PorousZone
LaminarZone
SourceTerms
FixedValues
臣Reaction
MotionPorous2oneRodionSourceTermsFixedValues
—勺
UpdateFromLineTool
Diredion-1Vector
0K
Cancel
Help
Figure7.19.1:
TheFluidPanelforaPorousZone
定义多孔介质区域
就像7.1节描述的那样,多孔介质区域就像一个特殊的流体区域。
点击fluid
panel的porouszone选项就将这部分流体区域设定为了多孔介质区域。
这时界面就被展开如图7.19.1。
定义多孔介质的孔隙速度方程
在多孔介质模型进行模拟的时候,求解面板有多孔介质速度方程区域,可以选择指导fluent使用虚假速度或者物理速度来进行求解。
默认状态速度是虚假速度。
详细情况见section7.19.7.
定义流过多孔介质的流体
选择fluidpanel的下拉菜单materialname中选择流过多孔介质的流体。
如果想检查或者修改所选材料的属性,点击edit,;这个面板只有选择了的材料属
性,而不像materials面板里面的所有材料属性。
如果模拟的是组分扩散方程或者多相流模型,materialname菜单将不会
出现在fluid面板里。
在组分扩散方程计算中,多孔介质区域和流体区域的混合材料就是定义在speciesmodel面板里面的材料。
多相流模型中,材料因相的不同而不同,详细建section23.10.3.
激活多孔介质的化学反应
如果相模拟组分扩散的化学反应,可以通过激活fluid面板里面的reaction选
项来激活在多孔介质中的化学反应。
如果化学反应中包含壁面化学反应,那么就需要定义surface-to—volume
ratio值。
它是单位体积的表面积(A/V),可以看作催化剂载荷的一种度量方式。
有了这个值,fluent就可以通过将它和网格里面的体积相乘来得到网格内化学反应发生的总表面积。
定义化学反应机理的详细部分见Section14.1.4,
壁面化学反应部分见Section14.2
包含相对速度阻力公式
PriortoFLUENT6.3,caseswithmovingrefereneeframesusedtheabsolutevelocitiesinthesourcecalculationsforinertialandviscousresistanee.
Thisapproachhasbeenenhancedsothatrelativevelocitiesareusedfortheporoussourcecalculations(Section7.19.2).UsingtheRelative
VelocityResistaneeFormulationoption(turnedonbydefault)allowsyoutobetterpredictthesourcetermsforcasesinvolvingmovingmeshesormovingrefereneeframes(MRF).Thisoptionworkswellincaseswithnon-movingandmovingporousmedia.NotethatFLUENTwillusetheappropriatevelocities(relativeorabsolute),dependingonyourcasesetup.
定义粘性和惯性阻力系数
粘性和惯性阻力系数在同一个面板里面定义。
定义这些系数的基本方法是定义一个方向矢量(二维)和两个方向矢量(三维),然后定义每个方向上的粘性或者、
和惯性系数。
二维情况下,第二个方向没显式表示出来而是由定义的这个方向与z方向矢量确定的平面的垂直矢量。
三维时,第三个方向矢量是由已经定义好的两个矢量确定的平面的垂直方向。
三维问题中,第二个矢量必须与第一个矢量垂直,如果不垂直,求解器为了确保这两个矢量垂直,fluent会将第一个矢量相
关的第二个矢量的组成忽略。
所以你必须保证第一个方向矢量正确。
也可以使用udf来定义粘性和惯性阻力系数。
书写和加载了fluent以后udf选
项就可以使用了。
注意必须使用DEFINEDPROFILE宏来定义系数。
详细的udf资料间udf帮助文件。
当使用轴对称的旋转流动时,可以为粘性和惯性阻力定义一个附加方向组成。
这个方向始终与定义的另两个矢量相切。
基于密度和基于压力求解器都可以使用这个模型。
三维时,也可以使用圆锥坐标和圆柱坐标系统来定义阻力系数,如下。
□注意粘性和惯性阻力系数都基于虚假速度定义。
定义阻力系数步骤如下:
1.定义方向矢量
•在笛卡尔坐标系中简单定义1-矢量,三维时候还要定义2-矢量,没有定义的那个矢量由上面说明的方式定义。
这些方向矢量对应于多孔介质的原始坐标轴。
有些问题中多孔介质的原始坐标轴与计算区域的坐标轴不一致,这时就有可能不知道多孔介质的先前方向矢量。
这种情况下,三维平面工作和二维线工具就能帮助来确定这些方向矢量。
(a)"Snap''theplanetool(orthelinetool)ontotheboundary
oftheporousregion.(FollowtheinstructionsinSection27.6.1or27.5.1forinitializingthetooltoapositiononanexisting
surface.)在多孔介质区域的边界上面使用快照来使用平面或者线工具。
(b)旋转工具轴到合适的位置。
(c)一旦位置合适以后,点击updatefromplanetool或者updatefromlinetool选项,fluent就会自动设定红色箭头防线为1-矢量,绿色箭头方向为二矢量。
*Touseaconicalcoordinatesystem(e.g.,foranannular,conicalfilterelement),followthestepsbelow.Thisoptionisavailable
onlyin3Dcases.
(b)SpecifytheConeAxisVectorandPointonConeAxis.TheconeaxisisspecifiedasbeinginthedirectionoftheConeAxisVector(unitvector),andpassingthroughthePointonConeAxis.Theconeaxismayormaynotpassthroughtheoriginofthecoordinatesystem.
(c)SettheConeHalfAngle(theanglebetweenthecone'saxis
anditssurface,showninFigure7.19.2).Touseacylindrical
coordinatesystem,settheConeHalfAngleto0.
Forsomeproblemsinwhichtheaxisoftheconicalfilterelementisnotalignedwiththecoordinateaxesofthedomain,youmaynotknowapriorithedirectionvectoroftheconeaxisandcoordinatesofapointontheconeaxis.Insuchcases,theplanetoolcanhelpyoutodeterminetheconeaxisvectorandpointcoordinates.Onemethodisasfollows:
(a)Selectaboundaryzoneoftheconicalfilterelementthatis
normaltotheconeaxisvectorinthedrop-downlistnexttothe
SnaptoZonebutton.
(b)ClickontheSnaptoZonebutton.FLUENTwillautomatically
"snap''theplanetoolontotheboundary.ItwillalsosettheConeAxisVectorandthePointonConeAxis.(Notethatyouwillstill
havetosettheConeHalfAngleyourself.)
Analternatemethodisasfollows:
(a)"Snap''theplanetoolontotheboundaryoftheporousregion.
(FollowtheinstructionsinSection27.6.1forinitializingthe
tooltoapositiononanexistingsurface.)
(b)Rotateandtranslatetheaxesofthetoolappropriatelyuntil
theredarrowofthetoolispointinginthedirectionofthecone
axisvectorandtheoriginofthetoolisontheconeaxis.
(c)Oncetheaxesandoriginofthetoolarealigned,clickon
theUpdateFromPlaneToolbuttonintheFluidpanel.FLUENTvillautomaticallysettheConeAxisVectorandthePointonConeAxis.(NotethatyouwillstillhavetosettheConeHalfAngleyourself.)
2.在粘性阻力下面定义粘性每个方向的阻力系数。
在惯性阻力现面定义每个方向的惯性阻力。
对于高非均质多孔介质惯性阻力,激活alternativeformulation选项,这个选项为非均项多孔介质计算时候提供了一个更好的稳定性。
多孔介质的压降与速度的大小成正比。
使用方程7.19-6推到出如下表达式:
乙-沙莎诃泌⑺19-10)
是否使用alternativeformulation依赖于计算结果时候与实验结果吻合情况
如果流场与网格曲线平行,那么是否使用这个选项将没有区别。
Section7.19-8有详细的涉及到高非均相多孔介质模拟。
0注意alternativeformulation选项只能使用在压力求解器中。
IfyouareusingtheConicalspecificationmethod,Direction-1istheconeaxisdirection,Direction-2isthenormaltotheconesurface(radial(■)directionforacylinder),andDirection-3isthecircumferential()
direction.
三维有三种可能的系数分类,二维有两种:
•各向同性情况时,所有方向的阻力系数都相同,这时需要明确的设定所有阻力系数相同。
•三维时如果两个方向系数相同,第三个方向不同,或者二维情况时其它两个方向都不同,必须小心设定每个方向的阻力系数。
例如,如圆柱形杆在
流动方向有许多小孔,这样流动就会很轻松低通过圆杆,但是其它方向的
就会很少。
如果一个垂直与流动方向的平板,流动就会过不去,而往其它
方向流动。
・三维情况时,有可能三个方向的系数都不相同。
如多孔介质包含一系列不规则的空间物体,物体间的流动三个方向都不相同。
你就需要在各个方向定义不同的系数。
取得粘性系数和惯性损失系数的方法如下:
在已知压损情况下取得多孔介质速度基于虚假速度。
当使用多孔介质模型时,必须记住多孔介质必须百分百开放,粘性阻力和或惯性
阻力系数定义都必须基于这个假设。
以下算例就是如何来计算惯性阻力系数的算例。
假设一个多孔板的空口面积为损失系数如下:
25%,流场经过后的压损为动力水头的0.5倍,
S=;乜:
卜轨仁汐:
(7.19-11)
是0.5,基于孔板的实际流速。
例如开口面积为25%,计算一个合适的惯性阻力系数如下:
1.通过空办的速度假设孔口为100%。
2.压损系数必须被转化为单位多孔介质长度动力水头损失。
注意第一项,第一步就是计算调整以后的压力系数,是基于开口面积为100%情
况下:
切%0JW1=4X地聊弼怦
■11,-"r■,⑺19-12)
or,notingthatforthesameflowrate,
=吠X(了
=8(7.19-13)
调整以后的压力损失系数为&注意第二项,必须将其转化为单位厚度孔板的压
力损失系数。
假设孔板厚度为
1.0mm这时的惯性损失就如下:
Ga=
thickness
(7.19-14)
..
注意,对于非均质多孔介质,这些参数必须在其它几个方向也进行计算。
使用ergun方程来提取填充床多孔介质参数
第二个例子是填充床模型计算,在湍流流动中,填充床模型中包含渗透和惯性阻力。
其中能用来提取合适的常数的方法是ergun方程,半经验相关系数能够运用与很广范围的雷诺数和多种流化床:
|Ap|_i50M(l-e)a„T一■—吒十Dp
当流体流动为层流的时候,根据blake—kozeny方程以上方程的第二项就可以忽略:
方程中丿'是粘性系数,「是平均颗粒粒径,是床厚度,-是孔隙率,定义
为孔隙体积除以流化床区域体积.比较Equations7.19-4and7.19-6与
7.19-15,各个方向的渗透阻力系数和惯性阻力系数如下:
and
使用经验公式来提取孔板湍流多孔介质模型输入参数
第三个例子采用vanwinkle方程来计算矩形孔板压降的多孔介质模型输入。
文献中作者采用的上矩形孔等边三角形板表达式如下:
式中
孔的总面积
板的总面积
孔直径与板厚度的
=thefreeareaortotalareaoftheholes
=theareaoftheplate(solidandholes)
=acoefficientthathasbeentabulated
ranges不
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 多孔 介质 边界条件
![提示](https://static.bdocx.com/images/bang_tan.gif)