null - 程序员宅基地

解变饱和条件下的修正Picard迭代方法的新收敛准则_picard迭代失效-程序员宅基地

技术标签：算法

摘要

变饱和多孔介质中水流的Richards方程解对水资源评估和环境管理越来越有用。

除了解的准确性，对所需的计算工作量也颇为关注，特别是在涉及高度非线性土壤水力特性和干燥初始条件时。

本文评估了在使用修正的Picard迭代方法解混合形式Richards方程时不同收敛准则的性能。

结果以计算机处理（CPU）时间和迭代次数进行比较。

通过使用含水量的Taylor级数展开导出的新的非线性收敛准则被实施到混合形式的数值算法中。

新准则的计算效率针对不同的土壤类型、边界条件、初始条件和分层土壤进行了评估。

尽管三个准则在计算的含水量、压力水头和水通量分布方面产生了几乎相同的结果，且质量平衡误差微乎其微，但所需的CPU时间明显不同。

总体而言，新的非线性收敛准则在计算效率上远远优于其他两个准则。

在高度非线性流问题中，其他两个收敛准则失败时，新的准则也更为健壮（即解保持收敛）。

本研究结果表明，当新的收敛准则应用于混合形式Richards方程的修正Picard解时，可以产生一种非常高效和准确的方法，用于模拟土壤中的变饱和水流。

1. 引言

在环境研究中，数值模型是评估非饱和带释放化学物质可能导致地下水污染风险的重要工具

地质材料的极端变异性和复杂性、干燥的初始条件以及不断变化的边界条件可能使得在可接受的准确性和计算工作量范围内解决流动和传输问题变得困难。

当前大多数可用的数值方法在准确性或计算效率之间做出了妥协。

例如，对于相对粗纹理的土壤的渗透通常难以建模，因为其具有高度非线性的水力特性。

干燥的初始条件和高度非线性的水力特性通常需要使用非常细致的空间和时间离散化来避免数值不稳定性。

这些条件导致数值算法变得CPU密集型，特别是当需要模拟长期和/或多维问题时。

因此，追求最佳数值算法的努力应该涉及优化方案的准确性和健壮性，以及最小化所需的计算时间。

本文的目标是提高Celia等人（1990）的混合形式算法用于解决变饱和流问题的计算效率。

通过引入新的非线性收敛准则，改进了算法的收敛速度。

所提出的准则将通过与其他两个广泛使用的收敛准则进行比较来评估其性能。

2. 背景

在等温条件下，可变饱和刚性土壤中的一维垂直水流通常用Richards方程（Richards，1931）描述，该方程可以用压力头或水含量的形式表示，即

$C\frac{\partial h}{\partial t}=\frac{\partial }{\partial t}(K\frac{\partial h}{\partial z}-K)-S\qquad(1)$

或

$\frac{\partial \theta}{\partial t}=\frac{\partial }{\partial z}(D\frac{\partial \theta}{\partial z})-\frac {\partial K}{\partial z}-S \qquad(2)$

其中， $h$ 为压力头L，

$\theta$ 为体积含水量（ $L^3 L^{-3}$ ），

$t$ 为时间（ $T$ ），

$z$ 表示距离土壤表面向下的垂直距离[L]，

$C = d \theta/dh$ 为比水含量容积（ $L^{-1}$ ），

$K$ 为水力导度（ $L T^{-1}$ ），

$D = K/C$ 为土壤水非饱和扩散率（ $L^2 T^{-1}$ ），

$S$ 为源/汇项（ $T^{-1}$ ）。

直到最近，大多数数值研究使用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方程的混合形式表达为：

$\frac{\partial \theta}{\partial t}=\frac{\partial }{\partial z}(K \frac{\partial h}{\partial z})-\frac {\partial K}{\partial z}-S \qquad(3)$

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）时间逼近，具体如下：

本文链接：https://blog.csdn.net/weixin_44259522/article/details/134960535

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

SaltStack数据系统_comment: no top file or master_tops data matches f-程序员宅基地

文章浏览阅读69次。文章目录1SaltStack数据系统2 SaltStack数据系统组件2.1 SaltStack组件之Grains2.2 SaltStack组件之Pillar2.3 Grains与Pillar的区别1SaltStack数据系统SaltStack有两大数据系统，分别是GrainsPillar2 SaltStack数据系统组件2.1 SaltStack组件之GrainsGrains是SaltStack的一个组件，其存放着minion启动时收集到的信息。Grains是SaltStack组件中非_comment: no top file or master_tops data matches found. please see master lo

error C2668: “swap”: 对重载函数的调用不明确_swap对重载函数的调用不明确-程序员宅基地

文章浏览阅读3.4k次，点赞6次，收藏3次。当出现这个错误提示的时候，很有可能代码并没有问题，而是编译器的设定问题。测试代码如下：#include using namespace std;//函数模板templatevoid swap(T &a,T &b){ T temp; temp = a; a = b; b = temp;}int main(int ,char **){ int a = 1, b _swap对重载函数的调用不明确

ZYNQ PCIe EP实现DMA+Linux交互，非常简洁的程序_c++ 往pcie发送数据-程序员宅基地

文章浏览阅读1w次，点赞8次，收藏120次。ZYNQ PCIe-DMA源码例程 PS-PL交互 linux/裸机 verilog C/C++ZYNQ PCIe-DMA的实现过程一、概述二、基础知识三、系统总框架四、工作原理与工作模式五、接口时序六、资源使用情况七、PS-PL交互以及测试程序ZYNQ PCIe-DMA的实现过程近期在网上淘来个源码，看了之后觉得还不错。完全刷新我对ZYNQ的认知啊，原来ZYNQ也可以这么玩的。PS-PL交..._c++ 往pcie发送数据

Python-Django毕业设计员工宿舍管理系统（程序+Lw)_用python设计一个宿舍管理系统,要求能实现寝室信息的增加、删除、修改、查询以及-程序员宅基地

文章浏览阅读173次。该项目含有源码、文档、程序数据库、配套开发软件、软件安装教程项目运行环境配置：Pychram社区版py项目技术：django + python+ Vue 等等组成，B/S模式 +pychram管理等等。环境需要1.运行环境：最好是python3.7.7，我们在这个版本上开发的。其他版本理论上也可以。2.pycharm环境：pycharm都可以。推荐pycharm社区版;3.mysql环境：建议是用5.7版本均可。_用python设计一个宿舍管理系统,要求能实现寝室信息的增加、删除、修改、查询以及

linux使用apache搭建http服务器(文件服务器)_linux怎么搭http apache服务器-程序员宅基地

文章浏览阅读8.4k次，点赞5次，收藏41次。一、安装Apache$ sudo apt-get install apache2二、修改服务器访问端口Apache2的默认访问端口为80,可修改为其他端口（当端口被占用时需要更改其访问端口）进入apache2的安装目录 /etc/apache2/，修改ports.conf文件$ cd /etc/apache2/$ sudo chmod 775 ports.conf$ vim po..._linux怎么搭http apache服务器

逆转广义表_将广义表(a,( (b, c),( ) ),( ( ( d ), e ),f) 逆转为 ((f,(e-程序员宅基地

文章浏览阅读898次。题目请编写递归算法，逆转广义表中的数据元素。例如：将广义表：(a,((b,c),()),(((d),e),f))逆转为：((f,(e,(d))),((),(c,b)),a)。代码#include <iostream>#include <cstring>#include <algorithm>using namespace std;void reverse_str(char *s, int l, int r){ int n = r-l; _将广义表(a,( (b, c),( ) ),( ( ( d ), e ),f) 逆转为 ((f,(e ,( d)),( ( ),(c,b), a)的代码

随便推点

HTC硬件介绍_htc是什么硬件-程序员宅基地

文章浏览阅读399次。其中最重要的是头盔中两个屏幕，用的是比较高端的三星AMOLED屏幕，屏幕尺寸对角线为91.8mm，PPI（每英寸集成的像素）达到了447（像素密度越高表现出来的画面就更细腻）。体验者可以在一个小范围的空间中行走，体验制作好的内容。是通过摄像头不停的采集周围环境的图像来达到设备的姿态估计（比如噪音等一些没用的信息也会影响到它的准确性）。是通过计算机视觉，以及设备自身自带的传感器比如陀螺仪、加速器来实现位置的追踪。不需要外部硬件来定位。由LED阵列和横向纵向的两个马达，马达以非常高的频率扫描场景中的定位。_htc是什么硬件

基于 Python 深度学习的车辆特征分析系统，附源码_pycharm实时检测车辆及源码-程序员宅基地

文章浏览阅读1.7k次，点赞25次，收藏29次。而在机动车的自动识别过程中，通过利用深度学习的算法来让计算机通过不断地获取信息要素形成信息库，可以更好的提升计算机对于车辆的识别能力。本次就是通过利用了深度学习技术结合Python开发工具来设计一款能够在线通过图片分析来识别车辆的品牌的软件。_pycharm实时检测车辆及源码

matlab采样点数傅里叶变换,【 MATLAB 】模拟信号采样及离散时间傅里叶变换（DTFT）案例分析...-程序员宅基地

文章浏览阅读1.1k次。中使用的模拟信号：为了研究在频域数量上的采样效果，对该信号使用两种不同的采样频率采样。a. 在 fs = 5000 对信号进行采样，求出并画出其离散时间傅里叶变换；b. 在 fs = 1000 对信号采样，求出并画出其离散时间傅里叶变换。题解：上篇博文也分析了，信号的带宽为2kHz，奈奎斯特频率就为 4000 样本/s，它小于第一问给出的采样频率，所以频谱混叠几乎不存在。我们通过MATLAB验证..._采样点时间轴变换matlab

程序员的选择-程序员宅基地

文章浏览阅读117次。程序员的15种选择前端程序员后端程序员全栈工程师运维工程师移动端开发工程师自由职业程序员测试工程师图像处理工程师游戏开发工程师交互体验工程师量化交易工程师数据科学家& 数据工程师研究型工程师创业公司程序员持续学习程序员转载于:https://www.cnblogs.com/huameixiao/p/11571598.html..._c程序员刚入门接触闭源和开源

vue 图片:src字符串拼接路径无效问题_vue 图片字符串拼接-程序员宅基地

文章浏览阅读2.1k次。无效代码：<img :src="'../../../../static/Hongkong1/img/commodity/' + oneArticleData.icon_name" />改为使用require获取图片编码即可<img :src="require('../../../../static/Hongkong1/img/commodity/' + oneArticl..._vue 图片字符串拼接

Notes Sixth day-渗透攻击-红队-打入内网_打入目标内网-程序员宅基地

文章浏览阅读3.5k次，点赞4次，收藏23次。**Notes Sixth day-渗透攻击-红队-信息收集(dayu)**作者：大余时间：2020-09-22请注意：对于所有笔记中复现的这些终端或者服务器，都是自行搭建的环境进行渗透的。我将使用Kali Linux作为此次学习的攻击者机器。这里使用的技术仅用于学习教育目的，如果列出的技术用于其他任何目标，我概不负责。我必须再重申一遍：务必不要做未授权测试！不要未经授权在真实网络环境中复现任何本书中描述的攻击。即使是出于好奇而不是恶意，你仍然会因未授权测试行为而陷入很多麻烦。为了个人能更好的继_打入目标内网