关于三维Navier—Stokes方程的粘性项计算

借助于张量分析和张量计算，在贴体曲线坐标系下本文讨论了不同求解变量导致了粘性项个数上的重大差异和不同大小的计算量，并提出了便于粘性计算的最佳形式。文中借助于有限体积离散技术，通过引进两个对称辅助矩阵[A]和[B]，使粘性项的计算量大大减少，这对完成三维粘性流的数值计算具有重要的指导意义。借助于上述方法，本文完成了某型真实进气道两种工况的三维粘性Navier-Stokes方程计算（即M∞=3.0,a=0°和设计工况M∞=2.65,a=0°），获得了满意的全场结果；对于M∞=2.65的设计工况。同实验数据作了比较，符合良好。由于本文的方法明显的减少了粘性项的计算量且节省了大量内存，以致于使三维流场的N-S求解能在普通微机上进行。     

基于模拟退火算法与响应面模型的三维气动优化设计方法

提出了一套适用于高耗时三维气动设计问题的优化设计体系.其主要思想是采用改进拉丁超立方体试验设计选取样本点,应用自行开发的三维粘性流场求解程序进行流场计算建立数据库,采用二次响应面方法建立近似模型,再应用高效模拟退火算法进行全局寻优.以NASA rotor57为对象,在详细进行流场计算基础上采用所提出的优化体系对其进行了三维积叠优化设计.在对流量、效率加以严格约束的条件下,总压比可提高1.8%,对流场结构进行了分析.优化结果表明本优化方法省时,适于三维气动设计的特点.     

用粘性体积力方法计算高速离心压气机叶轮内部流场

高速离心压气机是用途非常广泛的一种动力机械,充分了解其内部复杂的三维流场对设计高压比、高效率离心压气机具有很重大的意义。本文采用一种简化的Navier-Stokes方程的计算方法-粘性体积力方法,计算了Krain叶轮内部流场。使用H型网格离散微分方程,并采用多重网格方法和当地时间步长方法加速计算收敛。粘性应力通过Baldwin-Lomax模型计算。为了证明计算结果的有效性,文中也给出了相应的激光测量结果。      

一种双斜切双压缩面气道性能的数值模拟

采用LU－SGS隐式推进方法及ROE的通量差分分裂格式，结合带非线性限量因子的MUSCL插值方法编写了雷诺平均的三维N－S方程的数值求解程序，并供助该流场计算程序研究了一种双斜切双压缩面进气道在未经流场控制情况下的流动结构和气动性能（M＝1.9,H=11Km）。计算结果表明：所编制的流场计算程序能够准确快速地求解三维粘性流场，并对流场间断具有较高的分辨率；所研究的进气道即使在未经流场控制的情况下也具有较高的总压恢复系数，这体现了双斜切双压缩面进气道总压恢复系数高的优点，但由于流场中存在分离流动，应采取抽吸附面层等流场控制措施。     

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

三级涡流器环形燃烧室化学反应流场的数值研究

在三维曲线坐标系下采用多区域耦合法对包括三级涡流器在内的环形燃烧室三维两相反应整体流场进行了数值模拟。计算中采用区域法和偏微分方程法生成三维贴体网格，采用标准k双方程紊流模型，EBU-Ar-rhenius紊流燃烧模型和六通量热辐射模型。在非交错网格体系下，气相用SIMPLE法求解，液相采用颗粒群轨道模型，并用PSIC算法对其进行数值求解。计算结果表明，程序编制可靠，建立的三维网格生成和流场计算程序可为燃烧室优化设计和研制提供有用的数据。     

液氢诱导轮内部流场数值研究

为有效地预测液氢诱导轮在额定工况下的工作性能,采用了有限体积差分格式、显式时间推进法和多重网络法求解清华大学热能工程系设计的液氢诱导轮内部三维不可压粘性流场。计算表明,该液氢泵设计合理,可以满足工程要求。     

低速流场全方位校准五孔探针测量数据处理方法

论述了低速不可压流场全方位校准五孔探针测量数据处理的计算机程序,并对计算结果进行了对比分析。研究表明,采用速度迭代法构筑二阶拟合多项式对校准和测量数据进行处理,角度偏差不大于1°,速度偏差不大于1％,完全可以满足工程上定量测量三维流场的要求,具有较高的工程实用价值。     

燃烧室三维湍流的数字模拟

文中以SIMPLE方法为基础,采用k-ε两方程湍流模型,发展了一种模拟燃烧室内三维湍流内流场的数值计算方法和计算程序。通过对四种冷流状态进行数字模拟,初步结果表明,理论计算与已有实验结果一致性良好。     

燃气排导箱三维场的数值计算

对燃气排导箱内的典型三维稳态流动过程进行了粘性不可压数值计算,采用较为成熟的ＳＩＭＰＬＥ方法。所发展的计算程序可求解质量、动量和能量守恒方程,ｋ－ε双方程湍流模型方程,以及热辐射通量方程,并进行了程序的算例验证。三维燃气排导箱内的计算结果也得到了一些实验测量的数据的验证,得出一些有工程实用价值的结果和结论。     

静止叶栅三维粘性数值模拟

应用Ｎ－ Ｓ方程压力修正算法和ｋ－ ε两方程湍流模型,对端弯叶片和基准（无端弯）叶片的平面叶栅进行了三维粘性数值模拟。在采用简单交错网格的情况下,通过对计算控制体通量的速度进行修正的方法解决了Ｎ－ Ｓ方程压力修法算中遇到的压力场和速度场呈锯齿状分布的问题。用实验测量值对本数值模拟方法的精度和可靠性进行了检验     

超音进气道流场的数值计算

用Jameson的有限体积法,对某型号飞机在不同飞行状态下的进气道流场进行了数值模拟。求解的控制方程为三维非定常N-S方程。在对流项的计算上用二阶中心差分有限体积法,粘性项的计算则采用作者在［1］中发表的一种有限体积框架下的新的离散化技术,时间方向用多步Runge-Kutta方法进行推进,以增大CFL条件数。湍流模型则采用Baldwin-Lomax代数模型。求解过程中还采用了当地时间步长,隐式残量平均、人工粘性等加速收敛技巧。最后给出了各种飞行状态下进气逼出口截面的平均总压恢复与总压畸变,得出了比较满意的结果。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

分流叶片位置对高转速离心压气机性能的影响

运用三维粘性流动数值计算程序FineTM/Turbo对叶顶间隙泄漏存在时带分流叶片的高转速离心压气机模型级内部流动进行了数值模拟,重点分析了分流叶片不同起始位置及不同周向位置对压气机级内三维粘性流场及整级性能的影响.计算中采用Jameson的中心差分格式结合Baldwin-Lomax代数模型使用时间推进法求解雷诺平均N-S方程,计算模拟了模型级内部复杂的三维粘性流动过程及气体参数分布的详细结构和规律.计算结果表明:采用分流叶片在进口段会减少叶片阻塞,从而使更高的质量流量可以流经叶轮;分流叶片起始位置位于Ⅲ时,两个通道叶轮出口处速度分布最均匀;分流叶片越短,长叶片压力面无量纲静压载荷越大.当分流叶片长度达到某一数值后,长叶片载荷变化趋于平缓.分流叶片位于不同周向位置时,IBSA叶轮的模型级效率最高,压气机性能最好.     

区域分裂的SIMPLE算法用于管路粘性流场计算

为了研究气动系统中节流管的特性和计算节流管内部及尾部射流区的粘性流场,利用ＳＩＭＰＬＥ算法对两端带有节流口的管路粘性流场进行了数值计算,根据管路的几何形状和流动特点,对计算区域共划分了５个区,并给出了相应的分区算法和各种边界条件的处理方法,通过计算得到节流管内各项参数。计算结果表明这种算法不仅可以计算出管路内部的流场,而且可以计算出管路尾部射流区的流场,并可以较好地描述通量的守恒性。     

CISA 3D Navier—Stokes程序在压气机中的鉴别和应用

用压气机叶栅试验数据、多级压气机的第4级静子计算结果和跨声速转子试验结果对CISA 3D N-S程序的计算精度进行了鉴别研究。结果表明,该程序对设计点的落后角的计算精度较高,而对损失的计算还存在一定的分散度。采用 CISA 3D N-S程序和S_2-S_1程序对一台多级压气机试验件的第3级静子进行三维造型修改设计,取得了预期的效果。     

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

三维跨声速粘性-无粘相互作用的半逆方法

本文为计算三维跨声速粘性-无粘相互作用提供一半逆方法版本。湍流边界层采用积分反方法。证明了对于跨声速流,选取外部流线角α和等价的不可压缩形参作为输入量,边界层积分方程总是双曲型的。无粘外流采用FL027程序。通过粘性-无粘迭代求得小分离区情况下的与实验符合的收敛解。     

基于遗传算法与近似模型的全局气动优化方法

开发了基于精英保留策略与小生境技术的求解有约束优化问题的改进遗传算法,并用一多峰值有约束函数对其寻优能力进行了测试。为了降低气动优化设计时间,提出了集成实验设计、三维粘性流场求解程序、二次多项式近似模型、改进遗传算法的具有全局寻优能力的三维气动优化设计体系。采用该方法对NASA Rotor37动叶片叶型进行了以绝热效率最大为目标的优化设计,优化后叶片绝热效率提高1.1%。该结果表明了本优化设计体系省时、高效的特点。