土壤水流与溶质耦合运移问题的混合元-迎风广义差分法研究

本文研究了一维非饱和土壤水流与溶质耦合运移问题的数学模型, 建立了求其数值解的守恒混合元-迎风广义差分格式. 对非线性土壤水分入渗方程, 采用守恒混合元法进行离散模拟, 同时得到了土壤含水量和水分通量; 而对对流-扩散形式的溶质运移方程, 利用迎风的广义差分法离散求解. 且分析了解的存在唯一性, 并讨论了误差估计. 最后给出数值算例, 模拟结果表明利用本文格式来求解非饱和土壤水流与溶质耦合运移问题是可靠的, 且该格式具有稳定性和可实用性.     

蒸汽沉淀化学反应方程混合元法的离散格式研究

蒸汽沉淀化学反应过程有着极其广泛的应用,其数学模型归结为一个包含流速场,温度场,压力场和气体溶质场的非线性偏微分方程组．用混合有限元方法研究蒸汽沉淀化学反应方程组,导出其半离散化和全离散化的混合元格式,并证明这些格式的解的存在性和收敛性(误差估计)．用混合元法处理究蒸汽沉淀化学反应方程组,可以同时求出流速场,温度场,压力场和气体溶质场的数值解. 因此该研究既具有重要的理论意义,又具有广泛的应用前景．     

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

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

二维土壤溶质输运方程基于POD方法的降阶CN有限元外推算法

首先给出二维土壤溶质输运方程时间二阶精度的Crank-Nicolson(CN)时间半离散化格式和时间二阶精度的全离散化CN有限元格式及其误差分析.然后利用特征投影分解(proper orthogonal decomposition,简记为POD)方法对二维土壤溶质输运方程的经典CN有限元格式做降阶处理,建立一种具有足够高精度、自由度很少的降阶CN有限元外推格式,并给出这种降阶CN有限元解的误差估计和外推算法的实现.最后用数值例子说明数值结果与理论结果是相吻合的.     

污染物在非饱和带内运移的流固耦合数学模型及其渐近解

污染物在非饱和带中运移过程是多组分多相渗流问题.在考虑气相的存在对水相影响的前提下,基于流固耦合力学理论,建立了污染物在非饱和带内运移的流固耦合数学模型.对该强非线性数学模型采用摄动法及积分变换法进行拟解析求解,得出了解析表达式.对非饱和带内的孔隙压力分布、孔隙水流速以及污染物的浓度在耦合与非耦合气相条件下的分布规律进行解析计算.对该渐近解与Faust模型的计算结果进行了对比分析,结果表明:该模型解与Faust解基本吻合,且气相作用以及介质的变形对溶质的输运过程产生较大的影响,从而验证了解析表达式的正确性和实用性.这为定量化预报预测污染物在非饱和带中迁移转化和实验室确定压力-饱和度-渗透率三者之间的关系提供了可靠的理论依据.     

中空纤维透析器全血传质的双室模型

在本文中,我们建立了中空纤维透析器中全血传质的双室模型,并给出该模型的解析解;也给出了用尿素、肌酐和尿酸三种溶质进行实验验证的结果,从中发现传统的单室模型仅适于描述尿素在全血中的传质动力学,而对肌酐和尿酸则宜用本文的双室模型来描述。     

二维土壤溶质输运问题的CN有限元法

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

一类非线性奇摄动方程的匹配解法和解的精度估计

利用匹配渐近展开法,研究了一类非线性奇异摄动方程.在适当的条件下,得出了该类问题解的渐近展开式.并将结果应用于例子,对渐近解与精确解和用两变量方法求得的解进行比较,可知所得到的渐近解达到了较高精度.     

一类广义非线性强阻尼扰动发展方程的行波解

研究了一类非线性强阻尼广义扰动发展方程问题.它们在数学、力学、物理学等领域中广泛出现.首先,引入一个行波变换,把相应的偏微分方程问题转化为行波方程问题并求出原典型问题的精确解.再用小参数方法和引入伸长变量构造了问题的渐近解.最后, 用泛函分析的不动点理论证明了原非线性强阻尼广义扰动发展方程初值问题渐近行波解的存在性,并证明渐近解具有较高的精度和一致有效性.该文求得的渐近解是一个解析展开式, 所以它还可继续进行解析运算, 而单纯用数值模拟的方法是不行的.     

用区域分解法求不可压N-S方程的差分解

§1.引言 对不可压小粘性流的数值解,[1]和[2]用奇异摄动观点提出了一个区域分解法.从常微分方程(组)的奇异摄动问题出发,解分解为外部解加边界修正解(以下简称为修正解).外部解的边界条件有:给定(原边界条件)、待定(用原边界条件和修正解)和延拓类.修正解的边界条件有:给定(用原边界条件和外部解延拓)渐近(在边界层外缘)和待定     

The equimaterial approach for the numerical solution of the one dimensional transient diffusion equation with zero flux boundary conditions

A new approach is proposed for the grid motion for the numerical solution of a general transient diffusion equation in one spatial dimension with zero flux boundary conditions. The new criterion for grid motion is that the solute amount contained in each discretization section should be a pre-described fraction of the total solute amount at each time step. This requirement is not explicitly enforced to the solution technique but it is implicitly included in the equation through the appropriate variable transformation. The results showed that although the technique leads to the required grid motion the numerical results are of pure quality due to the appearance of singularities during the variable transformation procedure. Nevertheless, it is shown that by appropriate numerical handling of the solution at the singularity region the technique can lead to accurate results and potentially can replace the existing moving grid algorithms at least for the particular problem at hand.     

Effect of chemical reaction on the dispersion of a solute in a porous medium

A mathematical model is presented in this paper which describes the dispersion of a chemically active solute in the laminar flow in a sparsely packed porous medium. The validity of time-dependent dispersion coefficient is widened by using a generalized dispersion coefficient. The effect of porous parameter and chemical reaction on the dispersion coefficient is studied. The exact solution for the mean concentration distribution of a chemically active solute is obtained as a function of downwind distance and time. Results are also obtained for pure convection.     

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.     

Modeling of solute transport into sub-aqueous sediments

A model for investigating the solute transport into a sub-aqueous sediment bed, under an imposed standing water surface wave, is developed. Under the assumption of Darcy flow in the bed, a model based on a two-dimensional, unsteady advection–diffusion equation is derived; the relative roles of the advective and diffusive transport are characterized by a Peclet number, Pe. Two solutions for the equation are developed. The first is a basic control volume method using the power-law scheme. The second is a smear-free, modified upwind solution for the special case of Pe → ∞. Results, at a given time step, are reported in terms of a laterally averaged solute verse depth profile. The main result of the paper is to demonstrate that the one-dimensional solute concentration verse depth profile is essentially independent of any numerical dissipation present in the solute field predictions. This demonstration is achieved by (i) using an extensive grid refinement study, and (ii) by comparing Pe → ∞ predictions obtained with the basic and smear-free solutions.     

Numerical stability of the BEM for advection‐diffusion problems

A boundary element method (BEM) approach has been developed to solve the time‐dependent 1D advection‐diffusion equation. The 1D solution is part of a 3D numerical scheme for solving advection‐diffusion (AD) problems in fractured porous media. The full 3D scheme includes a 3D solution for the porous matrix, which is coupled with a 2D solution for fractures and a 1D solution for fracture intersections. As the hydraulic conductivity of the fracture intersections is usually higher than the hydraulic conductivity of the fractures and by at least one order of magnitude higher than the hydraulic conductivity of the porous matrix, the fastest flow and solute transport occurs in the fracture intersections. Therefore it is important to have an accurate and stable 1D solution of the transient AD problems. This article presents two different 1D BEM formulations for solution of the AD problems. The particular advantage of these formulations is that they provide one of the most straightforward and simplest ways to couple multiple intersecting 2D Boundary Element problems discretized with linear discontinuous elements. Both formulations are tested and compared for accuracy, stability, and consistency. The analysis helps to select the more suitable formulations according to the properties of the problem under consideration. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Wavelet based multiscale scheme for two-dimensional advection–dispersion equation

In this paper, wavelet based adaptive solver is developed for two dimensional advection dominating solute problem which generates sharp concentration front in the solution. In order to handle simultaneously smooth and shock-like behavior, the framework uses finite element discretization followed by wavelets for multiscale decomposition. Daubechies wavelet filter is incorporated to eliminate spurious oscillations at very high Peclet number. The developed solution is compared with the analytical solution to assess the accuracy and robustness. The advantages of the present method over the commonly used methods such as FDM and FEM for solving the problems which show non-physical oscillation in the numerical solution are demonstrated.     

ONE‐DIMENSIONAL TEMPORALLY DEPENDENT ADVECTION–DISPERSION EQUATION IN POROUS MEDIA: ANALYTICAL SOLUTION

Abstract Analytical solutions of one‐dimensional advection–dispersion equation in semi‐infinite longitudinal porous domain are obtained in this work. The solute dispersion parameter is considered temporally dependent along uniform flow. The first‐order decay term, which is inversely proportional to the dispersion coefficient, is also considered. Initially, the space domain is not solute free. Analytical solutions are obtained for uniform and varying pulse‐type input. A new time variable is introduced. The Laplace transform technique is used to get the analytical solutions.     

A two-dimensional finite volume method for transient simulation of time- and scale-dependent transport in heterogeneous aquifer systems

In this paper, solute transport in heterogeneous aquifers using a modified Fokker-Plauck equation (MFPE) is investigated. This newly developed mathematical model is characterised with a time-, scale-dependent dispersivity. A two-dimensional finite volume quadrilateral mesh method (FVQMM) based on a quadrilateral background interpolation mesh is developed for analysing the model. The FVQMM transforms the coupled non-linear partial differential equations into a system of differential equations, which is solved using backward differentiation formulae of order one through five in order to advance the solution in time. Three examples are presented to demonstrate the model verification and utility. Henry's classic benchmark problem is used to show that the MFPE captures significant features of transport phenomena in heterogeneous porous media including enhanced transport of salt in the upper layer due to its parameters that represent the dependence of transport processes on scale and time. The time and scale effects are investigated. Numerical results are compared with published results on the same problems.     

Modeling and analysis of reactive solute transport in deformable channels with wall adsorption–desorption

We show well posedness for a model of nonlinear reactive transport of chemical in a deformable channel. The channel walls deform due to fluid–structure interaction between an unsteady flow of an incompressible, viscous fluid inside the channel and elastic channel walls. Chemical solutes, which are dissolved in the viscous, incompressible fluid, satisfy a convection–diffusion equation in the bulk fluid, while on the deforming walls, the solutes undergo nonlinear adsorption–desorption physico‐chemical reactions. The problem addresses scenarios that arise, for example, in studies of drug transport in blood vessels. We show the existence of a unique weak solution with solute concentrations that are non‐negative for all times. The analysis of the problem is carried out in the context of semi‐linear parabolic PDEs on moving domains. The arbitrary Lagrangian–Eulerian approach is used to address the domain movement, and the Galerkin method with the Picard–Lindelöf theorem is used to prove existence and uniqueness of approximate solutions. Energy estimates combined with the compactness arguments based on the Aubin–Lions lemma are used to prove convergence of the approximating sequences to the unique weak solution of the problem. It is shown that the solution satisfies the positivity property, that is, that the density of the solute remains non‐negative at all times, as long as the prescribed fluid domain motion is ‘reasonable’. This is the first well‐posedness result for reactive transport problems defined on moving domains of this type. Copyright © 2015 John Wiley & Sons, Ltd.