CN110245385B - 一种考虑气液固三相耦合的曝气法离散元数值模拟方法 - Google Patents
一种考虑气液固三相耦合的曝气法离散元数值模拟方法 Download PDFInfo
- Publication number
- CN110245385B CN110245385B CN201910410513.2A CN201910410513A CN110245385B CN 110245385 B CN110245385 B CN 110245385B CN 201910410513 A CN201910410513 A CN 201910410513A CN 110245385 B CN110245385 B CN 110245385B
- Authority
- CN
- China
- Prior art keywords
- particles
- soil
- bubble
- particle
- liquid
- 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
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
-
- 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
- Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
- Y02W10/00—Technologies for wastewater treatment
- Y02W10/10—Biological treatment of water, waste water, or sewage
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑气液固三相耦合的曝气法离散元数值模拟方法。通过在颗粒介质模拟的土体模型中,实时注入气泡颗粒,同时以施加力的形式将液体以浮力、拖拽力的形式作用在浸泡于液体中的土体颗粒及气泡颗粒上,实现了固液气三相耦合,从而达到模拟曝气过程中气泡颗粒运动过程的目的。采用本发明技术方案,可以有效研究曝气法的运移规律以及曝气过程中的作用范围,为进一步从微观角度以统计角度研究分析曝气法中气泡颗粒运行的模式及机理提供了可能。
Description
技术领域
本发明属于岩土学科数值模拟技术领域,尤其涉及一种考虑气液固三相耦合的曝气法离散元数值模拟方法。
背景技术
随着社会经济的发展以及人类城市规模的扩大,污染土壤的修复受到越来越多地关注。污染场地修复不仅利于环保,同时可以扩大土地资源的利用率,符合国家可持续发展的战略。土壤修复技术大体分为原位修复与非原位修复两类。目前我国大力发展污染土壤原位修复技术。原位修复技术可以直接对受污染土体进行修复,不需要对土体开挖,从而大大减少了工程量,节约了修复费用的开支。其中,属于原位修复技术的曝气法因其可以高效地去除污染土中的有机污染物而得到大力的发展。该方法通过在土体深处注射气体,通过气体向地面运动进而挥发的过程中带出土体中的有机污染物,从而达到治理有机污染土的目的。曝气法因其成本低,修复时间短,成效明显,因此开发前景远大。
为了研究和推广曝气法,进一步提高曝气法的效率,科研人员发现仅仅从现场或室内试验的基础上无法满足对曝气法机理进行进一步的需求,特别是通过试验手段无法追踪气泡在土体中的运移路径及曝气挥发面上气体的挥发效果。因此,数值模拟是另一种研究曝气法机理的有效手段。
发明内容
发明目的:针对以上问题,本发明公开了一种考虑气液固三相耦合的曝气法离散元数值模拟方法。通过在颗粒介质模拟的土体模型中,实时注入气泡颗粒,同时以施加力的形式将液体以浮力、拖拽力的形式作用在浸泡于液体中的土体颗粒及气泡颗粒上,实现了固液气三相耦合,从而达到模拟曝气过程中气泡颗粒运动过程的目的。采用本发明技术方案,可以有效研究曝气法的运移规律以及曝气过程中的影响范围,为进一步从微观角度与统计角度研究分析曝气法中气泡颗粒运行的模式及机理提供了可能。
技术方案:为实现本发明的目的,本发明所采用的技术方案是:一种考虑气液固三相耦合的曝气法离散元数值模拟方法,包括如下步骤:
步骤A:根据曝气法应用场景尺寸建立相应的长宽高分别为x,y,z的立方体边界,在该边界内生成互不接触的土体颗粒,直到土体颗粒填满边界内部空间,然后进入步骤B;
步骤B:在已生成的土体颗粒上施加重力,土体颗粒自由落体在边界底部堆积,通过采用振动或压实的方式对土体颗粒进行处理,直到土体颗粒系统稳定,即土体颗粒系统的最大不平衡力小于设定值,得到满足预设孔隙率的土体试样,然后进入步骤C;
步骤C:根据岩土体离散材料属性,针对试样中土体颗粒与土体颗粒之间、土体颗粒与边界内壁之间的各个接触位置,分别设置接触模型,接触模型包括线性接触模型、Hertz模型、平行粘结模型,然后进入步骤D;
步骤D:根据实际场景设立液面高度,液面以上为空气,液面以下为液体;在模拟过程中,所有液面以下的土体颗粒都被施加浮力与拖拽力,液面以上的土体颗粒与气泡颗粒不受浮力与拖拽力的影响,将液面设置在边界顶部;设定模拟总时长T与单位时间步长Δt;将初始时长设为ti=0,然后进入步骤E;
步骤E:设置气体喷嘴的中心坐标及尺寸大小,喷嘴默认位于边界底部区域正中间,进入曝气循环计算,直到当前时刻的运行时长达到模拟总时长,从而实现曝气法整个进行过程的数值模拟。
进一步,步骤A所述生成土体颗粒互不接触,方法如下:
步骤A-1:在边界内随机生成一个初始土体颗粒,即随机生成土体颗粒的中心坐标及半径r,判断该中心点坐标是否位于边界之内,若不在边界之内则删除当前土体颗粒并生成新的土体颗粒,直到初始土体颗粒位于边界之内;
步骤A-2:在边界内随机生成新的土体颗粒,即随机生成土体颗粒的中心坐标及半径r,同时对此中心点坐标根据判定条件进行判定;判定条件一,中心点坐标应位于边界之内,若不满足则删除当前土体颗粒并生成新的土体颗粒;判定条件二,若该中心点坐标与相邻土体颗粒的中心点坐标距离小于两土体颗粒的半径之和,则可判定两土体颗粒接触,并删除当前土体颗粒,同时重新生成一个新的土体颗粒中心点坐标,直到生成的土体颗粒与相邻土体颗粒不接触为止;
步骤A-3:重复步骤A-2,直至土体颗粒充满整个边界空间。
进一步,步骤D中,施加在每个土体颗粒或气泡颗粒上的浮力计算公式为:
Fb=ρwVg
其中,Fb为浮力,ρw为液体密度,V为浸泡在液体中的颗粒体积,g为重力;
颗粒为球形颗粒时,若颗粒完全浸泡在液体中,浸泡体积V计算公式如下:
其中,r为球形颗粒半径;若颗粒没有完全浸泡在液体中,则分为两种情况:
情况一:液面hz在球形颗粒形心bz坐标上方,则浸泡部分体积V计算公式如下:
Δh=r-(bz-hz)
情况二:液面hz在球形颗粒形心bz坐标下方,则浸泡部分体积V计算公式如下:
Δh=r-(hz-bz)
其中,Δh为球面的上水平切面至液面的距离与下水平切面至液面的距离之中相对较小的距离,r为球形颗粒半径;
施加在每个土体颗粒或气泡颗粒上的拖拽力计算公式为:
进一步,步骤E所述曝气循环计算,步骤如下:
步骤E-1:在喷嘴区域内生成气泡颗粒,气泡颗粒受到浮力,拖拽力及重力的作用,三力合力向上,气泡颗粒向上运动;同时,搜索所有与液面接触的气泡颗粒,即气泡颗粒中心坐标至液面的距离,若该距离小于该气泡颗粒半径,删除该气泡颗粒,同时更新气泡颗粒总数,并进入步骤E-2;
步骤E-2:在当前时间步ti计算每个气泡颗粒受到的实时浮力与拖拽力,以及气泡颗粒与土体颗粒的接触力,并求得每个气泡颗粒所受的合力,根据牛顿第二定律,计算气泡颗粒的速度、加速度与位移,并进入步骤E-3;
步骤E-3:更新计算时长,即当前时刻的运行时长ti等于上一时刻的运行时长ti-1加上单位时间步长Δt,并进入步骤E-4;
步骤E-4:判断当前时刻的运行时长ti是否达到模拟总时长T,若当前的运行时长ti小于模拟总时长T,则返回步骤E-1;若当前的运行时长ti等于或大于模拟总时长T,则模拟终止。
有益效果:与现有技术相比,本发明的技术方案具有以下有益的技术效果:实现了采用数值模拟手段模拟曝气法的目标。目前曝气法仍停留在实验研究阶段,没有文献或专利发表成熟的曝气法的数值模拟的技术手段。采用本发明技术方案,考虑气液固三相耦合,实现了曝气法的数值模拟,从微观角度及统计学的角度对进一步研究分析曝气法中气泡运行的模式及机理提供了可能。
附图说明
图1是本发明的程序流程图;
图2是本发明数值模拟方法获取的曝气过程中气泡颗粒运行示意图;
图3是本发明数值模拟方法获取的气泡颗粒在土体颗粒中运移整体效果图及曝气法扩散角。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步的说明。
本发明所述的一种考虑气液固三相耦合的曝气法离散元数值模拟方法,图1展示了程序流程图。
包括如下步骤:
步骤A:根据曝气法应用场景尺寸建立相应的长宽高分别为x,y,z的立方体边界,在该边界内生成互不接触的土体颗粒,直到土体颗粒填满边界内部空间,然后进入步骤B;
步骤B:在已生成的土体颗粒上施加重力,土体颗粒自由落体在边界底部堆积,通过采用振动或压实的方式对土体颗粒进行处理,直到土体颗粒系统稳定,即土体颗粒系统的最大不平衡力小于设定值,得到满足预设孔隙率的土体试样,然后进入步骤C;
步骤C:根据岩土体离散材料属性,针对试样中土体颗粒与土体颗粒之间、土体颗粒与边界内壁之间的各个接触位置,分别设置接触模型,接触模型包括线性接触模型、Hertz模型、平行粘结模型,然后进入步骤D;
步骤D:根据实际场景设立液面高度,液面以上为空气,液面以下为液体;在模拟过程中,所有液面以下的土体颗粒都被施加浮力与拖拽力,液面以上的土体颗粒与气泡颗粒不受浮力与拖拽力的影响,本发明将液面设置在边界顶部;设定模拟总时长T与单位时间步长Δt;假设模拟总时长为1s,单位时间步长10-7s,则曝气循环将运行107次;将初始时长设为ti=0,然后进入步骤E;
步骤E:设置气体喷嘴的中心坐标及尺寸大小,喷嘴默认位于边界底部区域正中间,进入曝气循环计算,直到当前时刻的运行时长达到模拟总时长,从而实现曝气法整个进行过程的数值模拟。
步骤A所述生成土体颗粒互不接触,方法如下:
步骤A-1:在边界内随机生成一个初始土体颗粒,即随机生成土体颗粒的中心坐标及半径r,判断该中心点坐标是否位于边界之内,若不在边界之内则删除当前土体颗粒并生成新的土体颗粒,直到初始土体颗粒位于边界之内;
步骤A-2:在边界内随机生成新的土体颗粒,即随机生成土体颗粒的中心坐标及半径r,同时对此中心点坐标根据判定条件进行判定;判定条件一,中心点坐标应位于边界之内,若不满足则删除当前土体颗粒并生成新的土体颗粒;判定条件二,若该中心点坐标与相邻土体颗粒的中心点坐标距离小于两土体颗粒的半径之和,则可判定两土体颗粒接触,并删除当前土体颗粒,同时重新生成一个新的土体颗粒中心点坐标,直到生成的土体颗粒与相邻土体颗粒不接触为止;
步骤A-3:重复步骤A-2,直至土体颗粒充满整个边界空间。
步骤D中,施加在每个土体颗粒或气泡颗粒上的浮力计算公式为:
Fb=ρwVg
其中,Fb为浮力,ρw为液体密度,V为浸泡在液体中的颗粒体积,g为重力;
颗粒为球形颗粒时,若颗粒完全浸泡在液体中,浸泡体积V计算公式如下:
其中,r为球形颗粒半径;若颗粒没有完全浸泡在液体中,则分为两种情况:
情况一:液面hz在球形颗粒形心bz坐标上方,则浸泡部分体积V计算公式如下:
Δh=r-(bz-hz)
情况二:液面hz在球形颗粒形心bz坐标下方,则浸泡部分体积V计算公式如下:
Δh=r-(hz-bz)
其中,Δh为球面的上水平切面至液面的距离与下水平切面至液面的距离之中相对较小的距离,r为球形颗粒半径;
施加在每个土体颗粒或气泡颗粒上的拖拽力计算公式为:
步骤E所述曝气循环计算,步骤如下:
步骤E-1:在喷嘴区域内生成气泡颗粒,气泡颗粒受到浮力,拖拽力及重力的作用,三力合力向上,气泡颗粒向上运动;同时,搜索所有与液面接触的气泡颗粒,即气泡颗粒中心坐标至液面的距离,若该距离小于该气泡颗粒半径,删除该气泡颗粒,同时更新气泡颗粒总数,并进入步骤E-2;
步骤E-2:在当前时间步ti计算每个气泡颗粒受到的实时浮力与拖拽力,以及气泡颗粒与土体颗粒的接触力,并求得每个气泡颗粒所受的合力,根据牛顿第二定律,计算气泡颗粒的速度、加速度与位移,并进入步骤E-3;
步骤E-3:更新计算时长,即当前时刻的运行时长ti等于上一时刻的运行时长ti-1加上单位时间步长Δt,并进入步骤E-4;
步骤E-4:判断当前时刻的运行时长ti是否达到模拟总时长T,若当前的运行时长ti小于模拟总时长T,则返回步骤E-1;若当前的运行时长ti等于或大于模拟总时长T,则模拟终止。
图2展示了本发明数值模拟方法模拟取得的曝气进程中气泡颗粒运动过程;其中,图2的(a)、(b)、(c)、(d)分别展示了随着运行时间的进行,气泡颗粒逐渐扩散的过程。
图3是本发明数值模拟方法获取的气泡颗粒在土体颗粒中运移整体效果图及曝气法扩散角;其中,图3的(a)展示了气泡颗粒在土体中的扩散状态,图3的(b)展示了气泡颗粒的运行状态,以及气泡颗粒扩散角。
Claims (1)
1.一种考虑气液固三相耦合的曝气法离散元数值模拟方法,其特征在于:该方法包括如下步骤:
步骤A:根据曝气法应用场景尺寸建立相应的长宽高分别为x,y,z的立方体边界,在该边界内生成互不接触的土体颗粒,直到土体颗粒填满边界内部空间,然后进入步骤B;所述生成互不接触的土体颗粒,方法如下:
步骤A-1:在边界内随机生成一个初始土体颗粒,即随机生成土体颗粒的中心点坐标及半径r,判断该中心点坐标是否位于边界之内,若不在边界之内则删除当前土体颗粒并生成新的土体颗粒,直到初始土体颗粒位于边界之内;
步骤A-2:在边界内随机生成新的土体颗粒,即随机生成土体颗粒的中心点坐标及半径r,同时对此中心点坐标根据判定条件进行判定;判定条件一,中心点坐标应位于边界之内,若不满足则删除当前土体颗粒并生成新的土体颗粒;判定条件二,若该中心点坐标与相邻土体颗粒的中心点坐标距离小于两土体颗粒的半径之和,则可判定两土体颗粒接触,并删除当前土体颗粒,同时重新生成一个新的土体颗粒中心点坐标,直到生成的土体颗粒与相邻土体颗粒不接触为止;
步骤A-3:重复步骤A-2,直至土体颗粒充满整个边界空间;
步骤B:在已生成的土体颗粒上施加重力,土体颗粒自由落体在边界底部堆积,通过采用振动或压实的方式对土体颗粒进行处理,直到土体颗粒系统稳定,即土体颗粒系统的最大不平衡力小于设定值,得到满足预设孔隙率的土体试样,然后进入步骤C;
步骤C:根据岩土体离散材料属性,针对试样中土体颗粒与土体颗粒之间、土体颗粒与边界内壁之间的各个接触位置,分别设置接触模型,接触模型包括线性接触模型、Hertz模型、平行粘结模型,然后进入步骤D;
步骤D:根据实际场景设立液面高度,液面以上为空气,液面以下为液体;在模拟过程中,所有液面以下的土体颗粒都被施加浮力与拖拽力,液面以上的土体颗粒与气泡颗粒不受浮力与拖拽力的影响,将液面设置在边界顶部;设定模拟总时长T与单位时间步长Δt,将初始时长设为ti=0,然后进入步骤E;
其中,施加在每个土体颗粒或气泡颗粒上的浮力计算公式为:
Fb=ρwVg
其中,Fb为浮力,ρw为液体密度,V为浸泡在液体中的颗粒体积,g为重力;
颗粒为球形颗粒时,若颗粒完全浸泡在液体中,浸泡体积V计算公式如下:
其中,r为球形颗粒半径;若颗粒没有完全浸泡在液体中,则分为两种情况:
情况一:液面hz在球形颗粒形心bz坐标上方,则浸泡部分体积V计算公式如下:
Δh=r-(bz-hz)
情况二:液面hz在球形颗粒形心bz坐标下方,则浸泡部分体积V计算公式如下:
Δh=r-(hz-bz)
其中,Δh为球面上水平切面至液面的距离与下水平切面至液面的距离之中相对较小的距离,r为球形颗粒半径;
施加在每个土体颗粒或气泡颗粒上的拖拽力计算公式为:
步骤E:设置气体喷嘴的中心坐标及尺寸大小,喷嘴默认位于边界底部区域正中间,进入曝气循环计算,直到当前时刻的运行时长达到模拟总时长,从而实现曝气法整个进行过程的数值模拟;
所述曝气循环计算,步骤如下:
步骤E-1:在喷嘴区域内生成气泡颗粒,气泡颗粒受到浮力,拖拽力及重力的作用,三力合力向上,气泡颗粒向上运动;同时,搜索所有与液面接触的气泡颗粒,即气泡颗粒中心坐标至液面的距离,若该距离小于该气泡颗粒半径,删除该气泡颗粒,同时更新气泡颗粒总数,并进入步骤E-2;
步骤E-2:在当前时长ti内计算每个气泡颗粒受到的实时浮力与拖拽力,以及气泡颗粒与土体颗粒的接触力,并求得每个气泡颗粒所受的合力,根据牛顿第二定律,计算气泡颗粒的速度、加速度与位移,并进入步骤E-3;
步骤E-3:更新计算时长,即当前时刻的运行时长ti等于上一时刻的运行时长ti-1加上单位时间步长Δt,并进入步骤E-4;
步骤E-4:判断当前时刻的运行时长ti是否达到模拟总时长T,若当前的运行时长ti小于模拟总时长T,则返回步骤E-1;若当前的运行时长ti等于或大于模拟总时长T,则模拟终止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910410513.2A CN110245385B (zh) | 2019-05-17 | 2019-05-17 | 一种考虑气液固三相耦合的曝气法离散元数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910410513.2A CN110245385B (zh) | 2019-05-17 | 2019-05-17 | 一种考虑气液固三相耦合的曝气法离散元数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110245385A CN110245385A (zh) | 2019-09-17 |
CN110245385B true CN110245385B (zh) | 2023-05-26 |
Family
ID=67884556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910410513.2A Active CN110245385B (zh) | 2019-05-17 | 2019-05-17 | 一种考虑气液固三相耦合的曝气法离散元数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110245385B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112069719B (zh) * | 2020-09-02 | 2024-02-20 | 中国人民解放军陆军工程大学 | 落石冲击垫层的三维细观离散元分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529087A (zh) * | 2013-10-22 | 2014-01-22 | 东南大学 | 一种地下水曝气修复二维模型试验成像方法 |
CN103529190A (zh) * | 2013-10-22 | 2014-01-22 | 东南大学 | 一种曝气联合气相抽提二维试验装置 |
-
2019
- 2019-05-17 CN CN201910410513.2A patent/CN110245385B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529087A (zh) * | 2013-10-22 | 2014-01-22 | 东南大学 | 一种地下水曝气修复二维模型试验成像方法 |
CN103529190A (zh) * | 2013-10-22 | 2014-01-22 | 东南大学 | 一种曝气联合气相抽提二维试验装置 |
Non-Patent Citations (2)
Title |
---|
Study of Flow Transitions during Air Sparging Using the Geotechnical Centrifuge;Hu, Liming等;《Journal of Environmental Engineering》;20151231;全文 * |
地下水曝气法的两相渗流数值模拟;王建等;《第一次全国岩石力学与工程学术大会论文集》;20101018;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110245385A (zh) | 2019-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yao et al. | Numerical simulation of hydraulic performance with free overfall flow | |
CN204458973U (zh) | 一种简易液体阻尼tmd | |
CN203705311U (zh) | 一种研究变水头作用下双向管涌的实验装置 | |
CN110245385B (zh) | 一种考虑气液固三相耦合的曝气法离散元数值模拟方法 | |
CN107766978B (zh) | 不规则井网的智能优化方法 | |
CN107338758A (zh) | 河网交汇区复杂河床下水沙运动的实验系统及其模拟方法 | |
CN112257313B (zh) | 一种基于gpu加速的污染物输移高分辨率数值模拟方法 | |
CN111441306B (zh) | 一种水动力提升方法及其性能测试方法 | |
CN111767684B (zh) | 一种优化的摩阻源项隐式格式二维浅水方程建模方法 | |
CN112007943B (zh) | 原位注入高压旋喷注射修复地下水污染注入药剂扩散半径的确定方法 | |
CN107939372B (zh) | 针对小断块油藏的最优井位部署方法和装置 | |
CN110990926B (zh) | 一种基于面积修正率的城市地表建筑水动力学仿真方法 | |
JP2022542132A (ja) | 杭端に空洞が存在することによる杭基礎の軸受性能への影響をシミュレーションする試験装置及び試験方法 | |
CN112749475A (zh) | 一种确定梯级水库群连续溃坝风险分析方法 | |
CN108731892B (zh) | 一种水下振动台地形模拟假底和分区结构 | |
CN216209117U (zh) | 一种实验室裂隙含水层模拟装置 | |
CN115078237A (zh) | 一种研究物质循环的潮汐模拟系统 | |
Park et al. | Numerical analysis on energy dissipation of tsunami through vegetation zone | |
Davies et al. | " SUSPENSION OF SAND IN OSCILLATORY FLOW ABOVE RIPPLES: DISCRETE-VORTEX MODEL AND LABORATORY EXPERIMENTS | |
Yang et al. | One-and two-dimensional coupled hydrodynamics model for dam break flow | |
Xiao et al. | PHYSICAL AND ANALYTICAL MODELING OF RHYTHMIC KARST SPRINGS. | |
CN105956282A (zh) | 一种提高饱和土地基强夯施工效率的计算方法 | |
CN111324926B (zh) | 一种“以面代孔”地下工程排水孔幕的模拟方法 | |
Giri et al. | Slope stability of bioreactor landfill with leachate recirculation using horizontal trench system | |
CN212514573U (zh) | 一种岩土工程滑坡模型试验箱 |
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 |