- #386

- 475

- 26

Unfortunately the molar flow error didnt fix the output, although it did change it somewhat.

The molar flow out of a tank should have been defined as the below. It was previously this so I must have changed it somewhere along the way:

\begin{equation}

\dot{m}_j = \dot{m}_{j-1} + \frac{\rho_m \Delta z A_c}{T_g}\frac{dT}{dt} - (M_{CO_2}''*A_s)

\end{equation}

Further up in the code ##m_j## is defined as:

\begin{equation}

m_j = \frac{P}{RT_g}\epsilon A_C\Delta z

\end{equation}

So I now have:

\begin{equation}

\dot{m}_j = \dot{m}_{j-1} + \frac{\rho_m \Delta z A_c}{T_g}\frac{dT}{dt} - (M_{CO_2}''*A_s)

\end{equation}

which I've implemented as:

\begin{equation}

\dot{m}_j = \dot{m}_{j-1} + \frac{m_j \epsilon}{T_g}\frac{dT}{dt} - (M_{CO_2}''*A_s)

\end{equation}

Just posting this here in case I've done anything crazy there