Trouble with understanding "storage" concentrations in channels without flow (MCT model)


I have encountered some issues when fitting my data and it led me down a rabbit hole of confusion and what I think are mistakes in some general assumptions I had.

I will start in the beginning:

I am using the MCT unit to simulate tracer transport in plants and my data comes from PET scans containing 3-D data of tracer activity over time.

I have a setup in CADET-process with two MCT units (MCT1, MCT2) connected to each other via port 0 (channel 0). Both have an exchange parameter that allows storage in channel 1. There is no flow in channel 1.


When I choose the same storage parameter for both units, I will get something like this for the bulk concentration at 8 spatial regions of interest, 4 of which are in the first unit (MCT1) and the other 4 in the second unit (MCT2). The concentrations of all channels are summed up, as this matches the data I observe:


When I choose a higher storage parameter for the second unit (MCT2) I now get something like this:


Now this is something I would never see in my data, as I am observing proton decay directly relating to the number of labeled C-atoms in my tracer and that cannot become suddenly more at later times.

I tried to get deeper into this issue and had a look at the outlet (which is just the last spacial cell) of the MCT1 unit:

Channel 0:

Channel 1:

And I assume the problem has something to do with me not understanding, what the concentration in Channel 1 exactly means. I understand that, since there is no flow-rate there is kind of an accumulation happening. But still I don’t see, why for the upper example the concentration should rise to 2 mM/L when the original signal I put in was at 1 mM/L.

So now I am on the one hand trying to understand this issue more and on the other hand trying to find a solution to “equalize” several units behind each other with differing storage parameters.

Poplar_velocity_change_one_func_single_obj_problem.ipynb (299.1 KB)

Hi Hannah,

the max. concentration in channel 1 is not limited by the amplitude of the input signal. Since you’ve specified the exchange from channel 1 to channel 0 as 0.0, the concentration in channel 1 will rise as long as theres some concentration in channel 0.

Iam not sure what you mean by equalizing the units? Do you want ROI 1-4 and ROI 5-8 to be the same?


~ topic can be closed ~