用蒙特卡罗方法和区域分解算法求解一类对流扩散方程.pdfVIP

用蒙特卡罗方法和区域分解算法求解一类对流扩散方程 刘播 刘芳芳 刘春光 吉林大学数学研究所，130012，长春 liubom钊.en,liufang8265@126.com 摘 要 对对流扩散方程的求解域作矩形剖分，对时间层分层求解，每 一层用区域分解方法分割成若干个子域，子域内边界点用蒙特卡罗方法求解， 并给出了误差估计.这样每个子域的求解是独立的，可以进行并行计算.对一 维方程做了数值实验. 关键词 蒙特卡罗方法 区域分解 并行计算 对流扩散方程 1 弓1言 蒙特卡罗(MonteCarlo)方法，简称S:IC法，又称随机抽样方法，是一 种与一般数值计算方法有着本质区别的计算方法，它利用随机数进行统计试 验，以求得的统计特征值 (如均值等)作为待解间题的数值解。 流体力学问题的微分方程，用隐式格式求解是绝对稳定的，但是直接用传 统的方法 (如Gauss消去法，超松弛迭代法等)求相应差分方程组的解，运算 量很大，计算时间很长，对于高维情况，问题更显得突出.本文介绍了一种可 以用多处理机对差分方程进行并行计算的方法.先用区域分解算法建立并行 的差分算法，再利用NIC法计算各子域内边界点的函数值，这样各个子域的求 解是独立的，可以并行计算. 2 区域分解算法 区域分解算法包括不重叠区域分解法和重叠区域分解法，我们这里采用 不重叠区域分解法.它的基本思想是:将计算区域G分解为若干个子域，G= U=77=11G=，子域G;的形状尽可能规则.这样，可将原间题的求解转化为在各个 子域上求解.只要知道各子域边界上的函数值，就可用迭代法求解区域内各网 格点函数值.子域边界上的函数值可用mC法求解. 对于一维情形，我们只需用MC法求解p个单点值，就可以将求解区域 分成P十1个子区域，在每个子域上分别用迭代法求解，其差分方程组系数矩 阵的维数不超过NP--tl，其中N为空间域上的网点数，这样并行度大为提高· 考虑对流扩散方程 . . . . 了 U 毛 户 川 优 处 会·+hx(t)dx+x(,，，‘0x，‘0‘T, 夕 1 t t 冬 ) 0 二 沪(二)， (1) U 几 ) t 、 二 /l(t), u(l,t)二u2(t) 取空间步长h=IIN，和时间步长二=T/L，其中N,L都是自然数，用两族平 行直线x=x.f二jh(j二0,1.…，N)和t二tk=k,r(k=0,1,2,…、L)将矩形域 G={0x1,0￡T}分割成矩形网格，网格点分为两类，一类是内部网 格点集，记为Gi,，它的四个相邻网格点均位于区域内部和边界上;另一类是边 界网格点，我们把二=0，二二l上的网格点集记为r大，t=0上的网格点集记为 r;, j+1 图(1) 与方程 (1)相应的隐式差分方程组为 _1) 、十i 十“畏uk十‘ +2(uk.7十，一，tTli 一u(k.i++ll-2uhkj}+1+uak_+1l)+bk3 9 一 十马 L一1. .l二1.2，二 ，N一1,k二0.1,… 。 心 =};(xi), 姑 =P1(tk), u矢=P2(tk)- 假设对每一层已求出ni,n2,...,pn这P个点的值，即可将求解域分成x!p,XnJ. X[ni,Xni+1I,(y.二1,2,…,P一1),x[np,XNJ这;+1个子域，可在这些子域上求解 阶数小得多的差分方程. 下面我们用蒙特