Search the Community

Hi, We have been working with specifying custom rate laws for microbial respiration in React, and have been following the example specified on Page 67 of the reaction modelling guide. This integrates the differential equation for the time rate of change of the biomass so that the biomass is then updated according to: rpave = (1.0-Theta)*rprime0 + Theta*rprime biom = biomass0*exp((growth_yield*rpave-decay_con)*Deltat) We were wondering about the numerical justification for this. It seems like it assumes that the chemical concentrations change slower over time than the biomass does, because the concentration appears in "rpave" (actually in rprime). Basically this integral seems to have ignored fact that the chemical concentration is a function of time, which would be fine if the chemical concentration were varying in time much slower than the biomass, but we were just wondering about the exact numerical justification for this. Does anyone know about the numerical justification for this integral? Also we are a bit confused by the line rpave = (1.0-Theta)*rprime0 + Theta*rprime and what the variable Theta actually does. So if someone could shed some light on this question as well then that would be great. Thanks very much. Regards, Tim Bush