鱼类“钻石”状群游流动干扰及其减阻机理的数值研究

观察研究发现,在鱼类的群游中,鱼体可以有效利用相互之间有益的干扰,减少游动阻力,增强自己的推力和游动效率.本文基于作者以往建立的动态混合网格生成方法以及不可压缩流非定常计算方法,对四鱼呈"钻石"阵型排列的群游进行了数值模拟,通过不同横向间距、不同流向间距、不同相位差等游动状态下流场结构、力学特性的比较,分析了鱼体之间的横向干扰和流向干扰,对其降低阻力、减少功耗的流动机理进行了探讨.     

基于Delaunay背景网格插值和局部网格重构的变形体动态混合网格生成技术

建立了一种基于Delaunay背景网格插值方法和局部网格重构方法相结合的变形体动态混合网格生成方法.首先利用作者发展的定态混合网格生成技术生成变形体初始网格,即变形体附近利用结构化的四边形网格,外场采用自适应Cartesian网格,中间由三角形网格过渡.当物体运动或变形时, 首先利用Liu和Qin提出的基于Delaunay背景网格的插值方法对变形体附近的网格进行变形;当物体运动或变形位移很大导致局部网格质量急剧下降、甚至网格相交时,则在局部重新生成网格.利用该方法生成了多种变形体外形的动态混合网格,如单个鱼体、双鱼串列巡游,昆虫单翼、双翼扑动,多段翼型变形等,并在鱼体巡游的非定常数值模拟中得到了应用.     

二维不可压内流湍流附面层分离流的预测

本文采用适合于内流的准同时粘性/无粘性流匹配模型,快速求解二维通道内的不可压湍流附面层分离流动。无粘流采用由复变函数论的Plemeli公式推导的面元法;粘性流动区采用Head积分法求解。文中给出了两个算例。结果表明:所采用的准同时粘性/无粘性流匹配模型,不但可有效地克服分离形成的奇异性,而且收敛速度快,稳定性好。在湍流分离阶段,其预测结果与实验值基本吻合。文中还讨论了形状因子及表面摩擦系数关系式对预测结果的影响,形状因子经验关系对预测有比较明显的影响,而表面摩擦系数经验关系对预测无显著影响。文中还提出了改进模型的方向。     

双重互易边界元法在求解N—S方程中的应用

提出了用双互易边界元法求解Ｎ－Ｓ方程的新思路。用Ｐｅａｃｅｍａｎ－Ｒａｃｈｆｏｒｄ算子分理解将时间相依的Ｎ－Ｓ方程分解为线性和非线性子问题，线性问题用共轭梯度法解除压强－速度的耦合；非线性问题进行局部线化。对所得的Ｐａｓｓｉｏｎ方程及相关类型方程采用双重互易的边界元法求解，消除了传统Ｎ－Ｓ方程边界元解法的区域积分问题。     

不可压湍流平均流与拟序流的各向异性模式理论及其数值验证

使用侧偏统计平均方法理性推导了不可压湍流的封闭方程组。侧偏统计平均提供了湍流脉动量的一阶统计信息,在引入加权漂移速度对称性及正交各向异性后,导出了漂移流的连续方程、动量方程及能量方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组。方程组不含任何经验系数,保留了N- S方程的均化的非线性特性。其级数形式的能量方程与非线性现象多尺度层次现象相对应,具备了描述湍流统计平均流动及拟序结构流动的双重功能。      

Fourier伪谱方法在不可压缩平板边界层研究中的应用

以三维不可压缩平板边界层为研究对象,扰动形式N-S方程为控制方程,从空间模式的角度,直接数值模拟了三维不稳定T-S波传播的过程.时间离散采用三阶精度混合显隐分裂格式,空间离散则结合Fourier伪谱方法及高精度紧致有限差分逼近,法向采用非等间距网格坐标变换,出口边界条件采用嵌边函数法,程序采用MPI(Message passing interface)并行方法编写.实例验证,该方法计算结果与流动稳定性分析的结果一致.      

翼型低雷诺数层流分离泡数值研究

翼型在不可压和低雷诺数条件下工作时的性能越来越引起了人们的研究兴趣。本文采用Rogers发展的双时间步Roe格式,求解拟压缩性修正不可压N-S方程。数值模拟了低雷诺数条件下(Re=6.0×104,1.0×105,2.0×106),E387翼型在不同攻角条件下(α=0°、4°、7°)翼上表面后缘部分的层流分离现象,对应于长分离泡。时均化的结果同McGhee[12]的试验结果进行比较,验证了计算结果的正确性。研究了低雷诺数层流分离的非定常和时均化特性,并对漩涡脱落中主涡、二次涡以及出现的漩涡对并的周期性过程和对气动力脉动造成的影响进行了较为细致的研究。非定常计算结果表明,低雷诺数条件下的层流分离现象,是周期性的旋涡脱落过程。所谓长层流分离泡是其时均化积分的结果。     

毛细管束中层流粘性振荡流强化传热的若干问题

本文是作者从事振荡管流换热研究的总结。研究结果表明:流体在管内作轴向振荡运动时,可以极大地强化轴向的传热。本文给出了计及管壁效应的振荡管流换热的精确解析解,公布了变Pr数的实验结果。在此基础上,讨论了影响振荡管流轴向换热的6个无因次参数的物理含义,并指出其中的某些无因次参数在物理本质上就是傅里叶数的倒数的平方根,这将有利于理解振荡管流换热这种非定常换热过程。由本文的精确解析解,经过简化,可以得到国内外已公布的近似解。实验数据和解析解预言的基本规律符合很好。     

固体火箭发动机药柱不可压和近似不可压三维分析

基于不可压和近似不可压粘弹性有限元法,对固体火箭发动机药柱进行了三维分析。首先从Herrm ann泛函出发并结合对应原理,导出了适用于不可压和近似不可压粘弹性材料的本构关系,然后根据虚功原理,建立了三维问题的粘弹性增量有限元列式,最后对弹性约束的圆柱形中孔药柱进行了不可压和近似不可压分析。结果表明,该方法可用于泊松比接近甚至等于0.5的粘弹性问题计算,尤其适用于固体推进剂药柱的结构分析。     

第二代低阶面元法及在动力短舱上的应用

在气动力设计和计算中,面元法是主要实用计算方法之一。新近发展的第二代低阶面元法实用性强和计算精度较高。本文介绍了该面元法的理论,通过对比早期面元法,着重探讨了第二代低阶面元法的特点及其优越性。运用自行研制的计算程序,对低速动力短舱进行了计算和分析研究,结果表明与国外资料所提供的计算结果吻合得较好。