Upper limit or saturation
1 次查看(过去 30 天)
显示 更早的评论
Hi. I am currently trying to model a chemical equilibrium where there is an in flow of matter using the simbiology add-on.
There is an upper limit of one of the species, so my question is how do i specify an upper limit or saturation of a species in simbiology?
0 个评论
采纳的回答
Ingrid Tigges
2013-9-9
Should the upper limit come into play when running a simulation? In this case, you could set the upper limit using an event. Have a look at http://www.mathworks.de/de/help/simbio/ug/event-object.html
2 个评论
Ingrid Tigges
2013-9-10
Personally, I would contact technical support so that someone can look into your model.
更多回答(1 个)
Arthur Goldsipe
2013-9-10
编辑:Arthur Goldsipe
2013-9-10
Hi Peter,
The "right" way to implement an upper limit or saturation really depends on how you want to physically model the effect. What should happen to any reactions that produce this species after it reaches saturation? Should these reactions "turn off"?
As Ingrid said, the most natural way to model saturation is probably using an event. From your comment to her answer, you may be happy creating an event that simply ensures that your species never exceeds the upper limit. However, because of the way that events are detected, this event will not do what you want unless you also "turn off" all reaction that produce it.
It also sounds like your model uses dimensional analysis and/or unit conversion. In that case, here's how you might implement the event:
Let's assume that the species you want to limit is named D, has units of mole/liter, and is in a compartment named C. Let's also assume that D is the product of a reaction with a rate constant K.
First, add a parameter (with units) to the model that sets the saturation limit, for example, DMax = 3.9e-6 mole/liter. Next, add an event with the trigger "D > DMax" and with TWO event functions "C.D = DMax" and "K = 0*K". The first function ensures that D is set back to saturation, while the second function ensure that the reaction producing D is turned off. Note that the trigger must be written using an INEQUALITY not an EQUALITY. Also, you'll ned to set the ConstantValue property of K to false so that you can change its value by an event.
Any more specific suggestions would require knowledge of your model, and as Ingrid said, would best be handled by contacting technical support.
-Arthur
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Biotech and Pharmaceutical 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!