topic Re: Need help. Mathcad prime 7.0 : Parameter Estimation in Chemical Kinetics. in PTC Mathcad
https://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719627#M195017
<P>To start:</P>
<P>Your first of 14 expressions regarding C0(t) is simple. C0(t)=C_0*e^(-k0*t), where C_0 is a constant that can be found with the initial condition C0(0)=0.013, if t=0 then C0(t)=C_0*e^0 = C_0, so C_0 = 0.013. Thus:</P>
<P>C0(t)=0.013*e^(-k0*t)</P>
<P>Your third expression says that C2'(t)=k0*C0(t) = k0*0.013*e^(-k0*t). So you can find C2(t) by integrating k0*0.013*e^(-k0*t), and the integration constant is found through the initial condition C2(0)=0.</P>
<P> </P>
<P>Problems start to occur with your initial condition for C6(t). You state that C6(0)=C6e0, which is found from the data array. Now that makes C6(0) an array, where it can/should have only a single value. Besides it has a unit (mol/l) whereas none of the other initial conditions have a unit. That again may not be a problem, as long as they are 0, but the non-zero initial conditions might get you into unbalanced units.</P>
<P>Going further, I see that C9 and C10 are exactly equal. They're both defined with the same rhs, and their initial conditions are both 0. This will not help to solve the set of functions.</P>
<P> </P>
<P>So: make sure your units are balanced, sort out the initial condition for C6(t), such that it is a single value. You may want to make it a parameter. And check your equations.</P>
<P> </P>
<P>Success!<BR />Luc</P>
<P> </P>Thu, 18 Mar 2021 21:20:45 GMTLucMeekes2021-03-18T21:20:45ZNeed help. Mathcad prime 7.0 : Parameter Estimation in Chemical Kinetics.
https://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719487#M194997
<P>Help to sort out the problems. I have been using the program recently; I need to solve the estimation of parameters in chemical kinetics. Help to understand what the error is and how to solve it.</P>Thu, 18 Mar 2021 10:05:16 GMThttps://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719487#M194997K4L3KT0R2021-03-18T10:05:16ZRe: Need help. Mathcad prime 7.0 : Parameter Estimation in Chemical Kinetics.
https://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719512#M195004
<P>What error are you concerned with?</P>
<P> </P>
<P>Success!<BR />Luc</P>Thu, 18 Mar 2021 11:23:51 GMThttps://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719512#M195004LucMeekes2021-03-18T11:23:51ZRe: Need help. Mathcad prime 7.0 : Parameter Estimation in Chemical Kinetics.
https://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719529#M195007
<P> </P>
<P>This </P>
<P><span class="lia-inline-image-display-wrapper lia-image-align-inline" image-alt="K4L3KT0R_2-1616069236905.png" style="width: 400px;"><img src="https://community.ptc.com/t5/image/serverpage/image-id/39061i8D9B59B162EB47F1/image-size/medium?v=v2&px=400" role="button" title="K4L3KT0R_2-1616069236905.png" alt="K4L3KT0R_2-1616069236905.png" /></span></P>
<P> </P>
<P> </P>Thu, 18 Mar 2021 12:08:11 GMThttps://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719529#M195007K4L3KT0R2021-03-18T12:08:11ZRe: Need help. Mathcad prime 7.0 : Parameter Estimation in Chemical Kinetics.
https://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719530#M195008
<P>And do you see any other errors?If yes, please tell me which ones, I will be very grateful!</P>Thu, 18 Mar 2021 12:10:15 GMThttps://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719530#M195008K4L3KT0R2021-03-18T12:10:15ZRe: Need help. Mathcad prime 7.0 : Parameter Estimation in Chemical Kinetics.
https://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719627#M195017
<P>To start:</P>
<P>Your first of 14 expressions regarding C0(t) is simple. C0(t)=C_0*e^(-k0*t), where C_0 is a constant that can be found with the initial condition C0(0)=0.013, if t=0 then C0(t)=C_0*e^0 = C_0, so C_0 = 0.013. Thus:</P>
<P>C0(t)=0.013*e^(-k0*t)</P>
<P>Your third expression says that C2'(t)=k0*C0(t) = k0*0.013*e^(-k0*t). So you can find C2(t) by integrating k0*0.013*e^(-k0*t), and the integration constant is found through the initial condition C2(0)=0.</P>
<P> </P>
<P>Problems start to occur with your initial condition for C6(t). You state that C6(0)=C6e0, which is found from the data array. Now that makes C6(0) an array, where it can/should have only a single value. Besides it has a unit (mol/l) whereas none of the other initial conditions have a unit. That again may not be a problem, as long as they are 0, but the non-zero initial conditions might get you into unbalanced units.</P>
<P>Going further, I see that C9 and C10 are exactly equal. They're both defined with the same rhs, and their initial conditions are both 0. This will not help to solve the set of functions.</P>
<P> </P>
<P>So: make sure your units are balanced, sort out the initial condition for C6(t), such that it is a single value. You may want to make it a parameter. And check your equations.</P>
<P> </P>
<P>Success!<BR />Luc</P>
<P> </P>Thu, 18 Mar 2021 21:20:45 GMThttps://community.ptc.com/t5/PTC-Mathcad/Need-help-Mathcad-prime-7-0-Parameter-Estimation-in-Chemical/m-p/719627#M195017LucMeekes2021-03-18T21:20:45Z