Hey guys,
I'm quite new here and hope i'm in the right section. In trying to solve a PDE system with NDSolve the plots I obtain in the end are empty. I read through quite a few posts regarding this issue but haven't found a solution that I could apply to my code.
The PDE system:
deqs = {D[m[x, t], t] == D[vm[x]/am[x] D[vm[x] m[x, t], x], x] + f[{m[x, t], c[x, t]}],
D[c[x, t], t] == D[vc[x]/ac[x]*D[vc[x]*c[x, t], x], x] - f[{m[x, t], c[x, t]}],
m[x, 0] == M, c[x, 0] == C0, m[0, t] == m[L, t] , c[0, t] == c[L, t]};
Most important code:
sol1 = NDSolve[deqs, {m, c}, {x, 0, L}, {t, 0, 10}]
Plot3D[m[x, t] /. sol1, {x, 0, L}, {t, 0, 10}]
I assume my mistake already is somewhere in these lines, f is a polynome, vc and vm are constants and ac and am are arctan(x) like functions.
To be precise:
f[{a_, b_}] = g*a^2*b - d*a
am[x_]:= (A(ArcTan[B(x-L/2)]+ Pi/2))^(-1);
ac[x_]:=(A (ArcTan[B(-x+L/2)]+Pi/2))^(-1);
vm[x_] := vM;
vc[x_] := vC;
Does anyone have an Idea why this doesn't work?
Thanks a lot in advance!
Attachments:
