Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'Overdetermined NDSolve' topicStudent Support Forum > General > Archives > "Overdetermined NDSolve"

Help | Reply To Topic
Author Comment/Response
Burcu Ugur
03/21/13 5:27pm

Hello,

I am trying to model an adsorption process by doing mass and energy balances. I have 3 mass and 2 energy balances, and 5 dependent variables. However, mathematica tells me that my system is overdetermined.

Can you please help me what I did wrong?
I know equations does not look well therefore I attached the notebook as well.

Thank you so much.

NumericalSolution=NDSolve[{(1-ep)*rs*D[q[z,r,t],t]=ka*(3/rp)*(cp[z,r,t]-cpe), ec*Dz*D[cg[z,r,t],{z,2}]-vg*D[cg[z,r,t],z]=ec*D[cg[z,r,t],t]+ka*(3*(1-c)/rp)*(cg[z,r,t]-cp[z,rp,t]), ep*D[cp[z,r,t],t]==ep*Deff*(2/r)*D[cp[z,r,t],{r,2}]-(1-ep)*rs*D[q[z,r,t],t],k*D[Tg[z,r,t],{z,2}]-ec*vg*rg*Cpg*D[Tg[z,r,t],z]-(2*hfd/rc)*(Tg[z,r,t]-Tw[z,r,t])-n*D[Tg[z,r,t],t]-Hads*rp*D[q[z,r,t],t]=0, ((ro^2)-(rc^2))*rw*Cpw*D[Tw[z,r,t],t]=2*rc*hfd*(Tg[z,r,t]-Tw[z,r,t])-2*ro*ho*(Tw[z,r,t]-T0), cg[z,r,0]=0, cg[0,r,t]=cginlet*(1-Exp[-(t/tao)]), Derivative[1,0,0][cg][h,r,t]=0, q[z,r,0]=0, cp[z,r,0]=0, Derivative[0,1,0][cp][z,0,t]=0, -Deff*Derivative[0,1,0][cp][z,rp,t]=ka*(3/rp)*(cg[z,r,t]-cp[z,rp,t]), Tg[z,r,0]=T0, Tg[0,r,t]=T0, Derivative[1,0,0][Tg][h,r,t]=0, Tw[z,r,0]=T0}, {cg,cp,q,Tg,Tw}, {z,0,h}, {r,0,rp}, {t,0,8000}, Method{"MethodOfLines", "SpatialDiscretization"{"TensorProductGrid", "MinPoints"150,PrecisionGoal3}}]

Attachment: Support NDSolve.nb, URL: ,
Help | Reply To Topic