Inputs: rho = 0.414; c = 2029; m = 0.024; U = 013.811; A = 0.00405; Tsol = 300.6; Tin = 358; Tout = 354; L = 0.334; t200 = 1200; X1 = .01;
pde = (rhocD[V[t, X], t] == mD[V[t, X], X, X] + UA*(V[t, X] - Tsol));
sol = NDSolve[{pde, V[0, X] == 300, V[t, 0] == Tin, D[V[t, L]] == Tout}, V, {t, 0, 2000}, {X, 0.001, L}];
V[t200, X1] /. sol
plot[V[t, X] /. sol, {t, 0, 1000}, {X, 0.001, L}]
Outputs:
{357.737}
plot[{InterpolatingFunction[{{0., 2000.}, {0., 0.334}}, <>][t,
X]}, {t, 0, 1000}, {X, 0.001, 0.334}]
It is not giving me a plot. why is it so? Can you help me finding the error in my coding or conditions.
Thank You