对流扩散方程的QC3无网格法

对流扩散方程作为偏微分运动方程的分支,在流体力学、气体动力学等领域有着重要应用.为解决对流扩散方程难以通过解析法得到解析解的难题,采用二阶一致3点积分(Quadratically Consistent 3-Point Integration,简称QC3)提高无网格法的计算效率,通过对积分点上形函数导数的修正,改善无网格法的精度和收敛性.本文将QC3无网格法拓展到对流扩散方程问题中,时域离散采用广义特征线Galerkin法,空间离散采用QC3法.数值结果表明,应用QC3无网格法得到的对流扩散问题数值解与解析解十分接近,验证了QC3无网格法解决对流扩散问题的可行性.     

高熵合金辐照硬化与力学性能研究

高熵合金由于多主元元素混合引起高熵结构效应，使其具有优异的物理、力学和化学特性，如高强度、高耐磨性、耐蚀性、热稳定性、优异的抗辐照性能等。然而，辐照诱发高熵合金材料的硬化行为和力学性能预测仍缺少相关研究，严重地限制了对其长期服役后材料性能的评估。基于晶体塑性理论结合实验结果，研究了空洞形状依赖的硬化行为、位错环诱发的硬化行为以及氧化物弥散增强的高熵合金力学性能。研究发现，考虑多面体空洞与位错的概率依赖的空间交互作用，能够更加准确地预测辐照金属的屈服应力；晶格畸变对屈服强度，有着重要的贡献；氧化物弥散相对位错运动起强烈钉扎的作用，从而对强度产生影响，直接决定抗辐照性能。高熵合金作为一种具有综合优异力学性能的新型结构材料，在先进核能系统中有望被广泛应用，比如核反应堆的核燃料包壳管。     

竖直恒温面自然对流边界层不稳定性研究

在无扰动、随机式扰动以及正弦式扰动下,通过对竖直恒温面处状态Ra为1.328×10^9、Pr为6.24的自然对流进行模拟,探索了热边界层的不稳定性和共振强化自然对流换热。结果表明:(1)竖直自然对流边界层上游位置的随机式扰动对热边界层的影响主要体现在稳定阶段;(2)该状态下的竖直自然对流边界层的特征频率为15 067,且相比于无扰动状态,频率为15 067的正弦式扰动能在竖直恒温面处提高5.15%的换热量;(3)在竖直自然对流边界层上游位置加入特征频率的正弦式扰动,竖直恒温面处的局部努塞尔数Nu均出现明显波动,且波动随着边界层高度的增加而增大。     

对流扩散方程的缩减基有限元法的后验误差

将缩减基(RB)方法和有限元方法相结合,在保证偏微分方程的有限元离散格式具有足够高精确度前提下,能够大幅度地降低有限元离散格式的维数,从而大大降低计算中内存容量和计算时间的消耗.针对对流扩散方程建立基于RB方法的Crank-Nicolson有限元离散格式,并给出后验误差估计结果.     

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

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

倾斜层中的对流斑图及其临界条件

通过二维流体力学基本方程的数值模拟，探讨了Prandtl(普朗特)数Pr=6.99时，倾斜矩形腔体中的对流斑图和斑图转换的临界条件.根据倾角θ和相对Rayleigh(瑞利)数Ra_r的变化，倾斜矩形腔体中的对流斑图可以分为:单滚动圈对流斑图、充满腔体的多滚动圈对流斑图和过渡阶段的多滚动圈对流斑图.当θ一定时，随着Ra_r的减小，系统由充满腔体的多滚动圈对流斑图过渡到单滚动圈对流斑图.这时，对流振幅A和Nusselt(努塞尔)数Nu随着Ra_r的增加而增加.当Ra_r=9时，随着θ的增加，系统由充满腔体的多滚动圈对流斑图过渡到单滚动圈对流斑图，这时对流振幅A随着θ的增加而减小，Nusselt数Nu随着θ的增加而增加.在θ_c-Ra_r平面上对多滚动圈到单滚动圈对流斑图过渡的模拟结果表明， 在Ra_r=2时, 腔体中没有发现多滚动圈对流斑图.在Ra_r为2.5左右时，腔体中出现多滚动圈到单滚动圈对流斑图的过渡.当多滚动圈到单滚动圈对流斑图过渡的临界倾角θ_c<10°时，θ_{c随着Rar的减小而增加.当θc>10°时，θc随着Rar的增加而增加，在Rar≤5时，θc随着Rar的增加而迅速增加；当Rar>5时，θc随着Rar的增加而缓慢增加.θc与Rarθ的关系与Rar类似}     

求解一维对流方程的高精度紧致差分格式

本文提出数值求解一维对流方程的一种两层隐式紧致差分格式,采用泰勒级数展开法以及对截断误差余项中的三阶导数进行修正的方法对时间和空间导数进行离散.格式的截断误差为O(τ~4+τ~2h~2+h~4),即该格式在时间和空间上均可达到四阶精度.利用von Neumann方法分析得到该格式是无条件稳定的.通过数值实验验证了本文格式的精确性和稳定性.     

入口节流微通道换热器内相变传热特性

设计加工一种带有入口节流结构的铜基微通道换热器,理论分析其传热模型、实验测量微通道换热器内相变换热的传热特性和压力特性。结果表明:换热器内部的热传递过程为其主要换热模式;换热器表面温度随加热热流密度的增大而增大;微通道入口流速对表面温度影响较小;入口工质过冷度线性影响换热器的表面温度。热流密度在不同阶段对换热系数有不同影响,热流密度为360 W/cm~2时,换热器换热系数出现最大值;换热器压降随热流密度和系统流速的增加而增大。     

求解一维定常对流扩散反应方程的混合型紧致差分格式

借助显式紧致格式和隐式紧致格式的思想,基于截断误差余项修正,并结合原方程本身,构造出了一种求解一维定常对流扩散反应方程的高精度混合型紧致差分格式.格式仅用到三个点上的未知函数值及一阶导数值,而一阶导数值利用四阶Pade格式进行计算,格式整体具有四阶精度.数值实验结果验证了格式的精确性和可靠性.