网格间距及正交性控制在翼—身—尾组合体三维绕流中的研究.docVIP

网格间距及正交性控制在翼—身—尾组合体三维绕流中的研究 7一2^ 19973CHINESEJOU 计 RNALO 算 FCOM 物 PUTION札PHYSICS997年月ATAL:1.. 网格间距及正交性控制在翼一身一尾 组合体三维绕流中的研究 薹流Euk翼/,r,关键词塑;宴揸生虞间距控制正交性控制E1ller方程—砭旗爿髦/相舟娄号璧7刍.州 0弓I言. 1二维方法 1.1Hilgenstock2Dl方法 采用的椭圆型方程为 {:笔,c,【‰Q(e,) 在计算域求解的等价方程为 口(苷+:r{)+y(仰+)一2却=0(2) 其中 = ,P/口,妒=,Qlr.口=f,f,=?ry=ff ,=a(z.)/a(},,?),,=xi+ 由于源项P,Q的值分别对等}线和等线起牵拉作用(作用相同),且正的P,Q值将导 致等}线和等线分别向},增大的方向移动,可以利用这一特点,不断改变边界上的P,Q 值,逐渐将网格线牵拉到所期望的位置. Hilgenstock2D1方法是将网格线与边界夹角的实际值与期望值之差作为修正量,在= 常数边界上修正的值.在e=常数边界修正的值(外层迭代),比如在=边界上的 收穑日期:1996.01—24;恬回日期:1996-09—16 计算物理第14卷 修正方式如下 f:一?tanh(0一)… 【0=arccos[r{?r/(IrlIrI)]… 其中0和分别是网格线与边界夹角的实际值和期望值.在=一边界,(3)式中的tanh 前取反号.=常数边界上源项值作类似的修正.边界上修正后的源项值被线性均匀地插入 内场,求解网格生成方程(2)至收敛(内层迭代).重复这样的过程直至生成的网格满足边界上 的正交性要求.修正公式中采用衰减园子(一个正的小量)和阻尼函数tanh是为了避免修 正量过大引起内层迭代的不稳定.用该方法可以生成与四边正交或成任意指定夹角的网格. 1.2Hilgenstock2D2方法 采用如下的椭圆型方程 f+=gP l7l=玑=gQ 在计算域的等价形式为 g姐(链+Pr})+gn(卯+Q)一2g121.卸)=0(5) 其中和分别为曲线坐标系的协变和逆变度量张量 g:g=/g.g..=gtl/g.g:一g2t/g=一gl2/g gl1=?r车,g22:r-r,g12=?r,g=[a(x,y)/O(,)] Hilgenstock2D2方法将网格线与边界夹角和间距的实际值与期望值之差作为修正量分别修正 两个源项的值,两个源项在同一族边界上进行修正(外层迭代).比如在=常数边界上同时 修正P,Q,以1=1…为例,则采用公式 P()=P(一?tanh(or一0)(6) Q)=Q+-tanh(d一d)(7) 其中0和的意义同前,d和d,分别是离开边界的第一层阿格点与边界间距的实际值和期望 值.在=一边界,tanh函数前应取反号.边界上修正后的源项值被线性均匀地插入内场, 求解阿格生成方程(5)至收敛(内层迭代).重复这样的过程,直至生成的网格满足边界上的正 交性及与边界间距的要求.a和tanh的意义同前.用该方法可以生成与一族边界正交或成 任意夹角,且与其间距为指定值的网格. 2三维Hilgenstock方法 依照二维Hilgenstock方法的思路,构造一种三维源项修正方法--t-Iilgenstock3D1方 法[.所用的椭圆型方程为 f++=P(,,) ‰‰=Q(,,)(8) 【++=R(,,}) 在计算域的等价形式为 口1(蒋+Pr{)+2(卯+c,)+3(*+f)+2(卢12却+fl23r+f131强):0(9) 第2期张正科等:网格间距及正交性控制在翼一身一尾组合体三维绕流中的研究249 其中=j2P/口1,Q-jQ/n2, n1 n2 n3 卢12=(ref)(,)一(?,)ll 23=(?)【?re)一(r)l【 31=(f)(re)一(re?re)ll 方程中的三个源项卯,P.,分别对,=常数网格面起牵拉作用,并且正的源项值将导 致相应的网格面向曲线坐标增大的方向移动.根据这样的作用特点,在r/,=常数边界面上 修正卯,在,=常数边界面上修正PQ,在,=常数边界面上修正,其修正量就取网格 线与边界面夹角的实际值和期望值之差.以:边界面为例,设,分别为该边界面上 任一网格点处的网格线和网格线与过该点穿出该面的网格线的夹角(正向切向矢量的 夹角),,分别为这两个夹角的期望值(希望网格线与边界面正交时取为/2),那么在该 边界面上PQ,似可采用如下的修正公式. fP窖一d?tanh(0;一口)…, lP=P一d?tanh(《一f) 式中=cos一[,?re/(Ir【II)],:c0sI1[r?re/(1}I1I)].在=一边界面,公式 (10)中的tanh函数前取反号.其余边界面上的源项值均作类似的