Environmental Engineering Reference
In-Depth Information
The output parameters can be adopted from one of the other transport M-files,
'simpletrans.m', 'analtrans.m'
or
'pdepetrans.m'
. The discretization parameters are
also taken analogously as demonstrated above. The execution part consists mainly
of the call of the
pdepe
-module. The function names and parameters are of course
new for the slow sorption application:
The solution is again written in matrix form. As there are two unknown
concentrations at each position and for each output time level, there are three
indices: one speaks of a tensor of rank 3. Three indices are necessary to mark a
single element of
c
. The third index must be 1 or 2, depending on whether to denote
the concentration of the dissolved or the adsorbed species. In all output commands,
the third index needs to be specified in order to determine which concentration has
to be plotted. Where in former codes the variable c was sufficient, now
c(:, :, 1)
has to be inserted. With this command only the fluid phase concentration is plotted.
In order to obtain the solid phase concentration use
c(:, :, 2)
.
