技术标签: 算法
摘要
变饱和多孔介质中水流的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)时间逼近,具体如下:
文章浏览阅读2.5w次。另一方面,随着对数字化转型和敏捷性的更加关注,项目需要比以前更快地交付,但仍需要与更大的核心系统相互连接,这些系统的变化速度比其数字系统要慢得多。在这种情况下,一个人会有重复的功能(使用上述三种模式中的任何一种),但不是将功能流量定向到特定系统,而是将特定流量或条件流量定向到一个系统,并将剩余流量定向到另一个系统.微服务由不同的团队构建,并依赖于它们之间的交互。因此,这意味着团队之间需要非常密切的合作,以确保规范的设计和构建是正确的,不仅是在设计期间,而且在以后进行更改时也是如此,以避免互操作性问题。_decoupled architecture:
文章浏览阅读233次。(以安装DBI和DBD::mysql为例)方法一:cpan(这个是最简单的方法,首次使用时请按照方法二安装CPAN(v2.05)这个模块)>cpan>installDBI#使用该命令,DBI安装成功,但是DBD-mysql没有安装成功,报错如下methods to determineprerequisitesFailed during this command:CAPTTOFU/DBD..._[dbd-mysql-4.023] 错误 1
文章浏览阅读516次。在Oracle中的树形操作在Oracle中的树形操作1.取子节点及路径(正树):select t.id ,t.code, t.name ,t.pid,SYS_CONNECT_BY_PATH(t.id,'.')||'.' as IdPathfrom tas_catalog t--where id!=110start with id=110connect by pid = prior idorder..._oracle树轴
文章浏览阅读2.7k次,点赞5次,收藏16次。之前用tornado实现tcp server与硬件的client的连接,以便下发指令。基于django/tornado与wifi模块,lora设备的TCP/UDP即时通讯现在考虑到天猫精灵可以语音接入。也就用tornado来实现天猫精灵的OAuth2.0.发了两天空闲的时间去摸索,终于实现成功了。我的参考文章。Java实现Aligenie天猫精灵OAuth2.0认证授权流程重要的还是..._天猫精灵 python
文章浏览阅读4.7w次,点赞7次,收藏49次。系统整体对比对比说明/文件系统TFSFastDFSMogileFSMooseFSGlusterFSCeph开发语言C++CPerlCCC++开源协议GPL V2GPL _开源分布式文件存储系统应用举例和产品
文章浏览阅读384次。1、禁用缩放<meta name="viewport" content="width=device-width, initial-scale=1.0, user-scalable=no"/>这一方案在 safari 上并不起作用,而且大大降低了移动端页面的可用性和可访问性。2、使用指针事件 (IE10+)a, button{ -ms-touch-actio..._safari移动端下拉调分页接口延迟的问题
文章浏览阅读911次。安装以上包版本必须对应,不然会报错本人python 3.8,cpu通过参考各方大神,最后成功安装1)torch1.4,对应的torchvision版本是0.5.0对应版本参考:https://blog.csdn.net/AugustMe/article/details/109028989torchvision下载链接:https://pypi.org/project/torchvision/0.5.0/#filestorch1.4下载链接:https://download.pytorch.._torch_scatter下载
文章浏览阅读1.1k次。转自:http://blog.csdn.net/dapeng0802/article/details/50526365[Python标准库]decimal——定点数和浮点数的数学运算 作用:使用定点数和浮点数的小数运算。 Python 版本:2.4 及以后版本 decimal 模块实现了定点和浮点算术运算符,使用的是大多数人所熟悉的模型,而不是程序员熟悉的...
文章浏览阅读1.3k次,点赞28次,收藏29次。免费的 PDF 查看器很容易找到。除了 Adobe Acrobat Reader 之外,许多程序都可以显示 PDF 文件,但找到一个好的免费 PDF 编辑器(允许您编辑和保存现有 PDF 的程序)是一个更具挑战性的提议。问题在于,几乎所有能够编辑 PDF 文件的程序都试图向您推销附加功能,无论是按月订阅还是免费提供基本功能但对高级功能收费的“免费增值”模式。_奇客pdf
文章浏览阅读3.7k次,点赞3次,收藏2次。PostgreSQL时区SQL 标准通过"+"或者"-"是否存在来区分timestamp without time zone和timestamp with time zone文本。因此,TIMESTAMP '2021-03-06 18:02:00'是一个timestamp without time zone,而TIMESTAMP '2021-03-06 18:02:00 +08'是一个timestamp with time zone。+08:表示 时区与全球统一时间 ..._postgreql to_char 时区
文章浏览阅读1.3w次。在MyBatisPlus中添加注解信息即可//实体类(此处需要添加对应的MyBatisPlus主键信息)public class Notice{ private int noticeId; private String release_date; private String content;}//方案一:自定义SQL(自增主键) @Insert("insert into notice(release_date,content)values(#{releaseDate},#_mybatisplus新增返回主键
文章浏览阅读132次。因为 Android P 是默认禁止访问http的API的 所以想要跳过https的限制需要在。在Android O 升 P (8.0升9.0)的时候,会报以下异常。的Application中添加以下内容。_http failed: java.net.unknownserviceexception: cleartext communication to 19