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

Community Tip - Did you know you can set a signature that will be added to all your posts? Set it here! X

Newbie problem solving system of 3 ODE's

ddenholm
6-Contributor

Newbie problem solving system of 3 ODE's

 

Hello,

 

I was able (with some help from the community 🙂 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

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.

View solution in original post

10 REPLIES 10

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.

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.

 

 

ddenholm
6-Contributor
(To:Werner_E)

 

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.

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.

ddenholm
6-Contributor
(To:Werner_E)

 

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

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

Glad you are up and running!

 

ME_10255445
6-Contributor
(To:ddenholm)

Hi guys,

How are you doing?

Can anyone please send this in a pdf or whatever else bcs I can't read mcdx files due to issues with my laptop

 

Thx in advance

Here is a pic showing the content of Dougs latest worksheet (with some regions rearranged for compactness):

Werner_E_0-1652288411659.png

 

ME_10255445
6-Contributor
(To:Werner_E)

Thanks alot

ddenholm
6-Contributor
(To:Werner_E)

 

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

Top Tags