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

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

压缩性对细长体涡流非对称发展的影响

通过数值方法对大迎角细长体湍流流场的模拟,探讨压缩性对细长体非对称绕流发展的影响。结果表明细长体顶端的极小扰动诱发显著的非对称绕流,非对称的涡系结构沿轴向是逐步演化的;在亚临界横流马赫数区间,马赫数越高非对称越显著;在超临界横流马赫数区间,细长体两侧出现横流激波,非对称的发展被抑制,马赫数越高非对称越弱。     

空间高阶差分近似对三维ADI-FDTD数值色散的影响

研究采用空间高阶差分近似的三维交替方向隐式时域有限差分法(ADI-FDTD,Alternating-Direction Implicit Finite-Difference Time-Domain)数值色散问题,首先推导了采用空间高阶差分近似的数值色散迭代公式,分别对最小相速方向和最大相速方向进行二阶、四阶、六阶、十阶空间差分数值色散误差的数值计算和比较,并分析数值色散对空间差分阶数的依赖关系,最后发现在均匀网格中采用四阶空间差分能得到较好的效果.     

钢管混凝土柱极限承载力分析的分段合成法

在计算钢管混凝土受压构件弹塑性极限承载力的合成法基础上,提出了一种新的数值计算方法--分段合成法,计算假设仅包括平截面假设和考虑紧箍力影响的钢管混凝土的应力-应变关系,取消了原合成法中对柱变形曲线的假设,积分时不仅需要在截面上划分单元而且还需沿杆长划分单元,因而可考虑柱多个截面的平衡条件.与原合成法相比,该方法不仅适用于轴压构件和两端偏心距相等的偏压构件,而且还适用于两端偏心距不等的偏压构件的弹塑性极限承载力计算.利用提出的方法和程序,计算了圆形截面钢管混凝土轴压长柱、两端偏心距相等的偏压长柱以及两端偏心距不等的偏压长柱的弹塑性极限承载力,理论分析结果与文献中的理论结果和试验结果均吻合良好.     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

内外函混合加力燃烧室扩压器内液雾湍流两相流场和浓度场数值模拟

 本文用κ-ε湍流模型和液雾轨道模型,对二维内外函混合加力燃烧室扩压器内有蒸发的液雾两相流动和混合进行了数值预报。其中对气相速度场和温度场的计算结果与实验符合令人满意,液雾两相流的计算结果在趋势上也是合理的,表明本文的模拟方法可用于加力燃烧室性能估算及优化设计。     

GH4169合金反挤压成形模具磨损

基于修正 Archard 磨损模型,采用数值模拟方法系统分析了 GH4169合金反挤压成形过程中各挤压工艺参数对模具磨损的影响规律。结果表明：在选取的参数范围内,挤压凸模最易产生磨损失效的区域为凸模圆角处,模具最大磨损深度随凸模圆角半径及坯料预热温度的增大而降低,随摩擦系数的增大而增大;当挤压速率小于100 mm ／s时,模具最大磨损深度随挤压速率的增大而减小,当挤压速率大于100mm ／s 时,模具最大磨损深度随挤压速率的增大先增大后减小。最佳工艺参数坯料预热温度1020℃,摩擦系数0．05,变形速率100mm ／s,模具预热温度300℃时模具磨损量最小,为9．28×10－3 mm。     

高超声速飞行器跨流域气动力／热预测技术研究

分析了跨流域流动特点、工程背景需求、国内外研究现状,针对高超声速飞行器跨流域气动力／热预测存在的主要问题,从粘性干扰参数与克努森数搭接的跨流域模拟准则、微量天平结构优化及测力技术、红外大面积中低量值热流测量技术、流场显示与测量技术和 N-S／DSMC 自适应紧耦合数值模拟方法等方面,总结了近年来取得的研究进展,并讨论了下一步研究方向。     

油膜干涉测量翼型壁面摩阻低速风洞试验技术

油膜干涉试验可以定量测量得到模型壁面的摩擦应力,是开展摩阻特性研究的有力工具。通过对油膜干涉试验技术中的试验原理、试验方法、误差修正、数据结果分析等关键问题的研究,初步建立了试验系统,具备了开展模型壁面摩阻特性研究的试验能力。综合分析表面热膜试验结果与数值计算结果,在以模型弦长为特征长度,试验雷诺数 Re ＝1．3×106的条件下,翼型表面摩擦应力沿模型弦向的分布规律具有较好的一致性。     

气泡和液滴在上升或下降时大变形的数值模拟

本文在考虑界面张力的条件下用涡片方法导出了二维无粘气泡和液滴在另一种静止无粘流体中上升或下降时运动和变形的支配方程，提出了求解支配方程的数值方法，并在数值模拟的例子中研究了空气泡和煤油滴在水中上升时及水滴在煤油中下降时的变形规律，研究了数值模拟中界面张力对界面形状的影响，气泡上升的数值模拟结果与相应的试验结果基本一致。数值模拟结果表明本方法可正确模拟饣泡和液滴在上升或下降过程中的大变形，直到界面发