How to imply conditional boundary condition in PDEPE?

3 次查看（过去 30 天）

显示更早的评论

Sait Mutlu Karahan 2022-5-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1730280-how-to-imply-conditional-boundary-condition-in-pdepe

评论： Torsten 2022-5-31

Hi all,

I have problem my left side boundary condition is time dependent and it is

when I try to apply pdepe in matlab I couldn't write my pl values according the time values.

P.S. I'm a new learner please be kind when you are replying. Thanks :)

clear all;

clc;

close all;

L = 1000;

x = linspace(0,L,1001);

t = [linspace(0,5,20), linspace(5.1,50,10)];

m = 0;

sol = pdepe(m,@heatpde,@heatic,@heatbc,x,t);

colormap hot

imagesc(x,t,sol)

colorbar

xlabel('Distance x','interpreter','latex')

ylabel('Time t','interpreter','latex')

title('Heat Equation for $0 \le x \le 1000$ and $0 \le t \le 50$','interpreter','latex')

function [c,f,s] = heatpde(x,t,u,dudx)

c = 0.02;

f = dudx;

s = 0;

end

function u0 = heatic(x)

u0 = 50-(0.01*x);

end

function [pl,ql,pr,qr] = heatbc(xl,ul,xr,ur,t)

g=50-0.1*(t-1);

if (t >= 0) && (t <= 20)

pl = ul-50;

ql = 0;

elseif (t > 20) && (t <= 50)

pl = ul - g;

ql = 0;

end

pr = ur-40;

qr = 0;

end

8 个评论
显示 6更早的评论隐藏 6更早的评论

Sait Mutlu Karahan 2022-5-31

yes I'm talking about that but I couldn't find how to do that

Torsten 2022-5-31

在 MATLAB Online 中打开

t = [t1;t2];
sol = [sol1(:,:,1);sol2(:,:,1)];

请先登录，再进行评论。

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

回答（0 个）

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

类别

MATLAB Mathematics

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

产品

MATLAB

版本

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by