cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Newbie problem solving system of 3 ODE's

SOLVED
Highlighted
Level 7

Newbie problem solving system of 3 ODE's

 

Hello,

 

I was able (with some help from the community Smiley Happy to solve a single ODE using a MathCAD 5.0 Solve Block but now I want to solve a system of 3 ODE's but I am getting an error (see attached file).

 

I suspect I am making some newbie error. If someone can enlighten me as to what that error is and tell me how to do the calculation correctly, it would be much appreciated.

 

TIA,

 

--Doug

 

1 ACCEPTED SOLUTION

Accepted Solutions

Re: Newbie problem solving system of 3 ODE's

The problem is the equation for px'(t).

If you replace a for something dimensionless it works.

If k1 has the unit L^2/(min*mol^2) it works.

So something with that equation seems to be wrong.

7 REPLIES 7

Re: Newbie problem solving system of 3 ODE's

You have to apply units, even if the value is zero.

So to be on the save side you should write px'(0 s)=...

 

But you have a unit problem/mismatch in your equations:

You have an expression (1-px(t)) which means px(t) is dimensionless!

So your equation MeOH'(t)=-px'(t) meanst that MeOH(t) must be dimensionless, too.

But one of your IC says MeOH(0 s)=a*rx and this would mean unit mol/Litre which is NOT dimensionless.

So that's a discrepancy.

Maybe the "1" in (1-px(t)) is not unitless? But then there would be other discrepancies, I think.

Re: Newbie problem solving system of 3 ODE's

OK, I found another error which caused the "unknown error" you experience. You typed a capital "P" instead of a lower case "p" as argument of odesolve.

The solve block still fails as of the aforementioned unit mismatches.

I tried your solveblock without any units at all and it works OK.

 

 

Re: Newbie problem solving system of 3 ODE's

 

Thanks for looking at this for me.

 

The "px" variable is a dimensionless fractional conversion.

 

You are correct that both the MeOH and VOH variables have units of mole/L (i.e. concentration) which meant that the MeOH' derivative was not correctly defined in my original file. I also added "min" to the IC specifications.

 

But I still get an error for the first IC spec.

 

I have uploaded a revised version with MeOH' defined correctly with units of mole/L and "min" added to the IC spec's.

 

Thanks again for looking at this for me.

 

--Doug

Re: Newbie problem solving system of 3 ODE's

 

I did fix the typo with the capital "P" but then hit the new units of measure issue.

 

Thanks for point that out.

 

I am trying to solve this system for a colleague. I need to get back to him and get the units of measure issues resolved.

Re: Newbie problem solving system of 3 ODE's

The problem is the equation for px'(t).

If you replace a for something dimensionless it works.

If k1 has the unit L^2/(min*mol^2) it works.

So something with that equation seems to be wrong.

Re: Newbie problem solving system of 3 ODE's

 

Thanks very much for your help.

 

I spent some more time with my colleague (whose problem this is) until we got all the units of measures issues resolved.

 

Once we did that (and assigned units for all the IC's), MathCAD solved it properly (see attached file).

 

Thanks again!

 

--Doug

Re: Newbie problem solving system of 3 ODE's

Ah, I see the change in the equation (a was wrongly taken as factor originally).

Glad you are up and running!