CN107357759A - 基于渗流边界和运动微分方程条件的渗流求解方法 - Google Patents
基于渗流边界和运动微分方程条件的渗流求解方法 Download PDFInfo
- Publication number
- CN107357759A CN107357759A CN201710496148.2A CN201710496148A CN107357759A CN 107357759 A CN107357759 A CN 107357759A CN 201710496148 A CN201710496148 A CN 201710496148A CN 107357759 A CN107357759 A CN 107357759A
- Authority
- CN
- China
- Prior art keywords
- mrow
- equation
- boundary
- head
- seepage
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 4
- 230000008595 infiltration Effects 0.000 claims description 11
- 238000001764 infiltration Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000009189 diving Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种基于渗流边界和运动微分方程条件的渗流求解方法,其特征在于,包括以下步骤:步骤1. 对研究对象的宏观几何特征进行准确测量,并建立相对应的几何特征描述方程;步骤2. 研究研究对象的渗流系数分布特征,建立研究对象在研究区域的渗流系数分布方程;步骤3. 研究研究对象的边界条件流量、水头特征,并建立相对应的边界条件流量、水头的表示方程;步骤4. 选取水头的表示方程,使其满足相对应的运动微分方程、流量和水头的边界条件方程,并求解相对应的各常系数。该方法的提出将对大坝、边坡、路基、隧道、巷道、涵洞等各种营建构筑物和自然体的渗流运动研究和应用具有推动作用。
Description
技术领域
本发明属于渗流、土木工程和地质工程等与渗流和变形相关的技术领域,具体涉及一种基于渗流边界和运动微分方程条件的渗流求解方法。
技术背景
现行渗流求解往往建立在有限单元等数值计算方法上,有限单元法采用的是以点代面(注:对于二维问题)或以点代体(注:对于三维问题)),因而不同大小的单元,其计算结果并不相同;再则,数值计算时采用的是用线性方法求解非线性问题,当然不同的收敛标准其结果也是不一样的。但是,对于一个形状固定的承压水和潜水研究对象,其流态应该是确定的,目前缺乏对其进行渗流求解的方法。
发明内容
本发明是为了解决上述问题而进行的,目的在于提供一种基于渗流边界和运动微分方程条件的渗流求解方法,该方法基于一个形状固定的研究对象,其对应的流态也应是确定的事实,在假设水头满足运动微分方程和流量、水头的边界条件方程下,获得相对应的流量、水头理论解。这种方法将对现行渗流求解方法向前推进一大步。
本发明为了实现上述目的,采用了以下方案:
本发明提供一种基于渗流边界和运动微分方程条件的渗流求解方法,其特征在于,包括以下步骤:步骤1.对研究对象的宏观几何特征进行准确测量,并建立相对应的几何特征描述方程;步骤2.研究研究对象的渗流系数分布特征,建立研究对象在研究区域的渗流系数分布方程;步骤3.研究研究对象的边界条件流量、水头特征,并建立相对应的边界条件流量、水头的表示方程;步骤4.选取水头的表示方程,使其满足相对应的运动微分方程、流量和水头的边界条件方程,并求解相对应的各常系数。
本发明提供的基于渗流边界和运动微分方程条件的渗流求解方法,还可以具有这样的特征:在步骤1中,在研究对象的边界为直线的情况下,该边界的几何特征描述方程为直线方程;在研究对象的边界为曲线的情况下,该边界的几何特征描述方程为曲线方程。
本发明提供的基于渗流边界和运动微分方程条件的渗流求解方法,还可以具有这样的特征:在步骤3中,在研究对象为二维对象的情况下,边界条件流量、水头表达式为:
(I)第一类边界条件:如果在某一部分边界S1上,各点在每一时刻的水头都是已知的,这部分边界称为第一类边界,相应的边界条件表示为:
式中,H:水头,已知函数,x,y,z:坐标,t:时间;
(II)第二类边界条件:当知道某一部分边界S2单位面积上流入的流量q1时,这部分边界称为第二类边界,相应的边界条件表示为:
式中,n为S2的外法线方向,q1为已知函数,表示S2的侧向补给量,K为渗透系数;
(III)第三类边界条件:某边界上H和的线性组合是已知的,这部分边界称为第三类边界,相应的边界条件表示为:
式中,α,β为已知函数。
本发明提供的基于渗流边界和运动微分方程条件的渗流求解方法,还可以具有这样的特征:在步骤3中,在进行流量、水头求解时,对于四面体、六面体在已知两个面、三个面的情况下求解,对于三角形、四边形、五边形在已知两个边、三个边的情况下求解,相对应的其它面或边的边界条件流量、水头按求解的特征加以计算。
本发明提供的基于渗流边界和运动微分方程条件的渗流求解方法,还可以具有这样的特征:在步骤4中,针对二维稳定渗流研究对象:
假设水头表达式为:
H=a0+a1x+a2y+a3x2+a4xy+a5y2+a6x3+a7x2y+a8xy2+… (4)
假设对应的渗透系数方程为:
Kxx=b1,0+b1,1x+b1,2y+b1,3x2+b1,4xy+b1,5y2+b1,6x3+b1,7x2y+b1,8xy2+… (5)
Kyy=b2,0+b2,1x+b2,2y+b2,3x2+b2,4xy+b2,5y2+b2,6x3+b2,7x2y+b2,8xy2+… (6)
式中:ai,b1,i,b2,i为常系数,i属于零和整数;
满足力的运动微分方程为:
在任意坐标条件下,运动微分方程满足的必要条件为相对应的各项系数为零,在假设渗透系数为常数情况下,由方程(7)可得一系列常系数方程,由此可以决定一定数量的常系数,进而可以获得研究对象的流量、水头理论解。
发明的作用与效果
本发明提出了一种基于渗流边界和运动微分方程条件的渗流求解方法,该方法针对现有承压水和潜水分析技术的缺陷,在假设研究对象流量、水头满足渗流边界和运动微分方程条件下,根据研究对象的几何特征,求解研究对象的流量、水头等分布,可以求解任意几何形状(如平面问题:三边形、四边形、五边形、六边形等,三维问题:四面体、六面体、八面体等)的流量、水头等分布;在研究对象的边界流量、水头与边界条件流量、水头不一致的条件下,可以获得水头不连续解。该方法的提出将对大坝、边坡、路基、隧道、巷道、涵洞等各种营建构筑物和自然体的渗流运动研究和应用具有推动作用。
附图说明
图1为本发明实施例所涉及的第i研究对象的边界流量、水头特征示意图。
图中,AB、BC、CO和OA为相应的边界,为对应边界OA的水头条件;S1,S2分别为对应边界OA、BC。
具体实施方式
以下结合附图对本发明涉及的基于渗流边界和运动微分方程条件的渗流求解方法的具体实施方案进行详细地说明。
(1)在对研究对象准确测量研究的基础上,建立相对应的几何特征描述方程,如图1,相对应AO、AB、BC和CO的方程可以表征为:y=kx+b(如是曲线等形式,则可以以曲线等方程加以表示)。
(2)在对研究对象的渗透分布特征研究的基础上,建立研究对象在研究区域的渗透系数分布方程,相对应的渗透系数为Kxx,Kyy;
(3)在对研究对象的边界条件流量、水头特征研究的基础上,并建立相对应的边界条件流量、水头方程;对于二维问题,在图1的基础上,边界条件流量、水头表达式为:
(1)第一类边界条件(Dirichlet条件):如果在某一部分边界(见图1中AO边界)上,各点在每一时刻的水头都是已知的,则这部分边界就称为第一类边界或给定水头的边界,表示为:
式中:H:水头,已知函数,x,y:坐标,t:时间。
(2)第二类边界条件(Neumam条件):当知道某一部分边界(如图1中BC边界)单位面积(二维空间为单位宽度)上流入(流出时用负值)的流量q1时,称为第二类边界或给定流量的边界。相应的边界条件表示为:
式中,n为边界BC的外法线方向。q1则为已知函数,分别表示S2上单位面积的侧向补给量。K为渗透系数。
(3)第三类边界条件:某边界上H和的线性组合是已知的,即有:
又称混合边界条件,α,β为已知函数。
在流量、水头连续的条件下该表达式必须能描述相对应的全部边界条件流量和水头(如图1);如果流量、水头不连续,则该表达式不成立。
对于边界条件流量、水头:在进行流量、水头求解时,可以在已知两个面、三个面等(对于四面体、六面体等)或两个边、三个边等(三角形和四边形、五边形等)边界条件流量、水头下,求解对应研究对象的流量、水头,相对应的其它面或边的边界条件流量、水头按求解解的特征加以计算。
(4),选取水头的表示方程,使其满足相对应的运动微分方程、流量和水头的边界条件方程,并求解相对应的各常系数;针对二维稳定渗流问题,表述如下:
假设水头表达式(注:也可以取其它表示形式)为:
H=a0+a1x+a2y+a3x2+a4xy+a5y2+a6x3+a7x2y+a8xy2+… (4)
假设对应的渗透系数方程为:
Kxx=b1,0+b1,1x+b1,2y+b1,3x2+b1,4xy+b1,5y2+b1,6x3+b1,7x2y+b1,8xy2+… (5)
Kyy=b2,0+b2,1x+b2,2y+b2,3x2+b2,4xy+b2,5y2+b2,6x3+b2,7x2y+b2,8xy2+… (6)
式中:ai,b1,i,b2,i为常系数,i:属于零和整数。
满足力的运动微分方程为:
在任意坐标条件下,运动微分方程满足的必要条件为相对应的各项系数为零,在假设渗透系数为常数(注:可以研究渗透系数满足方程(5、6)的情况)情况下,则有:
由方程(7)可得:
a3Kxx+a5Kyy=0 (8)
3a6Kxx+a8Kyy=0 (9)
a7Kxx+3a9Kyy=0 (10)
……
对于第i研究对象(如图1),在边界条件流量、水头作用下,在流量、水头连续条件下边界流量、水头和边界条件流量、水头必须相等,如存在部分水头不连续,则相对应的边界水头和边界条件水头不相等;但边界条件流量应该相等。
按照上述步骤,可以决定一定数量的常系数,因而可以获得研究对象的流量、水头理论解,当研究对象复杂时,可以将整个研究对象划分为几个不同的小对象加以求解,但解必须满足几个不同研究对象流量、水头等之间的关系。
按照上述步骤,可以获得整个研究对象流量、水头的理论解,并可以与研究对象的现场状态加以对比,从而修正理论的各种渗流参数。
以上事例仅为举证,而并非是对本发明的实施方式的限定。除上述事例外,本发明还有其他实施方式。凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围。
Claims (5)
1.一种基于渗流边界和运动微分方程条件的渗流求解方法,其特征在于,包括以下步骤:
步骤1.对研究对象的宏观几何特征进行准确测量,并建立相对应的几何特征描述方程;
步骤2.研究所述研究对象的渗流系数分布特征,建立所述研究对象在研究区域的渗流系数分布方程;
步骤3.研究所述研究对象的边界条件流量、水头特征,并建立相对应的边界条件流量、水头的表示方程;
步骤4.选取水头的表示方程,使其满足相对应的运动微分方程、流量和水头的边界条件方程,并求解相对应的各常系数。
2.根据权利要求1所述的基于渗流边界和运动微分方程条件的渗流求解方法,其特征在于:
在所述步骤1中,在所述研究对象的边界为直线的情况下,该边界的几何特征描述方程为直线方程;在所述研究对象的边界为曲线的情况下,该边界的几何特征描述方程为曲线方程。
3.根据权利要求1所述的基于渗流边界和运动微分方程条件的渗流求解方法,其特征在于:
在所述步骤3中,边界条件流量、水头表达式为:
(I)第一类边界条件:如果在某一部分边界S1上,各点在每一时刻的水头都是已知的,这部分边界称为第一类边界,相应的边界条件表示为:
式中,H:水头,已知函数,x,y,z:坐标,t:时间;
(II)第二类边界条件:当知道某一部分边界S2单位面积上流入的流量q1时,这部分边界称为第二类边界,相应的边界条件表示为:
<mrow>
<mi>K</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>H</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,n为边界S2的外法线方向,q1为已知函数,表示S2的侧向补给量,K为渗透系数;
(III)第三类边界条件:某边界上H和的线性组合是已知的,这部分边界称为第三类边界,相应的边界条件表示为:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>H</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>&alpha;</mi>
<mi>H</mi>
<mo>=</mo>
<mi>&beta;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,α,β为已知函数。
4.根据权利要求1所述的基于渗流边界和运动微分方程条件的渗流求解方法,其特征在于:
在所述步骤3中,在进行流量、水头求解时,对于四面体、六面体在已知两个面、三个面的情况下求解,对于三角形、四边形、五边形在已知两个边、三个边的情况下求解,相对应的其它面或边的边界条件流量、水头按求解的特征加以计算。
5.根据权利要求3所述的基于渗流边界和运动微分方程条件的渗流求解方法,其特征在于:
在所述步骤4中,针对二维稳定渗流研究对象:
假设水头表达式为:
H=a0+a1x+a2y+a3x2+a4xy+a5y2+a6x3+a7x2y+a8xy2+… (4)
假设对应的渗透系数方程为:
Kxx=b1,0+b1,1x+b1,2y+b1,3x2+b1,4xy+b1,5y2+b1,6x3+b1,7x2y+b1,8xy2+… (5)
Kyy=b2,0+b2,1x+b2,2y+b2,3x2+b2,4xy+b2,5y2+b2,6x3+b2,7x2y+b2,8xy2+… (6)
式中:ai,b1,i,b2,i为常系数,i属于零和整数;
满足力的运动微分方程为:
<mrow>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>H</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>H</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
在任意坐标条件下,运动微分方程满足的必要条件为相对应的各项系数为零,在假设渗透系数为常数情况下,由方程(7)可得一系列常系数方程,由此可以决定一定数量的常系数,进而可以获得研究对象的流量、水头理论解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710496148.2A CN107357759B (zh) | 2017-06-26 | 2017-06-26 | 基于渗流边界和运动微分方程条件的渗流求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710496148.2A CN107357759B (zh) | 2017-06-26 | 2017-06-26 | 基于渗流边界和运动微分方程条件的渗流求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107357759A true CN107357759A (zh) | 2017-11-17 |
CN107357759B CN107357759B (zh) | 2020-10-30 |
Family
ID=60273581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710496148.2A Expired - Fee Related CN107357759B (zh) | 2017-06-26 | 2017-06-26 | 基于渗流边界和运动微分方程条件的渗流求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107357759B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110472309A (zh) * | 2019-07-29 | 2019-11-19 | 中国科学院力学研究所 | 一种二维渗流计算模型中考虑三维绕渗效应的方法、装置 |
CN112031754A (zh) * | 2020-08-20 | 2020-12-04 | 中国石油大学(北京) | 渗流边界的确定方法、装置和设备 |
CN114896837A (zh) * | 2022-04-27 | 2022-08-12 | 中铁城建集团有限公司 | 一种考虑支护结构厚度的基坑稳态渗流解析方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100082142A1 (en) * | 2005-11-22 | 2010-04-01 | Usadi Adam K | Simulation System and Method |
CN103452550A (zh) * | 2013-09-12 | 2013-12-18 | 中国石油大学(华东) | 一种地层渗流条件下稠油降粘剂降粘效果评价方法及装置 |
CN104504755A (zh) * | 2014-12-30 | 2015-04-08 | 华中科技大学 | 一种山体中分布式地下设施温度场仿真方法 |
CN105335607A (zh) * | 2015-10-12 | 2016-02-17 | 湖北工业大学 | 一种边坡渐进破坏潜在滑动面的计算方法 |
CN106503421A (zh) * | 2016-09-28 | 2017-03-15 | 西南交通大学 | 一种边坡渐进破坏全过程计算方法 |
-
2017
- 2017-06-26 CN CN201710496148.2A patent/CN107357759B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100082142A1 (en) * | 2005-11-22 | 2010-04-01 | Usadi Adam K | Simulation System and Method |
CN103452550A (zh) * | 2013-09-12 | 2013-12-18 | 中国石油大学(华东) | 一种地层渗流条件下稠油降粘剂降粘效果评价方法及装置 |
CN104504755A (zh) * | 2014-12-30 | 2015-04-08 | 华中科技大学 | 一种山体中分布式地下设施温度场仿真方法 |
CN105335607A (zh) * | 2015-10-12 | 2016-02-17 | 湖北工业大学 | 一种边坡渐进破坏潜在滑动面的计算方法 |
CN106503421A (zh) * | 2016-09-28 | 2017-03-15 | 西南交通大学 | 一种边坡渐进破坏全过程计算方法 |
Non-Patent Citations (1)
Title |
---|
卢应发等: "水–岩石相互作用力学参数的探讨", 《岩石力学与工程学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110472309A (zh) * | 2019-07-29 | 2019-11-19 | 中国科学院力学研究所 | 一种二维渗流计算模型中考虑三维绕渗效应的方法、装置 |
CN112031754A (zh) * | 2020-08-20 | 2020-12-04 | 中国石油大学(北京) | 渗流边界的确定方法、装置和设备 |
CN114896837A (zh) * | 2022-04-27 | 2022-08-12 | 中铁城建集团有限公司 | 一种考虑支护结构厚度的基坑稳态渗流解析方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107357759B (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Castro-Orgaz et al. | Non-hydrostatic free surface flows | |
Ern et al. | Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures | |
Pratt et al. | Rotating hydraulics: nonlinear topographic effects in the ocean and atmosphere | |
CN107357759A (zh) | 基于渗流边界和运动微分方程条件的渗流求解方法 | |
Nelson | In-place determination of permeability distribution for heterogeneous porous media through analysis of energy dissipation | |
CN106202746B (zh) | 模拟多孔介质中水流达西速度的Yeh-多尺度有限元方法 | |
Ostad-Ali-Askari et al. | Usage of rockfill dams in the HEC-RAS software for the purpose of controlling floods | |
CN108256193A (zh) | 一种河道溃决洪水一、二维数学模型耦合分析方法 | |
Hu et al. | Computationally efficient modeling of hydro-sediment-morphodynamic processes using a hybrid local time step/global maximum time step | |
CN105354362A (zh) | 模拟二维水流运动的三次样条多尺度有限元方法 | |
Birjukova Canelas et al. | Three-dimensional flow structure at fixed 70 open-channel confluence with bed discordance | |
CN103469780A (zh) | 一种重力坝坝基深层抗滑稳定临界滑动面的计算方法 | |
He et al. | Boundary integral representations of steady flow around a ship | |
Rebholz et al. | On the high accuracy NS-alpha-deconvolution turbulence model | |
Olson et al. | New finite element results for the square cavity | |
Snow | On inertial instability as related to, the multiple-vortex phenomenon | |
Kao | The phenomenon of blocking in stratified flows | |
CN109101463A (zh) | 一种广域层状大地格林函数的多精度求解方法 | |
El-Molla et al. | Seepage losses from trapezoidal earth canals with an impervious layer under the bed | |
Dong et al. | Adaptive moving grid methods for two-phase flow in porous media | |
Dong et al. | A semi-implicit method for incompressible three-phase flow in porous media | |
Bowen et al. | The linear limit of the dipole problem for the thin film equation | |
Shen et al. | Instability of Riemann solutions to a scalar conservation law with discontinuous flux | |
Mazo et al. | Absolute permeability upscaling for superelement modeling of petroleum reservoir | |
Morell et al. | Flow pattern at a river diffluence at the alluvial system of the Paraná River |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201030 |