平面突扩管道中的湍流流动及预混燃烧计算

应用随机涡法(Random Vortex Method)对平面后台阶突扩管道中的湍流流动及预混燃烧过程进行了数值计算。与平均流场所显示的气流结构不同,瞬时流场表现出大尺度的多涡结构。由于这些大尺度旋涡之间的相互作用,回流区长度在平均值附近发生拟周期性变化。预混燃烧计算结果表明,大尺度旋涡对火焰面在流场中的分布有重要影响。     

管道强制电流阴极保护电位计算的等效电阻法

针对受一个阴极保护站保护的埋地钢质管道,提出了管道电位和电流计算的等效电阻方法,并用解析方法推导出管道覆盖层破损后电位和电流的计算公式。等效电阻方法将所研究管道划分成若干管段,用一个电阻表示管段对轴向电流的阻碍作用,一个电阻表示对径向电流的阻碍作用,从而将管道强制电流阴极保护系统表示成一个等效电阻电路。利用网孔电流法求解该电路的网孔电流,并计算出管道各节点上的电位。计算结果表明,等效电阻方法得到的结果与解析法的计算结果完全吻合。该方法特别适合于管径和（或）过渡电阻发生变化的情况。     

基于流体网络法的喷管冷却气流量分配计算方法

流体网络法对喷管气膜冷却问题具有独特的优势,但难以反映出入流对冷却气流道静压变化的影响。针对此引入了可反映出入流导致静压变化的管道元件,并基于孔口速度落后角对静压变化量进行了建模。验证算例表明,这一改进使沿程静压变化趋势与实际吻合较好,使冷却气流量分配计算误差降低到0.5%以下。将改进后的流体网络法与主流通道的计算流体力学模拟耦合用于分析喷管冷却气流量分配,与全计算流体力学模拟相比结果差异小于1.2%,且大幅降低计算量;在某型喷管设计中的应用表明,可以有效反映冷却结构改变对冷却气流量分配的影响。      

飞机环控系统管道气流温度的测量

根据实际使用中出现的问题，对各种实用的管道气流温度传感器的误差进行了计算，对电阻式和热电偶式温度传感器的误差原因作了定量分析。由于带密封罩的温度传感器在实际使用中存在着明显的缺点，建议今后使用热电偶测量飞机环控系统管道气流温度，并提出了一种实用热电偶温度传感器的结构形式，同时也介绍了如何减小由于热传导和热辐射引起的误差。     

推进剂供应管道的集中参数近似模型研究——模态近似模型

提出了用模态近似法建立分布参数湍流管道的集中参数近似模型的新方法,发展了适用于水力和绝热气体湍流管道的两种模态近似模型：Ｉ－Ｒ－Ｃ模态和等效层流模型。模型中考虑了频率相关摩擦损失。利用单根管道阀门关断问题的仿真计算对本文的模型进行了验证,计算结果与考虑频率相关摩擦损失的特征线方法的计算结果十分吻合。     

用数值模拟方法确定喷气Z箍缩的负载特性

本文运用工程方法设计了喷气Z箍缩装置的超声速喷管,通过数值求解Navier-Stokes方程模拟了气流从储气室到喷管出口的非定常流动过程,最后用DSMC方法计算了喷管出口气体向真空中的膨胀,得到了喷气负载的气流位形和质量线密度.计算结果表明,除了储气室的压力、喷管的喉道和出口尺寸等因素外,负载的特性参数较强地依赖于阀门开启和关闭时间、阀门开启的速度、开启的口径大小以及阀门后真空管道的几何形状.采用定常流动模型比非定常流动模型预测的负载质量线密度大10倍左右,而后者与实验结果较为一致.     

气流温度测量的动态响应问题

本文是用试验和计算的方法来解决气流温度的动态响应问题，在计算中避开了考虑气流温度敏感头的材料、形式，形状和质量；被修正的气流温度应该等于气流的真实温度以及先拟合后修正的方法，使求得的气流温度敏感元件的时间常数增加了可靠性。     

超声速流场的流线-特征线坐标变换与应用

为了求解二维平面弯曲激波波后流场,讨论并发展了一系列基于流线-特征线坐标系变换的流场代数计算方法。该系列方法依据二维平面流场中气流角、静压沿特征线近似不变的特点,可快速求解平面弯曲激波波后流场。其中,等气流角近似法适合模拟流线轨迹,等静压近似法适合求解波后流场参数。在此基础上,又提出了一种改进的等气流角-等静压混合方法用于计算弯曲激波波后流场。和特征线法对比,等气流角-等静压混合法计算得到的流场基本特征与特征线法得到的结果相同,在来流马赫数Ma=6和Ma=4情况下误差分别仅为0.5%和0.15%,证实了该方法在求解二维平面弯曲激波流场中的适用性。     

静载下含平面缺陷的压力管道的失效率

参照 SAPV-95规范,用第 2水平法对静载下含平面缺陷的压力管道进行了可靠性分析,得出了功能函数。提出了计算失效率的快速收敛的迭代算法,并分析了失效率与广义强度和广义应力的均值及变异系数的关系。得出了用失效率可对静载下含平面缺陷的压力管道的安全状况进行评估和预测的结论     

超声速气流中横向煤油射流的数值模拟

为了对超声速气流中液体横向射流进行数值模拟,采用Eulerian Lagrangian方法描述气液两相流动现象,用随机轨道模型追踪液滴的运动历程,并考虑气流可压缩性、流场不均匀性及液滴变形对液滴运动的影响,开发了相应的计算程序。计算结果发现,液体射流与气流之间存在着强烈的相互作用,液滴在进入气流中不久就破碎成很小的子液滴,受煤油液雾的影响气流速度、温度急剧下降,同时在喷孔上游出现一道弓形激波。与液雾结构的纹影图像和液雾穿透的实验结果对比显示,数值计算结果基本合理。     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围。在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术。通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式。研究表明：在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值。     

叶栅三维流场的实验研究

本文通过流动显示和详细测量,研究了叶栅中的三维流动。在平面叶栅风洞中,对非平直端壁的大弯度压气机叶栅所完成的测量,包括沿叶高的叶型表面上气流参数分布和叶栅通道中的整个流场。为了进行对比和详细研究端壁效应,还对二维叶栅完成了同样的测量内容。对所观察到的复杂流动现象进行了分析。实验结果还与叶栅全三维无粘流计算进行了对比。本文所获得的完整的三维流场实验数据,为建立叶栅流动模型,研究小展弦比叶栅中气流展向掺混的机理,以及为检验全三维流动计算方法提供了可靠的依据。      

在激波边界层干扰中后掠柱形区的独立性原理研究

在后掠压缩角引起的三维激波边界层干扰中,柱形的上游影响区表现出准二维的流动特性。实验结果表明,在柱形区中通常的边界层独立性原理不再适用。但是,如果考虑由后掠引起的横流效应修正,柱形区中的上游影响距离的尺度和相应的二维干扰上游影响区尺度之间可以得到很好的相关关系。     

三维可压缩非定常流的壁面分离判据及其分离线附近的流动形态

通过引入依赖于密度的物面法向速度变换wr=-1/ρh1h2∫ozh1h2(e)ρ/(e)tdz,描述物理速度空间(ua,va,wa=u,v,w+wr)具有无滑移壁面条件的三维可压缩非定常连续方程可转换成变换速度空间(u,v,w)内具有无滑移条件定常连续方程.因此,采用定常壁面分离的分析方法和结论,再通过变换和研究wr的贡献,给出了三维可压缩非定常壁面分离的判则以及分离线附近的流动形态.研究指出,二维和三维情况下,都出现伴有壁外附着的壁面分离情况.数值模拟证实了理论和结论.     

二维超音速激波风洞及其应用

简述了二维超音速激波风洞的设计要点和性能,并给出在该风洞中低凸台诱导激波和湍流边界层相互作用的实验结果。实践表明：在被驱动段和喷管之间安装一个前缘光滑的矩形剖面短管道;并将二维喷管精加工成型,就能获得实验所需的均匀超音速气流。     

高精度数值方法在DNS中的应用

利用不同格式模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算结果的影响。研究了不同计算格式对自由剪切层扰动波增长的影响,结果表明:对于二维扰动波模型,采用二阶精度格式,扰动波在平板边界层中只能出现衰减的结果;在自由剪切层中,虽然能模拟出扰动波在传播过程中幅值的增长,但却在相当长的计算范围内看不见三维波的出现。但对于同样的扰动模型,如果采用高精度格式,在边界层中能出现增长的T-S波,并与线性理论非常吻合;在自由剪切层中,扰动波不仅迅速增长,并能在不很长的距离内,明显的看到三维波的出现,并激发出了三维结构。      

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

基于DES的二维和三维空腔流动特性研究

空腔流动广泛存在于航空航天工程中,对其流动特性的研究具有十分重要的工程意义。采用数值方法对比模拟了二维和三维空腔流动,控制方程采用N—S方程,空间离散采用有限体积方法,对流通量计算采用Roe格式,非定常时间离散采用双时间步长方法,湍流粘性计算采用基于SA模型的DES方法。数值计算所得的三维空腔底面压强分布与实验结果一致,所得的二维和三维空腔内的流动结构与相关文献中的分析相吻合。将二维和三维的计算结果进行比较,发现空腔底部压强分布、空腔内的压强振荡、空腔内部流线图及声压级等均有所不同,得出在空腔宽度有限时,横向流动简化了流动结构．削弱了振荡幅值．但使振荡过程复杂化．其三维效应不能忽略．     

湍流脉动作用理论在轴对称流动计算中的应用

本文在湍流压力脉动远距离作用理论及其平面问题计算基础上,将其应用于轴对称流动的计算,并与轴对称湍流射流的实验结果进行了比较,表明轴对称流动计算比用平面计算更与实验一致。