Solving advection diffusion pde

16 次查看（过去 30 天）

Deepa Maheshvare 2018-12-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437177-solving-advection-diffusion-pde

评论： Deepa Maheshvare 2019-1-2

I want to solve the above pde with the given boundary and initial conditions. I came across the pdepe function in MATLAB.

I had a chance to look at the example given here . I couldn't understand how pdex1pde function has to be defined for my case.

Could someone help?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Deepa Maheshvare 2018-12-27

在 MATLAB Online 中打开

Hi All,

Is the following implementation correct?

function DiffusionConvection
m = 0;
x = linspace(0,62,10);
t = linspace(0,10,100);
sol = pdepe(m,@pdefun,@icfun,@bcfun,x,t)
function [g,f,s] = pdefun(x,t,c,DcDx)
D = 900;
v = 10;
g = 1;
f = D*DcDx;
s = -v*DcDx;
end
function c0 = icfun(x)
c0 = 80;
end
function [pl,ql,pr,qr] = bcfun(xl,cl,xr,cr,t)
pl = cl -10;
ql = 1;
pr = cr;
qr = 1;
end
end

请先登录，再进行评论。

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

回答（1 个）

Bill Greene 2018-12-27

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437177-solving-advection-diffusion-pde#answer_354140

在 MATLAB Online 中打开

The only problem I see with your code is in the boundary conditions. I corrected your bcfun function and have attached my version of your code below.

function DiffusionConvection
m = 0;
x = linspace(0,62,10);
t = linspace(0,10,100);
D = 900;
sol = pdepe(m,@pdefun,@icfun,@bcfun,x,t);
    function [g,f,s] = pdefun(x,t,c,DcDx)
    v = 10;
    g = 1;
    f = D*DcDx;
    s = -v*DcDx;
    end
    
    function c0 = icfun(x)
    c0 = 80;
    end
    
    function [pl,ql,pr,qr] = bcfun(xl,cl,xr,cr,t)
    pl = -10*D;
    ql = 1;
    pr = 0;
    qr = 1;
    end
end