溶质运移模型的有限元数值解
对溶质运移模型的并行化处理
第32卷第1期吉首大学学报(自然科学版)Vol.32No .12011年1月Journ al of Ji shou Universit y (Nat ural Science Edit ion)J an.2011文章编号:1007-2985(2011)01-0026-04对溶质运移模型的并行化处理*邓永辉(湖南财政经济学院,湖南长沙410205)摘要:在了解和广泛收集前人研究资料的基础上,将并行算法引入到求解首采区卤水动态二维模型中关于溶质运移的问题中.把溶质运移方程按时间分裂法分成5个子方程,针对这5个子方程讨论了全域精细积分和子域精细积分的并行算法,给出了对流和扩散方程的子域精细积分并行算法.子域精细积分考虑了精细积分法高精度的特点,又避免了全域积分的大矩阵运算,其精度优于单点精细积分法.关键词:对流扩散方程;并行算法;偏微分方程中图分类号:O29文献标志码:A溶质运移模型的数值求解法主要包括有限元和有限差分法2种.由于有限差分法计算简单,物理意义直观,在溶质运移问题求解中得到普遍应用.但是,当对流项占优,即Peclet 数较高时,则出现数值波动和弥散.采用高阶的时间导数和空间导数离散逼近,并用极短的时间和空间步长求解,上游加权法可同时消除大部分的数值弥散,但求解效率却大大降低.总之,对差分法的改进,或消除了数值波动,或减少了数值弥散,仍然未能从根本上消除数值困难.与有限差分法相比,有限单元法比较灵活[1],求解精度高,因而成为早期求解溶质运移问题应用最广泛的数值法.1溶质运移数学模型通过几种简单模型的推导和分析,再结合首采区能引起浓度变化的各种因素,如对流、扩散、排泄、补给等等,可以导出最终目标模型如下:M(C t+V 1C x+V 2C y )=M (C x(D 11C x)+y (D 22C y))+Q(C q -C)+M f c.(1)(1)式只表明了浓度随时间变化的规律,要求出微分方程的解还需要一些定解条件,求出在特定条件下浓度的值,再计算区内边界条件如下:C(x,y,t)|t=0=C 0(x ,y),C(x,y,t)|1=C 1(x,y),M D n Cn |n=q c 2(x,y),12=,12=.其中:C 为浓度(mL -3);V 1,V 2为渗透速度(LT -1);D 11,D 22为弥散系数(L 2T -1);C q 为汇源项的浓度(mL -3);n 为外法线方向;f 为固液转化系数(T -1),若忽略固液转化作用,则f =0;C 0(x,y)为初始浓度分布函数(mL -3);C 1(x,y)为第1类边界1上的浓度分布函数(mL -3);q c 2(x,y)为第2类边界2上弥散通量分布函数(mT -1L -1);为渗流区;M 为含水层厚度(L);Q 为汇源补给项(LT -1);自由孔隙度率,即给水度,也就是从单位体积含水层中在重力作用下能够释放给出的水量所占该含水层体积的份数.*收稿日期基金项目国家科技攻关项目(B 6ZB )作者简介邓永辉(),女,湖南常德人,湖南财政经济学院讲师,硕士,主要从事应用数学领域研究:2010-11-24:2001A 0-09-1:1979-.2并行前的预处理多维的对流扩散方程,尤其当它为非线性和非定常的时候很难直接求解,而且在差分格式的选择上还有很大的局限性,即使得出结果,解的精度也不高.为了剔除这些不利因素,Yanenko 提出了一种利用时间分裂的概念[2],分步求解多维问题的方法.所谓时间分裂,笔者是这样理解的,即在某一固定的时间段内,一个物理或是化学过程中的若干反应所引起的效力,并不是同时完成的,而是存在主次先后之分.换句话说,在某一固定的时间段里,可以将一个物理或化学过程的每个子反应看成是互步相交的集合,而它们的并形成了这个时间段内完整的物理或化学过程的演变.因此,任意多维的问题都可以分解成若干个一维的偏微分方程,依次求解.对于方程(1),不步考虑排泄补给项和溶质的化学转换等,方程化为(C t +V 1C x +V 2C y)=(C x (D 11C x)+y(D 22C y))+S.(2)根据分步解法,方程(2)可以分裂成5个方程:C t=R r =1Q r (z -z r )(y -y r )(x -x r ),(3)C t =-V 1C x ,(4)C t =-V 2C x ,(5)C t =D 112C x 2,(6)Ct=D 112C y 2.(7)其中:()为Dirac 函数;Q r 为第r 个源头的源强;x r ,y r 和z r 为这个源的坐标.方程(3)为R 个源头对浓局度部变化C t 的贡献,(4)和(5)式为对流项对浓度局部变化Ct 的贡献,而(6)和(7)式为水平方向的扩散项对浓度局部变化C t的贡献.在对对流扩散方程采用时间分裂法,将多维的问题转换为5个一维的方程求解时,很大程度上增强了已知算法的可操作性.下面来研究一维扩散方程和一维对流扩散方程的并行化处理.3并行处理偏微分方程的精细时程积分法先在空间坐标上按差分法离散[3],而得到对时间的常微分方程组,然后采用全域精细积分法计算此方程组,可得到方程组的高精度解.这相当于给出了偏微分方程的半解析解.精细积分中指数矩阵T 的计算量很大,而T 的求解主要是矩阵乘法的运算.全域精细积分中,需要对矩阵运算进行并行化.文中使用的T ransputer 系统是分布式存储的多处理器系统,采用分块的矩阵乘法.变量较多时,全域精细积分的指数矩阵的存储量和计算量都很庞大.子域精细积分不需形成总体的指数矩阵,各子域的求解是相对独立的,计算量较大的指数矩阵求解可在各处理器中分别计算,它具有良好的并行计算的性质.3.1一维扩散方程的并行算法(1)全域精细积分的并行算法.考虑一维扩散方程C t =D 2Cx 2,当两端边界条件为0,用差分法在x 坐标方向离散后,可得到k +1=exp(H t)=exp(H t)k=Tk,(8)其中是常数矩阵,0是初值文献[]中给出了指数矩阵T 的精细计算方法T =x ()=[x ()]=[x (T)]=[I +T ],其中可选用=N 若是一不大的时间区段,则当N 并不很大时,=也将是一个非常小的时间27第1期邓永辉:对溶质运移模型的并行化处理H .4:e p H t e p H t/m m e p H m a m m 2.t t/m段.T 的计算相当于做N 次循环:T a =2T a +T a T a .(9)最后做T =I +T a ,这样可得到高精度的指数矩阵T =exp(Ht).Transputer 系统是分布式存储的并行多处理器系统,各处理器的数据由相互间的通讯得到.Fortr an 语言中数组的数据是按列存放的.为了通讯的方便和高效率,(9)式中矩阵的计算按列分块T a =(T a1,T a 2,,T an ).n 是处理器的数目.处理器i 中的矩阵计算为T ai =2T a i +T aT ai .由于相乘和相加是相同的T a 矩阵,因此计算时先将T a 阵分别送到各个处理器,各处理器计算T a i ;然后进行通讯得到新的完整的T a 矩阵;再将T a 阵送给各处理器,依次循环后,加上单位阵I 得最终的矩阵T.(8)式的时间步计算中,将T 阵按行分块T =(T 1,T 2,,T n ),k +1则分为n 段k +1=((k +11)(k +12)(k +1n)),k +1i=T ik,i =1,2,,n.由通讯收集得到k +1,再送给各处理器进行下一时间段的计算.(2)子域精细积分的并行算法.当未知数增多时,全域精细积分法指数矩阵的存储量和计算量都将迅速增加,这会影响到其解题的规模和效率.当时间步长t 很小时,指数矩阵T 是近似窄带宽的.因此,文献[5]提出了子域精细积分的方法,这可以使存储量和计算量大幅度下降.考虑子域的非齐次微分方程组,其向前精细积分格式为v n +1j=Tv n j +(T -I )H -1f njj =1,2,,J.(10)蛙跳格式为T(2t)=exp(2Ht),v n -1j=Tv n j +(T -I )H -1f n j j =1,2,,J.(11)格式(10)是一阶时间精度,格式(11)则是二阶时间精度,而两者的计算量几乎是一样的.当各子域划分相同时,只需计算1组T 和H -1阵.T 阵的并行算法与全域精细积分的相同,H -1阵的计算采用分块的Gauss Jordon 算法.将H 按列分块送到各处理器:H =(H 1,H 2,,H n).设有m 个内点.从H 阵的第1列到最后的第m 列依次计算.计算第j 列时,j 列所在的第I 个处理器从第j列选主元.第I 个处理器将选出的主元行号k 及第j 列向量a 送到每个处理器.各处理器交换H 阵中的第j,k 2行:tmp =H ijt ,H ijt =H ik l ,H ik l =tmpl =1,,m;i =1,,n.选出的主元行号存储于一维数组M 中.然后进行消去运算.设H 阵的第j 列在H i 阵中为第p 列,则第I 个处理器中的算法为:H Ijp =1/a j ,H It p =-a l /a j l =1,,j -1,j +1,,m;H I jl =H I jl /a jl =1,,p -1,p +1,,m I ;H Iq l =H Iq l -H Ijl a qq =1,,j -1,j +1,,m,l =1,,p -1,p +1,,m I .其他处理器中的算法为:H I jl =H I jl /a jl =1,,m i ,i =1,,n;H Iq l =H Iq l -H Ijl a qq =1,,j -1,j +1,,m,l =1,,m i ,i =1,,n.H 阵是按列存在各处理器中,交换的通讯量很大,可将各处理器中的H i 送到主控程序[4],由主控程序进行最后的列交换,得到逆阵H -1.再将T 和H -1阵送回各处理器.向量按照子域划分情况送到各处理器,按(10)或(11)式计算各子域的v n +1j .将v n +1j 中与其他子域重合部分的分量送出,同时接收其他处理器计算的这几个分量,并取平均值作为计算结果.3.2对流扩散方程的并行算法考虑更一般的一维对流扩散方程=D x x,()为改进差分法的数值精度和稳定性,文献[]给出了对流扩散方程的单点精细积分算法由方程()的28吉首大学学报(自然科学版)第32卷ut2u2-v u 124.122个通解u =常数和u =x -v t ,得到在空间坐标上离散的微分方程[5].u 0j =c 1u j-1-c 2u j +c 3u j+1.(13)(13)式的单点蛙跳格式为u n +1j =[u n -1j -(c 3u n j+1+c 1u n j -1)/c 2]exp(-2c 2t)+(c 3u n j+1+c 1u nj -1)/c 2.(14)此格式具有二阶时间精度,并且是无条件稳定的.为提高精度,可以采用多内点的子域精细积分算法.利用蛙跳格式(14),考虑3个内点的子域,其非齐次微分方程组是u 0j -1u 0j u 0j +1=-c 2c 30c 1-c 2c 30c 1-c 2u j -1u jj j+1+c 2u j-2c 3u j+2j =1,2,,J.(15)u j-2和u j +2作为常数,在t 时间段内是不变的,可以从相邻子域中得到这2点的值.采用与(10),(11)式相同的格式计算方程(15).当在空间坐标作等距离的离散,并且各子域内点相同时,T j 和H -1j 是相同的,则只需计算1组[6];算法与前相同.使用不同划分的子域时,先将向量x 按划分的子域分段地送到各处理器,各处理器计算各系数,并组装得到H j 矩阵;然后计算H -1j和T j .这一步计算量较大,且不需通讯,是大粒度的.4结语精细积分法有很高的计算精度,不论是全域积分还是子域积分都有很好的并行性,在并行机系统上可充分发挥其效能.其中子域积分法有更大的优点,它克服了全域积分矩阵过大的弱点,又可以提高计算的精度.在并行计算中能方便地将各子域放到各处理器中,算例表明此算法是高效率的.增加内点数加大了指数矩阵和每个时间步的计算量,但能够提高精度,所以应适当地选择子域的大小.参考文献:[1]王文科.地下水有限分析数值模拟的理论与方法[M].西安:陕西科学技术出版社,1996:102-126.[2]ROACH E P putational Fluid Dynam ics [M].H ermosa,Albuquer que,New Mcxico,1976:446.[3]孙纳正.地下水污染数学模型和数值方法[M].北京:北京地质出版社,1989.[4]钟万勰.子域精细积分及偏微分方程数值解[J].计算结构力学及其应用,1995,12(3):253-260.[5]罗焕炎,陈雨孙.地下水运动的数值模拟[M].北京:中国建筑工业出版社,2001.[6]DURBIN F.Numer ical Inver sion of Laplace Tr ansfor m:An Efficient Improvement to Dur bner and Aba re s Method [J].Comp.J,1973,17:371-376.Parallel Algorithms for a Solute Transport ModelDENG Yong hui(H unan Univer sity of Finance and Economics,Changsha 410205,China)Abstr act:This paper pr esents the parallel algorithms of precise integration methods.The subdomain par allel algorithms of pr ecise integration of convection diffusion equations are given.T he method take advan tage of the high precision of precise integration methods,and they also avoid a large amount of matr ices calculation of whole domain methods.The precision of the methods is better than that of the single point integration.Key words:convection diffusion equations;parallel algorithm;partial differential equations(责任编辑向阳洁)29第1期邓永辉:对溶质运移模型的并行化处理。
一维变密度溶质运移实验及参数推求
第24卷第3期2008年5月水资源保护W ATER RES OURCES PROTECTI ON V ol.24N o.3May 2008 基金项目:国家自然科学基金(50679025);教育部新世纪优秀人才支持计划(NCET Ο04Ο0492);高等学校学科创新引智计划(B08048)作者简介:马建良(1982—),男,山东乐陵人,硕士研究生,研究方向为地下水数值模拟。
E 2mail :maliang @ 一维变密度溶质运移实验及参数推求马建良1,陈 喜1,程勤波2,宋 轩2,鲍振鑫2(1.河海大学水文水资源与水利工程科学国家重点实验室,江苏南京 210098;2.河海大学水文水资源学院,江苏南京 210098)摘要:通过室内土柱注水实验,观测沙质土壤中氯离子浓度的变化过程。
以变密度水流连续性方程、溶质运移方程和达西方程为基础,运用有限单元法和差分法对这3个方程进行联立求解,建立了一维变密度水流和溶质运移数值模型。
利用实验数据反求变密度渗透系数、弥散系数等水动力参数。
关键词:变密度水流;一维对流—弥散方程;土柱实验;海水入侵中图分类号:O351.2 文献标识码:A 文章编号:1004Ο6933(2008)03Ο0008Ο04Identification of hydrodynamic parameters based on one 2dimensional variable density and solute transport numerical modelMA Jian 2liang 1,CHEN Xi 1,CHENG Q ing 2bo 1,SONG Xuan 2,BAO Zhen 2xin 2(1.State K ey Laboratory o f Hydrology 2Water Resources and Hydraulic Engineering ,Hohai Univer sity ,Nanjing 210098,China ;2.College o f Hydrology and Water Resources ,Hohai Univer sity ,Nanjing 210098,China )Abstract :The concentration of chloride ions in sandy s oil was observed through water 2filling s oil column tests in the laboratory.A one 2dimensional numerical m odel for groundwater f1ow of variable density and s olute transport was developed on the basis of equations of variable density groundwater flow and s olute transport as well as Darcy ’s Law.The equations were s olved using the finite element and finite difference methods.The hydrodynamic parameters for variable density flow ,such as infiltration coefficients and dispersion coefficients ,were calibrated against the observed data.K ey w ords :variable density groundwater f1ow ;one 2dimensional convection 2diffusion equation ;s oil column experiment ;seawater intrusion 20世纪70年代以来我国沿海地区陆续出现海水入侵。
地下水溶质运移数值模型
地下水溶质运移数值模型(资料性附录)水是溶质运移的载体，地下水溶质运移数值模拟宜在地下水流场模拟基础上，因此地下水溶质运移数值模型包括水流模型和溶质运移模型两部分。
DJ 地下水水流模型非均质、各向异性、空间三维结构、非稳定地下水流系统:1)控制方程σ∂h ∂hy 3(“∂h}∂(∂h ∖S,—=—K v —+—K Y —+—K ——+/∂t 3xI ∂x)为('∂y JAzI ~∂z)式中：SS 一一给水度［I/］；h --- 水位［1］；Kχf Ky,Kz ——分别为X,y,Z 方向上的渗透系数［EΓ∣］;T 一一时间［T ］；Qs 一一源汇项m注：方括号［］中的符号为量纲，以下同。
2)初始条件h(x y y 9z y t)=Zz 0(x,y,z)(x,y,z)∈Ω,/=O 式中：4*，y ，z)——已知水位分布：Q ——模型模拟区。
3)边界条件：第一类边界： 〃(x,y,z√)∣「=Λ(x,y,z√)(x,y,z)∈Γ1,r≥O式中：r '一一类边界； h(x,y,z,t)一一类边界上的己知水位函数。
第二类边界：式中：「2 --- 二类边界;∂nq(x,y,Z) (x,y,z)∈Γ2κ——三维空间上的渗透系数张量;nn——边界r2的外法线方向；q(x,y t z)——二类边界上已知流量函数。
第三类边界：r(k(h-z)-+ah)=q(x,y,z)加r3式中：0一一系数；「3一—二类边界；k一一三维空间上的渗透系数张量；n——边界G的外法线方向;q(x,y f z)——三类边界上已知流量函数。
D.2地下水水质模型1)控制方程R啜喘［吗(他C)Fe—/〜元式中：R——迟滞系数，无量纲Pb SC~Θ~∂Cph——介质密度IM1-3］；θ——介质孔隙度，无量纲；C——组分的浓度［M1，］；亍一一介质骨架吸附的溶质浓度［M1，］；t——时间［T］；X,y,Z一—空间位置坐标［1］；Dij——水动力弥散系数张量［1?T」］；Vi——地下水渗流速度张量［EΓ∣］;q s——源和汇［T∣］;CJ一一源或汇水流中组分的浓度［M1"；4一一溶解相一级反应速率［T」］；4一一吸附相一级反应速率［Tj］。
变异Henry问题检验网格Peclet数在数值计算中的重要性
变异Henry问题检验网格Peclet数在数值计算中的重要性马倩;李海龙【摘要】有限元法是求解地下水流和溶质运移对流-弥散方程的常用数值方法,它可以精确高效地处理以弥散为主的问题,但求解以对流为主的问题易引起显著的数值振荡.通过Galerkin有限元法对变异Henry问题进行模拟求解,得到了用不同的剖分网格及水动力弥散系数时,在特选节点处的浓度穿透曲线,分析并找到了浓度振荡的原因及合适的消除方法,即若出现浓度数值解在某值附近振荡,可以通过加密网格或增加水动力弥散系数将其消除.模拟结果及其分析表明:即使是研究区域相同,不同的边界条件、不同的水动力弥散系数对网格精度的要求不同;换言之,同一网格对不同模型参数的有效性也不同.网格Peclet数能够有效地判定给定的网格剖分是否会引起浓度振荡,对有限元法数值计算的网格剖分具有指导意义.【期刊名称】《水文地质工程地质》【年(卷),期】2014(041)003【总页数】5页(P1-5)【关键词】浓度振荡;有限元法;变异Henry问题;Peclet数;网格剖分【作者】马倩;李海龙【作者单位】生物地质与环境地质国家重点实验室/中国地质大学(北京)水资源与环境学院,北京 100083;地下水循环与环境演化教育部重点实验室/中国地质大学(北京),北京 100083;生物地质与环境地质国家重点实验室/中国地质大学(北京)水资源与环境学院,北京 100083;地下水循环与环境演化教育部重点实验室/中国地质大学(北京),北京 100083【正文语种】中文【中图分类】P641.2有限元法作为求解地下水流和溶质运移对流-弥散方程的常用数值方法，在水流模拟中有很多优点。
它有固定的网格，通常满足质量守恒定律;可以精确、高效地处理以弥散为主的问题;便于编写程序并得以实现等［1］。
但是，对于许多野外常见的以对流为主的问题，求解易引起显著的数值振荡。
为此，许多学者进行了大量的计算实践和理论分析指出［2］:对于以对流占优的溶质运移问题，有限差分法和有限元法计算时皆不可避免地出现数值振荡。
溶质运移模型数值分析
＋
ａ『
的各种定解问题来描述（ 参见文献 ［，］ ， 中 Ａ 为河床断面， 为流量）Ｄ 为弥散系数 ， １２ 等）其 Ｑ ， Ｃ是某组分子 在 ．断面处 时刻的浓度， ２ ｇ Ｓ表示各种源 、 汇项的代数和 在复杂情况下 ， Ｓ是 Ｃ的线性或非线性函数 ， 但在
河
南
科
学
第２卷 第 ３ Ｏ 期
本文在§１ 试图利用积分变换方法， 中， 寻求理想化条件下的简化模型（．） ０２的近似精确解（ 即理论上的近似解 析解） 的表达式。在 § 中， ２ 首先建立模型（．） ０１的古典差分格式， 给出模型（． 和（．） ０１ ０２的数值算例。 ）
１ 模型 《 ． ） ０ ２ 的积 分变换方法和近似解 析解
上式化 为
ｆ “ Ｄ嚣， ＋ ＝
｛（户 。（ｔ－ ｂ ， ０ ） ｃ ，ｅ出＝／ ， Ｊ ｏ）ｐ ｔ ｐ （３ １） ・
｝。户 ｉ ｃ。） ・ 。 ） （， ｔ （， Ｊ 。ｔ ０ ｄ
常微 分方程
Ｄ
嚣一 ＝ “ 一 ｏ
的通解 为
一
）＝
其中 ｅ ｃ ＝１ ， ｝ ， ｒ（ ｆ ） 一 ｒ（ ）而为误差函数。 证 明： ０２ 两端 同乘以 ， 将（ ． ） 并在 ０ ｏ 到ｏ范围内对时间 ｔ 求积分 ， 得到
＝
等 “
（ 一“
．
ｃ（ ．
ｚ
由（ ．） １２得
户 令 ＝ Ｄ
（户＝ｏ￣ ）， ）ｆ Ｃ ， ， （ｄ ｅ
两类模型在生物化学基因控制 、 口动力学系统和多孔介质物理学等方面都有许多应用背景 。因此， 人 已有许 多文献对它们作了定性分析和数值方法上的研究（ 见文献 ［ —６ 及其相关引用文献 等） 然而， ４ ］ 作为溶质运 移模型， 当地下水水流速度 Ｉ Ｉ “ 很小（ 亦即弥散项起主要作用） ， 时 可采用通常的差分方法进行数值计算。
二维地下水溶质运移耦合方程组的差分方法
中国海洋大学硕士学位论文二维地下水溶质运移耦合方程组的差分方法姓名：***申请学位级别：硕士专业：计算数学指导教师：***201206二维地下水溶质运移耦合方程组的差分方法摘要本文分析二维地下水溶质运移数学模型的控制方程，在此基础上，建立了求解控制方程的交替方向隐式差分格式。
古典显格式由于稳定性，在时间上限制很大；隐式差分格式和中心差分法需要求解多个未知量，导致计算时占用内存大，求解困难；本文采用的交替方向差分法，每层分两步计算，每步只需求解一个三对角方程组，在计算规模上相当于一维的计算量，并且它与隐式格式有相同的稳定性。
其次控制方程含有非线性系数，为了高精度地模拟模型，本文对非线性系数用三次样条插值逼近，提高了计算精度。
通过数值算例比较证明本文方法可以对地下水溶质运移进行有效的模拟，对解决此类非线性问题有参考价值。
最后分析了线性的溶质运移方程的稳定性。
文章由三部分组成：第一部分为前言。
介绍了地下水模型和求解地下水问题的数值方法的背景和发展状况，以及本文的主要研究内容。
第二部分为地下水溶质运移模型分析。
在这一部分里，介绍了地下水基本假设和根据均衡单元体的质量守恒建立地下水运动的连续方程和溶质运移方程。
再根据二维潜水问题的Ｄｕｐｕｉｔ假设，最终形成二维问题的模型。
还对方程中涉及到的参数做了详细的分析。
第三部分为二维地下水模型的有限差分方法。
给出方程的交替方向差分格式和二维时的三次样条插值计算。
对于方程非线性系数我们介绍了两种方法进行逼近，通过数值算例证明两层平均的三次Ｈｅｒｍｉｔｅ插值逼近更有效。
最后分析了线性的溶质运移方程的稳定性。
关键词：地下水连续方程；溶质运移方程；交替方向差分方法；样条插值；稳定性；二维地下水溶质运移耦合方程组的差分方法ＡＤＩｄｉ如ｒｅｎｃｅｍｅｔｈｏｄｆｏｒＣＩｌｍｅｎＳｌｏｎａＩｇｒｏｕｎｄｗａｔｅｒｃｏｕｐｌｅｄｓｙｓｔｅｍｏｆｔｗｏａｎｄｓｏｌｕｔｅｔｒａｎｓｐｏｒｔＡｂｓｔｒａｃｔＴｈｉｓｐａｐｅｒａｎａｌｙｚｅｓｔｈｅｃｏｎｔｒｏｌｅｑｕａｔｉｏｎｓｏｆｔｗｏ－ｄｉｍｅｎｓｉｏｎａｌｇｒｏｕｎｄｗａｔｅｒｓｏｌｕｔｅｔｒａｎｓｐｏｒｔ．Ｏｎｔｈｉｓｂａｓｉｓ，ｗｅｅｓｔａｂｌｉｓｈａｌｔｅｒｎａｔｉｎｇｄｉｒｅｃｔｉｏｎｉｍｐｌｉｃｉｔｄｉｆｆｅｒ—ｅｎｃｅｓｃｈｅｍｅｏｆｅｑｕａｔｉｏｎｓ．Ｃｌａｓｓｉｃａｌｅｘｐｌｉｃｉｔｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｄｕｅｓｔｏｔｈｅｓｔａｂｉｌｉｔｙｏｆｔｈｅｔｉｍｅｌｉｍｉｔ；Ｔｈｅｉｍｐｌｉｃｉｔｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｏｒｃｅｎｔｒａｌｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｒｅ—ｑｕｉｒｅｓｔｈｅｓｏｌｕｔｉｏｎｏｆｑｕｉｔｅａｌｏｔｏｆｕｎｋｎｏｗｎｖａｒｉａｂｌｅｓ，ａｎｄｔｈｅｎｌｅａｄｓｔｏｓｏｌｖｉｎｇｄｉｆｆｉｃｕｌｔｌｙ；Ｉｎｔｈｉｓｐａｐｅｒ：ｔｈｅａｌｔｅｒｎａｔｉｎｇｄｉｒｅｃｔｉｏｎｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｉｓｕｓｅｄ．Ｅａｃｈｌａｖｅｒｉｓｄｉｖｉｄｅｄｉｎｔｏｔｗｏｓｔｅｐｓａｓｗｅｌｌａ８ｅｖｅｒｙｓｔｅｐｓｉｍｐｌｙｎｅｅｄｓｔｏｓｏｌｖｅｓａｔｒｉｄｉａｇｏｎａｌｅｑｕａｔｉｏｎｓ．Ｉｎｔｈｅｓｃａｌｅｏｆｔｈｅｃａｌｃｕｌａｔｉｏｎ，ｉｔｉｓｅｑｕｉｖａｌｅｎｔｔｏｏｎｅ—ｄｉｍｅｎｓｉｏｎａｌ．Ａｎｄｈａｓｔｈｅｓａｍｅｓｔａｂｉｌｉｔｙｗｉｔｈｔｈｅｉｍｐｌｉｃｉｔｓｃｈｅｍｅ．Ｓｅｃｏｎｄｌｙ，ｔｈｅｃｏｎｔｒｏｌｅｑｕａｔｉｏｎｃｏｎｔａｉｎｓｔｈｅｎｏｎｌｉｎｅａｒｃｏｅｆｆｉｃｉｅｎｔ，ｔｈｅｃｕｂｉｃｓｐｌｉｎｅｉｎｔｅｒｐｏｌａｔｉｏｎｉｓｕｓｅｄｔｏａｐｐｒｏｘｉｍａｔｅｉｔｉｎｏｒｄｅｒｔｏｉｍｐｒｏｖｅｔｈｅｃａｌｃｕｌａｔｉｏｎａｃｃｕｒａｃｙ．Ｔｈｅｎｕｍｅｒ—ｉｃａｌｅｘａｍｐｌｅｓｈａｓｐｒｏｖｅｄｔｈｉｓｍｅｔｈｏｄｃａｎｂｅｅｆｆｅｃｔｉｖｅｔｏｓｉｍｕｌａｔｅｔｈｅｅｑｕａｔｉｏｎｏｆｓｏｌｕｔｅｔｒａｎｓｐｏｒｔｉｎｇｒｏｕｎｄｗａｔｅｒ．Ｔｈｉｓａｒｔｉｃｌｅｈａｓｒｅｆｅｒｅｎｃｅｖａｌｕｅｔｏｓｏｌｖｅｓｕｃｈｎｏｎｌｉｎｅａｒｉｓｓｕｅ．Ｆｉｎａｌｌｙ，ｗｅａｎａｌｙｚｅｔｈｅｓｔａｂｉｌｉｔｙｏｆｔｈｅｌｉｎｅａｒｓｏｌｕｔｅｔｒａｎｓｐｏｒｔｅｑｕａｔｉｏｎ．Ｔｈｅａｒｔｉｃｌｅｃｏｎｓｉｓｔｓｏｆｔｈｒｅｅｐａｒｔｓ：Ｔｈｅｆｉｒｓｔｐａｒｔｉｓｔｈｅ－ｐｒｅｆａｃｅ．Ｉｔｄｅｓｃｒｉｂｅｓｔｈｅｂａｃｋｇｒｏｕｎｄａｎｄｄｅｖｅｌｏｐｍｅｎｔｏｆｎｕｍｅｒｉｃａｌｍｏｄｅｌｉｎｇｒｏｕｎｄｗａｔｅｒ．Ｉｎａｄｄｉｔｉｏｎ，ｔｈｉｓｓｅｃｔｉｏｎｇｉｖｅｓｔｈｅｍａｉｎｃｏｎｔｅｎｔｓｏｆｔｈｉｓａｒｔｉｃｌｅ．Ｔｈｅｓｅｃｏｎｄｐａｒｔｉｓｔｈｅｓｏｌｕｔｅｔｒａｎｓｐｏｒｔｍｏｄｅｌｉｎｇｒｏｕｎｄｗａｔｅｒ．Ｉｎｔｈｉｓｓｅｃ—ｔｉｏｎ，ｉｔｐｒｅｓｅｎｔｓｔｈｅｂａｓｉｃａｓｓｕｍｐｔｉｏｎｓｏｆｇｒｏｕｎｄｗａｔｅｒ．Ａｎｄｗｅｅｓｔａｂｌｉｓｈｔｈｅｅｑｕａｔｉｏｎｓｏｆｔｈｅｇｒｏｕｎｄｗａｔｅｒｆｌｏｗａｎｄｔｈｅｓｏｌｕｔｅｔｒａｎｓｐｏｒｔａｃｃｏｒｄｉｎｇｔｏｍａｓｓｃｏｎｓｅｒｖａｔｉｏｎｏｆｔｈｅｅｑｕｉｌｉｂｒｉｕｍｕｎｉｔｃｅｌｌ．ＯｎｔｈｅｂａｓｉｓｏｆＤｕｐｕｉｔａｓｓｕｍｐｔｉｏｎ，ｗｅｕｌｔｉｍａｔｅｌｙｆｏｒｍａｔｔｗｏ－ｄｉｍｅｎｓｉｏｎａｌｍｏｄｅｌ．Ｆｉｎａｌｌｙｗｅｍａｋｅａｄｅｔａｉｌｅｄａｎａｌｙｓｉｓｏｆｔｈｅｐａｒａｍｅｔｅｒｓｉｎｖｏｌｖｅｄｉｎｔｈｅｅｑｕａｔｉｏｎ．Ｔｈｅｔｈｉｒｄｐａｒｔｉｓｔｈｅｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｏｆｔｗｏ－ｄｉｍｅｎｓｉｏｎａｌｇｒｏｕｎｄ—ｗａｔｅｒｍｏｄｅｌ．ＩｔｇｉｖｅｓｂｏｔｈｔｈｅａｌｔｅｒｎａｔｉｎｇｄｉｒｅｃｔｉｏｎｄｉｆｆｅｒｅｎｔｉａｌｆｏｒｍａｔｏｆｔｈｅⅨｏｆｔｗｏ－ｄｉｍｅｎｓｉｏｎａｌ．Ｗｅｉｎｔｒｏｄｕｃｅｅｑｕａｔｉｏｎｓａｎｄｔｈｅｃｕｂｉｃｓｐｌｉｎｅｉｎｔｅｒｐｏｌａｔｉｏｎｔｗｏｍｅｔｈｏｄｓｔｏａｐｐｒｏｘｉｍａｔｉｏｎｅｑｕａｔｉｏｎｏｆｎｏｎｌｉｎｅａｒｃｏｅｆｆｉｃｉｅｎｔａｎｄｎｕｍｅｒｉｃａｌｅｘａｍｐｌｅｓｐｒｏｖｅｔｈａｔｔｈｅｓｅｃｏｎｄｍｅｔｈｏｄｉｓｍｏｒｅｅｆｆｅｃｔｉｖｅ．Ｆｉｎａｌｌｙ，ｔｈｉｓｐａｐｅｒａｎ—ａｌｙｚｅｓｔｈｅｓｔａｂｉｌｉｔｙｏｆｔｈｅｌｉｎｅａｒｓｏｌｕｔｅｔｒａｎｓｐｏｒｔｅｑｕａｔｉｏｎｓ．Ｋｅｙｗｏｒｄｓ：Ｇｒｏｕｎｄｗａｔｅｒｃｏｎｔｉｎｕｉｔｙｅｑｕａｔｉｏｎ；ｓｏｌｕｔｅｔｒａｎｓｐｏｒｔｅｑｕａ－ｔｉｏｎ；ａｌｔｅｒｎａｔｉｎｇｄｉｒｅｃｔｉｏｎｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄ；ｓｐｌｉｎｅｉｎｔｅｒｐｏｌａｔｉｏｎ；ｓｔａ－ｂｉｌｉｔｙ；Ｘ第１章前言１．１背景介绍地下水是全球水资源的重要组成部分，在社会经济发展中发挥着重要的作用。
非饱和带溶质运移问题的SUPG有限元解法
非饱和带溶质运移问题的SUPG有限元解法
高俊合;朱学愚;赵维炳
【期刊名称】《河海大学学报（自然科学版）》
【年(卷),期】1998(000)001
【摘要】ＳＵＰＧ有限元法应用于土体的非饱和溶质运移问题，导出了其数值表达式，编制了计算程序。
以一维入渗问题为例进行了数值计算，与其他方法及野外试验结果进行了对比，表明该方法是正确的、稳定的，可用于非饱和带溶质运移问题的实际计算。
【总页数】5页(P54-58)
【作者】高俊合;朱学愚;赵维炳
【作者单位】河海大学岩土工程研究所;南京大学地球科学系
【正文语种】中文
【中图分类】P641.2
【相关文献】
1.非饱和流动问题的SUPG有限元数值解法 [J], 朱学愚;谢春红;钱孝星
2.超音速粘性流动的SUPG有限元数值解法 [J], 徐国群;张国富
3.饱和水流溶质运移问题数值解法综述 [J], 成建梅;胡进武
4.N-S方程SUPG有限元解法 [J], 熊杰;俞守勤
5.非饱和土壤中溶质运移对流占优问题的通量校正运移解法 [J], 徐绍辉;张佳宝;刘建立
因版权原因，仅展示原文概要，查看原文内容请购买。
二维地下水溶质运移耦合方程组的差分方法的开题报告
二维地下水溶质运移耦合方程组的差分方法的开题报告一、研究背景地下水对于很多领域的人来说是一个至关重要的资源，其中包括农业、工业和城市居民供水等。
越来越多的人开始关注和研究地下水的管理和保护。
地下水流和溶质运移是地下水研究中的两个重要方面。
为了更好地理解这些过程，建立数学模型是必不可少的。
二维地下水溶质运移耦合方程组是研究地下水流和地下水溶质运移相互作用的数学模型，对于解决地下水污染和地下水资源管理等问题至关重要。
目前，该方程组的差分方法是最常用的解决方法，因为它是一种简单而有效的数值求解技术。
因此，本研究旨在对二维地下水溶质运移耦合方程组的差分方法进行深入研究和探讨，以更好地理解这个模型和提高数值计算的精度和速度。
二、研究内容本研究将主要涉及以下几个方面：1. 理论分析：通过对二维地下水溶质运移耦合方程组的理论分析，深入了解该模型的特点和解决方法，为后面的数值计算提供基础。
2. 差分方法探究：对现有的差分方法进行分析和比较，选取最优的差分方法用于数值计算。
同时，也将尝试开发一种新的差分方法，以提高计算精度和效率。
3. 数值计算的实现：通过编写程序实现差分方法，对数值结果进行验证和影响分析，以对差分方法的优缺点和适用范围进行评估。
4. 应用案例分析：选取一些典型应用案例进行分析，以验证差分方法的适用性和精度。
同时，也将探讨如何在实际应用中进行该模型的参数确定。
三、研究意义本研究的意义主要有以下几个方面：1. 通过数学建模和数值方法开发，可以更好地理解和研究地下水溶质运移的过程，为地下水污染控制和地下水资源管理提供基础数据和技术支持。
2. 探究差分方法，可以提高数值计算精度和效率，为地下水模型的研究和应用提供更加准确的数据和预测。
3. 应用案例分析，可以对模型参数的确定和模型应用的可行性进行评估，为实际应用提供指导和支持。
四、研究方法本研究将应用理论分析和数值计算相结合的方法来完成研究任务。
主要研究方法包括：1. 理论分析：通过文献调研，研究二维地下水溶质运移耦合方程组的数学模型，探讨不同差分方法的数学原理和优缺点。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
关 键 词 ：对 流 扩 散 方 程 ；拉 普 拉 斯 变 换 ；混合 法
中 图分 类 号 ：Ｏ ３ ５ 文 献标 志 码 ：Ａ
察 尔汗 盐湖首 采 区含水 介 质为石 盐 层 ， 高浓 度 卤水 溶 液 （ 化 度 ３ ０～ ５ Ｌ ） 矿 ０ ４ ０ｇ・ 在这 种 易溶 岩 介 质 中 的溶 质 运移 涉及 到一个 突ｍ 的问题 是 卤水 与介质 问的 物相 转化 即 固液转 化 问题 ． 对 于在 开 采 而 过 程 中卤水 与易溶 介 质之 问 的同液转 化 问题 ， 由于涉 及 到矿 床地 质 、 文地 质 、 理 化学 、 球化 学 、 孑 水 物 地 多 Ｌ 介 质水 动力 弥散理 论 及水盐 体 系相 图原 理等 多方 面 的知识 ， 目前 国 内外 尚无 可借鉴 的 经验 和理论 依据 ．
文章 编 号 ：０ ４—１ ２ （ ０ １ ０ — ０ ５一 ４ １０ ７ ９ ２ １ ） １ ０ ２ ｏ
溶 质 运 移 模 型 的 有 限 元 数 值 解
邓 永 辉 ， 永 红 邓
（ 南财政经济学院 ， 湖 湖南 长 沙 ４０ ０ ） １２ ５ 摘 要 ： 混 合 拉 普 拉 斯 变换 有 限 单 元 法 引 入 到 求 解 首 采 区 卤水 动 态 二 维 模 型 的 溶 质 运 移 问 题 巾 ， 够 有 将 能
ｄ ）从 以上分 析可 知 ， 单位时 间 、 在 单位 体积 内浓 度随 时问 的变 化而 变化 ， 可表示 为
ｃ ）输入输 出可 以包括 弥 （ ） 和对流 引起 的现 象．式 （ ） （ ） 截 面积和 △ 去 除后 ， 扩 散 ３和 ４被 ￡ 使后 二项 系 数 为零 ， 即除以 Ａ Ａ Ａ 分别 可得 到 以下 方程 ｘｙｔ
［ （ Ｉ ＭＵ Ｏ Ｄａａ 】 ｌ￣ Ｃ ｘ ＼ 一ｘ Ｏ＝ ， ｃ Ｃ 【 ＼￣ ａ＝ Ｏ Ｄａａ 】 ｙ，ｃ 一 ’ （ ｃ Ｃ ．
隙度率 ， 即给 水度 ， 也就是 从单位 体积 含水 层 中在重 力作 用 下 能够 释放 给 出的 水量 所 占该 含水 层 体 积 的
份数 ．
ｂ ）顺 Ｙ轴方 向 由渗 流引起 物质迁 移 即渗流扩 散模 型 ， 用质量 守衡方 程求 出微分 方程解 等方 法来 确 可 定． 长为 △ 在 的单 元体 内物 质迁移 可表示 为
限地 消 除在 求 解 埘 流 ｒ优 的 地 下 水 溶 质 运 移 问 题 时 产 的数 值 扩 散 和 过量 的 现 象 ， 有 一 步 到 位 、 部 求 解 ｑ － 具 局 的优 点 ， 后 将 该方 法 应 到具 有 卒 间一 阶 导 数项 的 对 流 弥 散 方 程 ， 最 以检 验 该 方 法 的 数 值 有 效 性 和 求 解 溶 质
（ ２ ）
其中，， Ｄ 为扩 散系数 （ 表示 扩散 方 向 ） 为 真实 渗透 速度 ， 液体 浓度． Ｚ ， Ｃ为 １２ 二 维对 流 扩散 方程 ． 从 二维 流的 渗流 场 中 割离 出一 个 微 分单 元 体， 该单 元 的宽 度 为 ， 长度 为 Ａ ， ｙ 厚度 为
水 沿 ．轴 方 向从左 面 流人 单 元 体 ； Ｙ轴 方 向从 前 流人 ． 时段 ３ ｇ 沿 经 △ 后 分别 经 Ａ ， ｙ从 对 面流 出． ｆ ｘＡ
第２ 卷 第 １ ９ 期
２ １ 年 ３月 ０ １
海 南 大 学 学 报 自 然 科 学 版
ＮＡＴＵＲＡＬ Ｓ ＥＮＣＥ Ｊ ＣＩ ＯＵＲＮＡＬ ＯＦ ＨＡＩ ＮＡＮ ＵＮＩ ＶＥＲＳ ＴＹ Ｉ
Ｖ ｌ ２ Ｎ １ １ Ｉ． ９ （ １ ． Ｍａ ． ０ｌ ｒ２ １
海 南 大 学 学 报 自 然 科 学 版
［ ＤＯ＋ ） ＝］ （ ｘ ｖ）】 － ｃ Ｃ
ｘ＋Ａ ｘ
ＡｙＡｃ ，
输 入 ＝ 出 ＋积 累 ， 输 其 中 ，一Ｄ
（） ３
为扩散 通量 ， 其方 向与浓度 梯度 方 向相反 ， 即负号 的含 义 ； 为含 水层厚 度 （ ； Ｍ ）Ｉ 自由孑 ｘ为 Ｌ
［ Ｄａ＋ ） （ ｙ ｃ 州 － Ｃ △
，
［ ＤＯ＋ｙＡ （ ｙ Ｕ） － Ｃ引 Ｃ
ｍ
Ａ ＋Ｍ ｙ Ｃ， ｚ ｇＡ Ａ △
输 入 ＝ 出 ＋积 累 ， 输 其中， 一Ｄ 为 扩散通 量 ， 方 向与浓度 梯度方 向相 反 ， 其 即负号 的含义 ．
ｏ ｙ
（） ４
考 察渗 流 引起 的浓度 变化 ：
ａ ）顺 轴方 向由渗流 引起 物质 迁移 即渗 流扩 散模 型 ， 可用 质量 守
衡 方程 求 出微分 方程 解等 方 法 来 确 定． 长 为 在
移 可表 示 为
的单 元 体 内物 质迁
收 稿 日期 ： ００—１ ２１ ０一Ｉ ｌ 基 金 项 目 ：国家 攻 关 项 目 （０ １ Ａ ０ Ｂ一 ９—１ ２０ Ｂ ６ Ｚ ０ ） 作 者 简 介 ： 永 辉 （９ ９一） 女 ， 南 常 德 人 ， 南 财 政 经 济 学 院 讲 师 ， 士 邓 １７ ， 湖 湖 硕
１ 模 型 的 建 立
１ １ 对流 扩 散方程 ． 在 建立 目标 模 型之前 先 了解一 下一 维单 向的流场 的对 流扩 散方 程 ¨ ．
对 于一维 的对 流扩 散模 型 ， 以用质 量守 衡方 程求 出微分 方 程解 等方 法 来 确定 ． 期 一个 长 为 圆 可 周
柱形 内物 质迁 移表 示 为
（Ｄ ＋Ｃｒ］ ［ －ｒ ） ２ ＝（ ＇ ｒ Ａ Ｒ
Ｄ
＋Ｃｒ］ ＋ Ａ ｃ ） ２ ｘ・ ｘ一 Ｒ Ａ
（ １ ）
式 （ ） 以包 括扩 散 和对流 引起 的现 象 ． （ ） 面积 ， 和 △ 去 除后 ， １可 式 １被 ｆ 使二 项式 为零 ，
磐一ａ ， ０ Ｏ ＝ ｏ Ｃ