Multi-compartment model

From Wikipedia, the free encyclopedia

A multi-compartment model is a type of mathematical model used to describe the way materials or energies are transmitted among the compartments of a system. Each compartment is assumed to be a homogenous entity within which the entities being modelled are equivalent. For instance, in a pharmacokinetic model, the compartments may represent different sections of a body within which the concentration of a drug is assumed to be constant.

Hence a multi-compartment model is a lumped parameters model.

Multi-compartment models are used in many fields including pharmacokinetics, epidemiology, biomedicine, systems theory, complexity theory, engineering, physics, information science and social science. The circuits systems can be viewed as a multi-compartment model as well.

In systems theory, it involves the description of a network whose components are compartments that represent a population of elements that are equivalent with respect to the manner in which they process input signals to the compartment.

1 Assumptions
2 Single compartment model
3 Multi-compartment model
4 Model topologies
5 See also

[edit] Assumptions

Multi-compartment model assumes adopts several assumptions such that systems in physical existence can be modeled mathematically:

Instant homogeneous distribution of materials or energies within a "compartment".
The exchange rate of materials or energies among the compartments is related to the densities of these compartments.
Usually, it is desirable that the materials do not undergo chemical reactions while transmitting among the compartments.
When concentration of the cell is of interest, typically the volume is assumed to be constant over time, though this may not be totally true in reality.

Most of the time, multi-compartment model simulates only one substance, eg. the concentration of solutions, to simplify the mathematics.

[edit] Single compartment model

Possibly the simplest application of multi-compartment model is in the single cell concentration monitoring (see the figure above). If the volume of a cell is $V,$ the mass of solute is $q,$ the input is $u (t)$ and the secretion of the solution is proportional to the density of it within the cell, then the concentration of the solution $C$ within the cell over time is given by:

$\frac{\mathrm{d}q}{\mathrm{d}t}=u(t)-kq$

$C=\frac{q}{V}$

where $k$ is the proportionality.

[edit] Multi-compartment model

As the number of compartments increases, the model can be very complex and the solutions usually beyond ordinary calculation. Below shows a three-cell model with interlinks among each other.

The formula for n-cell multi-compartment models become:

$\begin{align} \dot{q}_1=q_1 k_{11}+q_2 k_{12}+\cdots+q_n k_{1n}+u_1(t) \\ \dot{q}_2=q_1 k_{21}+q_2 k_{22}+\cdots+q_n k_{2n}+u_2(t) \\ \cdots\\ \dot{q}_n=q_1 k_{n1}+q_2 k_{n2}+\cdots+q_n k_{nn}+u_n(t) \end{align}$

where

$k_{ii}=\sum^n_{j=0;i\neq j}{k_{ji}}$ for $i=1,2,\dots,n.$

Or in matrix forms:

$\mathbf{\dot{q}}=\mathbf{Kq}+\mathbf{u}$

where

$\mathbf{K}=\begin{bmatrix} k_{11}& k_{12} &\cdots &k_{1n}\\ k_{21}& k_{22} & \cdots&k_{2n}\\ \vdots&\vdots&\ddots&\vdots \\ k_{n1}& k_{n2} &\cdots &k_{nn}\\ \end{bmatrix}$ and

$\mathbf{q}=\begin{bmatrix} q_1&q_2&\cdots&q_n\end{bmatrix}^T$ and $\mathbf{u}=\begin{bmatrix} u_1(t)&u_2(t)&\cdots&u_n(t)\end{bmatrix}^T$

[edit] Model topologies

Generally speaking, as the number of compartment increases, it is challenging both to find the algebraic and numerical solutions of the model. However, there are speacial cases of models, which rally exist in nature, when the topologies exihit certain regularities that the solutions become easier to find. The model can be classified according to the interconnection of cells and input/output characteristics:

Closed model: No sinks or source, lit. all k_oi=0 and u_i=0;
Open model: There're sinks or/and sources among cells.
Catenary model: All compartments arranges in a chain, with each pool connecting only to its neighbors. This model has 2 or more cells.
Cyclic model: It's a special case of the catenary model, with 3 or more cells, in which the first and last cell are connected, i.e. $k_{1n}\neq 0$ or/and $k_{n1}\neq 0$ .
Mammillary model: Consists of a central compartment with peripheral compartments connecting to it. There're no interconnections among other compartments.
Reducible model: It's a set of un-connected models. It bears great resemblance to the computer concept of forest as against trees.