环境流体力学课件.pptxVIP

环境流体力学Matlab编程题目分析由于有害物质在池底均匀分布，池底及池壁对有害物质不完全吸收，则可按瞬时平面源一维扩散处理。已知条件：池底面积200m×200m水深50m有害物质质量4000kg扩散系数1.0cm2/s 1.解析法1.1 计算公式因底部不完全吸收，由底部沿垂向的扩散浓度应为式中x为自池底量起的距离， kg/m 扩散系数D=1.0cm2/s=8.64m2/天。2.有限差分法一维扩散的扩散方程为：2.1 网格剖分 令 x=xi=ih，h=1m，0≤i≤50 t=tk=k??，??=1天，0≤k≤365 2.2 古典隐格式——稳定性和收敛性差分格式为（与解析法区分，差分法用u表示C）：令步长比B表示成矩阵形式A3. Matlab程序clc;clear;close all; %清除屏幕上的内容，清除工作区里的内容，关闭画的图x = -8.64 * ones(1,48); A = 18.28 * eye(49) + diag(x,-1) + diag(x,1); InvA = inv(A);%求矩阵A的逆矩阵M = 0.1;%参数赋值D = 8.64; r = 8.64; l = 50;B = zeros(49,365);for t = 1:365 B(1,t) = r * 2 * M / sqrt(4 * pi * D * t); B(49,t) = r *