Skip to main content
1-Visitor
March 11, 2018
Solved

Problem with plotting results from odesolve

  • March 11, 2018
  • 1 reply
  • 4250 views

Hello,

 

I am working in mathcad prime 3.0 on an assignment and I am modeling a Trickle Bed Reactor (plug flow behavior). I have set up the mole balances, from which I know by now that they are correct. I use the odesolve function to solve a set of differential (and algebraic) equations. However when I plot the results I do not get the results I want. Attached is the file I made. I made a plot with the liquid concentrations of the components. I see only a difference in the concentration of A and I am wondering what I am doing wrong, because i should also a difference in the concentration of B and C.

 

Below a brief description on the case:

I have three reactions. The third reaction does not take place due to a reaction rate constant which is zero. 

  1. A + H2 -> B
  2. B + H2 -> C
  3. A -> G (This reaction is not taking place)

 

I am feeding hydrogen in the gas phase to the reactor together with liquid component A. I have set up 9 mole balances. The solid phase is in steady state, so, these are algebraic equations:

  1. Gas phase hydrogen
  2. Liquid phase hydrogen
  3. Solid phase hydrogen
  4. Liquid phase A
  5. Solid phase A
  6. Liquid phase B
  7. Solid phase B
  8. Liquid phase C
  9. Solid phase C

Maybe I am forgetting something trivial, but i'm am stuck for a while now. Any help would be appreciated. To be clear: I am wondering why i don't see a concentration difference for component B and C in the plot.

 

Thanks in advance

 

Kind regards,

Xafie

 

pc. unfortunately i can not attach the file. If needed i can also send the file via email

Best answer by Fred_Kohlhepp

Replace

Capture.PNG

with

CaptureA.PNG

and get

CaptureB.PNG

1 reply

Xafie1-VisitorAuthor
1-Visitor
March 20, 2018

Update:

 

Attached you'll find the mathcad file

 

24-Ruby IV
March 20, 2018

@Xafie wrote:

Update:

 

Attached you'll find the mathcad file

 


or see here http://twt.mpei.ac.ru/TTHB/New-Chem-Kin/En-ChemKin.html

23-Emerald I
March 22, 2018

Replace

Capture.PNG

with

CaptureA.PNG

and get

CaptureB.PNG