生物医学信号检测与处理实验2——复杂管道流动的三维模拟

复杂管道流动的三维模拟的几何模型的建模思路

首先建立一个通道截面的二维模型；然后通过“拉伸”操作形成一个三维模型。其中，通道截面的二维模型可以采用多个矩形进行布尔运算得到。

建模步骤

本模型采用流体流动模块，采用导航流程【新建】$\to$【模型向导】$\to$【三维】$\to$【流体流动】$\to$【单相流】$\to$【层流】$\to$【稳态求解器】$\to$【完成】，进入模型开发界面，开始进行建模。保存此工程文件为ComplexFlow3D.mph。进入的初始模型开发界面如下图。

第一步：设置参数与变量

参数在【模型开发器】$\to$【ComplexFlow3D.mph(root)】$\to$【全局定义】$\to$【参数1】中进行定义。参数与变量如下表所示

名称	表达式	值	描述
PIn	10[Pa]	10 Pa	入口压强
Pout	0[Pa]	0 Pa	出口压强
Rhos	1000[kg/m^3]	1000 kg/m3	流体密度
Mus	0.001[Pa*s]	0.001 Pa.s	流体动力粘度

其中“名称”、“表达式”和“描述”手动输入，“值”会根据“表达式”的输入自动显示。在输入参数数值（即“表达式”）时，用中括号[ ]来标识物理单位，物理单位之间通过数学运算符连接。参数设置时，还可以借助参数设置栏下面的功能图标，采用导入导出的方式批量操作。

第二步：构建几何模型

采用7个矩形进行布尔运算获得模型的纵截面。

首先要建立一个可供拉伸的二维模型，所以要先建立一个工作平面。建立流程为【模型开发器】【ComplexFlow mph(root)】→【组件1】→【几何1】，在【参数设置区】选择【长度单位】为um。鼠标右键点击【几何1】，在弹出菜单选择【工作平面】即可建立一个工作平面对象【工作平面1】。在【参数设置区】选择【平面】为【xy平面】，【z坐标】默认为0um。左键单击【工作平面1】下的【平面几何】，【模型可视化区】即进入几何设计模式。

建立7个矩形：右键点击【平面几何】，在弹出菜单中选择【矩形】，这样的操作重复7次。然后按照下表中的数据，在【参数设置区】设置矩形的尺寸和位置，其中xw表示矩形左下角x方向坐标，yw表示矩形左下角y方向坐标。设置完成后点击【参数设置区】顶端的【全部构建】。

矩形编号	宽	高	xw	yw
1	400	150	0	250
2	400	400	400	0
3	400	150	800	0
4	125	125	450	225
5	125	125	625	225
6	125	125	625	50
7	125	125	450	50

为了建立一个整体的二维通道模型，还需要对7个矩形进行布尔运算，分为下面两步进行。

(1)建立差集

点击“模型树”中【几何1】→【工作平面1】下的【平面几何】。点击【主菜单功能选择区】$\to$【工作平面】$\to$【布尔操作和分割】，点击其下拉菜单选择【差集】，即可在【平面几何】下建立对象【差集1】。
对【差集1】进行设置，在【参数设置区】的【差集】设置中，点击【要添加的对象】旁边的【激活】，使之切换为“on” 状态，之后选择r2，(在【模型 可视化区】中点击矩形r2）即可。再点击【要减去的对象】旁边的【激活】，之后选择r4、r5、r6、r7。再点击【参数设置区】顶端的【构建选定对象】即可。
在【模型可视化辅助菜单区】点击【缩放到窗口大小】，可得如下图结果。

(2)建立并集

点击“模型树”中【几何1】→【工作平面1】下的【平面几何】。点击【主菜单功能选择区】$\to$【工作平面】$\to$【布尔操作和分割】，，点击其下拉菜单选择【并集】，即可在【平面几何】下建立对象[并集1]。
对【并集1】进行设置，在【参数设置区】的【并集】设置中，点击【输入对象】旁边的【激活】，使之切换为“on” 状态，选择r1，dif1和r3，之后清除【保留内部边界】前面的复选框。再点击【参数设置区】顶端的【构建选定对象】即可。

右键点击【工作平面1】，点击【主菜单功能选择区】$\to$【拉伸】。在【参数设置区】，在【距离】→【指定:到面的距离】，设置为150，再点击【参数设置区】顶端的【构建选定对象】即可获得下图所示几何模型。

第三步：设置物理场属性

本例中，流质的设置不采用材料库中的【Water】，直接采用“空材料”定义动力粘度和密度。定义方法为:右键点击【模型开发器】→【ComplexF low. mph(root)】→【组件1】→【材料】，在弹出菜单中选择【空材料】。在【参数设置区】的【几何实体选择】中选择所有域；设置【材料属性明细】中的【密度】(变量rho) 为“1000” ，设置【粘度【 (变量mu) 为“0. 001”。

第四步：设置边界条件

【模型开发器】→【ComplexF low. mph(root)】→【组件1】→【层流】，右键点击【层流】，分别点击【入口】、【出口】添加边界条件，建立的对象为【入口1】和【出口1】。

点击【入口1】，在【参数设置区】可对入口边界进行设置。在【模型可视化区】在【边界选择】窗口中点击“x=0”处端面(编号为1)，添加入口边界。在【边界条件】的下拉菜单中选择【压力】，并在下面的【压力条件】中输入Pin。类似地，点击【出口1】，然后在【模型可视化区】，点击x=1200处的边界端面，该端面的序号为28。出口边界默认为压力边界条件，数值处输入“Pout"不用再进行设置。

第五步：网格划分

点击“模型树”中的【网格1】，在【参数设置区】进行网格大小设置。如果按照上例的方式，采用自动的网格设置，会出现一些问题。比如采用小划分的网格，在边界和拐角处的网格会比较密，而在通道内部的网格又会比较粗显然，默认的网格划分模式不能较好的满足需求。在这种情况下，需要在【参数设置区】的【网格设置】→【序列类型】中选择【用户控制网格】，选定后发现在“模型树”的【网格1】下面出现了新的对象，主要对其中的【大小】进行设置。

点击【大小】，进入【参数设置区】，在【单元大小】中选择【定制】，将【最大单元大小】设置为15 (通道宽度的10%)，【最小单元大小】设置为5。其他设置不变。最后点击【参数设置区】顶端的【全部构建】即可完成网格划分。

第六步：配置求解器

此例不涉及多场耦合，属于求解层流流动的基本情况，求解器的选择及相关参数采用默认值即可。

第七步：求解

求解方法: 【模型开发器】→【ComplexFlow.mph(ot)】→【研究1】，右键单击其下选项【步骤1:稳态】，在弹出的菜单中选择【计算选定步骤】，即可开始计算。

怕了怕了，算了好久还是这样。

第八步：后处理

本例中，后处理的任务有两个，第一是显示出通道纵截面的速度分布;第二是找到通道中速度最大值出现的位置和数值大小。

任务1

【结果】→ 【速度】→【切面1】，单击【切面1】进入【参数设置区】。在【平面数据】中设置【平面类型】为【快速】，设置【平面】为【xy平面】，设置【定义方法】为【平面数】，设置【平面数】为1，表示只选择-一个 显示平面。设置完毕后点击【参数设置区】上端的【绘制】，即可得到在“z=75” 处(通道中央)的速度场。

任务2

点击【结果】→【速度spf】，再到主菜单功能区找到【更多绘图】，在下拉菜单中找到【体图】，选择【体最大值/最小值】，即可在【结果】→【速度spf】下建立【体最大值/最小值1】的对象。点击该对象进入【参数设置区】，在【高级】做如下设置: 【单元细化】设为4 (进一步细化内部网格，提高最大值的定位精度) ；【显示精度】设为3 (保留3位有效数字)；【显示】选择【最大值】；勾选【包含单位】。在【着色和样式】中，【颜色】选择【白色】；【背景色】选择【黑色】。最后点击【参数设置区】上端的【绘制】，即可在原速度场上增加关于最大值的显示。