24小时热门版块排行榜    

查看: 1663  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 321求调剂 +4 大米饭! 2026-03-15 4/200 2026-03-16 08:41 by Linda Hu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 调剂 +8 调剂的考研学生 2026-03-09 8/400 2026-03-15 22:14 by Winj1e
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见