论文部分内容阅读
多孔介质内的流体流动和传热传质普遍存在于工农业生产及环境工程中,如核废弃物处理、热管、干燥过程、太阳能集热管、燃料电池及石油热采工艺等,其中部分填充多孔介质复合系统内的流动及热质传递问题具有广泛的应用背景。复合系统内的对流及传热传质机理复杂,尤其是交界面区域的流体动力学特性即滑移效应对腔体内的流动和传热传质机理有着重要影响。在深入研究单独的多孔介质和流体空间传递问题基础上,对复合系统的流动和传热传质研究也日趋深入。其中对多孔介质和流体空间交界面处滑移特性的正确把握是认识复合系统内流动及传热传质规律的基本前提。本文以理论分析、数值模拟与实验研究相结合的方法探讨了多孔介质与流体空间交界面处的滑移效应及其对复合腔体内的流动及传热的影响机理。采用多尺度方法对复合腔体内的自然对流传热过程进行理论分析,利用体积平均法建立两区域自然对流传热数学模型,并推导出界面处的滑移边界条件。本文从孔隙尺度下的流动分析入手,建立质点的连续性及动量传递方程。利用相平均算子,对连续性及动量方程进行相平均计算,建立基于介观尺度下表征体元的连续性及动量守恒方程。在流体空间内质点动量守恒方程等同于不可压缩粘性流体的纳维-斯托克斯(Navier-Stokes)方程,在均质多孔介质区域内,动量守恒方程可表示为达西-布瑞克曼-福希海默(Darcy-Brinkman-Forchherimer)方程。为获得界面处的边界条件,首先将不受任何尺度因子限制的表征体元动量守恒方程在整个复合腔体区域内积分,然后将Darcy-Brinkman-Forchheimer方程和Navier-Stokes方程分别在均质多孔介质区域与和流体空间内进行积分。利用前者与后者积分方程之差,可推导出界面处表示滑移效应的应力滑移效应边界条件,同时获得应力滑移系数的一般表达式。利用前述方法,对传热过程进行分析。采用相平均算子,对孔隙尺度下的流体相及固体相的质点能量传递方程进行计算,建立介观尺度下基于表征体元的能量守恒方程,并分别获得均质多孔介质区域和流体空间内的能量守恒方程。由局部热平衡假设,获得界面处温度和热通量连续的边界条件。通过连续性方程、动量传递及能量传递方程的确定建立同时包含多孔介质区域和流体空间的两区域自然对流传热数学模型。在表征体元的基础上可获得参数平均值,并将其代入标准传递方程中,获得宏观变量的传输规律。采用有限元数值方法,辅之以界面处的弱约束处理,求解宏观尺度下复合腔体内的自然对流传热数学模型,着重探讨了界面滑移效应及其对复合腔体内自然对流流动及传热的影响。界面条件及相关应力滑移系数的确定是获得流场与温度场分布的前提,本文利用介观尺度下界面区域的结构参数对粘性应力滑移系数进行求解。选用三种孔隙率与渗透率的经典经验关联式计算得到渗透率及界面区域内的单位表面积,进而求解出该条件下的粘性应力滑移系数。以界面处速度及其通量均连续的连续性条件为基准,通过中截面处的垂直速度变化曲线与腔体左壁面处的平均传热速率变化分析在粘性应力滑移效应和惯性应力滑移效应条件下对腔体内的流动及传热变化的影响规律。对流动进行数值分析的结果表明,在前两种经验关联式下,粘性应力滑移效应的存在使界面两侧速度差值增大,速度明显增加,且随孔隙率的增加,该趋势也明显增强。但第三种经验关联式下,粘性应力滑移效应的影响影响使界面处出现反向流动,仅在孔隙率大于0.9时反向突变点才逐渐消失。惯性应力滑移效应使界面两侧的速度差值减小,曲线变化更为缓和,但仅在孔隙率大于0.9时有较明显影响,且影响范围限于界面附近狭小区域,与连续性条件下的流场相差不大。对传热进行分析的数值结果表明,粘性应力滑移效应对传热的影响与孔隙率和孔径尺寸相关,当孔隙率大于0.7时,随孔径增加,粘性应力滑移效应使腔体内的平均传热速率增加,当孔隙率小于0.7时,平均速率可能增加或减小。在粘性应力滑移效应条件下,存在一个渗透率临界值使传热的影响出现逆转,该值因选用经验关联式不同而异。当渗透率大于该值时,传热速率大于连续性条件下的值,反之则小于其值。与之不同,在任何孔隙率下,随着孔径增加,惯性应力滑移效应均使腔体内的平均传热速率增加,且随着孔隙率的增加,增速加快。不论在何种条件下,惯性应力滑移效应的存在都略微削弱了腔体内的平均传热速率,但其影响很小。搭建可视化部分填充多孔介质复合腔体实验台,利用粒子成像测速系统(PIV)对复合腔体内的自然对流流场进行了测试,重点关注了界面处的流动特性。通过在两竖直壁面间设定不同的温差实现多种瑞利数工况,利用PIV对各工况下的稳态自然对流流场进行拍摄并记录图片,通过互相关计算及图片后处理获得流线分布以及流动速度和剪切应力等参数值。利用前述三种经验关联式计算得出复合腔体内多孔介质区域的渗透率,结合实测的速度及剪切应力值对粘性应力滑移系数进行求解。结果表明在实验测试孔隙率下,粘性应力滑移系数随渗透率、界面高度及瑞利数的不同而改变。考虑到实际应用,将沿界面高度方向粘性应力滑移系数的平均值作为腔体内的平均粘性应力滑移系数,并利用二次多项式对不同瑞利数工况下的三种粘性应力滑移系数分布分别进行了拟合,获得了基于瑞利数的拟合公式。将数值模拟与实验结果进行对比,使二者得到相互验证。利用实验获得的平均粘性应力滑移系数值代入到两区域模型的边界条件中,对数学模型进行求解。实验测试与数值模拟所得流线对比的结果表明,当瑞利数处于104～106数量级变化时,数值模拟与实验测试结果比较吻合,但在瑞利数变化至107时二者的差别较大。垂直速度分量比较的结果表明,数值模拟结果与实验测试结果比较温和。与数值模拟相比,实验测试所得多孔介质区域内的速度值偏小。其原因可归因于几个方面,其一,实验腔体的内壁面不够光滑,其二构成多孔介质的圆柱阻力与数值模拟多孔介质阻力不同,其三在界面处,圆柱对流动方向有较大影响,故可能导致所测速度值偏小。多孔介质区域与流体空间交界面处的滑移效应影响因素众多,包括流体物性、流动状况及界面结构参数等。本文首先分析了界面处的滑移效应,通过数值模拟探讨了滑移效应对腔体内自然对流传热的影响规律,然后利用实验测试了腔体内尤其是界面处的流动特性,求解出粘性应力滑移系数,并将实验测试与数值模拟结果进行了对比,为复合腔体内多孔介质与流体空间交界面处的滑移效应及其影响机理的研究奠定了基础。