Skip to main content
P
ptc-62020951-Visitor
1-Visitor
June 19, 2014
Question

Odesolve

  • June 19, 2014
  • 1 reply
  • 1737 views

problem1.png

Hello eveyone..

Can someone please help me of solving this problem by using Odesolve.

I have to solve numerically the dimensionless mathematical model of unsteady pipe reactor where the reaction takes place the first order.

Best regards, Milica

P.S. Sorry for bad English

1 reply

H
HarveyHensley12-Amethyst
12-Amethyst
June 19, 2014

You have a PDE not an ODE so you need to use pdesolve. Pdesolve doesn't exist yet in Prime so you will have to use Mathcad 15 or earlier. Numol should also work if you change the boundary conditions as I discuss below.

I'm not sure, but I think pdesolve won't be able to handle the Dankwert's boundary conditions. You probably will need to set the inlet concentration = 1 and dc/dx=0 at x=0. The exit derivative = 0 is not realistic for all reactors. It is also not realistic to believe that a reactor would somehow know how to perform upstream in order to satisfy the exit condition. Thus, the boundary conditions need to be set at the inlet.

You also have an error in the first line: you have dcdt twice.

P
ptc-62020951-VisitorAuthor
1-Visitor
June 19, 2014

Thank you very much for your reply. You are right, it is PDE. I am trying to solve problem with Pdesolve but MathCad displays an error. Can someone please help me?pde.jpg

H
HarveyHensley12-Amethyst
12-Amethyst
June 19, 2014

The argument in the Pdesolve list should only be c, not c(x,y) or c(x,t). Also, your PDE still has too many dc/dt terms.

I have discovered that Pdesolve does allow the mixed boundary condition. In fact, that is the only boundary condition that you need. The exit boundary condition you have will lead to a non-realistic result. It will force the reaction to complete regardless of the reactor length. Try it.

Update: The c variable being assigned to the pdesolve result should also not have arguments.

c:= pdesolve(c,x,....

Terms & ConditionsAccessibility statement