大气表面层气-固颗粒两相流的大涡模拟

发布时间：2020-06-20 05:42

【摘要】：本文基于大涡模拟和点力模型研究湍流边界层中颗粒与湍流的相互作用,首先,本文提出了一种近壁平均速度梯度改进涡粘计算的简单的修正方法,并基于此方法模拟了Re_??3?1 0~7的高雷诺数大气边界层壁湍流,其中亚格子模型采用尺度依赖动力涡粘模式。模拟结果表明:所提出的修正方法能够有效改善湍流平均速度廓线的预测结果,并合理捕捉湍流相干结构的主要特征。同时采用湍流大涡模拟与颗粒点力模型数值模拟了可侵蚀颗粒床面上高雷诺数水平槽道湍流颗粒两相流动,以分析颗粒对湍流的调制作用机理。结果表明:可侵蚀颗粒床面上的颗粒分布在流向高、低速流体中间。平均三维颗粒浓度结构近壁大远壁小,呈hill-like状。除非常靠近近壁(y~((10))(27)3)外颗粒平均速度均小于湍流平均速度;整个计算高度内,颗粒流中的流体流向平均速度相比单相流减小。近壁颗粒导致流向速度脉动、雷诺应力相比单相流衰减,远离壁面增强,衰减与增强的转折点在约y~((10))?40,随雷诺数增大而降低。垂向和展向速度脉动在整个计算高度内均增强。颗粒流中,由速度相关确定的湍流结构的流向与展向尺度在近壁内区增大,外区湍流结构尺度减小。雷诺应力的象限分析显示,第二象限(Q2)和第四象限(Q4)的应力幅值和发生频率在近壁减弱远壁增强,第一象限(Q1)和第三象限(Q3)象限应力幅值近壁几乎不变,但频率增大,远壁幅值增强频率不变。条件分析揭示,跃移层颗粒对外层流体表现为活跃粗糙元,是导致内层流向速度脉动和雷诺应力减弱而外层流体速度脉动与雷诺应力增强的主要原因。此外,颗粒上升和下降过程中的条件平均流场表明,颗粒由于重力作用导致的下降过程是使湍流结构受影响的主要原因。另外,本文采用将颗粒拉格朗日轨迹模拟与湍流含壁模型大涡模拟(WMLES)相结合的方法,成功再现了区域为600m×12m×75m的近地表风沙流输沙的沙条带现象。以浓度相关系数提取条带的特征尺度发现随风速增大减小、随粒径增大而增大。在典型的0.2mm粒径下,浓度相关系数0.05为0.2mm时,沙条带的流向特征长度尺度为24~30m,宽度可达3m,在粒径为0.5mm时,流向特征长度尺度为30~90m,展向宽度可达6~12m,为与已有野外观测结果吻合。本文分析表明:平坦、均匀可侵蚀地表沙条带现象与风场存在湍流超大/大尺度结构有关,高浓度输沙出现在条件平均的超大/大尺度结构扫掠地表时的尾迹区,可以看作是外区湍流结构在近壁跃移层的足迹;同时,沙条带现象受风沙流中沙粒和地表击溅过程的影响,在不可侵蚀的表面上,沙条带虽然可以向下游延伸很长的距离,但会变窄。
【学位授予单位】：兰州大学
【学位级别】：博士
【学位授予年份】：2019
【分类号】：O359
【图文】：

(c) (d)图 1-2 自然界中的风沙运动，（a）近地表风沙流运动，图片来自于 Baas,08(Geomorphology)，（b）沙尘暴，图片为 2012.12.19 德克萨斯州的沙尘暴，（c） 风雪图片来自于百度百科，（d）水沙运动，图片来自于百度百科。以风沙两相流为例，在沙尘暴中，由于气体的作用，不同大小的颗粒在运程中采用不同的运动方式运动。Bagnold (1941)将颗粒运动分为三类;悬移、、蠕移，如图所示。

(a) (b)(c) (d)自然界中的风沙运动，（a）近地表风沙流运动，图片来自ology)，（b）沙尘暴，图片为 2012.12.19 德克萨斯州的沙尘暴片来自于百度百科，（d）水沙运动，图片来自于百度百科相流为例，在沙尘暴中，由于气体的作用，不同大小同的运动方式运动。Bagnold (1941)将颗粒运动分为图所示。

【相似文献】

相关期刊论文 前10条

1 王杨;周丹杰;王超;;基于隐式大涡模拟的高超声速热流计算研究[J];航空科学技术;2015年10期

2 何子干，倪汉根;大涡模拟法的二维形式[J];水动力学研究与进展(A辑);1994年01期

3 苏铭德;弯曲槽道内湍流运动的大涡模拟[J];力学学报;1989年05期

4 苏铭德;大涡模拟的代数模型及其应用[J];中国科学(A辑 数学 物理学 天文学 技术科学);1989年07期

5 刘有军;平面突扩流动非稳定性的大涡模拟[J];计算力学学报;1998年02期

6 张井龙;刘崇江;许红丽;王尊策;;旋转变径管内螺旋流相干结构的大涡模拟研究[J];地震工程与工程振动;2019年06期

7 黄振宇;缪国平;;大涡模拟在水下航行体周围黏性流场计算中的初步应用[J];水动力学研究与进展(A辑);2006年02期

8 苏铭德;直方管内充分发展湍流的大涡模拟第一部分[J];空气动力学学报;1995年01期

9 张井龙;王尊策;吴卓;徐艳;;液-液旋流分离器内流动特性的大涡模拟[J];化工机械;2013年03期

10 黄顺洲;台阶后紊流场的大涡模拟[J];燃气涡轮试验与研究;1998年03期

相关会议论文 前10条

1 叶丰;包芸;;热液羽流混合与输运的大涡模拟研究[A];中国力学大会-2015论文摘要集[C];2015年

2 吴霆;何国威;张星;;周期山状流的大涡模拟[A];第八届全国流体力学学术会议论文摘要集[C];2014年

3 张健;何国威;陆利蓬;;条件滤波大涡模拟中亚格子模型的数值研究[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集（下）[C];2007年

4 夏振华;史一蓬;肖左利;陈十一;;约束大涡模拟方法在分离流中的应用[A];北京力学会第18届学术年会论文集[C];2012年

5 史一蓬;肖左利;夏振华;王建春;江舟;洪仁楷;张日葵;陈十一;;约束大涡模拟方法的进展[A];中国力学大会——2013论文摘要集[C];2013年

6 崔桂香;张兆顺;;城市大气微环境的多尺度大涡模拟研究[A];中国力学大会——2013论文摘要集[C];2013年

7 马明;黄伟希;许春晓;崔桂香;;一种基于修正壁模型的大涡模拟-浸没边界混合方法[A];第十届全国流体力学学术会议论文摘要集[C];2018年

8 朱辉;符松;石磊;王志坚;;基于高精度通量重构格式的新型RANS/ILES方法研究[A];2016海峡两岸流体力学研讨会论文摘要集[C];2016年

9 王锁柱;李椿萱;雷娟棉;;基于空间模式的超声速边界层转捩的大涡模拟研究[A];中国力学大会——2013论文摘要集[C];2013年

10 龙云;徐默涵;程怀玉;季斌;龙新平;;水翼非定常空化大涡模拟的验证和确认分析[A];第二十九届全国水动力学研讨会论文集（上册）[C];2018年

相关博士学位论文 前10条

1 冯圣军;大气表面层气-固颗粒两相流的大涡模拟[D];兰州大学;2019年

2 李敬法;表面活性剂减阻流动的本构模型及大涡模拟研究[D];中国石油大学（北京）;2017年

3 于洲;基于详细化学建表方法的湍流气相燃烧大涡模拟研究[D];中国科学技术大学;2019年

4 李雪松;基于可压缩流方法的大涡模拟及其工程化应用[D];中国科学院研究生院（工程热物理研究所）;2006年

5 孙小波;可压缩剪切湍流的直接模拟和大涡模拟研究[D];中国科学技术大学;2007年

6 金晗辉;气粒两相湍流拟序结构的大涡模拟研究[D];浙江大学;2002年

7 许岩韦;结合扩展小火焰模型的甲烷及合成气射流火焰的大涡模拟与化学反应机理简化[D];浙江大学;2015年

8 李钊;RANS/LES嵌入耦合方法研究[D];清华大学;2015年

9 隋斌;明渠弯道交汇水流水力特性的大涡模拟与实验研究[D];武汉大学;2017年

10 王博;激波/湍流边界层相互作用流场组织结构研究[D];国防科学技术大学;2015年

相关硕士学位论文 前10条

1 孙可心;基于间断伽辽金法的隐式大涡模拟研究[D];南京航空航天大学;2019年

2 吕泽锋;基于大涡模拟的海豹胡须流场分析[D];哈尔滨工业大学;2018年

3 胡秋皓;制动线大涡模拟耦合方法下的潮流能水轮机数值模拟[D];上海交通大学;2018年

4 张春;双旋转入流湍流场流动特性的大涡模拟研究[D];江苏大学;2019年

5 陈柏宇;ASP连铸结晶器流场和温度场的模拟研究[D];辽宁科技大学;2018年

6 张风睿;弯曲激波压缩超—高超声速进气道流场快速预测与大涡模拟研究[D];南京航空航天大学;2018年

7 邓涛;钝体—平板角区复杂流动的大涡模拟研究[D];重庆大学;2018年

8 常牧;压水堆波动管热分层和温度震荡分析研究[D];华北电力大学(北京);2017年

9 陈武忠;基于非线性亚格子模型的Sandia火焰D和F的大涡模拟[D];华中科技大学;2018年

10 刘一秀;方槽内单颗粒存在时湍流场的大涡模拟[D];北京化工大学;2014年

本文编号：2721964