非饱和土壤水流问题的广义差分法及其数值模拟

本文利用广义差分法建立了一维非饱和水流问题的守恒形式的数值模型,讨论了半离散广义差分解和全离散广义差分解的存在唯一性,并给出了其误差估计．数值结果表明,该格式具有计算量小和稳定性等特点,且对求解非饱和水流问题有较好的适应性．     

非饱和水流问题的迎风差分法及其数值模拟

本文考察了非饱和水流问题模型方程的守恒型迎风差分法.我们基于有限体积方法建立的非饱和流动的守恒形式,分别提出了一阶和二阶迎风差分格式,并对差分格式进行了误差估计,给出了收敛性定理.最后,数值模拟验证了计算格式的有效性.     

非粘性土壤中溶质运移问题的守恒混合有限元法及其数值模拟

奉文建立了非粘性土壤水中溶质运移问题的守恒混合元格式,讨论了广义解和混合元解的存在唯一性,并给出了误差估计.数值模拟结果表叫,用该方法模拟溶质运移问题是合理有效的,不仅提高了通量的模拟精度,而且使计算稳定.     

守恒混合元法在模拟污染物运移过程中的应用

研究污染物在土壤中运移的时空规律,为土壤环境质量评价及污染预测和防治提供科学的根据与途径,具有重要的理论和实际意义.通过建立土壤中污染物运移问题的全离散守恒混合元格式,讨论了守恒混合元解的存在唯一性,并给出了误差估计.最后给出了数值算例,数值模拟结果表明,用该方法模拟污染物运移问题是合理有效的.     

饱和多孔弹性杆流固耦合动力响应的保结构算法

研究了不可压饱和多孔弹性杆的流固耦合动力响应问题.基于多孔介质理论,根据多孔介质流固混合物动量方程、孔隙流体动量方程及体积分数方程,建立流固耦合不可压饱和多孔弹性杆的轴向振动方程;引入正则变量,构造饱和多孔弹性杆轴向振动方程的广义多辛保结构形式、广义多辛守恒律及广义多辛局部动量误差;采用中点Box离散方法得到轴向振动方程的广义多辛离散格式、广义多辛守恒律数值误差及局部动量数值误差;数值模拟不可压饱和多孔弹性杆的轴向振动过程及流相渗流速度分布,考察了流固两相耦合系数对轴向振动过程及广义多辛守恒律误差和局部动量误差的影响.结果表明,已构造的广义多辛保结构算法具有很高的精确性和长时间的数值稳定性.     

二维非饱和土壤水分运动问题的半离散有限体积元模拟

本文利用基于重心对偶剖分的有限体积元法建立了二维非饱和土壤水分运动问题的数值逼近格式,讨论了离散有限体积元解的存在唯一性,并给出了最优误差估计的证明.最后给出数值算例,模拟结果表明,利用有限体积元格式来求解二维非饱和土壤水分运动问题是可靠的,且该格式具有稳定性和可实用性.     

非饱和土壤水流方程的CN广义差分法

首先给出二维非饱和土壤水流方程时间二阶精度的Crank-Nicolson(CN)时间半离散化格式,然后直接从CN时间半离散化格式出发,建立具有时间二阶精度的全离散化CN广义差分格式,并给出误差分析,最后用数值例子验证全离散化CN广义差分格式的优越性.这种方法能提高时间离散的精度,极大地减少时间方向的迭代步,从而减少实际计算中截断误差的积累,提高计算精度和计算效率.而且该方法可以绕开对空间变量的半离散化广义差分格式的讨论,使得理论研究更简便.     

二维非饱和水流的全离散广义差分格式分析及其数值模拟

建立二维非饱和水流问题的全离散广义差分格式,讨论了全离散广义差分解的存在唯一性,并给出最优误差估计的证明．最后给出数值算例,验证方法的有效性.     

广义Boussinesq方程的多辛方法

广义Boussinesq方程作为一类重要的非线性方程有着许多有趣的性质,基于Hamilton空间体系的多辛理论研究了广义Boussinesq方程的数值解法,构造了一种等价于多辛Box格式的新隐式多辛格式,该格式满足多辛守恒律、局部能量守恒律和局部动量守恒律．对广义Boussinesq方程孤子解的数值模拟结果表明,该多辛离散格式具有较好的长时间数值稳定性．     

以中心- 迎风数值格式研究具有外力项的p 系统

本文以半离散中心- 迎风数值格式研究具有外力项的p 系统. 中心型数值格式用来处理双曲型守恒律或系统的优势是快速且简单, 因为不需要使用近似Riemann 解, 也不需要做特征分解. 我们的数值模拟验证了理论研究结果: 具有外力项的p 系统的解的收敛及爆破行为, 同时也指出一些尚待理论研究的问题.     

Mathematical Modeling to Simulate the Movement of Contaminants in Groundwater

The objectives of this paper are twofold. Firstly, we formulate a system of partial differential equations that models the contamination of groundwater due to migration of dissolved contaminants through unsaturated to saturated zone. A closed form solution using the singular perturbation techniques for the flow and solute transport equations in the unsaturated zone is obtained. Indeed, the solution can be used as a tool to verify the accuracy of numerical models of water flow and solute transport. The second part of this paper, deals with how the water level in a water reserve drops due to pumping water out of a well that is some distance away.     

非饱和水流问题的混合元法及其数值模拟

1.引 言 均质土壤中的地下水流动可归结为非饱和土壤水的流动,是土壤水未完全充满孔隙时的流动,是多孔介质流体运动的一种重要形式.非饱和流动的预报在大气科学、土壤学、农业     

Error estimates for an Euler implicit-mixed finite element scheme for reactive transport in saturated/unsaturated soil

We present convergence results for a fully discrete scheme based on the mixed finite element (MFE) method and an one-step Euler implicit (EI) method for simulating reactive solute transport in saturated/unsaturated soil. The results considered the low regularity of the solution of the degenerate parabolic equation describing the water flow in porous media. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

非饱和土壤水流问题基于POD 方法的降阶有限体积元格式及外推算法实现

利用特征投影分解(proper orthogonal decomposition, 简记为POD) 方法对非饱和土壤水流问题的经典有限体积元格式做降阶处理, 建立一种具有足够高精度维数较低的降阶有限体积元格式, 并给出这种降阶有限体积元解的误差估计和外推算法的实现, 最后用数值例子说明数值结果与理论结果是相吻合的. 进一步表明了基于POD 方法的降阶有限体积元格式对求解非饱和土壤水流问题数值解是可靠和有效的.     

Uncertainty analysis of transport of water and pesticide in an unsaturated layered soil profile using fuzzy set theory

Incomplete information is notoriously common in planning soil and groundwater remediation. For making decisions groundwater flow and transport models are commonly used. However, uncertainty in prediction arises due to imprecise information on flow and transport parameters like saturated/unsaturated hydraulic conductivity, water retention curve parameters, precipitation and evapo-transpiration rates as well as factors governing the fate of pollutant in soil like dispersion, diffusion, degradation and chemical transformation. Different methods exist for quantifying uncertainty, e.g. first and second order Taylor’s Series and Monte-Carlo method. In this paper, a methodology based on fuzzy set theory is presented to express imprecision of input data, in terms of fuzzy number, to quantify the uncertainty in prediction. The application of the fuzzy set theory is demonstrated through pesticide (endosulfan) transport in an unsaturated layered soil profile. The governing partial differential equation along with fuzzy inputs, results in a non-linear optimization problem. The solution gives complete membership functions for flow (suction head) and pesticide concentration in soil column.     

Generalized alternating-direction collocation methods for parabolic equations. II. Transport equations with application to seawater intrusion problems

The alternating-direction collocation (ADC) method combines the attractive computational features of a collocation spatial approximation and an alternating-direction time marching algorithm. The result is a very efficient solution procedure for parabolic partial differential equations. To date, the methodology has been formulated and demonstrated for second-order parabolic equations with insignificant first-order derivatives. However, when solving transport equations, significant first-order advection components are likely to be present. Therefore, in this paper, the ADC method is formulated and analyzed for the transport equation. The presence of first-order spatial derivatives leads to restrictions that are not present when only second-order derivatives appear in the governing equation. However, the method still appears to be applicable to a wide variety of transport systems. A formulation of the ADC algorithm for the nonlinear system of equations that describes density-dependent fluid flow and solute transport in porous media demonstrates this point. An example of seawater intrusion into coastal aquifers is solved to illustrate the applicability of the method. An alternating-direction collocation solution algorithm has been developed for the general transport equation. The procedure is analogous to that for the model parabolic equations considered by Celia and Pinder [2]. However, the presence of first-order spatial derivatives requires special attention in the ADC formulation and application. With proper implementation, the ADC procedure effectively combines the efficient equation formulation inherent in the collocation method with the efficient equation solving characteristics of alternating-direction time marching algorithms. To demonstrate the viability of the method for problems with complex velocity fields, the procedure was applied to the problem of density-dependent flow and contaminant transport in groundwaters. A standard example of seawater intrusion into coastal aquifers was solved to illustrate the applicability of the method and to demonstrate its potential use in practical problems.