翼型外形高气动效率/低可探测性的优化

 采用ｖａｎＬｅｅｒ矢通量分裂格式求解Ｅｕｌｅｒ方程的方法计算了绕翼型的气动特性;采用矢通量分裂方法计算了绕翼型的时域电磁散射场特性及雷达散射截面积（ＲＣＳ）;采用一种简单而有效的数值优化方法对流场解和电磁场解进行了翼型外形高气动效率／低可探测性的优化计算。算例结果表明,本方法提供了一种对翼型既可作气动优化设计亦可进行多学科综合设计的有效工具。     

二维Euler方程的自适应网格矢通量算法

给出了二维轴对称欧拉方程的自适应网格隐式矢通量分裂差分格式和模拟结果,隐式矢通量分裂格式不需人工粘性项与人工耗散项,可以计算跨声速流场,也可以较准确地捕获激波,利用流场解的梯度信息对网格进行调整,可以提高激波的捕获精度,算例表明,该方法精度高稳定性好,不依赖于初始条件。     

欧拉方程隐式迭代矢通量分裂跨音流场计算

 本文用矢通量分裂格式的有限体积法求解时间相关欧拉方程组。文中证明,矢通量分裂的有限体积法不能再采用物面反射边界条件。二维问题用隐式迭代法解代数方程组,第三维方向上则用强隐格式往返扫描。此法CFL数可达100的数量级,收敛快,占内存少。跨音流场计算结果与实验符合良好,激波只占一两个网格。     

用高效矢通量分裂格式求解二维叶栅的跨音流动

提出了一种高效隐式矢通量分裂格式。该格式通过在显格式上添加合适的高阶小量, 避免了近似因子分解和系数矩阵运算, 大大地减少了计算量。同时该格式在收敛性及激波捕捉精度等方面均保持了矢通量分裂格式所具有的优势。将该格式用到了二维叶栅的无粘跨音流场的数值模拟中去, 结果令人满意, 证实了无矩阵运算、无近似因子分解的隐格式构造有一定的推广意义。      

矢通量分裂显式格式及二元喷管射流模拟

本文在较常用的显式两步格式的基础上,采用 Steger- Warming形式的矢通量分裂 (简称矢分或 FVS)的原理,提出了一种矢分二阶精度显式两步格式。在分析了格式的各种特性后,与Mac Cormack格式进行了实算对比。又给出了某二元喷管射流问题的模拟实例。计算结果表明,矢分与该显式格式的配合具有矢分法的基本优点和计算稳定、速度较快的优点。      

倒置喷管的湍流流场计算

数值模拟采用基于 Van Leer的矢通量分裂格式 ;应用线 Gauss- Seidel迭代法求解 ;并结合局部时间步长技术对收敛过程加速。对 Re=1 .66× 1 0 6 的倒置喷管流场进行数值计算的结果与实验符合     

基于Steger-Warming矢通分裂的NND格式

将Steger-Warming矢通分裂方法和NND格式相结合构成一种基于Steger-Warming矢通分裂且能自动捕捉到激波的差分算法.通过对激波管和对可燃气体DDT过程的数值分析,证明该算法解决了Steger-Warming矢通分裂方法在零特征点附近存在数值解震荡的问题,有着良好的计算精度和计算效率.文中采用的计算方法对各种有关矢通分裂计算方法的实现有一定的参考价值.     

平面叶栅跨音流动的数值计算及变工况性能的研究

提出一种隐式矢通量分裂差分格式并用来直接求解Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组。该方法避开了近似因子分解及矩阵运算，具有精度高、稳定性好、计算量少等优点。在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

近空间高超声速流场通量分裂型DSMC-IP方法

针对近空间高超声速流场的特点,采用Van Leer格式对直接模拟蒙特卡洛-信息保存(DSMC-IP)方法质量守恒方程中的计算格式进行改进.将局部马赫数作为分裂通量的标准,并重构单元分界面两侧的左右输运通量,使得计算格式具有通量分裂的特点,解决了DSMC-IP方法在高超声速流场计算中的应用问题.将改进后的通量分裂型DSMC-IP方法引入非结构网格中,对二维近空间高超声速流场进行数值模拟.计算结果表明:通量分裂型DSMC-IP方法所得出的数值结果与实验值及参考值符合较好,明显降低了直接模拟蒙特卡洛(DSMC)方法所带来的统计耗散.当来流气体的稀薄程度增加时,其非平衡效应也更加明显,而通量分裂型DSMC-IP方法的计算结果与参考值相差均在10%以内,良好地反映了非平衡条件下的流场特征,验证了通量分裂型DSMC-IP方法的可行性和有效性.      

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

扩压器内跨音速湍流的数值模拟

采用Johnson-King非平衡代数雷诺应力湍流模型(J-K模型)和Baldwin-Lomax零方程湍流模型(B-L模型),数值模拟较强激波/边界层相互作用时扩压器内的分离流动。计算结果与实验值进行了比较,表明J-K模型比B-L代数湍流模型可较好地计算出分离流动的再附点位置,并且可更好地计算出激波强度和沿流程的压力分布,仅增加很少的计算量,并更易推广应用于三维湍流问题的数值模拟。      

跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

扩压器进气道内流及激波／边界层干扰的数值研究

应用隐式矢通量分理解差分格式结合改进的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型直接求解时间平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，对扩压器跨声流动进行了数值计算，准确地捕获了激波，预测了由于强激波作用引起的湍流边界层分离现象，与实验数据比较，结果令人满意。数值实验表明：该方法具有精度高、稳定性好、计算量少、收敛快等优点，计算结果可用于高速飞行器气动特性的研究 。     

关于两方程湍流模型的考虑

通过对两方程湍流模型的分析,提出了在主流区或远壁区用两方程模型,近壁区用代数湍流模型的思路,即所谓的“双层”模型概念。通过对平板湍流附面层的计算表明,这种分层计算措施获得的速度分布、壁面剪切力较单纯应用低雷诺数两方程湍流模型时的结果准确。对非均匀来流条件下90°弯管内高湍流度流场的计算表明,双层模型可获得更加准确的速度分布。从而表明本文提出的分层计算措施是有效的,可满足工程计算要求。      

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

离心叶轮内稀释固液两相流动的数值模拟

 对固液两相流动进行了数值模拟。液相和固相分别采用 k- ε湍流模型和代数颗粒流( AP)湍流模型并考虑了液固两相间的相互作用。用此模型结合 SIMPLE- C算法对离心叶轮内的稀释固液两相流动进行了计算并研究了颗粒对流场的影响。计算结果同试验结果基本一致     

环形燃烧室冷态流场数值模拟中的数学方法

针对带双级径向旋流器环形燃烧室火焰筒,研究不同紊流模型(standardk-εmodel,RNG(renormalization group)k-εmodel和代数应力模型(algebraic stress model))及算法(SI MPLE(semi-i mplicitmethod for pressure-linked equations),PISO(pressure i mplicit split-operator))对复杂流场求解准确性及效率的影响.火焰筒三维贴体网格采用型线定点与偏微分方程相结合的方式生成,用壁面函数法处理近壁区域流动.计算结果表明,PISO算法提高求解效率;不同紊流模型预测结果均与试验数据符合较好,其中代数应力模型最为准确,更适用于复杂流场的数值预测.      

气固两相流k-E-kp-Ep双流体模型及应用

以流体k-E双方程模型为基础, 提出了一种包括颗粒湍动能输运方程和耗散率方程的k-E-k_p-E_p两相湍流模型。并对两个垂直上升管道内气固两相湍流进行了数值计算, 计算结果和已有实验结果吻合得很好。