三维复杂区域对流换热问题的控制方程及计算

研究并发展了一套计算复杂区域对流换热问题的通用方法并开发了相应的可视化计算程序。针对所建立的一组以直角坐标速度分量表示的非正交曲线坐标系下的三维流动及传热控制方程,采用有限体积方法离散,SIMPLER算法求解,并利用网格界面速度插值方法解决采用非交错网格时出现的压力波动问题。本文利用所开发的程序计算了低散热发动机冷却通道自然对流换热问题。      

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

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

方形截面导弹摇滚特性数值研究

采用隐式耦合求解RANS方程与滚转运动方程的方法研究了方形截面导弹的摇滚运动特性。耦合求解过程采用隐式处理提高了求解精度,采用强迫摇滚得到了自由滚转静导数和动导数;采用与摇滚方程的耦合求解得到了方形截面导弹模型在不同迎角下的自由摇滚历程。两种研究结果均显示随着迎角的增加,方形截面导弹将由滚转稳定变为不稳定,并形成极限环振荡,即摇滚;随着迎角进一步增加将形成滚转发散。两种分析方法所取得的结果与试验结果吻合较好。     

对流-扩散型方程的一种简单、有效的欧拉-拉格朗日分裂格式

本文提出了一种简单而有效的数值求解一维对流-扩散型方程的欧拉-拉格朗H(E.-L.)分裂格式。无论对于线性对流-扩散方程和非线性Burgers方程,或者是对流占优或非对流占优,它都能给出较好的结果。文中给出了四个典型而又较困难的计算例子,与精确解的比较表明,数值结果是相当令人满意的。     

非平衡流解耦方法及其计算激波诱导燃烧的应用验证

介绍并发展了一种改进的化学非平衡流动解耦方法.将流动控制方程与化学反应生成源项解耦处理,组分对流项的矢通量分裂基于流动方程密度对流项的分裂形式给出.对流项采用五阶WENO格式求解,化学反应源项的刚性采用简化的隐式方法进行处理.将该方法应用于预混H₂/O₂和H₂/Air激波诱导燃烧计算,得到的定常及非定常燃烧过程与实验观测相符,振荡频率与实验测量较吻合,结果表明该方法应用于多组分多步反应流计算是可行的.      

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

纤维增强复合材料币状I型裂纹桥联问题

利用平面等效方法，将轴对称币状裂纹问题转化平面裂纹问题，由此给出复合材料内币状Ⅰ型裂纹的桥联积分控制方程。对权函数法应用于求解币状裂纹问题的可行性提出质疑，并给出修正的建议。然后对币状Ⅰ型裂纹桥联问题给予了断裂力学分析。利用一种新的分段插值方法求解了桥联积分控制方程。针对该积分程提出了一种近似分析解法，其结果与数值结果符合得相当好。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

非正交界面下的插值耦合传热方法

针对复杂热防护结构,开发了一种求解非正交性流/固界面的耦合传热程序。流体和固体区域采用同一积分、守恒型的RANS(Reynolds-Averaged Navier-Stock)方程,通过有限体积法进行离散求解。为了保证非正交界面上的温度和热流密度连续,提出了一种结合网格周边信息计算界面热流密度的插值方法。利用该插值耦合方法模拟了双层复合平板和喷管的热传导。数值结果表明:该插值方法在空间上具有2阶精度;喷管壁面上的对流换热系数沿轴向先增大后减小,在喉部上游达到最大值,当喷管入口压力增加3.38倍时,对流换热系数的最大值相应增加了3.13倍;喉衬与壳体界面上温度的计算值和试验结果存在一定差异,这是由于本文数值计算未考虑接触热阻引起的。     

数值微分的小波最小二乘方法

给出了利用小波求解数值微分的方法.将微分问题转化为等价的第一类Fredholm积分方程,然后选取小波基,使用最小二乘求解该积分方程.我们给出的分析框架容许使用具紧支集的小波,利用正则化保证了计算的稳定性.数值计算实例表明,本方法稳定,运算速度快,计算精度高.     

多孔介质有效导热系数的计算方法

高孔隙率多孔介质如泡沫陶瓷在新型多孔介质燃烧器技术中应用日益广泛,其重要传热特性参数--有效导热系数反映了两相流气、固相导热、对流和辐射的综合效应,对其研究尚非常缺乏.本文基于实验测定的温度分布,给出多孔介质有效导热系数的初始估值,用有限体积法求解二维控制方程,采用二维寻优搜索的办法,确定使测定点上测量与计算温度均方根误差为最小的径向与轴向有效导热系数,是一种逆计算方法.对球粒子颗粒床进行的有效性试验证明了方法的可行性.     

扩缩管充分发展层流的强化传热场协同分析

采用二维轴对称不可压缩层流N-S方程,利用有限体积方法和SIMPLE算法,在雷诺数100～1300范围内,对轴对称周期性扩缩管的充分发展层流流动和传热特性进行了数值模拟,重点研究了不同雷诺数、几何外形与换热强度的关系,并利用场协同理论对计算的结果进行了分析,获取了层流态下扩缩管对流换热强度随雷诺数、几何尺寸变化的规律,对进一步研究扩缩管换热问题、探讨强化传热理论和发展换热技术均有实际意义.     

超声速对流推力矢量喷管的数值模拟

利用二次流对流剪切激励主射流方式来产生射流对流推力矢量方法是一种效率很高的矢量化喷流新方法,又称为对流推力矢量法。本文运用三阶精度MUSCL TVD有限体积格式求解二维非定常、可压缩、雷诺平均N-S方程和Spalart-Allmaras湍流模型,计算了对流推力矢量喷管的超声速射流问题,获得几种不同二次流吸气压强下对狭缝和喷管肩部壁面压强分布、剪切层诱导旋涡的结构特征和主射流矢量角等结果,计算结果与实验数据对比取得良好的一致性,并对计算结果进行了详细分析和讨论。     

改进节点积分的无单元Galerkin法及其在流动问题中的应用

无单元Galerkin法需要在背景网格上积分,计算量大,且在求解对流占优问题时会出现非物理的数值伪振荡现象。为此,基于局部Taylor展开思想,采用节点处的局部Taylor展开计算积分,建立了局部Taylor展开积分无单元Galerkin法。该方法同时解决了标准的无单元Galerkin法计算量大和对流占优时会出现数值伪振荡的问题。一维定常对流扩散方程和二维Burgers方程的求解说明了该方法的有效性,且计算效率远高于无单元Galerkin法。     

母线对扩缩管传热强化影响的计算研究

采用二维轴对称不可压缩N-S方程,利用有限体积方法和SIMPLE算法,对不同母线型式的扩缩管对流换热特性进行了数值模拟.在给定入口速度、温度型线的条件下,分别研究了层流和湍流状态下均匀、线性、抛物线等三种典型入口型线对对流换热性能的影响.探索了不同母线型式对传热强化的影响,得到了对流换热随母线形状的变化规律.在给定的状态中得到了一个比较适当的母线外形,使得扩缩管的传热性能有所优化.     

发动机经验故障方程的建立与求解

本文提出了根据发动机故障样本建立发动机故障方程并利用最优化方法(主特征量模型)求解的故障诊断方法。经验故障方程是不同于发动机小偏差方程的另一类故障方程,它为发动机故障诊断提供了一条简便易行的途径,并且使发动机故障诊断范围扩展到气路分析方法难以适用的场合。文申讨论了经验故障方程的建立方法,求解方法和应用特点,给出了利用经验故障方程进行发动机故障诊断的实例。     

长尾喷管传热及烧蚀特性研究

针对长尾喷管的热防护问题,本文采用实验和数值仿真结合的方法,研究了采用SQ-2推进剂固体火箭发动机复合结构长尾喷管的传热以及烧蚀特性。仿真程序采用基于格心的有限体积法,对流体域求解Navier-Stokes方程,对固体域求解热传导方程,湍流模型采用k-ω SST模型,通过保证热流密度大小相等、温度连续的方法实现耦合传热计算;在以上仿真计算方法的基础上建立了基于热解动力学的变热物性模型描述碳/酚醛材料的体积烧蚀,并采用简化的异相反应模型对喷管热化学烧蚀量进行计算。结果表明在收敛段和等直段区域压强变化较小,炭化层厚度和烧蚀量与温度分布相似, 并在距离等直段入口约1.6mm位置达到极小值。碳/酚醛材料炭化速率随时间减小,烧蚀速率随时间增长趋于稳定,而喉衬区域C/C复合材料烧蚀速率随时间增加。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。