版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、讹赤工章大每读书报告Matlab在求解扩散系统之浓度分布中的应用池雨一、问题的提出管中储放静止液体B,高度为L=10 cm,放置于充满A气体的 环境中。假设与B液体接触面之浓度为CA0=0.01mol/m3,且此浓度不 随时间改变而改变,即在操作时间内(h=10夭)维持定值。气体A在 液体B中之扩散系数为Dab=2X 10 9m2/s。试决定A与B不发生反应; 情况下,气体A溶于液体B中之流通量(flux)。参考如图所示的装 置。气体A在澈体B中的扩敝二、知识背景Fick第一定律:实验表明,在稳态扩散的条件下,单位时间内通过垂直于扩散方向的单位截面面积的扩散物质的通量与浓度梯度成 dc正比。数
2、学表达式为:J = -D dxFick第二定律:根据质量平衡关系即在微小体积中积存的物质(留入的物质量)J(留出的物质量)J得出d = J1-J212 dt dx因此Fick第二定律的数学表达为:竺=D竺dtdx 2三、问题求解根据题意不同时间t和距界面厚度不同处x的浓度C=f(z,t);因气 体A与液体B不发生反应,故其扩散现象的质量平衡方程根据Fick_n g第二定律。-.-依题意,其初始及边界条件为:I.C. CA0(Z,0)=0, Z0B.C. Ca(0,D=Ca0, t 涌;与 z=l= 河“z=0dC 在获得浓度分布后即可应用Fick第一定律求得流通量即NAz(t)= DABdZ四
3、、matlab程序设计偏微分方程(Partial Differential Equation,简称PDE)就是涉及 到两个自变量以上的微分方程。在化学工程领域,为了更好的进行过 程设计、优化和控制,经常需要了解化工设备(如反应器)中的温度、 浓度和速度在不同空间上的分布以及随时间的动态变化规律,因此涉 及到许多偏微分方程的问题。Matlab函数pdepe()可用于求解偏微分方程,模型为:du. du _ -m d . m du . du. cjl dt exexdx用以解含上述初始值及边界值条件的偏微分方程MATLAB命 令pdepe的用法如下:sol =icfun.bcfun.tspan.
4、options)若要获得特定位置及时间下的解,可配合以pdeval命令。使用格式如下:如血 tZwcnMdx = pdevaKm. xm.sh, ui.xuut)f (乙t,匕=0第一步:Fick化为标准形式后即8c1c(乙t,七云) = DAB8c第二步:编写偏微分方程的系数向量函数:function c,f,s=ex20_3_2pdefun(z,t,CA,dCAdt)c=1/Dab ; f=dCAdt; s=0;第三步:编写起始条件:function CA_i=ex20_3_2ic(z)CA_i=0;第四步:编写边界条件:function pl,ql,pr,qr=ex20_3_2bc(zl
5、,CAl,zr,CAr,t)global DAB k CA0pl=CAl-CA0;ql=0;pr=0;qr=1/DAB;第五步:取点。例如t=linspace(0,h,100);z=linspace(0,L,10);第六步利用pdepe求解。m=0;sol=pdepe(m,ex20_3_2pdefuna,ex20_3_2ic,ex20_3_2bc,z,t);Sol(j,k,i)第一维代表时间第二维代表空间位置z,第三维代表解向量u的第i个元素,即ui=sol (:,:, i)第七步显示结果。u=sol(:,:,1);surf(x,t,u)title(标题)xlabel(位置)ylabel(co
6、nc. (mol/mA3)zlabel(u)若要获得特定位置及时间下的解,可配合以pdeval命令。使用格式如下S挪皿站1 =河冲加珈绊即/蝴本题的程序编写为:for i=1:length(t)CA_i,dCAdz_i=pdeval(m,z,CA(i,:),0);NAz(i)=-dCAdz_i*DAB;End附录(matlab程序):function ex20_3_2%扩散系统之浓度分布clearclcglobal DAB k CA0%给定数据a=input(请输入要求通量处厚度a=)CA0=0.01;L=0.1;DAB=2e-9;h=10*24*3600;%取点t=linspace(0,h,
7、100);z=linspace(0,L,10);m=0;sol=pdepe(m,ex20_3_2pdefuna,ex20_3_2ic,ex20_3_2bc,z,t);CA=sol(:,:,1);for i=1:length(t)CA_i,dCAdz_i=pdeval(m,z,CA(i,:),a);NAz(i)=-dCAdz_i*DAB;end subplot(211)surf(z,t/(24*3600),CA) title(case (a) xlabel(length (m)ylabel(time (day)zlabel(conc. (mol/mA3) subplot(212)plot(t/(24*3600),NAz*24*3600)xlabel(time (day)ylabel(flux (mol/mA2.day) %* % PDE函数 %* function c,f,s=ex20_3_2pdefuna(z,t,CA,dCAdz)global DAB k CA0c=1/DAB;f=dCAdz;s=0;%* %初始条件函数%* function CA_i=ex20_3_2
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年电力因数校正控制器的设计与应用
- 2024全新房产过户协议文档下载
- 赞助合作合同范本
- 搭建loft 合同范本
- 涵管销售合同范本
- 2024家装施工与服务协议条款
- 2024年企业派遣人员服务协议样本
- 六年级数学上册 期中常考易错真题检测卷(试题)(苏教版)
- 2024年度铁矿石大宗买卖协议
- 交叉学科教学法模板
- 国开2024年《中国法律史》平时作业1-3答案
- 第八讲 地形图应用(二)
- 普铁避雷器检修作业指导书
- 下水管道施工合同通用版
- 铁合金生产工艺
- 钢结构策划书(范本)
- 焦化厂生产工序及工艺流程图
- 呕吐(急性胃肠炎)诊疗指南(制订)编制说明排版
- (外研版)初中英语语法汇总[新版]
- 李燕璇植树问题卡通版5
- 有砟轨道铺设的施工讲解
评论
0/150
提交评论