Mixed boundary condition in pdepe solver
4 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
I'm solving the following 1D convection diffusion system in pdepe solver.

I've to impose no-flux boundary condition at both left and rigth boundary.
function sol=main(Input)
format short
global D v
m = 0;
D = Input.D;
v = Input.v;
x = Input.xmesh;
t = Input.tspan;
sol = pdepe(m,@pdefun,@icfun,@bcfun,x,t)
function [g,f,s] = pdefun(x,t,c,DcDx)
g = 1;
f = D*DcDx;
s = -v*DcDx;
end
function c0 = icfun(x)
c0 = Input.c0;
end
I'm not sure how to implement the no-flux boundary condition according to the general form given here
Is it correct to re-write the BC: - vC + DdC/dx = 0
as
C -( D/v) dC/dx = 0
where p = C
q = -(1/v)
f = DdC/dx
Is the following correct?
function [pl,ql,pr,qr] = bcfun(xl,cl,xr,cr,t)
pl = cl;
ql = -1/v;
pr = cr;
qr = -1/v;
end
or
function [pl,ql,pr,qr] = bcfun(xl,cl,xr,cr,t)
pl = -cl;
ql = 1/v;
pr = -cr;
qr = 1/v;
end
0 Kommentare
Antworten (0)
Siehe auch
Kategorien
Mehr zu Mathematics and Optimization finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!