三维C-H网格下跨音压气机转子Euler流场数值模拟

介绍了三维C-H网格的生成和全三维转子流场的数值解。以法国宇航院设计的一个高负荷跨音转子为算例, 并与他们的计算结果作了比较。应用这套数值模拟技术得以计算带进口预旋的某跨音压气机转子流场, 并对进口预旋对其下游的转子流场的影响作了数值模拟, 得出了若干有意义的结论。      

基于伴随方程的翼型数值优化设计方法研究

对Ｅｕｌｅｒ方程为流动模型的翼型绕流,在微分方程生成的网格中,以Ｂ样条控制顶点为设计变量,用连续伴随方程法对翼型气动外形进行数值优化设计。给出了两个亚声速的反设计算例,验证了该方法的可行性。     

跨声速三维Euler方程的显式／隐式格式有限元解法

本文采用有限体积Ｇａｌｅｒｋｉｎ法和非结构网格数值求解了跨声速三维Ｅｕｌｅｒ方程,通量的求解采用ＶａｎＬｅｅｒ－阶逆风矢能量发裂方式,并通过外插使其上升为二阶精度,分别采用了显式及隐式两种格式求解。     

求解线性矛盾方程组的递推Householder变换法

 求解线性矛盾方程组是工程实践中经常遇到的问题。目前主要采用法方程组解法和基于各种正交变换的三角化解法。由于法方程组系数矩阵的病态条件数是原矛盾方程组系数矩阵病态条件数的平方,所以计算精度较差,有时甚至得到完全错误的结果。因此一些计算精度要求高的问题和高阶方程组的求解,多采用三角化法。但是,三角化法只能进行一次性处理,尚未出现随数据更新的实时递推算法。许多需要在线处理和数据存储最大的问题,只好采用一般递推最小二乘法。     

Navier—Stokes（NS）方程组差分计算中的物理和网格尺度效应及NS方程组的简化

对于高Re数流动计算，在通常二阶精度NS差分格式和网格数条件下，存在某些粘性项落入修正微分方程截断误差项的问题。这类NS方程组计算实际是计算某种简化NS方程组，而且重复计算误差物理粘性项既浪费机时和内存，误差积累又会对数值解产生不可预测的影响，避免外述缺陷的办法一个提高NS差分格式的精度，另一个是丢掉可能落入截断误差项的物理粘性项，把NS方程组简化为广义NS方程组，广义NS计算避免了误差物理粘性项误差积累对数值解的不可知影响，又可节省内存和机时，对高Re数流体工程计算很有好处。利用广义NS方程组计算超声速绕前向和后向台阶流动的结果表明：广义NS方程组与NS方程组的数值结果很好相符符。     

三维机翼气动结构多学科优化方法

采用Euler方程、连续伴随优化方法、自由形面变形技术（free-form deformation,FFD）以及一种四边形线性壳单元模型对三维机翼气动结构多学科优化方法开展了研究。目标函数为阻力系数和机翼翼盒结构质量的加权和,设计变量为FFD控制点以及10段蒙皮的厚度。气动方面升阻力系数对设计变量的梯度由伴随方法求得,结构方面翼盒质量对于蒙皮厚度的梯度可直接计算,应力的梯度采用有限差分方法获得。对跨声速大展弦比机翼进行气动结构多学科优化,并将该计算结果与机翼气动单学科优化设计后的静气弹计算结果进行对比,结果表明:该气动结构多学科优化方法能够在满足约束条件下,实现阻力系数和翼盒质量同时降低,保证优化结果的合理性。      

CFD技术在目标电磁特性计算中的应用

电磁学中麦克斯韦方程组和流体力学中无粘流动欧拉方程一样,都是具有实特征值的双曲型偏微分方程组,相同的数学特性使得计算流体力学(CFD)技术能够在计算电磁学(CEM)中得以应用。采用MUSCL格式结合Steger-Warming分裂计算电磁通量,采用4阶Runge-Kutta法计算非定常时间推进,借鉴CFD方法计算电磁场边界条件;采用时域有限体积法(FVTD)数值模拟宽带脉冲波、完全导电体、含吸波材料涂层介质/金属混合体以及复杂外形目标的电磁散射。结果表明:基于CFD的FVTD方法能够高精度地计算目标的电磁特性。     

钟形返回舱空气动力特性数值模拟

文章基于Euler方程及N－S方程的数值求解方法，对回舱亚、跨声速和高超声速的流场及气动特性进行了数值模拟，其中Euler方程数值求解采用二阶Godunov有限体积法；N－S方程数值求解采用二阶Harten-Yee格式的差分法，得到同实验值一致的物面压力、气动力系数和在不同速度范围出现的激波、流动分离及旋涡等流场特征。通过完全气体、平衡气体和非平衡气体的流场数值模拟结果分析比较，得出真实气体效应对返回舱气动力特性影响较小这一结论。计算结果表明数值模拟方法是预测返回舱气动特性的有效手段。     

石油钻井井身轨道设计的优化计算方法

在井眼轨道设计中．假设井眼轨道上的每一段都是空间平面圆弧或直线段，相互之间在连接处相切。根据这一假设，建立了描述井眼参数之间相互关系的非线性方程组，并采用数值优化理论求解非线性方程组的方法来进行井身轨道参数的数值计算。