MATLAB的答案

Pdepe:空间离散化已经失败了。离散化只支持抛物线和椭圆方程,万博1manbetx通量项涉及空间导数。

1视图(30天)
路易吉皮耶罗Di鲣鱼
路易吉皮耶罗Di鲣鱼 2020年10月3日
编辑: 路易吉皮耶罗Di鲣鱼2020年10月3日
你好,
我试图解决pd系统关于扩散的催化剂(失活模型建模),收益率:
错误使用pdepe(第293行)
空间离散化已经失败了。离散化只支持抛物线万博1manbetx
和椭圆方程,通量项涉及空间导数。
错误分辨率(30行)
索尔= pdepe (m, @eqn2, @initial2 @bc2 x, t);
在attchament阅读模式后,请参见下面的代码:
全球曹Def T Dp R k kd kO cW alphaOγ
%的参数
T = 333;% = K
Dp = 2;% = m
Def = 10 ^ 5;%
曹= 0.0001;% (=)kmol / m ^ 3开始索绪尔语言学的浓度
连续波= 0.0001;% (=)kmol / m ^ 3
R = Dp / 2;
k = (1.191 * 10 ^ 5) * exp (-7.544 / (0.001607460438947 * T));
kd = (4.299 * 10 ^ 3) * (exp (-6.86 / (0.002102114315754 * T)));
kO = k *连续波;
γ= 1;
alphaO = 1;%开始催化剂活动的价值
% PDE2: MATLAB脚本解决PDE的m文件
存储在eqn2 %。米,bc2。m, initial2.m
m = 2;%球状的
x = linspace (0, R, 10);
t = linspace (0, 1, 10);
索尔= pdepe (m, @eqn2, @initial2 @bc2 x, t);
u1 =索尔(:,:1);
u2 =索尔(:,:,2);
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% EQN2: MATLAB的m文件,其中包含的系数
%两个PDE的系统在时间和空间维度。
函数[c、b s] = eqn2 (x, t, u, DuDx)
全球Def kO kdγ
c = [1;1);
b = [Def;0]。* DuDx;
s = [- kO * u (1) * (u (2)); kd * u (1) * ((u(2)) ^(γ+ 1)));
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% INITIAL2: MATLAB函数m文件,它定义了初始条件
%的系统两个PDE在时间和空间变量。
函数值= initial2 (x);
全球曹alphaO
值=[曹;alphaO];
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% BC2: MATLAB函数m文件,它定义了边界条件
%的系统两个PDE在时间和空间维度。
函数[pl, ql,公关,qr] = bc2 (xl, ul, xr, ur, t)
全球曹R
pl = [0;0);
ql = [1;0);
公关= (ur (R)曹;0);
qr = [0;0);

在回答这个问题。

接受的答案

比尔•格林
比尔•格林 2020年10月3日
编辑:比尔•格林 2020年10月3日
你的第二个PDE的边界条件是无效的。这些工作:
pl = [0;0);
ql = [1;1);
公关= (ur(1)曹;0);
qr = [0;1);
3评论
显示2年长的评论

登录置评。

更多的答案(0)

类别

找到更多的在特征值问题在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!