24小时热门版块排行榜    

查看: 1666  |  回复: 5

ghw_nit

铁杆木虫 (正式写手)

[求助] 一个偏微分方程组的求解问题 已有1人参与

在一篇文献中看到的某系统的方程,及相应的一些条件如图所示,文章中说用pdepe就可以求解,但是我就是解不出来呀,请大家看看有没有什么问题。
一个偏微分方程组的求解问题
360截图-3516750.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

这是我编写的matlab程序,大家都可以调一下,看看怎么样,总是说步长不够小,我已经调到很小了,比文献中说的要小好多,都是不能求解的,不知道是怎么回事。
function [c,f,s]=pdefun(x,t,u,dudx)
Pe_C=1;
Pe_T=1;
B_C=1;
B_T=2;
gamma=20;
beta_T=1;
u_c=0;
y=u(2)/(1+u(2)/gamma);
R_11=B_T*B_C*(1-u(1))*exp(y);
R_12=-beta_T*(u_c-u(2));
R_13=-B_C*(1-u(1))*exp(y);
c=[1 1]';
f=[1/Pe_C 1/Pe_T]'.*dudx;
s=[R_13-dudx(1) R_11+R_12-dudx(2)]';

function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t)
pl=[-ul(1),-ul(2)]';
ql=[1 1]';
pr=[0 0]';
qr=[1 1]';

function u0=icfun(x)
u0=[0 0]';
m=0;
t=linspace(0,1,500);
x=linspace(0,1,100);
sol=pdepe(m,@pdefun,@icfun,@bcfun,x,t);
2楼2014-01-23 19:57:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

这个方程应该是没有什么问题的,这是一个温度与浓度的传播方程,好多篇文献中都以这个方程为例进行讨论,尽管写成了不同的形式,但都是这个方程,不管什么形式我都求解不出来呀,怎么回事,大家看看。如果大家有兴趣我也可以相关文章贴出来。
3楼2014-01-23 20:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

自己顶一下
4楼2014-01-25 07:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfg25aw

新虫 (初入文坛)

【答案】应助回帖

最好用迎风格式去解。
5楼2014-02-06 10:57:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)

有人给我解答了,模型本身存在问题,因此求解不出,只是我看到的这个模型在很多的文献中都被作为例子进行讨论,而且不同的论文我都找了一下,没有一篇我是能够仿出来,真是奇怪,不知道这些作者是怎么弄的,我是无论如何也搞不出来。
这些论文中的参数没有差别,只是在最后的常数项前面的正负号有差别,我都试过了,还是不行呀。哎放弃了。
6楼2014-02-08 14:58:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ghw_nit 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 材料工程,326分,求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
请填处理意见