Clear Filters
Clear Filters

Different dilution effect on the kinetics of a time course

5 views (last 30 days)
Hi to everyone!
To simplify the explanation and the problem, I simulated the kinetics of an irreversible first-order reaction, A -> B. I implemented it in two independent compartments, R and P. I simulated the effect of a dilution in R by doubling at t= 0,1 the R volume. I programmed in P that, at t = 0.1, the instantaneous concentration of A and B would be reduced by half. I am sending an attach with the implementation of these simulations in the Simbiology interface.
When the simulations of the two compartments are plotted, it can be seen that the responses are not equal. That is, from t = 0.1 s, the reaction follow an exponential function in R with half of the initial amplitude and half of the initial value of k1. That is, the relaxation time is doubled. Meanwhile, in P, from t = 0.1, the reaction follows exponential kinetics with half the amplitude value but maintaining the initial value of k = 10. Without a doubt, the correct simulation is the latter (compartment P) where only the effect is observed in the amplitude and not in the relaxation time. Could you tell me what the error is that makes these kinetics that should be equal not be?
Thank you in advance!
Luis B.
  2 Comments
Umar
Umar on 15 Jun 2024
To address the issue and ensure consistent simulation results between compartments R and P, you need to review the code that handles the dilution and concentration reduction processes. Check for any inconsistencies or errors in how these processes are applied to each compartment.
Luis B. Walter
Luis B. Walter on 17 Jun 2024
Thanks, Umar.
I didn't write any code on my part. I just implemented three Events (as can be see in the attachment) in the Simbiology interface for this elemental reaction:
  1. time >= 0.1, R = 2*R
  2. time >= 0.1, P.A = P.A/2
  3. time >= 0.1, P.B = P.B/2
Luis B.

Sign in to comment.

Accepted Answer

Arthur Goldsipe
Arthur Goldsipe on 17 Jun 2024
Edited: Arthur Goldsipe on 17 Jun 2024
Hi Luis,
Because you don't use units in your model, SimBiology interprets the reaction rates to have dimensions of amount/time. (This is documented here.) The easiest way to fix this is to add units of 1/second to your rate constants. Then, the two appoaches for dilution behave identically.
One debugging tip mentioned on the above page is to view the differential equations for your model. You can find more details on how to do that in the SimBiology Model Builder app here.
-Arthur

More Answers (0)

Communities

More Answers in the  SimBiology Community

Categories

Find more on Extend Modeling Environment in Help Center and File Exchange

Tags

Products


Release

R2023b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!