高速气流中激光加热平板数值模拟与分析

 采用流固耦合方法,数值模拟了高速流场中激光作用下来流速度对平板温度分布的影响。结果表明:无激光辐照时,高速气流中平板有较高的气动生热平衡温度,且平板-气流之间的换热系数随来流速度增大而增大;在平板前沿换热系数增长最快,沿平板长度方向增速趋于平缓。分析了激光辐照时高速气流中激光加热平板的温度分布情况,考察了来流速度不同时,气动生热、散热和激光辐照对平板温度的影响,给出了激光辐照后的温升情况和温度分布,分析了在不同速度来流下,对流散热、摩擦生热和激光加热之间的竞争关系,结果表明,平板温度具体分布主要是加热过程竞争的结果。     

Convection-Dominated Accretion Flows with Radiative Cooling

By numerically solving the set of basic equations describing black hole accretion flows with low accretion rates, we show that although the dynamical structure of these flows is essentially unaffected by radiative processes in comparison with the case in which the radiation is not considered, the radiative cooling can be more important than the advective cooling in the flow＇s convection-dominated zone, and this result may have implications to distinguish observationally convection-dominated accretion flows from advection-dominated accretion flows.     

基于蒙特卡罗随机有限元方法的随机多孔介质内流体自然对流不确定性研究

为分析孔隙率不确定性对多孔介质方腔内自然对流换热的影响,发展了一种基于KL（Karhunen-Loeve展开）-蒙特卡罗随机有限元算法的随机多孔介质内自然对流不确定性分析数理模型及有限元数值模拟程序框架。通过K-L展开及基于拉丁抽样法生成多孔介质孔隙率随机实现,并耦合多孔介质自然对流有限元程序,进行随机多孔介质内自然对流传热数值模拟,得出了多孔介质内流场与温度场平均值与标准偏差,并分析了孔隙率不确定性条件下Da数对Nu数的影响。结果表明,孔隙率不确定性对多孔介质方腔内自然对流有重要影响。随机多孔介质内流场及温度场与确定性条件下的流场及温度场存在一定偏差,Nu数标准偏差随着Da的增大先增大后减小。     

强约束球形装药反应裂纹传播和反应烈度表征实验

炸药燃烧的高温高压气体产物可以进入基体裂纹中引发炸药表面热传导燃烧,形成所谓的对流燃烧。在一定约束条件下,不断上升的气体压力反过来又使炸药基体产生更多的裂纹,为对流燃烧提供更多的通道和燃烧表面积,快速生成大量产物气体导致高烈度反应现象的产生。本文中设计了一种新型强约束球形装药中心点火实验,针对一种HMX为基的PBX炸药,对高烈度反应条件下燃烧裂纹传播和反应增长过程进行了观测,实验中采用测得的反应压力和壳体速度历程对反应烈度进行了量化表征。在带窗口结构中,早期炸药中的燃烧裂纹不可见;中期燃烧裂纹扩展到药球表面时,先形成4条沿经线方向近似对称的主裂纹,随后环向贯通并扩展到整个药球表面;最后的剧烈反应造成强烈发光。上述反应演化经历低压增长阶段约为100 μs,之后伴随着壳体变形膨胀产生剧烈的反应,此时产物压力在约10 μs时间内超过1 GPa,并形成约20%相对于裸炸药爆轰的超压输出。在全钢结构中,20 mm厚的壳体膨胀速度最大可达到500 m/s,此时壳体完全破裂。

     

薄膜沿加热平板下落的稳定性分布

本文分析了薄膜沿加热平板下落的稳定性。在时间模式下，发现流动的不稳定性是由表面波不稳定和加毛细不稳定构成的，同时当流体的热扩散越大以及界面热量损失越小时，热毛细不稳定越剧烈，在时空模式下，流动随着Marangoni数的增大。流动有可能从对流不稳定过渡到绝对不稳定，这一结论尚待实验验证。     

大尺寸液桥热毛细对流失稳性地面实验研究

开展大尺寸液桥浮力-热毛细对流地面实验, 探究流场转捩的临界条件及临界状态附近的流动情况. 通过粒子图像测速方法(PIV) 获得流体速度场, 研究液桥内部定常和转捩后的流场结构以及流体运动规律;并用红外热像仪测量液桥自由面温度分布, 研究流体流动的时空演化和温度振荡. 实验发现大尺寸半浮区液桥浮力-热毛细对流临界值与几何参数有关, 在大普朗特(Prandtl) 数情况下, 流场存在由稳定态向不稳定态再到混沌的转捩过程, 在临界马兰哥尼(Marangoni) 数附近, 流场内会出现行波现象, 流动模式也会随高径比的变化而发生变化;当继续增大马兰哥尼数, 流动会进入混沌状态.      

热毛细对流速度场测试研究

为了探索微重力状态下表面张力驱动流速度场的测试技术，本文在地面进行了模拟实验，对位于热壁下流体中气泡周围的表面张力驱动流进行了研究，并用ＰＩＶ技术测量了流场的速度分布．     

浸没边界-简化热格子Boltzmann方法研究及其应用

针对流固耦合传热问题,本文提出了一种基于浸没边界-简化热格子玻尔兹曼方法(immersed boundary method-simplified thermal lattice Boltzmann method,IB-STLBM)的耦合模型.不同于传统的格子玻尔兹曼方法使用分布函数演化流场和温度场,简化热格子玻尔兹曼方法(simplified thermal lattice Boltzmann method,STLBM)的演化过程不需要依赖分布函数,只涉及平衡态分布函数和非平衡态分布函数,能够直接演化宏观量,极大减小了计算过程中所占用的虚拟内存,简化了边界条件的实现方式,同时具有较高的稳定性.传统的浸没边界法对流场的计算采用欧拉网格,对固体边界采用拉格朗日网格,认为固体边界是对流场产生某种体积力.在应用浸没边界法时,汲取介观的思想,把固体的介入看作是对流场的干扰,打破了固体附近流体介观微团颗粒原始的平衡状态,这种干扰可以看作是在耦合边界上产生的一个非平衡项,可用非平衡态分布函数来表示.基于此,在模型中浸没边界法与简化热格子玻尔兹曼方法更紧密联系在一起,更大程度发挥二者的优点,整个计算过程更加简单直观,符合物理特性.通过对热圆柱绕流和内含热颗粒的封闭方腔自然对流问题的模拟以及对其结果的分析,验证了该算法在求解流固耦合传热问题的有效性和可行性.      

一种新的求解对流扩散边值问题的无网格方法

基于配点法和楔形基函数,提出了一种新的求解对流扩散边值问题的无网格方法。通过一维和二维的问题验证了该数值方法的可行性;并根据数值算例和分析,可以看到该数值方法能达到满意的收敛效果。该数值方法的隐格式形式能够有效地消除对流占优问题的数值振荡现象,是一种真正的无网格方法。     

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。