Different dilution effect on the kinetics of a time course
5 views (last 30 days)
Show older comments
Luis B. Walter
on 14 Jun 2024
Commented: Luis B. Walter
on 17 Jun 2024
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
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.
Accepted Answer
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
See Also
Categories
Find more on Extend Modeling Environment in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!