Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration     

Xiaoxian Zhang  A. Glyn Bengough  John W. Crawford  Iain M. Young 《Journal of Hydrology》2002,260(1-4):75-87

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.  相似文献   

贡嘎山东坡亚高山森林区蒸散力的估算   总被引：3，自引：1，他引：3  

周杨明  程根伟  杨清伟 《山地学报》2002,20(2):135-140

以海螺沟 30 0 0m气象站的观测资料为基础 ,运用Penman公式法　空气饱和差法和桑斯维特公式法 ,计算了贡嘎山东坡亚高山林林的年平均蒸散力 ,其结果为分别为 431 81mm、171 4mm和 44 6 4mm。通过分析蒸散力的影响因素 ,对这三种计算蒸散力的方法作了比较 ,对计算结果存在的差异作了比较合理的解释 ,认为用Penman公式法可以估计出本研究区的蒸散力。同时 ,对蒸散力及其影响因自进行了相关分析 ,指出湿度和风不是速制约研究区蒸散力的主导因子 ,并分析了蒸散力与水面蒸发的关系 ,由此推导出估算蒸散力的简便方程 :PE =6 6 77 0 6 91E60 1 0 75 1ITm 0 0 396Pm。  相似文献   

云南省金沙江流域土壤流失方程研究   总被引：17，自引：2，他引：17  

杨子生 《山地学报》2002,(Z1)

云南金沙江流域是长江中上游水土流失最严重的区域。本项研究以“通用土壤流失方程”(USLE)为蓝本,运用小区实验等手段,综合分析了各个侵蚀因子,建立了云南金沙江流域土壤流失方程A=R·K·LS·c·P,并确定了方程中诸因子的求算方法和数值,以及该流域土壤允许流失量,为方程的应用提供了基本的技术数据。同时,还进行了方程的检验,方程计算值与小区实测值的相对误差在6.3％以下,表明该方程在实际应用中是可靠的。该方程的建立,可为云南金沙江流域预测预报土壤侵蚀,制定土地合理利用规划方案、水土保持措施和土地生态安全格局提供了一套可靠的科学方法和依据。  相似文献   

云南省金沙江流域水土流失基本特征分析     

杨子生  梁洛辉  王云鹏 《山地学报》2002,(Z1)

在应用新近建立的云南金沙江流域土壤流失方程测算该流域各县(市、区)年均土壤流失总量和各地类年均 土壤流失量基础上,分析了该流域水土流失的总体特征、各地类水土流失特征和不同坡度级坡耕地的水土流失特 征,为因地制宜地防治水土流失提供了依据。  相似文献   

Determining the magnetic field in the core–mantle boundary zone by non-harmonic downward continuation     

L. Ballani  H. Greiner-Mai  D. Stromeyer 《Geophysical Journal International》2002,149(2):374-389

  相似文献   

行播冬小麦热辐射方向性的影响因子   总被引：1，自引：0，他引：1  

范闻捷  徐希孺 《地理与地理信息科学》2003,19(4):84-88

混合像元的热辐射方向性由于受到多种因素的影响而变得极为复杂。该文以行播冬小麦为例，建立矩阵模型，研究混合像元热辐射方向性问题。结果表明，混合像元热辐射方向性主要受到混合像元内部各组分在三维上分布的不均一、像元比辐射率的方向性以及由外界影响造成的各组分温度差异的影响。在这三个主要影响因子中，各组分的比辐射率是地物固有的特征，不随外界条件而发生变化。在同温条件下，像元热辐射方向性完全由像元比辐射率的方向性决定。而各组分温度虽然没有方向性，但可以通过组分空间分布对其产生调制作用。  相似文献   

挠曲线复位的微分方程解法求梁的位移   总被引：5，自引：0，他引：5  

喻晓今 《华东地质学院学报》2003,26(3):271-273

梁位移求解的方法主要有两种：积分法和叠加法。积分法的困难是当多个荷载同时作用时，要按控制面分段列弯矩方程，这样，两次积分时带来的积分常数较多，导致依边界条件、连续条件定这些常数时要解多个联立方程，计算繁杂。梁位移求解的叠加法是借用典型荷载下梁已知积分结果来进行对应位置位移的叠加，具体是查转角方程、挠曲线方程表格。从表中知，梁任一位置处的位移方程大多是多次多项式，很难记忆，故离开表格便无法采用。本文找到了一种新的方法求梁的位移，通过设置比拟梁，使比拟梁的挠曲线复位于原梁，从而建立了相应的微分方程，求解此方程后便可得梁的位移，避开了分段多次积分或查表的弱点。  相似文献   

机械卡盘设计中的计算误差分析与应用          下载免费PDF全文

陈宝杰  卜长根 《探矿工程》2003,30(3):30-32,34

通过对普通机构二维力因素影响的分析，利用修正公式计算求值，以修正单因素影响最大值计算所产生的误差。利用这一计算方法，对CD—3型钻机卡盘进行了计算，得到了正确的设计结果，说明用修正公式对多因素影响采用单因素影响最大值计算进行修正是可行的。  相似文献   

改进的最小二乘法在水文分析计算中的应用     

王淑英  高永胜 《水文》2003,23(5):5-9

在水分析计算中，经常涉及到变量之间的线性或非线性拟合，而在拟合各种特性曲线时，通常应用以实测资料与拟合曲线间的误差平方和最小作为目标函数的方法——最小二乘法，但这种方法忽视了所有实测点应与拟合曲线间的相对误差尽量不超过某一百分比的原则，为了达到上述要求，提出了非线性的加权最小二乘法及线性相关方程的最小距离平方和法，探讨改进了传统的最小二乘法达到优化的效果。最小距离平方和法与常用的图解法相比，本法所得成果较为客观；与传统的单方向(x或y方向)最小二乘回归法相比，所求线性方程不会因坐标系的选取而改变。最后应用算例进行了初步讨论。  相似文献   

MATLAB及其在重力基点网平差中的应用     

冯治汉 《物探化探计算技术》2003,25(4):336-339

MATLAB作为最优秀的数学软件之一，在国内外工程技术界产生愈来愈大的影响，其强大的数据计算和图形显示功能尤为突出。复杂的物探重力基点网平差计算通过用MATLAB的矩阵计算功能，变得简单易懂。而编制的MATLAB的M文件程序，计算快速准确。相信MATLAB在物探数据处理工作中将发挥重要的作用。  相似文献