塞式喷管三维流场的数值模拟

从曲线坐标下的三维平均雷诺的N-S方程出发,并用k-ε两方程湍流模型封闭方程组,采用二阶精度的精度无波动、无自由参数的耗散差分格式(NND格式),发展了模拟塞式喷管三维流场的数值程序,并对线性塞式喷管和瓦状塞式喷管在不同落压比下的三维流场进行了数值模拟研究.结果较好地预示了塞式喷管的三维流场特征,并且显示瓦状塞式喷管除了圆柱形结构比平板更能承载外,其塞锥压强分布可以进一步改善其受力情况.     

附属小圆柱对主圆柱绕流影响的数值模拟

从小马赫数下的可压缩Navier-Stokes方程出发,采用覆盖网格分区算法,对在主圆柱尾流场的适当位置放置附属小圆柱时的流动特性及其对主圆柱绕流问题的影响进行了详细的数值模拟研究.根据模拟结果,主要探讨在小雷诺数范围内,附属小圆柱位置、个数、大小以及雷诺数变化等因素对主圆柱绕流的流场结构及其非定常演化过程的影响.计算结果表明,在小雷诺数范围内,无论附属小圆柱的个数如何,都存在最佳抑制区域.当附属小圆柱放在这个区域时,涡街被完全抑制,整个流场达到准定常状态;同时,阻力系数也显著下降.在一定程度上,随着小圆柱直径的增加,主圆柱涡脱落更易被抑制.另外,相对于单个附加小圆柱而言,使用两个附加小圆柱能使主圆柱涡脱落的抑制提高到更高的雷诺数.      

绕三角翼流动中的非定常现象研究

利用流动显示技术、片光技术及PIV技术,对三角翼面上的各种非定常特性做了研究.分别对涡层中的小涡结构,二次涡同主涡相互作用所引起的二次涡的喷射效应,螺旋破裂扰动的非定常特性,涡破裂点沿轴方向振动的非定常特性,以及完全分离流中的非定常特征进行了较为细致的研究,给出了不同扰动的频率特性.      

脉冲射流强化圆柱喷流混合的三维数值模拟

对利用脉冲射流强化圆柱喷流混合的流场进行了数值模拟,这项工作将高频高强度信号射流射入亚音剪切层以强化喷流混合.利用有限体积法和重整化群(RNG)k-ε湍流模型求解N-S方程.分别进行了低亚音马赫数及高亚音马赫数喷流直接排入空气的数值模拟.主要计算3个状态:未受激励的状态,斯德鲁哈尔数St分别为0.2和0.4时的状态.数值模拟结果显示,脉冲激励引起喷流在与激励源平行的平面迅速扩展,在与之垂直的平面收缩,当脉冲射流激励的St在0.2附近时,相关的混合较好,而且少量的射流流量可以极大地增加速度的衰减率.      

接触角法测玄武岩及玻璃纤维表面能实验

通过Wilhelmy吊片法测试液体与玄武岩、玻璃纤维单丝的前进接触角,应用OWRK理论计算得到纤维表面能及其色散与极性分量.实验研究表明:未经表面处理的玄武岩纤维表面能高于S-2玻璃纤维,去除浸润剂后2种纤维表面能的极性分量均明显降低,并且测试液体的极性分量对纤维表面能的测试结果影响甚大.分析结果为纤维与树脂的表面张力合理匹配、形成良好界面浸润和粘接提供了理论依据.      

基于非Favre过滤的大涡模拟

将直接过滤的可压缩Navier-Stokes(NS)方程组作为可压缩湍流大涡模拟方程组.方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,得到可压缩张量扩散模型.与传统的采用质量加权过滤(Favre过滤)的可压缩大涡模拟方程组相比,由于这种方法没有密度信息的损失,所以能反映真实的密度脉动信息;方程所求量为流场物理量的过滤值,物理意义明确,使大涡模拟结果与直接数值模拟或实验结果更具有可比性.用这种大涡模拟方程组对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验、Baldwin-Lomax(BL)模型及Smagorinsky亚格子模型的结果进行了对比.      

涡街流量计信号处理方法与实现

设计一种基于变增益运算放大器和带通开关电容滤波器的新型涡街流量计信号处理电路.单片机利用模拟数字转换电路(A/D)实时检测涡街流量计传感器输出信号的幅度和频率,同时利用数字模拟转换电路(D/A)对变增益运放进行控制,为后续的信号检测电路提供幅度相对稳定的有效信号.带通开关电容滤波器在单片机的控制下实时跟踪变增益运放输出的信号频率,保证带通滤波器的中心频率接近信号频率,以达到良好的滤波效果.滤波后信号通过单片机进行频率信息的精确计算以及流量的解算.实验结果表明,变增益运算放大电路有效解决了涡街流量计传感器输出信号幅度变化范围大而造成的放大电路复杂,分段放大信号幅度不连续等问题;带通滤波器通过单片机实时调整中心频率,对信号中各种干扰起到较好滤波效果.      

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

航空发动机一维数值仿真

分析了航空发动机一维气动热力数值仿真模型,此模型能够用于航空发动机整机数值仿真.进行了整机性能数值仿真,求解过程考虑了粘性效应以及气体流动损失对于整机性能的影响;采用了具有TVD性质的高阶精度Godunov格式,此格式具有良好的数值稳定性和激波捕捉能力;应用了时间推进法,交替使用显/隐格式.计算了某型带加力的混排涡扇发动机的速度特性以及节流特性并且与设计数据进行了比较分析,数值仿真结果与设计数据比较吻合.分析了典型工作点的计算过程,模型收敛速度较快.