超声速化学反应流动的LU-SGS伪时间迭代空间推进求解

介绍了求解超声速多组分有限速率化学反应流动的伪时间迭代(lower-upper symmetric-Gauss-Seidel,LU-SGS)方法.空间推进求解多组分抛物化(parabolized Navier-Stokes,PNS)方程时,在一个推进面上采用修正的LU-SGS方法迭代至收敛,把得到的结果作为初值赋给下一个推进面.沿推进面依次迭代直至求解完整个流场.采用伪时间迭代LU-SGS方法求解化学反应PNS方程,计算结果的准确性和时间迭代求解完全Navier-Stokes(N-S)方程相当,求解效率提高一个数量级.      

配置点谱方法求解非均匀介质内辐射传热问题

发动机内的燃气等高温介质随着组分和浓度的变化，会引起折射率在空间上的非均匀分布，从而导致辐射能束沿着曲线传播，其相应的辐射传热过程也更为复杂。为了避免射线追踪方法的复杂计算、提高计算效率，本文提出了配置点谱方法求解二维非均匀介质内辐射传热问题。在求解过程中，角向采用离散坐标法处理，空间采用配置点谱方法处理。通过将三种非均匀介质内辐射传热问题的配置点谱方法结果与文献结果进行对比分析，发现配置点谱方法可以在较少的节点数下，获得准确、有效地计算结果。并且，采用配置点谱方法求解三种算例的计算时间均消耗较少，均在20分钟以内。这将为进一步开展发动机复杂结构内高温燃气辐射快速仿真提供基础。     

气动热数值模拟中的网格相关性及收敛

针对气动热数值模拟中的网格相关性和收敛问题,采用计算流体力学(CFD)计算方法,以典型的钝头为算例,进行了数值模拟研究。研究结果表明:网格是气动热数值模拟中的关键因素,壁面附近法向网格间距最为敏感,热流结果随网格间距不同会出现数倍乃至数量级的差异,在得出正确的壁面压力结果的网格上不能保证得到正确的热流结果。分析认为壁面附近客观存在的巨大温度梯度是网格要求敏感的原因。计算表明:气动热计算的收敛过程比压力和流场收敛慢得多;在压力和流场完全收敛时,热流与最终的迭代收敛值偏差达20%以上;以方程残值的下降量级和压力收敛作为热流收敛的判别标准是不合适的。在气动热计算中,应直接观察热流数据的收敛,确保得到真正的收敛解。     

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响      

机翼电加热防冰表面内外传热的耦合计算

分析了机翼防冰表面的传热传质现象,对电加热防冰表面内外传热耦合计算进行了研究.采用附面层积分法对机翼表面传热系数进行计算.根据传热传质比拟得到了防冰表面蒸发量.建立表面微元的能量与质量守恒差分方程组,利用高斯-赛德尔迭代法求解,将外部防冰载荷、蒙皮弦向导热和内部加热热流三者耦合,采用亚松弛系数稳定迭代过程,最后得到平衡表面温度、溢流水和溢流结冰的范围.计算结果和文献中试验以及其他数值仿真结果比较,吻合较好,表明了模型及计算方法的正确性.      

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

TBCC燃烧室/喷管一体化壁面温度计算

基于Navier Stokes(N-S)方程组对包括隔热屏、隔热屏内外流、大气外流在内的涡轮基组合动力(TBCC)发动机燃烧室/喷管进行了一体化的气/热耦合数值模拟,考虑了燃气组分输运、辐射换热等影响,研究了其在某典型飞行状态下TBCC冲压发动机燃烧室/喷管筒体及隔热屏内外壁壁面温度、辐射换热热流及对流换热热流分布.结果表明:燃烧室/喷管筒体与对称面上下交线的壁面温度在轴向距离为0.5~2.6m内变化较小,在轴向距离为2.6~3.1m内急剧增加,在轴向距离为3.1~3.5m内急剧下降.之后,上交线筒体壁面温度沿流向减小,下交线筒体壁面温度先升高后降低.筒体壁面温度最高点在喷管下调节板收缩段,为1577K.隔热屏内壁面辐射热流在370~500kW/m2变化,上下交线处的辐射热流较外壁面的辐射热流约高300kW/m2,辐射热流沿流向先减小后增加.隔热屏外壁面辐射热流在50~200kW/m2范围内分布.      

高超声速压缩拐角湍流计算的收敛性分析

应用基于shear stress transport(SST)湍流模型的隐式紧耦合方法对高超声速压缩拐角流动进行了数值计算.考察了差分格式、限制器以及湍流时间步长等7个因素对湍流计算过程稳定性和收敛性的影响.结果表明:差分格式和限制器对收敛性有很大影响,总体而言耗散性越大的差分格式和限制器收敛性越好;增加湍流时间步长对收敛性影响很小;适当降低湍流变量最低限制值有利于湍流方程的收敛;增大Courant-Friedrichs-Lewy(CFL)数和隐式内迭代次数可显著加速收敛,但均存在最优值,超过最优值则加速效果不明显.      

某型飞机短舱辐射换热计算

针对某飞机地面慢车工况,进行了短舱内复合换热数值计算研究,分别运用球谐函数法、离散坐标法和蒙特卡洛法计算辐射换热.结果显示:不同计算方法所得结果总体上趋于一致,热辐射与自然对流耦合作用明显.与单纯自然对流的计算结果相比,舱内温度显著升高,各区域温度值上升7%~30%.考虑空气介质对不同波段热辐射的选择性吸收后,短舱内各计算点温度进一步上升.综合分析表明:地面条件下飞机短舱换热分析过程中热辐射的影响不可忽略.      

飞机整体瞬态热状况的数值仿真研究

 在分析飞机内外热环境的耦合传热机制基础上，建立了描述飞机整体瞬态热状况的物理数学模型。引入壁面热流函数，将飞机蒙皮外的气动对流与辐射换热作用转化为舱内热分析的浮动热边界条件，实现了飞机内外耦合热作用的解耦计算，避免了难以实现的直接耦合求解。采用区域分解技术，将飞机整体分解为特性不同的多个求解区域。采用蒙特卡罗法求解舱内设备之间的辐射换热，采用热网络法求解设备内部及各设备之间的辐射-导热-对流耦合换热。利用该方法对某型飞机飞行过程中的瞬态热状况进行了数值模拟，获得了飞机舱内的瞬态温度场，分析了相关因素的影响，为飞机设计及相关研究提供了重要依据。     

双排叉排孔气膜冷却特性研究 (二) 绝热有效温比

本文介绍双排叉排孔气膜冷却绝热有效温比η的测量结果和数值计算。根据η沿主流和横向的变化, 特别是在近孔区的分布规律分析了前后排孔二次射流的不同流动特点。数值计算采用求解椭圆方程的三维差分程序, 吹风比M等于0.6。      

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

求解非线性代数方程组的一种建议方法

飞行动力学研究中常遇到的求解非线性代数方程组的问题，现提出一种求解实函数非线性代数方程组的组合迭代方法，以多步交替使用的方式将梯度下降法和拟牛顿代法结合起来，综合两者的优点，组成一种对大范围内偏离精确解的任意初值均能收敛，且有一定收敛速度的迭代法，通过算例，对三种方法进行了对比和分析，计算结果证明，该方法是优越的。     

燃气透平叶型气动与传热的优化设计

建立了一个透平叶型气动与传热优化的两步模型.第1步采用解析多项式曲线构造叶片的参数化造型，基于CFD数值模拟，通过正交设计和敏感性分析，优选几何设计参数.第2步采用贝塞尔(Bezier)曲线参数化叶型型线，以总压损失和表面热流加权函数为优化目标，调整叶型型线，获得了气动性能与热负荷较优的叶型.设计案例的结果表明:增加安装角可以降低叶片压力面热负荷；增大叶片前缘直径可以降低前缘滞止点热负荷.调整叶型的表面曲率分布，有可能推迟吸力面的转捩，弱化叶片表面的传热.两步优化法计算量小，收敛性好，具有工程应用价值.      

横流对冲击射流换热特性的影响

用实验和数值计算的方法研究了横流对单孔短距冲击换热特性的影响规律.结果表明:①在实验研究范围内,冲击靶面平均努塞尔数Nu_av,is随横流的增大而增大;②冲击靶面前区努塞尔数Nu_fia随横流的增大而减小,这是因为壁面附近x负向流动的气流流速绝对值变小,换热削弱;③靶面冲击点附近及冲击靶面后区努塞尔数Nu_bia随横流的增大而增大,原因是在冲击点上游出现的马蹄涡的旋涡方向和冲击射流自身形成的卷吸涡拉伸涡对方向相同,这两种涡的叠加使该区域冲击靶面的换热得到增强.      

网络连接机群上CFD计算的一种负载平衡方法

讨论在网络连接的分布式存储并行机群上进行CFD计算的一种负载平衡方法和相关问题。当数值计算的方法是区域分解法时,各节点机上的负载与其处理的子区域的网格节点数成正比,所以同样大小的子区域可保证负载平衡。采用负载再分配即区域再剖分的方法来达到动态负载平衡,构造了一种对多块网格(每块内为结构网格但各块间的连接关系可以是非结构的)进行近似平均分划的方法,并用一种数据结构和算法保证了剖分出的多个子区之间正确的互联和边界信息传递关系。给出了用这种剖分算法对常见多块网格进行分划的例子,并对三维流场进行了实际计算。结果证明本文方法是正确的,能有效地均衡计算负载,达到较理想的并行效率。     

几何非线性机翼本征梁元素模型的高效化改进

 采用Hodges等提出的时间-空间离散化的几何精确非线性本征梁通用模型处理柔性机翼结构动力学问题时,当离散化的节点数增大时,该方法的未知数数量成倍地增长,而且方程组是严重病态的,因此数值模拟计算的速度非常缓慢。针对机翼中最常见的悬臂梁结构,根据空间离散化的边界条件,提出了空间缩聚法把空间离散差分方程缩聚为常系数矩阵格式,得到了只与时间相关的微分方程组,进一步推导得到了该方程组的雅可比矩阵,因而大大减少了方程组的数量以及求解过程的循环和迭代步数。采用Gear方法分别求解了原始的本征梁元素模型和本文提出的缩聚模型,结果表明空间缩聚模型在相同条件下可提高运算速度约5.1倍,而且对不同类型的外载荷都具有较好的通用性、稳定性和高效性。     

高超声速再入飞船气体热辐射计算边界虚网格方法

针对再入飞行器激波层辐射加热,应用有限体积法计算辐射传输方程,引入了一种边界处理方式——虚网格方法,根据化学非平衡流场数据,计算出飞行器表面辐射加热热流。通过分析发现:采用虚网格方法处理辐射传输边界,能够确保真正的边界条件得到应用,满足真实的物理意义,且辐射传输计算的空间网格、立体角网格、初始温度值的收敛性很好,计算结果精度误差小于2.5%。利用虚网格方法计算飞船再入返回时的辐射加热,结果表明计算所得壁面辐射热流与其他研究者辐射热流计算结果符合良好。      

解大型结构特征值问题的同时迭代的一个新算法

本文提出了同时迭代法的一种变型算法,用来求解大型结构动力分析中的广义特征值问题。使用了E_k子空间和特征方向的概念,证明了该算法的收敛率,并讨论了如何在计算机上实现的某些细节问题。通过统计每步迭代所需运算量和大量的实际计算,表明该算法与当前几种常用的同时迭代法或子空间迭代法相比,收敛率相同,但减少了每步迭代所需的计算量,提高了计算效率。     

管内高温介质流动入口段的辐射与对流耦合换热数值模拟

将求解辐射传递方程的离散坐标法与求解对流换热的控制容积法相结合 ,数值模拟了高温下圆管流动入口段参与性介质的辐射与对流耦合换热。考察了介质光学厚度及管壁温度对温度分布、壁面热流密度及局部努谢尔数分布的影响。结果表明 ,高温下介质辐射使介质的温度分布和对流换热状况发生较大改变 ,对管内流动入口段的换热起着重要作用