下载文档原格式
4.2顶盖驱动流4.2.1物理模型在一个正方形的二维空腔中充满等密度的空气，方腔每边长为0.1m，其顶板以0.1m/s 的速度向右移动，同时带动方腔内流体的流动，流场内的流体为层流。
计算区域示意图如图4-2-1所示。
图4-2-1 计算区域示意图4.2.2在Gambit中建立模型Step1：启动Gambit并选择求解器为Fluent5/6。
Step2：创建面操作：→→打开对话框如图4-2-2所示。
输入长度和宽度10，在Direction中选择XY Centered。
图4-2-2 创建面设置对话框Step3：划分面网格操作：→→打开对话框如图4-2-3所示，Shift+鼠标左键选择正方形面，Internal size=0.5，其它保留默认，点击Apply确认。
划分后的网格如图4-2-4所示。
图4-2-3 网格划分设置对话框图4-2-4 计算区域网格图Step4：设置边界类型操作：→●在Name栏输入边界名称wall-1，将Type栏选为Wall，在Entity栏选取Edges，并选中方腔顶部边线。
●在Name栏输入边界名称wall-2，将Type栏选为Wall，在Entity栏选取Edges，并选中方腔其它三条边线。
Step5：输出网格文件操作：Fil m→export→mesh打开对话框如图4-2-5所示，选中Export 2-D mesh 前面的复选框，输出网格文件。
图4-2-5 网格文件输出对话框4.2.3求解计算Step1：启动Fluent选择2d单精度求解器，点击Run，如图4-2-6所示。
图4-2-6 启动求解器图4-2-7 网格尺寸设置对话框Step2：导入并检查网格1．读入网格文件操作：Fil e→Read→Case...找到文件后，单击OK按键确认。
2．检查网格操作：Grid→Check网格读入后，一定要进行网格检查，注意最小体积不能为负值。
3．网格比例设置操作：Grid→Scale...在Gambit中，生成网格使用的单位是cm，在Grid Was Created In下拉菜单中，选取cm，如图4-2-7所示，然后单击Scale，关闭对话框。
水泥3D打印喷头内浆体流动的MRT-LBM分析吴伟伟;黄筱调;方成刚;李媛媛【摘要】水泥3D打印时,水泥浆体在喷头内的流动性对挤出成型有重要影响.以水泥、水、聚羧酸减水剂、改性棒土和纤维素醚作为原始打印材料,根据流变特性测试,采用Herschel-Bulkley模型作为水泥浆体的流变模型,并利用多松弛时间的格子Boltzmann方法(MRT-LBM)对浆体在喷头中的流动进行了分析,引入矩函数来描述松弛过程中的碰撞步,并利用泊肃叶流的理论解对MRT模型进行了验证.在进行实际仿真时,以雷诺数作为准则进行物理单位和格子单位之间的转换,获得相关的流线图和速度场分布图.结果表明:水泥浆体在螺槽截面的流动呈现环流形状,环流的中心在(0. 5W,0. 7h)位置处;根据速度场的分布,在螺槽横截面的左下角和右下角不存在流体流动,适当增大螺杆速度或螺槽宽度有助于水泥浆体在螺槽内的有效输送.【期刊名称】《南京工业大学学报（自然科学版）》【年(卷),期】2018(040)005【总页数】6页(P79-84)【关键词】水泥3D打印;Herschel-Bulkley流体;多松弛时间的LBM;数值模拟;流线;速度分布【作者】吴伟伟;黄筱调;方成刚;李媛媛【作者单位】南京工业大学机械与动力工程学院,江苏省工业装备数字制造及控制技术重点实验室,江苏南京 211800;南京工业大学机械与动力工程学院,江苏省工业装备数字制造及控制技术重点实验室,江苏南京 211800;南京工业大学机械与动力工程学院,江苏省工业装备数字制造及控制技术重点实验室,江苏南京 211800;南京工业大学机械与动力工程学院,江苏省工业装备数字制造及控制技术重点实验室,江苏南京 211800【正文语种】中文【中图分类】TH12对浆体在喷头内部的流动进行分析,可以很好揭示其流动规律,流线图可以描述浆体在腔体中的主要流动区域和环流的形成,速度场分布图可以描述各速度分量沿不同方向的分布情况,速度较大的区域更有利于浆体的输送,速度较小的区域往往受黏性系数和机械结构的影响,因此合理的机械结构有助于浆体在喷头内部的流动。
一、问题描述方腔顶盖驱动流动如图1所示的一个简化两维方腔（高，宽都等于L)，内部充满水分。
上表面为移动墙，非维化速度为ｕ／ｕ0 =1。
其他三面为固定墙。
试求方腔内水分流动状态。
u=1, v=0u=0, v=0 u=0,v=0u=0, v=0图1常微分方程理论只能求解极少一类常微分方程；实际中给定的问题不一定是解析表达式，而是函数表，无法用解析解法.二、离散格式数值解法:求解所有的常微分方程 计算解函数 y(x) 在一系列节点a = x 0< x 1<…<x n = b 处的近似值),...,1()(n i x y y i i =≈节点间距为步长，通常采用等距节点，即取 hi = h (常数)。
步进式：根据已知的或已求出的节点上的函数值计算当前节点上的函数值，一步一步向前推进。
因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。
欧拉方法1(,) 0,1,...n n n n y y h f x y n +=+=几何意义在假设 y n = y (x n )，即第 n 步计算是精确的前提下，考虑公式或方法本身带来的误差: R n = y (x n +1) y n +1 , 称为局部截断误差.截断误差: 实际上，y (x n ) ? y n ， y n 也有误差，它对y n +1的误差也有影响，见下图。
但这里不考虑此误差的影响，仅考虑方法或公式本身带来的误差，因此称为方法误差或截断误差。
局部截断误差的分析：由于假设y n = y (x n ) ，即y n 准确，因此分析局部截断误差时将y (x n +1) 和 y n +1都用点x n 上的信息来表示，工具：Taylor 展开。
显式欧拉公式一阶向前差商近似一阶导数推导如下：223111232()[()()()()][ (,)] ()()h n n n n n n n n n h n R y x y y x hy x y x O h y hf x y y x O h +++'''=-=+++-+''=+1()()()n n n y x y x y x h+-'≈111()()() ()()(,)n n n n nn n n n n y x y x hy x y x y y x y y h f x y +++'≈+↑≈≈=+隐式欧拉公式xn +1点向后差商近似导数 推导如下：几何意义设已知曲线上一点 P n (x n , y n ),过该点作弦线，斜率为(x n +1 , y n +1 ) 点的方向场f (x ,y )方向,若步长h 充分小，可用弦线和垂线x =x n +1的交点近似曲线与垂线的交点。
流函数- 涡量法的二维方腔流数值模拟基本方程：22v tξγξψξ∂=∇-∇∂-∇=⎧⎨⎩t ———时间步长, s;γ———流体的运动粘度, m2 /s;ξ———涡量, s- 1;ψ———流函数, m2 /s;v ———速度矢量, m/s;差分格式：采用FTCS 格式，对于ξ有：1,,1,1,,1,11,,1,,1,,1,2222()22()()n n n nn nn n nn n ni j i ji j i ji j i j i j i j i ji j i j i j nni ji juvtxyx y ξξξξξξξξξξξξυ++-+-+-+-+----+-+=--++∆∆∆∆∆采用FTCS 格式，对于ξ有：1,,1,,1,,1,,1,2222)()()n nn n n n n n i j i ji j i j i j i j i j i j i j tx y ψψψψψψψψξ++-+---+-+=++∆∆∆边界条件：速度的边界条件：固体壁面处设置无滑移边界条件, 即u=0, v=0; 流函数的边界条件：在固体壁面及平板驱动处的流函数ξ=0; 涡量的边界条件： 平板驱动处：2uxξ=∆ 左右壁面上:;,1,,22()()i j i j i j y ξξξ+-=∆网格划分：采用等距结构化网格划分（40*40）编程计算：本算例采用MATLAB进行编译，其主要优势是语言简单，可以方便地描绘出方腔环流的等值线图等。
主要语句：while norm(c1)>1e-4|norm(c2)>1e-4n=n+1;O1=O;E1=E;t=t+dt;for i=2:Ifor j=2:JO(i,j)=O(i,j)-dt*(u(i,j)*(O(i+1,j)-O(i-1,j))/(2*dx)+v(i,j)*(O(i,j+1)-O(i,j-1))/ (2*dy))+dt/re*((O(i+1,j)-2*O(i,j)+O(i-1,j))/dx^2+(O(i,j+1)-2*O(i,j)+O(i,j-1))/d y^2);endendfor i=2:Ifor j=2:JE(i,j)=E(i,j)+0.25*(E(i+1,j)+E(i-1,j)+E(i,j+1)+E(i,j-1)-4*E(i,j)+dx^2*O(i,j)); endend% 0为涡量，E流函数%计算效果图：分别设置雷诺数为200，500,1000雷诺数Re=200时的流函数图及速度矢量图雷诺数Re=500时的流函数图及速度矢量图雷诺数Re=1000时的流函数图及速度矢量图。
二维流体力学流函数方程解巴巴二维流体力学的流函数方程是用来描述二维流体流动的一种数学方法。
在流体力学中，流函数是一个与速度场相关的标量场，是通过将速度场的分量与流体的密度相乘得到的。
流函数方程可以通过求解一个偏微分方程得到，其解可以给出流体在平面上任意点的速度和流线的分布。
∇²Ψ=0其中，∇²表示拉普拉斯算子，Ψ表示流函数。
这个方程是一个二阶偏微分方程，可以通过数值或解析方法求解。
为了求解流函数方程，需要满足一些边界条件。
常见的边界条件有：1.无物体边界条件：在流体的边界上，流函数Ψ的值可以设定为常数。
这个常数可以根据物理问题的条件来确定。
2.固壁边界条件：如果存在凸起的固体物体，可以设定在固体物体表面上流函数Ψ的梯度为零。
这意味着流体在固体物体表面上流动时，流线是与固体表面垂直的。
3.对称边界条件：当存在对称性时，可以设定在对称线上流函数的梯度为零。
这意味着在对称线上流体的速度分量与对称线垂直。
对于简单的几何形状，可以使用解析方法求解流函数方程。
例如，对于圆柱体的流动，可以使用柱坐标系下的方程来求解。
而对于复杂的几何形状或非线性问题，通常需要使用数值方法来求解流函数方程。
常用的数值方法有有限差分法、有限元法和有限体积法等。
一旦求解了流函数方程，就可以通过流函数求解速度场和流线的分布。
流函数的等值线是与速度场的切线垂直的，在流线上速度矢量的切线方向与流函数梯度的方向相同。
因此，通过流函数可以方便地描绘流体流动的图像。
总之，二维流体力学的流函数方程是求解二维流体流动的一种数学方法。
通过数值或解析方法求解流函数方程，可以得到流体在平面上的速度和流线的分布。
这对于研究和理解流体流动的特性和行为具有重要意义。
MAC算法计算二维方腔顶盖流动李江飞;石兆东;段兴华;李岩芳;张康;逯国强;陈颖超;任亚东【摘要】二维方腔流动是不可压缩黏性的典型流动,可以用来检验各种数值算法计算精度和可靠性,目前尚不能求得它的解析解.基于Matlab编程,采用交错网格MAC算法求解二维方腔流动,计算采用控制容积积分法离散控制方程,对流项和扩散项采用中心差分格式,得到流动达到稳定状态时各物理量的分布.【期刊名称】《宜宾学院学报》【年(卷),期】2015(015)006【总页数】4页(P28-31)【关键词】数值模拟;方腔流动;控制容积积分法;MAC算法;离散【作者】李江飞;石兆东;段兴华;李岩芳;张康;逯国强;陈颖超;任亚东【作者单位】承德石油高等专科学校热能工程系,河北承德 067000;承德石油高等专科学校热能工程系,河北承德 067000;承德石油高等专科学校热能工程系,河北承德 067000;承德石油高等专科学校热能工程系,河北承德 067000;承德石油高等专科学校热能工程系,河北承德 067000;承德石油高等专科学校热能工程系,河北承德067000;承德石油高等专科学校热能工程系,河北承德 067000;承德石油高等专科学校热能工程系,河北承德 067000【正文语种】中文【中图分类】TB126Li JF,ShiZD,Duan XH,etal.Calculation of Two-dimensionalCavity Flow Based on MAC[J].Journal of Yibin Univer⁃sity,2015,15（6）：28-31.二维不可压缩黏性流体方腔流动顶盖拖动速度为utop，方腔的长度和高度均为H,流体密度为ρ、动力粘度为μ.边界条件：流动速度u、v采用无滑移边界条件，利用动量方程推导压力p的边界条件[1].流动与传热的控制方程如下：其中，p为压力，u、v分别为x、y方向速度分量.用高度H、流体密度ρ和拖动速度utop作为无量纲标尺，将控制方程无量纲化，流场初始状态为静止，Re=1000求流动达到稳定状态时，x方向中垂线（x=H/2）上的无量纲速度U，y 方向中垂线（y=H/2）上的无量纲速度V，绘制出速度分布曲线，并求出中垂线上||U、||V的平均值.1.1 涡量控制方程无量纲化以高度H、流体密度ρ和速度utop作为无量纲标尺，将控制方程无量纲化[2]：将上述无量纲量代入题中流动与传热的控制方程，得出如下的无量纲方程：1.2 边界条件边界条件为：流动速度采用无滑移边界条件，壁面处法向速度恒为0，切向速度也为零.顶盖u=1，v=0；其余u=v=0.采用均分网格，网格数80×80的交错网格来离散方程，将压力和速度放在不同位置，压力放在网格中心，以Pi,j为主节点，背离P点的U、V与P点有相同的编号，速度分量U与P在X方向位置相错半个网格，速度分量V与P在Y方向位置相错半个网格，具体如下[3]：P：X方向：0—81，左边点0，右边点81，边点与内点距离为，其余ΔX；Y方向：0—81，下边点0，上边点81，边点与内点距离，其余ΔY；U：X方向：0—80，左边点0，右边点80，相邻两点距离ΔX；Y方向：0—81，下边点0，上边点81，边点与内点距离，其余ΔY；V：X方向：0—81，左边点0，右边点81，边点与内点距离，其余ΔX；Y方向：0—80，下边点0，上边点80，相邻两点距离ΔY.对于MAC算法而言，采用交错网格，用控制容积积分法离散控制方程，对流项和扩散项采用中心差分格式.时间步长为Δτ，空间步长为ΔX、ΔY.对速度分量U进行离散[4-5]：内点处理：非稳态项：对流项：扩散项：压力项：边界点处理：对于上边界点，扩散项：对于下边界点，扩散项：内点离散后的动量方程为：对速度分量V进行离散：内点处理：非稳态项：对流项：扩散项：压力项：对于左边界点，其扩散项：对于右边界点，其扩散项：离散后的动量方程为：将，代入离散的连续性方程，如下：整理化简可得压力离散方程：其中：aP=aE+aW+aN+aS,aE=aW= aN=与上边界相邻的内点：与下边界相邻的内点：与左边界相邻的内点：与右边界相邻的内点：求解步骤如下[6-8]：①确定网格信息，如空间步长、时间步长：ΔX，ΔY，Δτ；②定义变量，给速度场和压力场赋初始值和边界值；③经过（1）、（2），可得完整的速度场离散结果，据公式求；④根据（3）求解压力泊松方程，采用Gauss-Seidel迭代求解，循环直至满足收敛条件；⑤用该时层满足收敛条件最新的压力场去更新速度场，得到下一时层的，；⑥用下一时层的，返回（4），直到稳态的解，求出速度场和压力场.程序流程如图1所示.U、V、P放置在相同的网格位置时，取中心差分，切断相邻点间关系，计算结果会出现非物理意义的振荡；通过交错网格，引入相邻两点的压差，在物理上保证了相邻点压力的相互影响，避免了振荡的压力解.稳定状态时，计算结果如图2、图3所示.【相关文献】[1] Peng Y F,Shiau Y H,Hwang R R.Transtion in a 2-D lid-driven cavity flow[J].Computer＆Fluids,2003,32（3）：337-352.[2] 陶文铨.数值传热学[M].第二版.西安：西安交通大学出版社，2010.[3] Abdallah S.Numerical solutions for the pressure poisson equation with neumann boundary conditions using a non-staggeredgrid[J].Journalofcomputationalphysics,1987,70（1）：182-192.[4] Hortmann M,PerićM,ScheuererG.Finite volumemultigrid predic⁃tion of laminar natural convection：Bench-mark solutions[J].Inter⁃national Journal for Numerical Methods in Fluids,1990,11（2）：189-207.[5] DemirdžićI,PerićM.Finite volumemethod for prediction of fluid flow in arbitrarily shaped domainswithmoving boundaries[J].Inter⁃national Journal for Numerical Methods in Fluids,1990,10（7）：771-790.[6] Brandt A.Multi-level adaptive technique（MLAT）for fast numeri⁃cal solution to boundary value problems[C].Proceedings of the Third International Conference on NumericalMethods in Fluid Me⁃chanics,Springer Berlin/Heidelberg,1973：82-89.[7] Wang J,Li JF,ChengW X,et parison of finite difference and finite volume method for numerical simulation of the incom⁃pressible viscous driven cavityflow[J].Advanced Materials Re⁃search,2013（732-733）：413-416.[8] Li JF,Long J,Yuan L A,parison of finite difference and finite volumemethod for numerical simulation of driven cavity flow based on MAC[C].Computational and Information Sciences（IC⁃CIS）,2013 Fifth International Conference on,Shiyang,2013：891-894.。
