风吹雪数值模拟的两方程模型方法 Two-equation model for numerical simulation on snow drifting.pdfVIP

第30卷第5期 空气动力学学报 Voi．30。No．5 2012年10月 ACTAAERODYNAMICASINICA 0ct．，2012 0640—06 文章编号：0258—1825(2012)05 风吹雪数值模拟的两方程模型方法 周喱毅1，李雪峰h2，顾明1，张艳萍1 (1．同济大学土木工程防灾国家重点实验室，上海200092；2．上海市城市建设设计研究总院，上海200125) 摘要：采用两相流理论，根据雪颗粒的运动机理，提出了新的两方程数值模型方法，认为跃移层传输方程中的源 项与壁面摩擦速度、雪层阀值摩擦速度的相对大小有关，并引入待定系数来确定源项的缩放比例。接着基于数值 试验的方法开展参数分析，采用将不同参数条件下的跃移质量浓度与经验公式结果进行比较的方法来确定待定系 数。为了验证该方法的适用性，对开阔场地条件下的风致积雪运动进行了模拟，并与实测数据进行了比较。结果 显示，数值模拟结果与实测数据吻合较好，能够定量地重现风致积雪运动的重要特性。模拟的跃移层雪质量浓度 与实测结果比较一致；而悬移雪质量浓度偏高。总质量传输率的模拟结果落在实测值范围之内。 关键词：风致积雪运动；两方程模型；雪质量浓度；质量传输率 中图分类号：V211．3 文献标识码：A 0 引 言 风致积雪运动的基本过程按雪颗粒离开地面的 单方程方法考虑非充分发展流场中的雪流量。在国 程度分为蠕移、跃移及悬移运动。研究风致积雪运动 的方法大致分为四类：理论分析、实地观测、风洞试验 大跨屋盖表面的风致积雪重分布时采用了单方程模 和数值模拟。随着计算机技术的迅速发展，数值模拟 型。2010年李雪峰进一步考虑了雪颗粒的浮力、重 方法越来越受到了雪工程研究人员的关注。 力和拖拽力(气动阻力)对空气相的影响，对单方程模 F1uid 风致积雪运动的CFD(ComputationalDy— 型进行了改进，并模拟了立方体周围的积雪分布口]。 namics)数值模拟一般采用欧拉一欧拉方法，利用两相 流理论，在考察风致积雪运动的力学机理基础上，通 程模型研究了风导向板在控制雪堆积方面的应用。 过在空气相的N—S方程中增加雪相控制方程进行求 解计算。由于不同运动类型对雪颗粒迁移的贡献大 型共同点是，跃移控制方程中都没有考虑雪的扩散效 小有所不同，一般仅模拟跃移及悬移运动。根据雪相 应；不同之处是对于源项的处理提出了各自的观点。 控制方程的个数，模拟方法大致可以分为丽类：单方 总体来讲，风雪运动的CFD数值模拟方法研究 程模型方法和两方程模型方法。单方程模型方法中 在国际上还属于发展阶段，各种数值模型也不尽相 对雪相的控制方程没有区分跃移运动和悬移运动，控 同。尽管前人的研究提供了开阔场地条件下(无障碍 制方程只有一个；两方程模型则考虑到雪跃移和悬移 物)跃移层、悬移层的实测结果，但将CFD数值模拟 运动的力学机理不同，分别对跃移和悬移运动建立控 与现场实测两者结合起来进行比较分析的文献很少。 制方程。在单方程模型方法方面，1994年I瓜ton建为此，本文提出了新的两方程模型方法，认为跃移层 立了风致积雪运动的数值模型模拟开阔地貌下风致 传输方程中的源顼与壁面摩擦速度、雪层阀值摩擦速 积雪运动，跃移运动采用了Pomeroy的经验公式口]，度的相对大小有关，并引入待定系数来确定源项的缩 悬移运动的模拟则采用了基于质量守恒的控制方 放比例。文中基于数值试验的方法开展参数分析，采 02—24 *收稿日期：2011一0921；修订日期：2012