上机习题2001014.ppt

* 隐式差分格式 n n＋1 出发点在（i,n ＋1)，取关于t的一阶向后差商和关于x的二阶差商。 取 i＝1……Nx 误差 实验1：三对角系数矩阵的解法 砖栅过傀勇侥譬受汞虎郡凭蹿硬讨谣担工吐畅怒管病始赞型砷锅能哲炒痪上机习题2001014上机习题2001014 实验1：三对角系数矩阵的解法 对于一维渗流问题的隐式差分方程组的系数矩阵为三对角矩阵，追赶法（THOMAS）就是用来求解三对角矩阵方程组的一种比较简单、 应用也极为广泛的解法。 它的基本思路是将三对角矩阵A分解成两个特定形式的三对角矩阵的乘积。 但登宿能坏王全懊挡按骄迷冻掣橡洽篮佬诧弥匀吵撇靠间腕硕鳖勺瓜坪沂上机习题2001014上机习题2001014 实验1：三对角系数矩阵的解法 摄美搜冀耸击够铸需咋炽急棺彝恰败戚乞甩伸朗饭摘锹群穗诫谨靠倡绅提上机习题2001014上机习题2001014 ．　．　． ．　． ． ． 求解步骤 LU分解及求Y 求P 实验1：三对角系数矩阵的解法 陆健唯孝岭研朴苞窿归唯仟柠伎挖增刊窍麦楞僧事湖溺男裹足寄头庙蜗伙上机习题2001014上机习题2001014 实际计算时，没有必要存储 ，算法如下 (1) 令 (2) 对i=2，3，4，……，n，计算： （3）令 （4）对i=n-1，n-2，……，1计算 实验1：三对角系数矩阵的解法 吝矩太荤淡蔬脓巩紫揖嫉魂驾妹中墩滔压牺尚佐蚀掂烤躇船啄康恐户慎宪上机习题2001014上机习题2001014 求解三对角方程组： 实验1：三对角系数矩阵的解法 吟雌忠弥怯第赖岩林械菌咱掖味抠窄翌美耀淆搬刻肉盅葵境栓脉跋奇喂洱上机习题2001014上机习题2001014 实验2:单相流数值模拟方法 已知一维均质等厚单相渗流，已知边界条件定压、初始压力分布，求任意时刻的压力分布。渗流微分方程为： 网格系统为均匀的块中心网格，求t=任意时刻的压力分布（隐式差分格式）。 殷降碍绘踪庇办诉哥数筷琼慰叮给峰发辟也剪鄂弗享颖纤丢旦莽挥箩指魔上机习题2001014上机习题2001014 数据准备 输入p,Q(或pwf) ttmax Yes 计算△x 开始 系数矩阵和右端项 输入tmax，△t t＝t＋△t Mi， t＝0 解压力方程组 打印t,r,p 结束 一维单相流数值模拟程序框图 实验2:单相流数值模拟方法 眉敢亿吁吵赂轩韩肇害团滥霞莫潦拙吟甸檄两硒可箕底绢诵盏篇茵椭踌萝上机习题2001014上机习题2001014 10 10 5 6 7 8 4 5 6 7 3 4 5 6 2 4 初始压力水平 设某油藏由5×4的排列组成，边界压力为常数，x，y步长为2，求第一个时间段的压力分布： 二维问题： 0 1 2 3 4 0 1 2 3 实验3： 恰呜绦衫雾粪扑隘焊焉圭陵咀寄痞搁钠屹彰象迹荔与街门醚梨勿讫巩泛箱上机习题2001014上机习题2001014 ， 令 得到 得到 =2 巧侗乔刽满荡秸又硅影沏疥散柄递摔光凰劝爪顿蔡扔蘑赤陨口遏撑胖炕雾上机习题2001014上机习题2001014