WOLFRAM|DEMONSTRATIONS PROJECT

Controlled Release of a Drug from a Hemispherical Matrix

​
τ
2
α
0.1
This Demonstration models the release of a drug from a hemispherical matrix.
Drug transfer occurs only though the hemispherical cavity at the bottom of the matrix. The governing dimensionless equation for this system is:
∂θ
∂τ
=
1
2
ξ
∂
∂ξ
2
ξ
∂θ
∂ξ
, with the initial and boundary conditions:
θ(ξ,0)=1
,
θ(α,τ)=0
and
∂θ
∂r
(1,τ)=0
.
The dimensionless variables are:
θ=
c
c
o
,
α=
r
i
r
o
,
ξ=
r
r
o
and
τ=
t
2
r
o
,
where
c
is the drug concentration,

is the drug diffusivity,
r
is the radial coordinate,
t
is time,
c
0
is the initial drug concentration, and
r
i
and
r
o
stand for the inner and outer radii of the matrix, respectively.
Siegel [1] solved this problem using separation of variables. The expression for the dimensionless drug concentration is:
θ(ξ,τ)=
2α
ξ
∞
∑
n=1
1
λ
n
2
λ
n
+1
2
λ
n
(1-α)-α
exp-
2
λ
n
τsin[
λ
n
(ξ-α)]
,
and the fractional drug release as a function of dimensionless time is given by
Φ(τ)=1-
6
2
α
1-
3
α
1
2
λ
n
2
λ
n
+1
2
λ
n
(1-α)-α
exp-
2
λ
n
τ
,
where
λ
n
is the
th
n
positive root of
λ-tan(λ(1-α))
.
In this system, the drug is released at a near-constant rate. This is a desired goal of all controlled-release drug-delivery mechanisms, because this pattern can guarantee maintenance of plasma drug concentrations within therapeutic levels. Diffusion-controlled monolithic systems such as slabs, cylinders and spheres usually do not yield constant release rates.