Hello,
I am trying to solve a system of differential equations to calculate the change of temperature and pressure in switchgear during an arc fault. The equations are as follows:
where C2 equals:
"T" is temperature and "p" is the pressure, PL is the power of the arc which is a function of time and is calculated earlier. The rest are constants.
I've attemped to calculate this in Mathcad Prime 8 using the rkfixed function, inputting the full equations with the C2 factored into the them:
I've input "T" as "y0" and "p" as "y1" in the equations using the matrix index operator. The first equation if for temperature and the second for pressure. The initial conditions are for the initial temperature of 300 and pressure of 101325.
I'm getting an unknown error which roughly translates to "calling object expection". It does seem to dissapear when i replace one of the "y0" or "y1" variables with a constant in the second equation, but not in the first. However the results are then incorrect.
I would like to ask if it's even possible to solve such a system of equations using this method, or is there something else entirely I should be doing?
If it is possible, then what might be causing this error?
I'm attaching the worksheet, the shown above equations are at the botton of page 2. On page one i've calulated the currents and power "PL".
Thank you for any help!
ACCEPTED SOLUTION
Accepted Solutions
This should be possible.
I'd have used Odesolve in a solve block, but rkfixed can also be used.
In an English version the error message is:
Your function D takes the root of a difference. That difference might become negative, to result in imaginary (non-real) roots...
You say that "The rest are constants.", but you've replaced c.v with an expression involving T (y0)....
Oh, I see now that your value tau a is imaginary:
What's this value 3.14 do in your expressions for iL1 and iL3? Shouldn't that be pi?
Attached is my Odesolve attempt. It does produce a result.... I may have made one or more mistakes.
Success!
Luc
This should be possible.
I'd have used Odesolve in a solve block, but rkfixed can also be used.
In an English version the error message is:
Your function D takes the root of a difference. That difference might become negative, to result in imaginary (non-real) roots...
You say that "The rest are constants.", but you've replaced c.v with an expression involving T (y0)....
Oh, I see now that your value tau a is imaginary:
What's this value 3.14 do in your expressions for iL1 and iL3? Shouldn't that be pi?
Attached is my Odesolve attempt. It does produce a result.... I may have made one or more mistakes.
Success!
Luc
Indeed, sorry I forgot to mention that "cv" is a function of T and equals to 0.094*T+682.06.
The 3.14 is meant to be an approx. of pi and for τa I just didnt think of using the determinant operator, but i got that fixed. I've also corrected some other simplifications i made earlier while calculating the currents.
Your solution is working after I put the "α" constant back in the equation, since I forgot I had it removed before while I was trying to fix the error:
It's meant to be a factor reducing the gas flow out the decompression flaps.
I also multiplied the T(t) in the second equation by 2:
The results are finally making some sense 🙂
The "A" is the area of the pressure relief flaps and "V" is the volume of the compartment in which the arc is ignited. Changing them seems to affect the results accordingly.
I'm attaching the worksheet where i seem to have the mistakes corrected.
Thank you very much for the help!
