基于Fluent风沙环境下路基风蚀破坏模拟的研究.pdfVIP

第3l卷增刊1 温州大学学报·自然科学版 2010年11月 Vol 1 JournalofWenzhou Sciences Nov,2010 31，Sppl University·Natural 基于Fluent风沙环境下路基风蚀破坏模拟研究 黄浩，李驰，孙兵兵 (内蒙古工业大学土木工程学院，内蒙古呼和浩特010051) 摘要：通过FLUENT气固两相流计算，运用k-e湍流模型，对挟沙风作用下沙漠路基的风蚀破坏规 律进行数值模拟研究．分析路基不同横断面下风沙流扰动、增速、减速、恢复的过程，以及路基障碍 物的存在对风沙流扰动后造成的风速减弱区和恢复区，总结路基断面特征点的风速变化规律，与现场 实测结果作比较，二者具有很好的一致性．研究结果表明，不同路基横断面下风速减弱区和恢复区的 分布对路基高度变化敏感，受边坡坡率的影响较小，不同路基高度、不同边坡坡率下路基沿程风速变 化明显不同．边坡坡率一定，随路基高度的增加，路基周围流场扰动被增强，迎风坡坡顶风速增大显 著，背风坡坡底风速下降愈明显；路基高度一定时，边坡坡率越小，路基沿程风速变化越平缓．当边 坡坡率为l：3时，路基模型高度由1m增加到3m时，迎风坡坡顶风速增加12％，背风坡坡底风速降至 80％．建议为避免沙漠路基的风蚀破坏，路基高度不宜太大，边坡坡率较小为宜． 关键词；FLUENT两相流；挟沙风；沙漠路基；风蚀破坏 长期以来，沙漠公路风蚀破坏是公路交通部门亟待解决的关键问题。由于风沙土是沙漠地区 最丰富、最廉价的路基填料。其颗粒粒径大都小于0．50ram，粒间无粘聚力，抗风蚀能力弱，属 于最易蚀性颗粒【11。当风沙土路基置于“风大沙多”的沙漠风沙环境下时，路基在挟沙风作用下被 掏蚀、吹蚀和磨蚀、堆蚀将不可避免。国内外关于风沙物理学的研究已取得一些研究成果，主要 集中在风沙流的整体结构、风沙层内气流速度、颗粒平均速度、颗粒浓度、输沙通量等物理量的 分布，以及着眼于风沙流中单个颗粒的运动状态，通过分析它的受力情况研究颗粒起跳的机理和 颗粒在空中运动状态的变化规律等方面12-7]。 目前的研究手段有两种：一是实验观测，包括风洞实验和野外观测，这是风沙物理学传统的研 究方法，文献【8，9】采用多点遥测风速仪、高速摄像仪、风洞试验等方法对沙漠公路扰动下风沙流 整体结构以及风速变化进行研究。二是数值模拟，即通过建立风沙运动的数学模型对风沙流进行 求解，这类工作近10年来发展很快。成为该领域内一个活跃的研究方向。模拟内容主要有风沙 流、沙波纹、沙丘和粉尘运动等方面[9-tl】，文献【9．10】采用FLUENT单相流对风沙流的整体结构 进行研究，为课题的深入研究提供了很好的借鉴一 本文采用FLUENT6计算软件【lM31，研究气固两相流复杂条件下，路基周围风速流场的变化， 以及路基障碍物的存在对风沙流扰动后造成的风速减弱区和风速恢复区，总结路基风蚀破坏规 律，为进一步揭示路基风蚀破坏机理提供理论依据。 基金项目：国家自然科学基金项目，内蒙自然科学基金项目(2009BS0704) 作者简介：黄浩(1978．)，男，安徽人，硕士研究生，研究方向：沙漠环境岩土工程 黄浩等：基于Fluent风沙环境下路基风蚀破坏模拟研究 169 1流体动力学与FLUENT两相流计算模型 1．1计算原理 运用流体动力学原理，主要考虑流体的连续方程、动量方程、动量矩方程和能量方程。由于 动量矩方程和能量方程和本文的研究内容无关，在此仅列出微分形式的连续方程和动量方程，见 公式(1)和公式(2)。 连续方程为： 娑+挈：0 (1) 靠d， 动量方程为： D +心 +p =船 + 动 + + = + D A p 船 讹一西抛一西 却i孙|奄 和～匆如一砂 监缸监砂 监匆堡玉 1．2 n^脚两相流计算模型 当