自适应直角坐标网格方法在耦合传热数值模拟中的应用

首先介绍了耦合传热计算区别于一般的导热和对流换热过程的特点,然后介绍了常用的计算方法应用于耦合传热计算遇到的困难及自适应直角坐标网格方法在耦合传热计算中的优势。最后用两个算例证明了采用自适应直角坐标网格方法数值模拟耦合传热问题是适用的。     

空域空间网格标识原理及应用

在对空域空间有限元的相关研究中,构建空域空间的网格标识系统,是将信息技术应用于空域管理、推进空域数字化的重要基础。在美军全球区域参考系统和通用地理位置参考系统基础上,提出了一种同时具备标识平面和立体空间位置能力的参考系统。通过在地球表面以上特定高度之下的空间内,建立一种递归剖分策略方法,将空域空间划分为若干空间网格,不同尺度下的网格既是一种地理空间位置参考系统,也是一种组成空域的基本空域体,以此实现空间位置参考与空间位置标定的混合表达,为开展空域数字化管理提供底层模型。      

超音速分离线喷管大摆角状态下化学烧蚀动态仿真

针对超音速分离线喷管大摆角状态下化学烧蚀导致的壁面退移,基于动网格技术建立了相应的动态仿真模型,实现了对不同燃烧室条件下喷管化学烧蚀率的预示。初步稳态计算得到喉部烧蚀率为0.048 6 mm/s,高出试验结果5.67%,验证了仿真设置的合理性。以此状态结果为瞬态计算的初场,进行相应的化学烧蚀动态仿真计算。该喷管的矢量角放大系数在0.5 s仿真时间内因壁面退移减小了0.42%,对称面下侧分离线结构附近因燃气流动受阻成为喷管烧蚀最严重的位置,烧蚀率为0.074 5 mm/s。增大燃烧室压强或温度,会导致同周向位置的分离线后侧与前侧壁面烧蚀率比值减小。对于分离线附近型面变化较小处,压强5.5 MPa增加到7.5 MPa,该比值减小了10%,温度3200 K增加到3600 K减小了15%。     

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

基于在轨组装维护的多载荷协同观测卫星平台构型设计

针对高轨遥感卫星多载荷协同观测和一致测量应用需求,在前期总体技术方案研究的基础上,通过开展现有在轨组装航天器和高轨卫星公用平台的构型调研,提出了一种适用于多载荷协同观测需求的平台构型设计方案。平台利用CZ-5运载火箭发射并在轨展开,采用“网格式中心承力筒+板+桁架”式的总体构型,开展了主传力路径、贮箱承载结构、网格式中心承力筒等的细化设计工作,并进行了结构产品的初步设计,建立了平台结构整星的有限元模型。仿真分析表明:平台结构整星的各向一阶频率均能满足运载火箭的频率要求。基于标准化、通用化的平台接口设计,使其具备可扩展及可维修的能力,大大提升了平台鲁棒性。基于全网格筒的主承力设计,增强平台结构承载能力,并改善平台内部热环境。该设计方案可为基于在轨组装维护的高轨遥感大平台结构的设计提供有价值的参考。     

飞机结冰中水滴撞击特性欧拉法的网格影响分析（英文）

过冷水滴撞击特性计算研究是飞机结冰预测与防除冰系统性能分析的基础,欧拉法是过冷水滴运动及撞击特性计算的常用方法之一,分析欧拉法在计算复杂表面水滴运动撞击时的准确性具有重要意义。以NACA0012翼型、冰风洞风道、S形进气道及三段翼为研究对象,利用欧拉法采用不同网格计算获得了水滴运动及局部水收集系数,结果表明：当水滴运动没有受到上游部件的导流或者遮挡时,欧拉法计算结果具有网格无关性;当水滴运动受到上游部件影响发生轨迹偏转时,欧拉法计算结果具有网格相关性,不同网格计算得到的液态水含量及局部水收集系数不一致,需要进一步考虑水滴运动受到上游效应影响时网格对欧拉法计算结果的影响。     

随机场TPDF在非结构网格平台的开发及检验

为实现对燃烧室内部燃烧过程的高精度数值模拟,增强航空发动机燃烧室研究手段,采用开源的非结构网格平台Saturne,开发可以耦合详细化学反应机理,高精度解析湍流燃烧过程的概率密度函数输运方程（TPDF）湍流燃烧模型。采用随机场方法求解TPDF方程,在原有程序基础上研发了TPDF程序模块及配合的加速算法等模拟单元,发展了针对燃烧室燃烧过程模拟的软件功能。采用射流火焰Flame D和旋流火焰TECFLAM对新软件进行了测试,结果表明：模拟结果与试验数据趋势一致,精度尚可。对某型燃烧室性能模拟验证中,计算的出口温度分布与试验结果一致,平均温度误差小于3.8%,新的软件可进一步应用于燃烧室性能的研究。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

多维对流通量重构方法的改进

为了研究多维对流通量重构格式中不同形式主方向对激波稳定性、压力场计算方面的影响,分别采用以笛卡尔方向、密度梯度及压力梯度作为主方向进行计算,对结果分析后发现,单纯依赖物理量梯度的主方向会造成压力等值线曲折的非物理解。为了消除这种非物理现象,采用一种基于压力函数的混合方式来确定主方向。单层及展向多层高超声速圆柱绕流问题模拟分析的结果表明:混合方式所确定的主方向不仅可以消除笛卡尔坐标系确定的主方向计算带来的红玉现象,同时还可以消除单纯的密度或压力梯度作为主方向所带来的压力场"zigzag"的非物理现象;通过与结构网格结果的比较表明:网格相关格式驻点处热流的相对误差为29.2%,而混合主方向确定方法所计算驻点热流相对误差只为4%,表明此改进方法计算得到的壁面热流是合理的。     

颗粒无序分布的堆积床内部流动与传热分析

为了详细分析单组元发动机催化床内部的流动与传热,针对颗粒无序分布的堆积床,提出了一种网格生成途径,适合于对大量颗粒自然堆积的管道流动进行数值模拟。根据该方法,对不同颗粒直径、不同进口速度的管道堆积床进行了液体介质的流动与非定常传热模拟,并与实验结果和Ergun方程计算进行了对比。结果表明,本文提出的方法可以有效模拟颗粒无序分布的堆积床内部流动与传热,适合于大量颗粒的堆积床仿真;网格尺度小于1/20颗粒直径时,模拟结果与实验结果符合很好;在相同雷诺数下,摩擦系数随颗粒直径增大而减小,Ergun方程在高雷诺数下计算的流阻偏大;非定常传热时,壁面效应会使壁面附近温度先接近流体温度。