Please help! pdepe does not stop!!

Question

Hüseyin Cagdas Yatkin 2019-12-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/497124-please-help-pdepe-does-not-stop

回答： Guru Mohanty 2020-1-23

Hello guys, I trying solve below question with pdepe and whe I run the code as mantioned in the below code, it never stop working.

v = 49 680 m/day

kd = 0.24 1/day

ks = 58,8 m/day

h = 5 m

C(t = 0,x) = 13.75 mg/L initial condition

C(t, x = 0) = 13.75 mg/L

C(t, x = 1000m) = 10.05 mg/L

clc
clear
x = linspace(0,1000,1000);
t = linspace(0, 7, 28);
m = 0;
sol = pdepe(m, @pde,@pdeic,@pdebc,x,t);
u = sol(:,:,1);
sol(1000,:);
plot(x, sol(1000,:))
function [c, f, s] = pde(x, t, C, dCdx)
c = 1/(0.575*24*3600);
f = C;
s = (0.24+2.45*24/5)/(0.575*24*3600)*C;
end
function u0 = pdeic(x)
u0 = 13.75;
end
function [pl,ql,pr,qr] = pdebc(xl,ul,xr,ur,t)
pl = ul - 13.75;
ql = 0;
pr = ur - 10.05;
qr = 0;
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Guru Mohanty 2020-1-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/497124-please-help-pdepe-does-not-stop#answer_411600

在 MATLAB Online 中打开

Hi, I understand you are trying to solve this partial differential equation. You can do this by using pdepe. Here is a sample code.

clc
clear all;
x = linspace(0,1000,1001);
t = linspace(0, 7, 28);
m = 0;
sol = pdepe(m, @pde,@pdeic,@pdebc,x,t);
val=sol(28,:);
figure
surf(t,x,sol');
figure
plot(x,val);
function [u, f, s] = pde(x,t,u,dudx)
u = 1/(0.575*24*3600);
f = u;
s = (0.24+2.45*24/5)/(0.575*24*3600)*u;
end
function u0 = pdeic(x)
u0 = 13.75;
end
function [pl,ql,pr,qr] = pdebc(xl,ul,xr,ur,t)
pl = ul - 13.75;
ql = 0;
pr = ur - 10.05;
qr = 0;
end