| 查看: 2296 | 回复: 13 | ||
小木虫zb木虫 (正式写手)
|
[求助]
求解一个三维的偏微分方程组已有1人参与
|
该模型为固定床动态吸附模型，模型和模型参数如下： 有两种求解方法： 1.对第一个方程经行有限差分离散，对第二个方程采用Crank-Nicolson方法离散，然后通过r=R处的边界条件联系起来求解，参见文献 《1Effect of Concentration-Dependent Surface Diffusivity on Simulation of Fixed Bed Sorption Systems》附录所示； 2.对第一个方程的床层方向z离散，对第二个方程粒径方向r离散，得到一系列关于时间的常微分方程组，然后通过r=R处的边界条件联系起来求解，参见文献 《1Adsorption dynamics of p-nitrophenol in structured fixed bed with microfibrous entrapped activated carbon》中第四页。 两种方法都尝试过，但是在两个方程的连接处不知道怎么处理（因为涉及到三维z,r,t，所以不知道怎么将两个方程通过边界条件联系起来求解），希望哪位大神帮帮忙，帮我用matlab编个程序，或者给我指点一二，谢谢了！[ Last edited by 小木虫zb on 2014-2-26 at 10:00 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : 1Effect_of_Concentration-Dependent_Surface_Diffusivity_on_Simulation_of_Fixed_Bed_Sorption_Systems.pdf
- 附件 2 : 1Adsorption_dynamics_of_p-nitrophenol_in_structured_fixed_bed_with_microfibrous_entrapped_activated_carbon.pdf
2014-02-25 22:51:19, 398.28 K
2014-02-25 23:04:34, 1 M
» 猜你喜欢
感谢小木虫，已完成目标
已经有0人回复
-
零点能
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有55人回复
-
科研求助
已经有0人回复
-
【25申博】博导推荐求助
已经有6人回复
-
青岛大学李山东教授招生
已经有2人回复
-
新泽西理工大学Prof. Junjie Yang招实验凝聚态物理博士
已经有1人回复
-
新泽西理工大学Prof. Junjie Yang招实验凝聚态物理博士
已经有1人回复
-
2024年度“中德博士后”项目启动申请，截止时间6月30号
已经有0人回复
-
新泽西理工招实验凝聚态物理博士
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有偿求助，偏微分方程求解编程~（急）
已经有3人回复
- 求助解复杂非线性方程组的好的方法
已经有24人回复
- 哪位大神给指导下如果求解一个偏微分方程组
已经有11人回复
- 三维非稳态导热微分方程的MATLAB解法
已经有10人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 非线性二阶常微分方程组求解--matlab
已经有16人回复
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 多项式方程次数与约束方程个数之间的关系
已经有4人回复
- fluent-udf自定义标量方程边界条件
已经有5人回复
- [原创]在Origin8.6 中用作图法求解超越方程
已经有45人回复
- 拜求，这种带参数的方程怎么解比较好？
已经有4人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab数值模拟求助（二元N-1次方程组）
已经有3人回复
- 请问如何求解二元一阶微分方程组
已经有24人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 怎么用MATLAB生成一个三维的Shepp-Logan头骨模型并保存下来，谢谢了
已经有5人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 一个微分方程组，求解
已经有18人回复
- 各位，在论文里写明数值模拟是用matlab进行的，会不会有版权问题？
已经有3人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】微分方程组拟合参数求助
已经有9人回复
- 【求助】请问怎么求函数方程的函数表达式？
已经有11人回复
Nonebull
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5663.1
- 红花: 6
- 帖子: 463
- 在线: 306.9小时
- 虫号: 571775
- 注册: 2008-06-10
- 专业: 化学工程及工业化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫zb: 金币+50, ★有帮助 2014-03-05 22:41:56
fegg7502: 金币+1, 鼓励交流 2014-03-07 09:36:29
小木虫zb: 金币+50, ★有帮助 2014-03-05 22:41:56
fegg7502: 金币+1, 鼓励交流 2014-03-07 09:36:29
|
楼主 首先谢谢你的参考文献了 很有启发 编程还是要你自己来解决 第一篇文献里的编程思路已经很清晰了，不过还是有点小问题，例如两片文献其实都没有把模拟时使用的具体参数值列完全，导致重复他们的结果比较困难，另外在1330页，cs(iz)=f(Q(iz))里面在实践过程中应该改成cs(iz,it+1)=f(Q(iz,it))，这样就好理解了； 具体其实只有q是三维的，c和cs是二维的，算c的方程很简单，主要是要列出算q的A和B矩阵，两个方程的边界条件是通过平衡假设来连接的，cs和最外层q(R)满足sips的吸附平衡关系，希望能有帮助！ |
2楼2014-03-01 23:51:22
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1558.5
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
3楼2014-03-02 00:07:30
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1558.5
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
★
fegg7502: 金币+1, 鼓励交流 2014-03-07 09:36:39
fegg7502: 金币+1, 鼓励交流 2014-03-07 09:36:39
|
您好，收到您的回复非常感谢，对我有些启发，其实我不是希望您给我完整的程序，我只是对1330页的算法过程有些地方不是很清楚，例如C(IZ)=[E1*C(IZ-1)+C(IZ)+E2*CS(IZ)]/E3，按照文章的意思该处完整的意思是C(IZ，it)=[E1*C(IZ-1,it)+C(IZ,it-1)+E2*CS(IZ,it)]/E3，而在第一次求解中C(1，1)求得的是包含CS(1，1)的代数式，将这个代数式带入B中的（C(n,i)-CS(n,i）+C(n-1,i)-CS(n-1,i）,然后用AX=B求解，不明白的地方是用此代数式能求出Q吗？希望能有机会和你交流 |
4楼2014-03-03 14:27:58
Nonebull
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5663.1
- 红花: 6
- 帖子: 463
- 在线: 306.9小时
- 虫号: 571775
- 注册: 2008-06-10
- 专业: 化学工程及工业化学
5楼2014-03-03 23:14:51
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1558.5
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
★
fegg7502: 金币+1, 鼓励交流 2014-03-07 09:36:54
fegg7502: 金币+1, 鼓励交流 2014-03-07 09:36:54
|
还以为是代数方程带入（AX=B），然后加上q（r=R）和cs的关系，联立解这个隐函数方程组求解； 用CS（iz,it-1)带入求得C（iz,it)的数值，然后带入求AX=B，得q(r=R),通过它求得CS（iz,it)，是不是此处还要检验一下CS（iz,it)是不是等于CS（iz,it-1)，我的意思是：您是不是用前一个时间的值做初始值，然后迭代求解，如果结果和初始值不一致，则改变初始值继续求解？还是考虑吸附平衡的非瞬间达到，而认为下一个时间的浓度CS（iz,it+1)和上一个时间的平衡吸附量q(R,it)平衡？ |
6楼2014-03-04 12:42:44
Nonebull
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5663.1
- 红花: 6
- 帖子: 463
- 在线: 306.9小时
- 虫号: 571775
- 注册: 2008-06-10
- 专业: 化学工程及工业化学
7楼2014-03-04 21:01:09
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1558.5
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
★
送红花一朵
fegg7502: 金币+1, 鼓励交流 2014-03-07 09:37:09
送红花一朵fegg7502: 金币+1, 鼓励交流 2014-03-07 09:37:09
|
前辈，受你的启发，我编了如下程序，不过q的求解不是采用文献中的AX=B追赶法，而是采用了对粒径方向向后差分化为关于时间的常微分方程组，然后调用matlab中的ode15s求解，其他地方不变，但是运行的结果全是零，我检查了很多遍，不知道问题出在哪里，有空您能帮我看下吗？万分感谢！ function HSDM_MOL_test_1 clear all;clc kf=200;n=1/0.52355;k=0.12456;%液膜扩散系数，Freundlich常数1/n和k a=0.363;u=0.361032/a;c0=7.03;%a为空隙率，u为空隙间流速，c0为初始浓度 L=8.31e-2;nz=20;dz=L/nz;nt=20;dt=20/nt;R=0.35e-3;nr=10;dr=R/nr; E1=u*dt/dz;E2=3*kf*dt*(1-a)/(a*R);E3=1+u*dt/dz+E2; cs=zeros(nt,nz);c=zeros(nt,nz);cs(1,1)=c0;cs(1,2:nz)=0; for it=1:nt waitbar(20/it) for iz=1:nz if (it==1&&iz==1) c(it,iz)=(E1*c0+E2*cs(it,iz))/E3; end if (it==1&&iz~=1) c(it,iz)=(E1*c(it,iz-1)+E2*cs(it,iz))/E3; end if (it~=1&&iz==1) c(it,iz)=(E1*c0+c(it-1,iz)+E2*cs(it,iz))/E3; end if (it~=1&&iz~=1) c(it,iz)=(E1*c(it,iz-1)+c(it-1,iz)+E2*cs(it,iz))/E3; end c_q=c(it,iz);cs_q=cs(it,iz); q0=zeros(1,10); options=odeset('relTol',1e-6);tspan=[0:20]; [t,q]=ode15s(@particle,tspan,q0,options,c_q,cs_q);%求得的c和cs带入方程求q cs(it+1,iz)=(q(it,nr)/k)^n;%通过颗粒r=R处的q（it,nr)求cs(it+1,iz)，迭代下一步求c(it+1,iz) result_c(it,iz)=c(it,iz); end end disp(result_c) function dqdt=particle(t,q,c_q,cs_q)%采用对粒径r向前差分，运用ode15s求解 Ds=10;R=0.35e-3;nr=10;dr=R/nr;p=436.8/(1-0.363);kf=200;%p为颗粒密度，Ds为表面扩散系数 dqdt(1)=3*Ds*(q(2)-q(1))/dr^2; for j=2:nr-1 dqdt(j)=Ds/dr^2*((2/j+1)*q(j+1)-(2/j-2)*q(j)+q(j-1)); end dqdt(nr)=Ds/dr*((2/nr+1)/Ds*kf/p*(c_q-cs_q)-(q(nr)-q(nr-1))/dr); dqdt=dqdt'; |
8楼2014-03-06 15:07:32
Nonebull
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5663.1
- 红花: 6
- 帖子: 463
- 在线: 306.9小时
- 虫号: 571775
- 注册: 2008-06-10
- 专业: 化学工程及工业化学
9楼2014-03-06 23:21:33
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1558.5
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
10楼2014-03-07 00:01:33
5