二维三维流固交界面网格生成新方法.docxVIP

二维三维流固交界面网格生成新方法 流场和流固耦合研究包括网络切割问题。网格生成与更新是对连续空间区域的剖分,将连续空间分为有限数量的子区域。流固耦合问题中,流固耦合交界面是两种不同介质间的边界,不同于固体力学或流体力学中的边界问题。流体与固体的相互作用主要体现在移动交界面处;当交界面发生运动变化,即产生流固耦合现象。流体域和流固交界面处附近区域的网格生成与更新质量直接影响流固耦合计算的精度和收敛性,交界面处理及其数据传递决定流固耦合问题计算的成败。 一般,网格生成可分为结构化网格方法和非结构化网格方法两类[2―3]。这两种方法各有优缺点,均广为应用:在结构化网格方法中,每一节点与其邻点间的连接关系固定不变,无需设置数据确认节点与相邻之间的联系,且网格数量容易控制[4―5];在非结构化网格方法中,单元与节点的编号无规定规则。 本文针对流场和流固耦合问题中的初始网格生成和网格更新,基于形式代数理论,运用位移映射协调流固交界面处的网格,提出虚拟结构化映射的新型方法。在虚拟规则多边形上建立虚拟网格,利用函数建立一一对应的映射关系,把虚拟规则多边形投影到不规则流固交界面上,把虚拟网格映射到真实的流固边界网格上。该方法简洁、精确,适用范围较宽。 1 基于理论和方法的基础 1.1 维虚拟网格的构建 形式代数是利用逻辑关系表达空间关系的逻辑数学,是空间关系和逻辑关系的转换工具。形式代数中的形式集合由离散的形元组成,而空间结构化网格同样由节点按照一定规律组成,两者相似。因此,形式代数适用于结构化网格的生成。 形元是网格图形的基本单元。结构化网格的一条边可看作一个逻辑形元。该形元由两个节点、一条边组成。所有逻辑运算均可建立在最小的形元运算基础之上。形式函数是对形元进行运算的工具。函数的形式有平移、镜像、放大等。通过逻辑运算,形元可按一定的逻辑规律排列。 构建平面网格形式前,应先建立基本形元,并通过形式函数运算,得到形式。图1为网格生成的示意图;其中,图1(a)示5×5二维虚拟网格的编号,可看出形元生成过程。其中:0号形元为水平基本形元,通过平移函数作用,先进行水平方向平移,后竖直方向平移,直至生成所有水平形元;30号形元为竖向基本形元,由相同的过程生成所有竖向形元。上述过程可表达为: 其中:F代表平面网格形式;F横杆代表横杆形元组合;F竖杆代表竖杆形元组合;#代表和运算;col和row分别代表网格总列数和总行数;basicnum(x,y)为基础形元的编号和坐标。 平移函数tranid(j,i,x,y)可写为: 其中:basicnum(x)、basicnum(y)分别为基础形元的x轴、y轴坐标;rowlength、collength分别为网格横向和纵向总长度。 1.2 两维网格线性映射 在生成虚拟网格基础上,对真实的计算域可运用映射算法进行网格生成[7―9]。兹以两维网格为例推导网格的映射函数。若真实的计算域边界为直线边界,可采用线性映射方法,映射关系写为: 若计算域边界为曲线边界,则可采用非线性映射方法,如二次非线性映射的映射关系为: 式中:ui代表真实计算域中边界上参考节点i的坐标;xi和yi分别代表虚拟网格中对应于节点i的虚拟坐标;n代表真实计算域网格边界上参考点的数目,参考点数目可依边界复杂程度选取。若边界为直线,则可仅取4个参考点。以两维网格线性映射为例推演映射列式,将式(3)扩展为: 其中,i、j、k、l为真实网格边界上的四个参考点。令: 得到: 其中: 定义节点i的转换函数为: 则映射函数写成: 式中:u代表真实计算域中任意节点的坐标,亦可代表其他标量;ui、uj、uk和ul分别代表真实计算域网格边界上i、j、k、l等4个参考点处与u相同的标量。 把虚拟网格边界参考点的坐标i(-1,-1)、j(1,-1)、k(1,1)和l(-1,-1)代入式(6),由式(7)可得: 同理,经推导并整理,两维网格非线性映射的节点转换函数为: 1.3 流体固体交界网格的网格更新 计算域中流体域和固体域交界面处,两种介质的网格节点共存并有部分共点节点。一般,对流体网格可划分得更细密些,即耦合交界面上,固体域的节点间常设有若干流体域的节点,如图2所示。首先对固体域进行网格剖分且基于此剖分,在交界面处把流场域划分为若干子域,并在每个子域中运用映射关系建立网格。这样,在流固交界面上的流体域网格和固体域网格的对应关系清晰,两种介质的计算数据交换简便。若每个子区域均足够小,可直接利用线性映射关系。 当固体域的节点发生移动,需及时调整流体域的网格分布,以避免流固交界处两种介质网格的脱离或纠缠。传统处理方法之一是将流体网格看成拟弹性介质,对流体域网格进行更新[10―11]。该方法不破坏流固交界面处的流体固体节点拓扑关系,但不能保证移动后的网格不发生扭曲,且调整后的网格