模拟地下水流运动的有限体积Yeh多尺度有限元法制造技术

本发明专利技术公开了一种模拟地下水流运动的有限体积Yeh多尺度有限元法，包括：确定粗、细单元的尺度；将研究区剖分为粗单元；将每一粗单元剖分为细单元；以粗网格上的每一未知节点为中心，连接与该节点相关粗单元的中心，将研究区剖分为互不重叠的体积元；在每一粗单元上构造基函数，构造x和y方向速度矩阵；在每一体积元上将水流方程积分，应用散度定理进行变换，通过速度矩阵将达西速度项应用水头的粗尺度解线性表示，获得该体积元粗尺度水头解的方程，结合差分格式，相加水头的总方程；采用有效的计算方法求解总方程，同时通过速度矩阵获得达西速度值。与多种经典方法相比，本发明专利技术的有限体积Yeh多尺度有限元法具有更高的计算效率。

Finite Volume Yeh Multiscale Finite Element Method for Simulating Groundwater Flow Movement

【技术实现步骤摘要】
模拟地下水流运动的有限体积Yeh多尺度有限元法
本专利技术属于水力学
，具体涉及一种模拟地下水流运动的有限体积Yeh多尺度有限元法。
技术介绍
21世纪以来，地下水污染、海水入侵等问题的重要性与日俱增，建立精确高效的水流模型综合分析水头、达西速度等参数具有重要意义。精确高效的水头、达西速度解法能够精确描述地下水位及运移状态，能够令地下水工作者能更好的了解地下水系统情况，具有重要的研究意义。地下水系统的非均质性常跨越很多尺度，应用传统有限元法直接进行水头求解非常困难，需要非常大的计算消耗。因此，科研工作者提出了多尺度有限单元法(HouandWu1997)改进了传统有限元法，通过在粗单元上构造多尺度基函数抓住介质的细尺度信息，并直接在粗尺度上求解水头，能够大幅降低计算消耗。随后，科学工作者在“Finitevolumemultiscalefiniteelementmethodforsolvingthegroundwaterflowproblemsinheterogeneousporousmedia”一文中提出了有限体积多尺度有限元法(HeandRen)，将有限体积法和多尺度有限单元法本文档来自技高网...

【技术保护点】
1.一种模拟地下水流运动的有限体积Yeh多尺度有限元法，其特征在于，包括步骤如下：（1）根据研究区域确定所要模拟的地下水问题的边界条件，设定粗尺度，对该研究区域进行网格剖分，得到粗单元；（2）设定细尺度，对上述每一粗单元进行网格剖分，得到细单元；（3）根据步骤（1）中的粗尺度，以步骤（1）中获得的研究区剖分网格上的每一未知节点为中心，连接与该节点相关粗单元的中心，获得该节点的体积元，将研究区剖分为互不重叠的矩形体积元；（4）在步骤（2）中的网格剖分下，根据渗透系数、粗单元顶点上的多尺度基函数值以及多尺度基函数边界条件公式，在粗单元上求解退化的椭圆型问题，获得多尺度基函数在粗单元所有节点上的值；...

【技术特征摘要】
1.一种模拟地下水流运动的有限体积Yeh多尺度有限元法，其特征在于，包括步骤如下：（1）根据研究区域确定所要模拟的地下水问题的边界条件，设定粗尺度，对该研究区域进行网格剖分，得到粗单元；（2）设定细尺度，对上述每一粗单元进行网格剖分，得到细单元；（3）根据步骤（1）中的粗尺度，以步骤（1）中获得的研究区剖分网格上的每一未知节点为中心，连接与该节点相关粗单元的中心，获得该节点的体积元，将研究区剖分为互不重叠的矩形体积元；（4）在步骤（2）中的网格剖分下，根据渗透系数、粗单元顶点上的多尺度基函数值以及多尺度基函数边界条件公式，在粗单元上求解退化的椭圆型问题，获得多尺度基函数在粗单元所有节点上的值；（5）在步骤（2）中的网格剖分下，根据渗透系数、粗单元上的基函数，结合Yeh的伽辽金有限元模型框架，在粗单元上求解x方向上的达西定律方程，获得粗单元上关于Vx的速度矩阵；（6）在步骤（2）中的网格剖分下，根据渗透系数、粗单元上的基函数，结合Yeh的伽辽金有限元模型框架，在粗单元上求解y方向上的达西定律方程，获得粗单元上关于Vy的速度矩阵；(7)在步骤（3）中每一个体积元上将水流方程积分，根据有限体积多尺度有限元法中水头的粗尺度解的定义，将水头的粗尺度解对时间偏微分导数代入水流方程，应用散度定理对水流方程进行变换；（8）将步骤（7）中的经过散度定理变换得到的水流方程离散到与该体积元相关的粗单元上，在每个粗单元上通过速度矩阵将水流方程中的达西速度项应用水头的粗尺度解线性表示，获得该体积元上关...

【专利技术属性】
技术研发人员：谢一凡，吴吉春，鲁春辉，叶逾，薛禹群，谢春红，
申请(专利权)人：河海大学，南京大学，
类型：发明
国别省市：江苏,32