利用fluent对全机的数值模拟lsh.docVIP

利用fluent对全机的数值模拟 绪言 很高兴能站在这儿给大家说说我最近的工作。从今年7月份开始，我利用fluent作为工具，开始着手歼－8教练机的全机模拟。下面我将谈到工作中遇到的一些问题和我的解决方法，肯请大家指正。 网格的工具，机器配置 gambit， 双奔3◎866mhz，2G内存 数据处理 原始数据文件是文本文件，给出了所有节点的笛卡儿座标，考虑用gambit里的journel功能进行批处理生成节点，简单的小程序可以实现文本文件到journel脚本的转换。 总的节点数目约250万。 做网格的基本思路 网格生成的基本顺序是点联线，线围面，面围体，然后切分网格，其中每个过程都有很多方法可供选择，在这次工作中用的最多的方法分别是用nurbs方法联线，用net surface(U-V)方法围面，用stitch方法围体 现在进行到具体的gambit处理部分，我简单说说在工作中遇到的问题和处理方法： 尾翼和机体交线过于复杂 交线过于复杂给围面带来了困难，同时也给将来的网格生成带来潜在的隐患，所以要进行合理的简化，重新作出一条简单的交线。飞机尾部和尾翼分别略去这条edge后围面后求交依然得不到理想的结果，原因是尾部形状在这儿有个弯度，两个面直接求交并不理想。我的处理方法是引进了第三个面。过这两点连一直线，水平往外平移得到一个矩形，和尾翼面相交得到这条线，认为它就是交线。 机翼与机体分离 从这个图可以看到，从所得数据绘出得机翼并没有插入得机身里面，必须人为得把机翼往里延伸才能得到封闭的体。处理方法是把这条edge拷贝一份，并适当放大后往里平移，然后uv围面，现在这两个面是相交的了，可以顺利找到两个面的交线，切完后机身是这个面，而机翼则去掉这条临时拷贝出的edge，换成交线重新围出。 围机翼时候U-V方法失败 在围机翼的面的时候有几个小问题，第一是在这个地方有个比较大的拐角，为了保持几何特征是分成两个面来围了，还有个问题是在翅膀靠里的这个面围成的时候，uv方法总是报错，后来去掉里边的某一个edge就通过了，由于这几个edge变化比较平稳，去掉其中一个几乎没有什么影响。 飞机的几何翻转，注意点，并线时候遇到问题及解决办法 前面所有的面都是半个飞机的，而所给的数据也只是半个飞机的，利用对称的特性可以得到另一半。最简单的做法是在网格做完后完成翻转，然而这样并不顺利，因为网格连续要求必须所有的线都是唯一的，而翻转网格的做法使得在对称面上的线都有两份，它们的几何位置完全一样，但因为上级有面，体的存在而难以合并。所以翻转的工作适合在这个时候完成，也就是在所有的面刚做完之后来完成。面翻转之后，把对称面上的所有的点，线合并。 计算区域 至此飞机的面全都围好了，考虑加上计算区域，计算区域应该有足够大以淡化边界效应，在这个例子中计算区域是一个矩形，长宽高分别为翼展的xxx，飞机在矩形的中心部 虚面的出现及影响 前面飞机表面的处理都是用实面来完成的，进行到围体阶段，不可避免的出现了虚面，有些操作，如将飞机的某些面切分成多个面，用实面操作很难完成，因为实面操作能接受的误差非常小，而用虚面却能轻易的做到，另外有些功能像edge的合并只能用虚的操作来实现，而虚的edge所在面必须也是虚的，虚的面所在体也必须是虚的，如此一来一个虚的edge就会导致多个体都是虚的，所以在这儿看这幅图大部分几何都是这种蓝色，它代表的是虚体。 先将所有体围好后做网格的坏处 第一次做网格的时候我是将所有的体都围完后才开始做网格的，结果是白费了许多时间，教训是在没有很好的把握下，最好还是做完一个体就画好一个体的网格，因为会有这样的问题产生，即有的体很难划分结构化网格，需要对几何进行改动，而一个体的改动会牵动相邻的体，进而波及更多的体，因此我觉得一个体一个体做网格比较稳妥，可避免浪费工作。 先从机翼处围起的好处 网格是从飞机的主翼根部开始围起的，好处是显而易见的，这儿的网格最为复杂，最容易产生问题，这儿处理好之后接着往两头推就比较容易了，而如果从两头往里做的话，一旦到这儿遇到问题难以避开的话，势必要返工前面的网格。 交错网格的问题 这是我在这次工作种遇到的最为头疼的问题，当这个体做完网格后检查时候，报网格inverted的错，网格方案的调整始终没有帮助 交错网格的解决办法 后来的处理方法是比较悲惨的，我改变了映射方法，从原来的像侧面投影改成向底面投影，代价是必须把所有的体都删了从来，一次大的返工 网格的完成 之后的工作比较简单了，依次做体做网格，这张图显示的就是整个的网格图，有将近400万的网格，试验室的机器并起来后内存依然不够，所以必须要简化。 网格的简化 简化是件轻松的工作，看这张描绘体的图，飞机所在的这几个