129: compute the mean transit time

Description Arguments Details Author(s) References


This method computes the mean transit time for the linear time invariant system that can be constructed from the given operator and input distribution.

It relies on the mehtod getTransitTimeDistributionDensity using the same arguments.




To compute the mean transit time for the distribution we have to compute the integral

\bar{T} = \int_0^∞ T \cdot S_r≤ft( \frac{\vec{I}}{I},0,T\right)\; dT

for the numerically computed density. To avoid issues with numerical integration we dont use as upper limit but cut off the integragion interval prematurely. For this purpose we calculate a maximum response time of the system as Lasaga

τ_{cycle} = \frac{1}{|\min(λ_i)|}

where λ_i are non-zero eigenvalues of the matrix {\bf A}.


Carlos A. Sierra, Markus Mueller


Lasaga, A.: The kinetic treatment of geochemical cycles, Geochimica et Cosmochimica Acta, 44, 815 – 828, doi10.1016/0016-7037(80)90263-X, 1980.

Search within the SoilR package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.