You're reading the documentation for a development version. For the latest released version, please have a look at v4.4.0.

Two Dimensional General rate model (GRM2D)

The general rate model as introduced in Section General rate model (GRM) assumes homogeneity in the cross sections of the column. This allows to consider transport along the axial dimension only. However, due to packing irregularity and inhomogeneous flow at the inlet (i.e., frits), this assumption may be a crude approximation. This model can be improved by introducing a radial coordinate \(\rho \in [0, R]\), where \(R\) is the column radius, in the interstitial volume Eq. 1:

(16)\[\begin{split} \varepsilon_c \frac{\partial c^l_i}{\partial t} = &-\varepsilon_c u \frac{\partial c^l_i}{\partial z} + \varepsilon_c D_{\text{ax},i} \frac{\partial^2 c^l_i}{\partial z^2} + \frac{1}{\rho} \frac{\partial}{\partial \rho} \left( \rho D_{\text{rad},i} \frac{\partial}{\partial \rho} \left( \varepsilon_c c^l_i \right) \right) \\ &- \left(1 - \varepsilon_c\right) \sum_j d_j \frac{ 3 k_{f,j,i} }{r_{p,j}} \left[ c^l_i - c^p_{j,i}(\cdot, \cdot, \cdot, r_{p,j}) \right] + \varepsilon_c f_{\text{react},i}^l\left(c^l\right).\end{split}\]


  • \(c^l_i\colon \left[0, T_{\text{end}}\right] \times [0, L] \times [0, R] \rightarrow \mathbb{R}^{\geq 0}\),

  • \(c^p_{j,i}\colon \left[0, T_{\text{end}}\right] \times [0, L] \times [0, R] \times [r_{c,j}, r_{p,j}] \rightarrow \mathbb{R}^{\geq 0}\), and

  • \(c^s_{j,i,m_{j,i}}\colon \left[0, T_{\text{end}}\right] \times [0, L] \times [0, R] \times [r_{c,j}, r_{p,j}] \rightarrow \mathbb{R}^{\geq 0}\)

depend on \(\rho\). Additionally, the porosity \(\varepsilon_c\), axial dispersion coefficient \(D_{\text{ax},i}\), radial dispersion coefficient \(D_{\text{rad},i}\), and interstitial velocity \(u\) may depend on \(\rho\).

The dependence of the parameters on \(\rho\) is not arbitrary. For simplicity, it is assumed that the parameters are piecewise constant, that is, the range \([0, R]\) is divided into disjoint zones in which all parameters are constant. These zones are used for radial discretization and can be supplied to the simulator. Continuous dependence of the parameters can be realized by piecewise constant approximation.

The Danckwerts boundary conditions at the column in- and outlet, Eq. 12 and 11, are modified to account for the radial coordinate:

(17)\[\begin{aligned} u(\rho) c_{\text{in},i}(t,\rho) &= u(\rho) c^l_i(t,0,\rho) - D_{\text{ax},i}(\rho) \frac{\partial c^l_i}{\partial z}(t, 0, \rho) & \forall t > 0, \rho \in (0,R), \end{aligned}\]
(18)\[\begin{aligned} \frac{\partial c^l_i}{\partial z}(t, L, \rho) &= 0 & \forall t > 0, \rho \in (0,R). \end{aligned}\]

Conditions for the radial direction are added:

(19)\[\begin{split}\begin{aligned} \frac{\partial{c^l_i}}{\partial \rho}(\cdot, \cdot, 0) &= 0, \\ \end{aligned}\end{split}\]
(20)\[ \begin{aligned} \frac{\partial{c^l_i}}{\partial \rho}(\cdot, \cdot, R) &= 0. \end{aligned}\]

While the inner condition Eq.19 represents symmetry at the column center, the outer condition Eq. 20 is a no-flux condition.

Using the inlet boundary condition Eq. 17, each radial zone is equipped with its own inlet and outlet port. That is, this unit operation has as many inlet and outlet ports as it has radial zones (parameter NRAD in the discretization group). This allows each radial zone to have its own inlet profile, which enables modeling of flow distribution in the frits by sending the feed through varying hold-up volumes before injecting it into a radial zone.

Specification of flow rate / velocity and direction

Since the column radius \(R\) and the zones \((\rho_k, \rho_{k+1})\), \(k = 0, \dots, N_{\text{rad}} - 1\), are known, the interstitial velocities \(u_k\) are inferred from the volumetric flow rates via

\[\begin{aligned} u_k = u_{\text{int},k} = \frac{F_{\text{in},k}}{\pi \left( \rho_{k+1}^2 - \rho_k^2 \right) \varepsilon_{c,k}}, \end{aligned}\]

where \(F_{\text{in},k}\) denotes the volumetric flow rate into zone \(k\).

The direction of flow inside the radial zone of the unit operation is governed by the sign of the interstitial velocity \(u_k\). A positive sign results in (standard) forward flow, whereas a negative sign reverses the flow direction. Note that in case of reversed flow, the chromatogram is returned at the unit operation’s inlet port, which may not be returned from simulation by default.

Note that, contrary to the standard general rate model as presented in Section General rate model (GRM), the interstitial flow rate is always given by the volumetric flow rate. The velocity parameter only determines the flow direction.

For information on model parameters see Two dimensional general rate model.