一种关于平面网格生成的通用数值技术的研究.pdfVIP

CSAA98．PC．005 一种关于平面网格生成的通用 数值技术研究 虢天总公司三院三十一所 施发树刘兴洲 f摘要】本文建立了一种应用于田四条不规则边界线封闭的平面域通用网格生成数值方 法，应用于某型小涡扇发动机环形燃烧童的CFD计算。实践证明：该方法较简单，计算 快速、准确，在匠{复杂形状边界组成的喷气发动机燃烧童流场计算中具有重要应用价值。 主题词：CFD，网格生成，燃烧 术语 乱坐标x沿X方向的增景或修正量 山 ．坐标Y沿Y方向的增萋或修正量 下标： XX 沿x坐标方向或i方向的二纷中心差分 YY 沿Y坐标方向或i方同的二价中心差分 ———————————j}——～ —夸—鳃为拜秘扛瑚肆瑚一 V 沿y坐标方向的一价中心差分 1．引言 网格生成技术是CFD研究的～个重要领域。如果物理区域边界规则且简单，就可以 用不太复杂的数学方法直接生成正交网格，照对于复杂边界，传统的方法就很难适用了。 现代涡喷发动机燃烧皇结构十分复杂。它主要田多个串联的火焰筒或～体兜环形燃 烧室组成。在其内部，燃油充分燃烧形成高温燃气流：在其外部，气流对高温壁面进行冷 却，最后进入内部与高温燃气流掺混形成合适的出口温度场。在对燃烧室进行介质数值模 拟时，要区分不同的物理区域。而它们都是由燃烧室内、外机匣和相邻的边界封闭的，一 收地，边界很不规则。这些边界可以分为内、外‘缘’和进、出‘口’。所谓内、外是依 边界离发动机轴线是近还是远．而进、出口则是依气流是进还是出该区域。这样，闯题就 转化为如何生成由四条不规ⅢU边界封闭域的网格。 28 2．网格生成方法 采用柱坐标系，分别沿与内外缘和进出口方向相对应的径同和轴向，利用线睦插值方 法初步划分出网格。这肘每个网格元是一四边形元，每个网格点周围分布有八个点。由于 边界形状极不规则，育不少’瘸态单元’，需要对网格进一步修正，建立与边界相适应的 接近于正交的网格元，量网硌线比较光滑。如何建立修正准则和数值模型是关键。 2．1初始网格 缓设燃烧宝中某～物理区域如图1所示。 domain X Fig．1physical 圈中CED、AFB分别为内、外边界，CA、DB分别为进、出边界，而E、F分别是 内、外缘上导数不连续点。首先，分别对径向、轴向边界，利用点分布函数方法＆l口求 搭密度：其次，秀别求径向、动向相对应边界点坐标的平均值L冠瓢)，同时计算各向平 均坐昧点构成的折线的累积长度与总长的比铡r／， f霓：05。(朋f，1)+豇f、M1)；√(否一萨(i一国3 辅向：i墨二。5×c r(i．t，十H，．M，，’2{：了i毫；；i焉‘滓2’3，…N～i： ∑托F丽盯乒可 虢嚣袅鼢；揣吨= j=2．3，…M．～1 ∑托F丽Aj4开*≯j-葡w ， 然后，先沿轴向利用r值线性插值计算出”中间”坐标，即 {琢f√：)：!。xf．V√)+(1—0×．rd√) xH吖√)+(1一t)×取i√) I琢i√)=t 最后沿径