技术标签: 算法
摘要
变饱和多孔介质中水流的Richards方程解对水资源评估和环境管理越来越有用。
除了解的准确性,对所需的计算工作量也颇为关注,特别是在涉及高度非线性土壤水力特性和干燥初始条件时。
本文评估了在使用修正的Picard迭代方法解混合形式Richards方程时不同收敛准则的性能。
结果以计算机处理(CPU)时间和迭代次数进行比较。
通过使用含水量的Taylor级数展开导出的新的非线性收敛准则被实施到混合形式的数值算法中。
新准则的计算效率针对不同的土壤类型、边界条件、初始条件和分层土壤进行了评估。
尽管三个准则在计算的含水量、压力水头和水通量分布方面产生了几乎相同的结果,且质量平衡误差微乎其微,但所需的CPU时间明显不同。
总体而言,新的非线性收敛准则在计算效率上远远优于其他两个准则。
在高度非线性流问题中,其他两个收敛准则失败时,新的准则也更为健壮(即解保持收敛)。
本研究结果表明,当新的收敛准则应用于混合形式Richards方程的修正Picard解时,可以产生一种非常高效和准确的方法,用于模拟土壤中的变饱和水流。
1. 引言
在环境研究中,数值模型是评估非饱和带释放化学物质可能导致地下水污染风险的重要工具
地质材料的极端变异性和复杂性、干燥的初始条件以及不断变化的边界条件可能使得在可接受的准确性和计算工作量范围内解决流动和传输问题变得困难。
当前大多数可用的数值方法在准确性或计算效率之间做出了妥协。
例如,对于相对粗纹理的土壤的渗透通常难以建模,因为其具有高度非线性的水力特性。
干燥的初始条件和高度非线性的水力特性通常需要使用非常细致的空间和时间离散化来避免数值不稳定性。
这些条件导致数值算法变得CPU密集型,特别是当需要模拟长期和/或多维问题时。
因此,追求最佳数值算法的努力应该涉及优化方案的准确性和健壮性,以及最小化所需的计算时间。
本文的目标是提高Celia等人(1990)的混合形式算法用于解决变饱和流问题的计算效率。
通过引入新的非线性收敛准则,改进了算法的收敛速度。
所提出的准则将通过与其他两个广泛使用的收敛准则进行比较来评估其性能。
2. 背景
在等温条件下,可变饱和刚性土壤中的一维垂直水流通常用Richards方程(Richards,1931)描述,该方程可以用压力头或水含量的形式表示,即
或
其中,为压力头L,
为体积含水量(),
为时间(),
表示距离土壤表面向下的垂直距离[L],
为比水含量容积(),
为水力导度(),
为土壤水非饱和扩散率(),
为源/汇项()。
直到最近,大多数数值研究使用Richards方程的压力水头或含水量形式来描述可变饱和土壤中的流动(例如,Davis和Neuman,1983;Huyakorn等,1983,1989;Hills等,1989;Kool和Van Genuchten,1991;Kirkland等,1992)。
使用方程(1)的基于含水量的方案。
(2) 可以被写成保守质量的形式,因此在大多数情况下,无论时间步长和网格间距如何,都应在计算域内保持质量(Hills等人,1989)。
Huyakorn和Pinder(1983)表明,对于最初是干燥均质土壤,使用这样的方案是有优势的。
基于含水量的公式的一个局限性是,这种形式不能用于描述饱和带的流动,而且分层土壤中的流动也不容易模拟。
此外,即使这种公式在流动系统内部准确保持质量,基于含水量的算法在边界处可能出现质量平衡错误。
基于压力头的公式(1)被认为对涉及分层或空间异质土壤的实际问题以及变饱和流问题更有用。
不幸的是,使用基于压力头的公式模拟干燥和/或高度非线性土壤的入渗通常在保持质量方面面临困难。
有关基于压力头和基于水含量的Richards方程形式的相对优势和劣势的更详细讨论已由Hills等人(1989)提供。
一些研究人员探索了解决基于压力头或基于水含量形式的Richards方程的替代数值技术。
这些研究的目标包括获得更稳定的数值算法,加速计算,最小化质量平衡误差,并针对不同的土壤类型或初始和边界条件实现更精确的解。
例如,Milly(1985)提出了一种质量保守的解决方案,其中使用了有效元素土壤水容量项。
这种方法,结合质量集中法(Neuman,1973),有效地确保了基于压力头的方程的整体质量平衡。
Rathfelder和Abriola(1994)通过扩展和离散化土壤水容量开发了类似基于压力头方程的质量保守解决方案。
其他人使用了Kirchhoff类型的转换(例如,Ross和Bristow,1990)或替代函数(Ross,1990;Pan和Wierenga,1995)以简化非线性流的描述。
在不同的方法中,Gottardi和Venutelli(1992)使用了一种移动有限元方法,其中网格点沿润湿前沿移动,从而允许较少的节点而不损害数值精度。
然而,移动网格法在应用于分层系统或用于时变边界条件时存在一些局限性。该方法也被发现不如传统的固定网格公式具有较差的质量保守性。
El-Kadi和Ling(1993)提出了Peclet和Courant数标准,用于空间和时间离散化,以描述解决Richards方程的数值方案的准确性和效率。
通过引入源项,Hills等人(1989)成功地用基于含水量的算法解决了涉及一维水流进入分层土壤的问题。
Kirkland等人(1992)随后开发了一个基于含水量的算法,通过变量转换模拟了二维变饱和流。
最近,Huang等人(1994)提出了一种基于特征的粒子跟踪技术,用于解决基于压力头的Richards方程,用于高度非线性的入渗问题。
察觉到现有基于压力头和基于水含量的Richards方程解决方案的缺点,许多人尝试将这两种方法的优势结合起来。
Richards方程的混合形式被认为保持了基于水含量方程固有的质量守恒特性,同时提供了以压力头h表示的解。
Richards方程的混合形式表达为:
Brutsaert(1971)是最早使用混合形式Richards方程来解决饱和-非饱和流问题的人之一。
他将混合形式方程的有限差分逼近与牛顿迭代方案结合起来,以有效处理陡峻的润湿前缘。
Allen和Murphy(1985,1986)在其共轭有限元算法中使用了Richards方程的混合形式,并采用了“拟牛顿”迭代方法。
最近,Celia等人(1987)和Celia等人(1990)提出了一种质量保守的数值方案,使用“修正的Picard迭代”解混合形式Richards方程(3)。
由于完美的质量平衡,修正的Picard迭代技术被证明是对早期Picard方法的重大改进。
修正的Picard迭代方法在模拟具有陡峭润湿前缘的非饱和流方面也表现出很大的潜力(Celia等人,1990;Celia和Bining,1992)。
Ray和Mohanty(1992)随后重新审视了混合形式算法,并通过几个数值实验展示了它相对于基于压力头的方案的优势。
与基于压力水头和基于含水量的算法类似,Celia等人(1990)的修正Picard迭代方案基于Richards方程(3)的全隐式(向后Euler)时间逼近,具体如下:
文章浏览阅读3.5k次,点赞2次,收藏13次。为了从FTP服务器下载文件,需要要实现一个简单的FTP客户端。FTP(文件传输协议) 是 TCP/IP 协议组中的应用层协议。FTP协议使用字符串格式命令字,每条命令都是一行字符串,以“\r\n”结尾。客户端发送格式是:命令+空格+参数+"\r\n"的格式服务器返回格式是以:状态码+空格+提示字符串+"\r\n"的格式,代码只要解析状态码就可以了。读写文件需要登陆服务器,特殊用..._ftp 登录返回230
文章浏览阅读648次。前提:systemctl stop firewalld 关闭防火墙关闭selinux查看getenforce临时关闭setenforce 0永久关闭sed-i'/SELINUX/s/enforcing/disabled/'/etc/selinux/configselinux的三种模式enforcing:强制模式,SELinux 运作中,且已经正确的开始限制..._centos7 安装rabbitmq3.6.5
文章浏览阅读5.8k次。满意答案s55f2avsx2017.09.05采纳率:46%等级:12已帮助:5646人新版Android Studio/IntelliJ IDEA可以直接导入eclipse项目,不再推荐使用eclipse导出gradle的方式2启动Android Studio/IntelliJ IDEA,选择 import project3选择eclipse 项目4选择 create project f..._android studio 项目导入idea 看不懂安卓项目
文章浏览阅读860次,点赞2次,收藏6次。AI大模型技术已经在自然语言处理、计算机视觉、多模态交互等领域取得了显著的进展和成果,同时也引发了一系列新的挑战和问题,如数据质量、计算效率、知识可解释性、安全可靠性等。城市运维涉及到多个方面,如交通管理、环境监测、公共安全、社会治理等,它们需要处理和分析大量的多模态数据,如图像、视频、语音、文本等,并根据不同的场景和需求,提供合适的决策和响应。知识搜索有多种形式,如语义搜索、对话搜索、图像搜索、视频搜索等,它们可以根据用户的输入和意图,从海量的数据源中检索出最相关的信息,并以友好的方式呈现给用户。_ai大模型应用开发
文章浏览阅读8.2k次,点赞12次,收藏121次。为什么要测量阻抗呢?阻抗能代表什么?阻抗测量的注意事项... ...很多人可能会带着一系列的问题来阅读本文。不管是数字电路工程师还是射频工程师,都在关注各类器件的阻抗,本文非常值得一读。全文13000多字,认真读完大概需要2小时。一、阻抗测试基本概念阻抗定义:阻抗是元器件或电路对周期的交流信号的总的反作用。AC 交流测试信号 (幅度和频率)。包括实部和虚部。图1 阻抗的定义阻抗是评测电路、元件以及制作元件材料的重要参数。那么什么是阻抗呢?让我们先来看一下阻抗的定义。首先阻抗是一个矢量。通常,阻抗是_阻抗实部和虚部
文章浏览阅读955次。前面章节分享试用了pyzero,pygame但随着想增加更丰富的游戏内容,好多还要进行自己编写类,从今天开始解绍一个新的python游戏库arcade模块。通过此次的《连连看》游戏实现,让我对swing的相关知识有了进一步的了解,对java这门语言也有了比以前更深刻的认识。java的一些基本语法,比如数据类型、运算符、程序流程控制和数组等,理解更加透彻。java最核心的核心就是面向对象思想,对于这一个概念,终于悟到了一些。_arcade语言 like
文章浏览阅读1.1k次。源码简介与安装说明:2021增强版短视频去水印源码 去水印微信小程序源码网站 去水印软件源码安装环境(需要材料):备案域名–服务器安装宝塔-安装 Nginx 或者 Apachephp5.6 以上-安装 sg11 插件小程序已自带解析接口,支持全网主流短视频平台,搭建好了就能用注:接口是公益的,那么多人用解析慢是肯定的,前段和后端源码已经打包,上传服务器之后在配置文件修改数据库密码。然后输入自己的域名,进入后台,创建小程序,输入自己的小程序配置即可安装说明:上传源码,修改data/_去水印机要增强版
文章浏览阅读557次。1. 触发器是FPGA存储数据的基本单元2. 触发器作为时序逻辑的基本元件,官方提供了丰富的配置方式,以适应各种可能的应用场景。_fdre #(.init(1'b0) // initial value of register (1'b0 or 1'b1) ) fdce_osc (
文章浏览阅读560次。本该是不同编译器结果不同,但是尝试了g++ msvc都是先计算c,再计算b,最后得到a+b+c是经过赋值以后的b和c参与计算而不是6。由上表可知,将q复制到p数组可以表示为:*p++=*q++,*优先级高,先取到对应q数组的值,然后两个++都是在后面,该行运算完后执行++。在电脑端编译完后会分为text data bss三种,其中text为可执行程序,data为初始化过的ro+rw变量,bss为未初始化或初始化为0变量。_嵌入式面试笔试c语言知识点
文章浏览阅读2.3k次。57 Things I've Learned Founding 3 Tech CompaniesJason Goldberg, Betashop | Oct. 29, 2010, 1:29 PMI’ve been founding andhelping run techn_mature
文章浏览阅读1.9k次。问题:先讲下需求,有若干个文本文件(txt或者csv文件等),每行代表一条数据,现在希望能合并成 1 个文本文件,且需要去除重复行。分析:一向奉行简单原则,如无必要,绝不复杂。如果数据量不大,那么如下两条命令就可以搞定合并:cat a.txt >> new.txtcat b.txt >> new.txt……去重:cat new...._python 超大文本合并
文章浏览阅读489次。这个过渡页是第一次打开小程序展示的,点击某个小程序前把手机的开发者->network link conditioner->enable & very bad network 就会在停在此页。比如《支付宝运动》这个小程序先看这个类的.h可以看到它继承于DTViewController点击左上角返回的方法- (void)back;#import "DTViewController.h"#import "APBaseLoadingV..._类似支付宝页面过度加载页