非达西渗流模型模拟青藏铁路抛石路基冻土层温度.pdfVIP

非达西渗流模型模拟青藏铁路抛石路基冻土层温度牢 侯彬彬1，2卢德唐1·2陈永辉1，2孔祥言1 (1．中国科学技术大学力学系渗流实验室‘n，合肥230027； 2．安徽省计算与通讯软件重点实验室‘孙，合肥230027) 摘要：青藏铁路穿越冻土带，路基下的冻土融化会威胁铁路运行的安全。铁路没计中一般单纯靠增加热阻的消 极保护冻土的思路(增加路堤高度，使用保温材料)，本文采用碎石冷却路基来保护冻土的方法，用多孔介质中的 非达西自然对流来描述碎石中的空气流动，并给出了数值计算和分析结果。 关键词：冻土冷却路基渗流有限差分温度场 多孔介质里面的自然对流在近些年来吸引了很多研究人员的重视，因为里面包括了复杂物质和热的输运。这 拟的结果，揭示被动冷却效果。但是，以上的工作都是用达西流的方程计算的。这个在流速非常低的情况下 是适应的，对于青藏铁路来说，就不是很适应了。下面的理论和计算分析会揭示这一点。本文采用非达西流 的假设进行数值模拟计算。 1数学和物理模型 1．1 雷诺数和渗流瑞利数分析 渗流雷诺数尺P定义为 Re=—pV—d (1) 舡 其中矿是渗流速度(或称比流量)，≯是多孔介质的孑L隙度，d是同体颗粒的直径，将有关数据代入可得，对于 块石区： 尺P：筚掣掣~90 (o．4)(1．75×10。5) 10，流动为湍流，就需要考虑惯性一湍流效应，特别是湍流效应。 2．2渗流瑞利(Rayleigh)数分析 渗流瑞利数尺口定义为 R口：—pflgK—HAT：—pflglt—3AT．姜 (2) ／2Czt ∥口H‘ 式中：∥为空气的热膨胀系数；q2‰)，为总的有效热扩散系数；日为路堤高譬；K为多孑L介质的渗透 率；△丁为上下温差。 等温和定压边界，临界瑞利数一般为12．20。就是说，当下边界温度高于上边界，且尺口大于这个数值时就 会产生对流。 +基金项目：973项目(G1999032805)，中国科学院知识创新重大项目(kzcxl一SW-04)联合资助 第一作者简介：侯彬彬，1982年生，男，在读研究生，主要研究方向：渗流力学计算。 ·495· 从上面的分析可以知道，本文采用的是非达西流的假设，用的方程是非达西流的动量方程。对于块石区， 流动是湍流，其非达西渗流的动量方程可写成 彪詈+即=一警矿一加m+见g[1一∥(丁一ro)l (3) 肚‰。 究表明与G边界壁面、流动区域的几何条件等因素有关。Beavers等(1973)给出： ，(4) Cr=0．55(1_5．s吾) 。， 和高度。我们近似取C≈0．4。所涉及的材料的物理参数见表1，计算区域图如图1所示。 表1材料物理参数 比热 夕因子B 材料 密度P(kg／m3) 热导率(W／m·K) 孔隙度中 渗透率g(m2) Cp(J／(kg‘K)) (L／m) 空气 0．64l l 004 O．02 抛石(多孔介质)