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

Community Tip - Want the oppurtunity to discuss enhancements to PTC products? Join a working group! X

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
24-Ruby V
(To:ddenholm)

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.

10 REPLIES 10
24-Ruby V
(To:ddenholm)

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.

24-Ruby V
(To:Werner_E)

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.

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.

24-Ruby V
(To:ddenholm)

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.

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

24-Ruby V
(To:ddenholm)

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

Glad you are up and running!

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

24-Ruby V
(To:ME_10255445)

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

6-Contributor
(To:Werner_E)

Thanks alot

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

Announcements
Top Tags