- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
研究流体的运动目前主要采用欧拉方法和拉格朗日方法。拉格朗日方法是从分析流体各粒子出发,分析粒子的密度、压强、速度等参数随时间的变化,从而研究流体的运动。
这里的SPH是拉格朗日描述下的一种无网格方法,本质上是一种插值法,它使用加权平均去近似各种物理性质及其导数。
SPH最早用于天文物理学,最近其在流体上的运用也大放异彩。这次研究中SPH模型建立最关键的便是正确模拟不同性质的流体和固体的自由面和交界面。
水和沙都可看作是弱可压缩流体,这里用的方法便是WCSPH(weakcompressiblesmoothedparticlehydrodynamics)。
水流的运动很简单,对泥沙运动的模拟是关键。这里首先定义泥沙颗粒间的强度阈值,当泥沙运动时水流的剪切力没有达到此阈值,此时只考虑泥沙颗粒间的剪切应力(刚性运动),而当水流剪切力超过此阈值,此时可看作泥沙按照水流的运动方程来运动的(拟流运动),此强度阈值可由宾汉-非牛顿本构模型得到。
研究者做了两个模型,一个是基于实验数据做了两相溃坝的模拟,另一个是基于实验数据做了土石坝漫顶的模拟。
WCSPH数值模拟的步骤和其他数值模拟方法类似。首先确定初始条件和边界条件。初始条件就是在初始时刻运动应该满足的初始状态,边界条件是指在求解区域边界上所求解的变量或其导数随时间和地点的变化规律。
然后列出流体运动的控制方程,并对这一系列的偏微分方程进行离散处理,不同于其它网格方法,这里采用SPH粒子近似,详细的初始条件、边界条件和控制方程可看参考文献[3]。接下来我们一起来看看这两种模型的模拟结果。
两相溃坝模型
第一个模型是两相(液-液)流体运动的模拟,初始几何参数如下图:
该模型流体粒子有21068个,粒子间初始间距为0.002m,时间步长0.0002s。上图中蓝色的是水,密度1000kg/m3,动力粘度为10-6m2/s,下游棕色液体为PEO(氧化聚乙烯),密度与水相似,动力粘度为0.935×10-6m2/s,两者均为牛顿流体。
实验时以1.5m/s的速度完全打开闸门,模拟结果如下图所示:
图中黑色的线为实验数据,可以看出WCSPH拟合效果和精度是非常好的,有一些小的差异也是因为实验中闸门的粗糙度和渠道宽度方向上的影响(实验比模拟多了一个维度)。渠道内的静止流体往往会阻碍流动,与运动的流体前缘碰撞产生向上推力形成蘑菇形状流体。
计算结果表明,用WCSPH进行多相流体的模拟准确度是足够的,所以接下来,研究者将此方法运用到土石坝漫顶溃决的研究中去。
土石坝漫顶模型
第二个模型采用的数据如下:土石坝剖面为梯形,高0.2m,宽0.2m,顶长为0.1m,上游/下游面坡度1:2,底宽0.9m,流量为定值11.31L/s。该坝建筑材料为无粘性的沙,粒径2mm,密度2650kg/m3,摩擦角37°。
该模拟假定坝前水库初始没有蓄水,而在几秒钟就充满水。初始粒子间距0.005m,总共采用了21000个粒子。模拟结果如下图所示:
从图中可以明显看出,在破坏过程的开始阶段,顶部削减是主导机制,然后逐渐形成破裂并开始扩张。当水流通过下游边坡时,形成了流场中的剪应力,当这种剪应力超过河床内的阻力时,侵蚀过程就开始了。
在洪水过顶初期,由于溢流流量小,下游冲刷带走的泥沙量较小,所以产生较小的侵蚀。随后,随着流量和流速的增加,形成了高应力导致了侵蚀速率的增加。溢流侵蚀土石坝的下游面,最终溢流面和下游面趋于平行。
下图给出了模拟值和实验数据的对比关系。模型与实验结果吻合较好,这同样也表明WCSPH多相模型成功地再现了泥沙的运动状态。
研究者还分析了WCSPH模型分辨率、泥沙的流变性质。
下图给出了模型分辨率的影响,粒子间距分别为0.01、0.005、0.0025m。
由图可以得出模型分辨率越高,即粒子间距越小,模拟结果越准确。
为研究泥沙流变性质对结果的影响,此处设三种情况:case1为稠度系数为0.6pa﹒s,摩擦角为37°;case2为稠度系数为0.6pa﹒s,摩擦角为0°;case3为稠度系数为0.06pa﹒s,摩擦角为35°。
由图可以得出若泥沙摩擦角为0°,泥沙运动就趋近于牛顿流体,溃决过程就会越快,而减小泥沙的稠度系数,土石坝的侵蚀面积也会随之减少。
二维WCSPH多相模型在非粘性土石坝因漫顶造成的破坏侵蚀的模拟中精确度是很高的。数值模拟结果与实验数据吻合较好,并且再现了流动的主要特征,包括湍流、混合过程、界面破碎和冲刷。
您可能关注的文档
- 道路医生的无损检测就不是X光那么简单了 - 副本.docx
- 用“膜”法向海要水 - 副本.docx
- 轨道交通行业网络空间安全现状与未来发展 - 副本.docx
- 工业控制系统安全防护技术发展研究 - 副本.docx
- 量子网络系统研究进展与关键技术.docx
- 新能源应用安全风险防控战略研究.docx
- 煤炭工业数字化发展战略研究.docx
- 航运安全保障智慧化管控发展研究.docx
- 欧洲空中交通管理的过去、现在与未来.docx
- 多式联运联接共性关键技术体系构建.docx
- 原电池电动势的测定实验报告.pdf
- 与业主、设计、总包、监理和他承包人的配合措施.pdf
- 公司管理流程.pptx
- 2024_2025学年新教材高中地理第1章地球的运动素养综合训练新人教版选择性必修1.doc
- 2024_2025学年新教材高中地理第3章大气的运动第1节常见天气系统第1课时锋与天气分层作业新人教版选择性必修1.doc
- 2024_2025学年新教材高中地理第1章地球的运动第2节地球运动的地理意义第4课时正午太阳高度的变化四季更替和五带划分分层作业课件新人教版选择性必修1.pptx
- 2024_2025学年新教材高中地理第2章地表形态的塑造第2节构造地貌的形成第1课时地质构造与地貌课件新人教版选择性必修1.pptx
- 2024_2025学年新教材高中地理第1章地球的运动问题研究人类是否需要人造月亮课件新人教版选择性必修1.pptx
- 五片小雪花课件.pdf
- 2024_2025学年新教材高中地理第3章大气的运动第2节气压带和风带第1课时气压带和风带的形成分层作业课件新人教版选择性必修1.pptx
文档评论(0)