公开(公告)号：CN106202746B

公开(公告)日：2019-04-16

申请号：CN201610556975.1

申请日：2016-07-14

Applicant: 南京大学

Inventor： 谢一凡 ,  吴吉春 ,  薛禹群 ,  谢春红

IPC: G06F17/50

Abstract: 本发明公开了一种模拟多孔介质中水流达西速度的Yeh‑多尺度有限元方法，运用伽辽金法将求解问题变分；将研究区剖分为粗网格单元，对所有粗网格单元进行剖分，得到细网格单元；在每一粗网格单元上求解退化椭圆方程，得到基函数；运用基函数求解变分形式，得到总刚度矩阵；根据研究区的源汇项与边界条件得到右端项；联立得到水头方程组；运用有效的数值方法求解该方程组得到研究区的节点水头；结合Yeh的伽辽金有限元模型，运用所构造的基函数和研究区的水头值，在研究区域直接求解达西方程，得到粗尺度节点上连续的达西渗透流速，再利用基函数线性表示细尺度达西渗透流速。与现有技术相比，该方法具有相近的精度和更高的效率。

公开(公告)号：CN107657075A

公开(公告)日：2018-02-02

申请号：CN201710732804.4

申请日：2017-08-24

Applicant: 南京大学

Inventor： 谢一凡 ,  吴吉春 ,  薛禹群 ,  谢春红

IPC: G06F17/50

Abstract: 本发明公开了一种模拟地下水介质交界面处达西速度的区域分解有限元法，包括：运用伽辽金法对地下水流问题变分，剖分研究区域，应用有限元法获得水头；根据研究区域的介质组成，应用不同介质的交界面将研究区域分解为若干单一介质子区域，应用区域分解法将研究区域上的地下水达西速度求解问题分解为子区域上的子问题；选取一个子问题，运用伽辽金法变分，应用Yeh的伽辽金模型获得达西速度，结合折射定律获得该子区域和其他子区域交界面上另一侧的达西速度，作为相邻子问题的第一类边界条件；选取下一个子问题求解，重复这一过程直到所有子问题求解完毕。该方法通过区域分解法降低了达西速度的计算消耗，应用折射定律保证交界面处达西速度的精度。

公开(公告)号：CN105354362A

公开(公告)日：2016-02-24

申请号：CN201510647088.0

申请日：2015-10-08

Applicant: 南京大学

Inventor： 谢一凡 ,  吴吉春 ,  薛禹群 ,  谢春红

IPC: G06F17/50

CPC classification number: G06F17/5018

Abstract: 本发明公开了一种模拟二维水流运动的三次样条多尺度有限元方法，该方法首先将需要求解的问题转换成变分形式；确定研究区边界条件，设定粗网格单元尺度h，剖分研究区域，得到粗网格单元，将粗网格单元剖分为细网格单元；在每一粗网格单元上求解退化椭圆型问题，构造基函数；采用有效的计算矩阵的方法求解总刚度矩阵和右端项联立方程组；求得研究区上每个节点(粗尺度)的水头；采用三次样条技术求得粗尺度水头导数，根据达西定律获得粗尺度达西渗流速度；通过基函数和插值公式得到粗单元上每个节点(细尺度)水头；采用三次样条法求得细尺度水头导数，根据达西定律获得细尺度达西渗流速度。与现有技术相比，该方法的精度与其相近，计算效率更高。