CN112784504B - 一种强耦合固液多相流数值模拟方法 - Google Patents
一种强耦合固液多相流数值模拟方法 Download PDFInfo
- Publication number
- CN112784504B CN112784504B CN202110116957.2A CN202110116957A CN112784504B CN 112784504 B CN112784504 B CN 112784504B CN 202110116957 A CN202110116957 A CN 202110116957A CN 112784504 B CN112784504 B CN 112784504B
- Authority
- CN
- China
- Prior art keywords
- solid
- phase
- liquid
- model
- force
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 239000007788 liquid Substances 0.000 title claims abstract description 35
- 238000004088 simulation Methods 0.000 title claims abstract description 31
- 238000010168 coupling process Methods 0.000 title claims abstract description 15
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 15
- 239000007791 liquid phase Substances 0.000 claims abstract description 71
- 239000007790 solid phase Substances 0.000 claims abstract description 66
- 239000013598 vector Substances 0.000 claims abstract description 33
- 230000016507 interphase Effects 0.000 claims abstract description 18
- 239000012530 fluid Substances 0.000 claims abstract description 17
- 239000012071 phase Substances 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000009825 accumulation Methods 0.000 claims abstract description 5
- 239000002245 particle Substances 0.000 claims description 26
- 230000001133 acceleration Effects 0.000 claims description 9
- 230000005484 gravity Effects 0.000 claims description 8
- 239000007787 solid Substances 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 6
- 230000002706 hydrostatic effect Effects 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 4
- 239000013618 particulate matter Substances 0.000 claims description 3
- 239000013612 plasmid Substances 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims 1
- 239000000126 substance Substances 0.000 abstract description 17
- 230000008878 coupling Effects 0.000 abstract description 9
- 230000007246 mechanism Effects 0.000 abstract description 7
- 238000005191 phase separation Methods 0.000 abstract description 2
- 230000015556 catabolic process Effects 0.000 abstract 1
- 238000006731 degradation reaction Methods 0.000 abstract 1
- 230000008569 process Effects 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 7
- 230000002265 prevention Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 101100001671 Emericella variicolor andF gene Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种强耦合固液多相流数值模拟方法,包括步骤:通过物理力学分析确定固相与液相受力特征,并将固相与液相的受力特征写入数值模型;获取待测山地的地形数据与物源数据,并通过网格坐标转化(x,y,z,h)四维网格数据;定义待测山地的相间粘性拖曳力、虚拟质量力,获取待测山地的沟道启动、运动、堆积位置信息,构建守恒向量模型、固液相间作用力模型、浮力模型、粘性托曳力模型、虚拟质量力模型、多相流模型、液相底床摩阻力模型,进行强耦合固液多相流数值模拟。该方法适用于不同物质多相复杂物质的强耦合力学机制,可以消除传统数字模拟方法中因相间作用力耦合不足导致的椭圆退化、相脱离、流态紊乱、流体失真等问题。
Description
技术领域
本发明属于山地灾害防治、环境治理技术领域,具体涉及一种强耦合固液多相流数值模拟方法。
背景技术
固液多相流是由固相颗粒和液相流体组成的混合流体,固相颗粒和液相流体在运动中互相影响、共存并有明确分界面,其涉及流动力学、热力学、传热传质学、燃烧学、化学和生物反应以及相关工业过程中的共性科学问题,对多相流科学的发展与进步对生态与环境的变迁、山地灾害防治分析、环境治理等具有极为重要的意义。
现今对固液多相流数值模拟不能精确实现山地灾害过程中涉及的多种多相复杂物质流体,目前对于山地灾害多相流中各种物理机制不明确,数值模拟方法过于简化,如:
(1)、专利申请号为201710958968.9的文献中公开一种泥石流全过程数值模拟及数值计算方法,该方法基于Savage-Hunter理论的单相NS方程,将要模拟的多相流介质通过简化为单相流复杂流体介质进行模拟,且用一套浅水波方程定义;虽然在宏观上可以起到很好的模拟及预测作用,但由其控制方程简化,无法考虑其中一些很深入的固液耦合物理机制(相间作用力,固相及液相本构模型等)故此技术无法更深入的分析多相流的各种物理受力机制及动力学信息,只能用于快速的宏观数值模拟及反演之效用;且该方法中通过人为定义颗粒体积浓度达到定义介质内固体和液体的目的,但这种比例是固定不变的。
(2)、在专利申请号为201811522183.8的文献中公开一种固液多相适用于泥石流的动力学数值模拟方法,该方法首次将固液两相方程独立定义,并且施用适当且普适的本构方程定义固相及液相的流态及受力特征,然而由于计算的简便性对本构方程做了如下简化,假设相间作用力约为粘性拖曳力的作用,因为其他应力相比于粘性拖曳力较小;但该方法经研究发现很多多相固液密集颗粒流中虚拟质量力的作用尤为明显,且该方法中使用μ(I)本构关系作为固相颗粒介质的本构关系,但模型μ(I)只考虑了颗粒在惯性作用下的摩阻力及运动特征,忽略了多相流中液相流体物质对固相的影响作用,模拟可靠性降低。
发明内容
有鉴于此,本发明的目的之一在于提供一种强耦合固液多相流数值模拟方法,该方法适用于不同物质多相复杂物质的强耦合力学机制,实现多物质、多流态的固液耦合多相的动力学数值模拟。
为实现上述目的,本发明的技术方案为:一种强耦合固液多相流数值模拟方法,包括以下步骤:
通过物理力学分析确定固相与液相受力特征,并将所述固相与所述液相的受力特征写入数值模型;
获取待测山地的地形数据与物源数据,并通过网格坐标转化(x,y,z,h)四维网格数据;其中,(x,y)为所述四维网格数据内各点相对于固定参考点的相对经纬距离,z为各点上的地形高程值,h为物源高程值;
定义所述待测山地的相间粘性拖曳力、虚拟质量力,获取所述待测山地的沟道启动、运动、堆积位置信息,构建守恒向量模型、固液相间作用力模型、浮力模型、粘性托曳力模型、虚拟质量力模型、多相流模型、液相底床摩阻力模型,进行强耦合固液多相流数值模拟;其中,
所述守恒向量模型为:
其中,g为重力加速度,hs、hf分别为固液两相的高度,us、vs分别为固相在x、y方向上分量,uf、vf分别为液相在x、y方向上分量,Tsx、Tsy分别为固相相底床摩阻力在x、y方向上的分量,Tfx、Tfy分别为液相相底床摩阻力在x、y方向上的分量,fix、fiy分别为固液相间作用力在x、y方向上的分量,U、F(U)、G(U)、S均为向量,为向量U对时间t作偏导数,为向量F(U)对x作偏导数,为向量G(U)对y作偏导数;
所述固液相间作用力模型为:
fi=fb+fd+fVM;
其中,fb为流体静力学浮力,fd、fVM分别代表粘性托曳力、虚拟质量力;
所述浮力模型为:
所述粘性托曳力模型为:
hd=hs+hf;
其中,Cd为固液相间动量交换系数,为固相速度矢量,为液相速度矢量,hd为多相流流体高度,hs、hf分别为固液相物质高度,UT为单个颗粒在环境流中的下沉速度,有其中g为重力加速度,d为颗粒粒径大小,为密度比,分别为固相和液相体积分数,γ为液相和固相密度之比,P为无量纲固定值,Rep为颗粒尺度雷诺数,F(Rep)为G(Rep)为Rep的相关函数、代表其固相和液相分别对拖曳力系数的影响;
所述虚拟质量力模型为:
其中,CVM为虚拟质量力系数,t为时间;
所述多相流模型为:
Ts=ρsμ(K)ghs;
其中,Ts为固相摩阻力,μ1、μ2分别为固相颗物质粒在准静止、高速剪切状态下的内摩擦角,K0为无量纲常数,K为惯性粘性数,μ(K)为K的函数,Is为颗粒惯性数,J为颗粒粘性数,∈为无量纲参数,为深度平均剪切率,Ps为固相颗粒物质底床静压,ηf为液相粘度;
所述液相底床摩阻力模型为:
进一步地,所述数值模型的控制方程为基于非可压纳维叶-斯托克斯方程通过深度积分法简化得到的浅水波方程。
与现有技术相比,本发明具有如下优点:
(1)本发明在现有单相或多相控制方程数值方法的基础上进一步明确了固相液相复杂的受力特征及相互影响反应机制,可以将泥石流等复杂多相流体固液相互作用下的复杂运动模式模拟出来,将不同性质固液物质通过相间作用力相联系,用于消除因相间作用力耦合不足导致的椭圆退化、相脱离、流态紊乱、流体失真等问题,其具体形式为固-液作用、固-固作用、液体间相互作用,提高在泥石流防治,环境保护等领域的科学性、可靠性;
(2)过数值模拟,研究不同地域,不同浓度,不同地形下泥石流动力学特征,将空间地域分布引入到泥石流研究中,为我国特定地区泥石流灾害防治提供技术支持;
(3)本发明中涉及的数值模拟方法及数值计算具有简单、高效等特点,适合提供给相关灾害领域工作者使用,提高泥石流防治,环境保护的科学化水平;
(4)本发明具体可变高精度,适用于泥石流动力学特征等优点,针对泥石流动力学启动-运动-堆积等过程,能较好地将整个动力学过程完整的模拟出来,提高泥石流动力学过程研究的科学性和可靠性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍。显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1为本发明中实验槽切面示意图;
图2为本发明中一种强耦合固液多相流数值模拟方法一实施例实验结果示意图
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
所举实施例是为了更好地对本发明进行说明,但并不是本发明的内容仅局限于所举实施例。所以熟悉本领域的技术人员根据上述发明内容对实施方案进行非本质的改进和调整,仍属于本发明的保护范围。
需要注意的是,本实施例中的各代表式的下标如无特别提出解释,则仅仅用于区分,无特殊含义。
实施例1
本实施例中公开一种强耦合固液多相流数值模拟方法,具体包括以下步骤:
S1:通过物理力学分析确定固相与液相受力特征,并将固相与液相的受力特征写入数值模型;
本步骤中,先通过物理力学分析分别确定固相及液相受力特征,获得本数值方法预处理阶段所有需要准备的参数,并将其应力张量、体力等受力特征以数学形式带入数值模型中;
优选地,数值模型的控制方程为基于非可压纳维叶-斯托克斯方程通过深度积分法简化得到的浅水波方程,其特点在相比于传统的纳维叶-斯托克斯方程其通过深度方向积分平均使得方程组合理简化,计算效率明显提高,在满足其深度尺度小于长度尺度假设下可以较为精确地进行大流域重力场山地灾害过程的数值模拟,通过不同性质物质所反映的本构关系及固液相间强耦合作用力,可以精确地对相关灾害过程进行快速精确地数值模拟,不仅对于微观尺度上多相流的研究分析,或是宏观尺度上的灾害过程的正反演、防灾减灾都有更为明确的意义。
本发明中采用的质量守恒方程为:
进一步地,将上述质量守恒方程过深度积分简化后得到固液耦合多相强耦合数值方法的浅水波方程格式质量守恒方程,即:
其中,hs、hf分别为固、液相物质高度,us、vs分别为固相在x、y方向上分量,uf、vf分别为液相在x、y方向上分量,上述方程右边恒等于0代表其固液相物质总量始终保持一致,不会发生相变等特殊物理化学变化;
S2:获取待测山地的地形数据与物源数据,并通过网格坐标转化(x,y,z,h)四维网格数据;
本实施例中,通过野外科学考察和物理及力学实验获取山地灾害沟道启动、运动、堆积位置信息;再通过地理信息系统Arcgis、高精度3D地形扫描仪、无人机等方式确定计算区域内高精度的地形数据及物源数据;通过网格坐标转化表达为(x,y,z,h)四维网格数据,分别为地形数据网格内各点相对于固定参考点的相对经纬距离(x,y)、各点上的地形高程值z和物源高程值h;
S3:定义待测山地的相间粘性拖曳力、虚拟质量力,获取待测山地的沟道启动、运动、堆积位置信息,构建守恒向量模型、固液相间作用力模型、浮力模型、粘性托曳力模型、虚拟质量力模型、多相流模型、液相底床摩阻力模型,进行强耦合固液多相流数值模拟;
本步骤中,通过界定不同的山地灾害体分别定义不同的相间粘性拖曳力、虚拟质量力;
具体地,参考步骤S1中的浅水波方程格式质量守恒方程,本发明的一种强耦合固液多相流数值模拟方法采用的动量守恒方程为:
然后将以上动量守恒方程通过深度积分简化后得到本发明中方法需要的浅水波方程格式的动量守恒方程,即:
其中,分别为固相、液相体积分数,且满足τs=(τsx,τsy),τf=(τfx,τfy)分别为固、液相底床摩阻力矢量,fi为强耦合相间作用力,其形式为fi=fb+fd+fVM;θ为斜面倾角,τsx和τsy为固相底床摩擦力在x、y方向上的分量,τfx和τfy为液相底床摩擦力在x、y方向上的分量,fix、fiy分别为固液相间作用力在x、y方向上的分量,ρs为固相密度,ρf为液相密度;
基于以上推到,获得守恒向量模型为:
其中,g为重力加速度,hs、hf分别为固液两相的高度,us、vs分别为固相在x、y方向上分量,uf、vf分别为液相在x、y方向上分量,Tsx、Tsy分别为固相相底床摩阻力在x、y方向上的分量,Tfx、Tfy分别为液相相底床摩阻力在x、y方向上的分量,U、F(U)、G(U)、S均为向量,为向量U对时间t作偏导数,为向量F(U)对x作偏导数,为向量G(U)对y作偏导数;
固液相间作用力模型为:
fi=fb+fd+fVM;
其中,fb为流体静力学浮力,fd、fVM分别代表粘性托曳力(Viscous drag)、虚拟质量力(Virtual mass),粘性拖曳力由相间相对速度引起,虚拟质量力由相间相对加速度引起;
本实施例中的浮力模型考虑减小固相物质有效重度进行耦合,设为:
优选地,粘性托曳力模型为:
hd=hs+hf;
其中,Cd为固液相间动量交换系数,为固相速度矢量,为液相速度矢量,hd为多相流流体高度,hs、hf分别为固液相物质高度,UT为单个颗粒在环境流中的下沉速度,有其中g为重力加速度,d为颗粒粒径大小,分别为固相和液相体积分数,γ为液相和固相密度之比,且P为无量纲固定值,一般取值0.5,用于调整固相和液相影响因子的相对大小,Rep为颗粒尺度雷诺数,且F(Rep)与G(Rep)为Rep的相关函数、代表其固相和液相分别对拖曳力系数的影响,M(Rep)为雷诺数的一个函数,其取值范围为(2,4);
虚拟质量力模型为:
再进一步地,本发明中的固相摩阻力Ts=(Tsx,Tsy)采用新型多相流模型μ(K):
Ts=ρsμ(K)ghs;
其中,Ts为固相摩阻力,μ1、μ2分别为固相颗物质粒在准静止、高速剪切状态下的内摩擦角,K0为无量纲常数,K为惯性粘性数,μ(K)为K的函数,其值为0.246±0.008,Is为颗粒惯性数,J为颗粒粘性数,∈为无量纲参数,为深度平均剪切率,Ps为固相颗粒物质底床静压,ηf为液相粘度;
本实施例中的液相底床摩阻力模型Tf=(Tfx,Tyf)采用牛顿体模型:
实施例2
基于实施例1中的方法,本实施例中进行相关实验操作,具体为模拟实验槽内复杂多项流体的运动,所有参数将取自物理实验所得,用于比较数值实验及物理时间的结果,用于揭示复杂多相流流体的受力模式及运动机理。
具体地,实验槽设置为下段坡脚为5度,上段坡脚为20度,每次实验方量0.05立方米,实验槽最大承载方量为0.17立方米,其实验槽的模拟切面图可参考图1;对于固体及流体的性质,本实施例中对于固体采用来描述其摩阻力特性,而对于液相则采用牛顿流体,对于相见作用力,我们采用本发明定义的相关应力模型,其它参数选取参见表格1:
表格1数值实验参数选取表
将相应的参数赋值及模型建模后参考实施例1中的方法,计算总时间为5秒,采用可变时间步计算方式,得到如图2所述的数字结果对照图。
从图2中可知,本发明中提出的方法使用μ(K)得到的模拟深度线条图相比于测量深度(Measured height),其流体龙头模拟的更准确,而传统模型的模拟深度(库伦摩阻力)不能很好的模拟出多相流的流态甚至是出现分离等现象。
本发明提车一种适用于多种山地灾害过程(携沙洪水、泥石流、颗粒流、滑坡等)、涉及多物质、多流态的固液耦合多相的动力学数值模拟方法,其采用一种新型固液耦合多相流动力过程的物理数值模型,可以适用于不同物质多相复杂物质的强耦合力学机制,将不同性质固液物质通过相间作用力相联系,用于消除因相间作用力耦合不足导致的椭圆退化、相脱离、流态紊乱、流体失真等问题。本方法先通过物理力学分析分别确定固相及液相受力特征,获得本数值方法预处理阶段所有需要准备的参数,并将其应力张量、体力等受力特征以数学形式带入数值模型中。本发明在现有方法上提高了对多流态山地灾害过程的动力过程数值模拟研究的科学性和可靠性,提高了灾害防治的针对性、增强预防效果为泥石流减灾提供技术支撑。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。
Claims (2)
1.一种强耦合固液多相流数值模拟方法,其特征在于,包括以下步骤:
通过物理力学分析确定固相与液相受力特征,并将所述固相与所述液相的受力特征写入数值模型;
获取待测山地的地形数据与物源数据,并通过网格坐标转化(x,y,z,h)四维网格数据;其中,(x,y)为所述四维网格数据内各点相对于固定参考点的相对经纬距离,z为各点上的地形高程值,h为物源高程值;
定义所述待测山地的相间粘性拖曳力、虚拟质量力,获取所述待测山地的沟道启动、运动、堆积位置信息,构建守恒向量模型、固液相间作用力模型、浮力模型、粘性托曳力模型、虚拟质量力模型、多相流模型、液相底床摩阻力模型,进行强耦合固液多相流数值模拟;其中,
所述守恒向量模型为:
其中,g为重力加速度,hs、hf分别为固液两相的高度,us、vs
分别为固相在x、y方向上分量,uf、vf分别为液相在x、y方向上分量,Tsx、Tsy分别为固相相底床摩阻力在x、y方向上的分量,Tfx、Tfy分别为液相相底床摩阻力在x、y方向上的分量,fix、fiy分别为固液相间作用力在x、y方向上的分量,U、F(U)、G(U)、S均为向量,为向量U对时间t作偏导数,为向量F(U)对x作偏导数,为向量G(U)对y作偏导数;
所述固液相间作用力模型为:
fi=fb+fd+fVM;
其中,fb为流体静力学浮力,fd、fVM分别代表粘性托曳力、虚拟质量力;
所述浮力模型为:
所述粘性托曳力模型为:
hd=hs+hf;
其中,Cd为固液相间动量交换系数,为固相速度矢量,为液相速度矢量,hd为多相流流体高度,hs、hf分别为固、液相物质高度,UT为单个颗粒在环境流中的下沉速度,分别为固相、液相体积分数,γ为液相和固相密度之比,P为(0,1)间无量纲固定值,用于衡量液相及固相对拖曳力的影响大小,Rep为颗粒尺度雷诺数,F(Rep)、G(Rep)为Rep的相关函数、代表其固相和液相分别对拖曳力系数的影响,M(Rep)为雷诺数的一个函数;
所述虚拟质量力模型为:
其中,CVM为虚拟质量力系数,t为时间;
所述多相流模型为:
Ts=ρsμ(K)ghs;
其中,Ts为固相摩阻力,μ1、μ2分别为固相颗物质粒在准静止、高速剪切状态下的内摩擦角,K0为无量纲常数,K为惯性粘性数,μ(K)为K的函数,Is为颗粒惯性数,J为颗粒粘性数,∈为无量纲参数,为深度平均剪切率,Ps为固相颗粒物质底床静压,ηf为液相粘度;
所述液相底床摩阻力模型为:
2.根据权利要求1所述的方法,其特征在于,所述数值模型的控制方程为基于非可压纳维叶-斯托克斯方程通过深度积分法简化得到的浅水波方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110116957.2A CN112784504B (zh) | 2021-01-28 | 2021-01-28 | 一种强耦合固液多相流数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110116957.2A CN112784504B (zh) | 2021-01-28 | 2021-01-28 | 一种强耦合固液多相流数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112784504A CN112784504A (zh) | 2021-05-11 |
CN112784504B true CN112784504B (zh) | 2022-08-30 |
Family
ID=75759303
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110116957.2A Active CN112784504B (zh) | 2021-01-28 | 2021-01-28 | 一种强耦合固液多相流数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112784504B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378493B (zh) * | 2021-06-25 | 2022-09-06 | 中国科学技术大学 | 一种全质量守恒的多相流数值模拟方法 |
CN114239352B (zh) * | 2021-12-14 | 2024-04-30 | 西南交通大学 | 一种深度积分流体模型和块体系统的流固耦合方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819650B (zh) * | 2012-08-16 | 2015-04-29 | 同济大学 | 一种岩土材料流滑灾变的计算模拟方法 |
NO337063B1 (no) * | 2014-01-24 | 2016-01-11 | Ledaflow Tech Da | Fremgangsmåte for transient kvasi-tredimensjonal simulering av flerfasefluidstrømning i rørledninger |
CN108520549B (zh) * | 2018-04-09 | 2021-10-22 | 华北电力大学(保定) | 一种基于物质点法的多尺度泥石流现象模拟方法 |
CN109657322B (zh) * | 2018-12-13 | 2020-08-21 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种固液多相适用于泥石流的动力学数值模拟方法 |
CN112100567B (zh) * | 2020-09-11 | 2023-06-06 | 中国科学院、水利部成都山地灾害与环境研究所 | 稀性泥石流垂向流速分布、表面流速、平均流速测量方法 |
-
2021
- 2021-01-28 CN CN202110116957.2A patent/CN112784504B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112784504A (zh) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112784504B (zh) | 一种强耦合固液多相流数值模拟方法 | |
Kaitna et al. | Experimental study on rheologic behaviour of debris flow material | |
Pinder et al. | Essentials of multiphase flow and transport in porous media | |
Hill | The computer simulation of dispersed two-phase flow | |
Jalaal et al. | Homotopy perturbation method for motion of a spherical solid particle in plane couette fluid flow | |
Van den Akker | Lattice Boltzmann simulations for multi-scale chemical engineering | |
Goodman et al. | Numerical simulations of marine hydrothermal plumes for Europa and other icy worlds | |
Liu et al. | Dissipative particle dynamics simulation of pore‐scale multiphase fluid flow | |
Feder et al. | Physics of flow in porous media | |
Otomo et al. | Studies of accurate multi-component lattice Boltzmann models on benchmark cases required for engineering applications | |
CN109657322A (zh) | 一种固液多相适用于泥石流的动力学数值模拟方法 | |
Gao et al. | A two‐phase, pore‐scale reactive transport model for the kinetic interface‐sensitive tracer | |
Shi et al. | Effects of topological disorder in unsaturated granular media via a pore-scale lattice Boltzmann investigation | |
Chen et al. | DEM study of the influences of the geometric and operational factors on the mechanical responses of an underwater mixing process | |
Cheng et al. | Numerical simulation of phosphorus release from resuspended sediment | |
Lee et al. | Particle dynamics at fluid interfaces studied by the color gradient lattice Boltzmann method coupled with the smoothed profile method | |
Schulze | Development of an Application-Oriented Approach for Two-Phase Modelling in Hydraulic Engineering | |
Gerlach et al. | Rivulets of finite height | |
Baysinger et al. | Particle interactions with permeable drops in shear flow | |
Shokrollahzadeh | Terminal settling velocity of a sphere in a non-Newtonian fluid | |
Tai et al. | A Simplex Multi-Phase Approach for Modelling Debris Flows in Smoothed-Terrain-Following Coordinate System | |
Gharbi et al. | Analysis of morphological changes in rivers | |
Zawadzki et al. | 9 Modeling of two-phase flows in rotating packed beds | |
Kijanski et al. | Modelling of Non-Spherical Particles in Dilute Non-Colloidal Suspensions Using SPH | |
Sadeghi et al. | Simulation of wetting tendency of fluids with high density ratios using RK Lattice Boltzmann method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |