Concentration Waves: Chromatographic Theory and Experimental Verification


Joachim Gruber
E-Mail: Joachim_Gruber
Address at the time this paper was finished:
Arbeitsbereich Umweltschutztechnik, Technische Universität,
Eissendorfer Str. 40, D 21073 Hamburg, Germany


ABSTRACT
Aqueous transport of negative ions (OAc-, Cl-, OH-) in a column filled with a positively charged ion exchange resin is described with a quasilinear set of hyperbolic conservation laws (mass balance equations). Analytical solutions of the Riemann problem (a single abrupt change of the ion concentrations at the column input) are found and discussed for the exact and asymptotic systems. An experiment performed 1984 by Barbara J. Bennett is interpreted with the presented formalism.


Introduction

There are two major research areas in the field of ion transport in aqueous porous media:

This paper will address the second area, assuming thermodynamic equilibrium in the flow channels. The mass action laws describing the chemical reactions are generally non-linear.

The chemical composition of the liquid or solid phases (speciation) in spatially homogenous (so called batch) systems change qualitatively when concentrations move through certain boundaries (thresholds) given by the mass action law constants describing the chemical interactions [see e.g. page 65, Morel and Hering, 1993]. The reason for this is the non-linearity introduced by the mass action laws. It is virtually impossible to predict the chemical composition on the other side of a concentration threshold when experimental or numerical results are available only on this side of the threshold and thus even the existence of the threshold is unknown. It's like walking around the edge of a building: You simply cannot know ahead of time what you will be facing behind the building.

After non-linear geochemistry has been added to transport models, the same became true for the resulting models. The Conclusions section will summarize the manifestations of non-linear processes in the simple exchange system presented here. This will exemplify, that by giving us access to the governing equations, i.e. the speciation laws and the hyperbolic conservation laws, respectively.

Finally, I will show to what extent the presented chromatographic solutions of the transport problem coincide with the observed concentration waves.

General Properties of Multicomponent Chromatography

The following theorems are basic to the theory of multicomponent chromatography and will be stated without proof (see Gruber, 1990, 1995, 1996 for a general introduction into the solutions of hyperbolic conservation laws using a notation consistent with the one in this paper).

Definition:
Riemann problem and centered rarefaction waves

The set of conservation laws will be solved for a single abrupt change (Riemann jump) of the chemical composition of the feed water (Riemann problem):

c(x ≤ 0, t=0) = c- (feed water composition),
c(x>0, t=0) = c+ (pre-equilibrant water composition).

The solution considered here is the centered k-rarefaction wave, because from a practical point of view this class of solutions will comprise also the shocks, the other relevant class of solutions. The sequence of k-waves c(x = l, t) is called breakthrough curve.

In a centered k-rarefaction wave

Theorem:
Riemann invariants

Properties and physical meaning of the Riemann invariant:

Theorem:
Grid of shocks and rarefaction-waves

When the separation factors are independent from the concentrations (as in (4)) and complex formation in water is negligible,

2-Component System: Na - OAc - OH

In the first example an aqueous solution of only sodium acetate migrates through the ion exchange column (for meaning of the symbols that follow see Notation Section in the back):

Definition 1:
Number of species and constraints


Definition 2:
Components and conservation laws



Theorem 1:
Waves in the Na - OAc - OH system



Solution of Riemann problem for Na-OAc-OH system
FIGURE 1: Solution of the Riemann problem for Na-OH-OAc system. Lines ("characteristics") are location of arbitrarily selected concentrations c(ξ) within the waves. Waves are bounded by heavy lines, the concentrations of the constant states. Constant state vectors are c- = {[OAc-]-, [OH-]-}, c = {[OAc-], [OH-]}, c+ = {[OAc-]+, [OH-]+}. ξ = x/t is the slope of the characteristic for c(ξ). Constant state concentrations in fan of waves are given for the Riemann invariants, [Na+] = [OAc-] + [OH-] - [H+] and [OAc]. 2-wave is non-retarded.

Waves in concentration space for Na-OAc-OH system
FIGURE 2. Waves in concentration space. 1-waves are the ones that are parallel and have a slope -1. 2-waves form the fan originating from the coordinate origin. They are (straight) lines because of 1/1 exchange process. [Na+]- = 0.02, 0.04, ..., 0.2 mol/L. Abscissa: concentration of OAc in solution, ordinate: concentration of OH in solution.


The set (9) and (11) of conservation laws has two spatially continuous solutions, the 1- and the 2-rarefaction wave. Fig. 1 gives a schematic of the waves and Riemann invariants in physical {x, t} space. Fig. 2 shows some 1- and 2-waves in concentration space.

The 1-rarefaction wave ("ion exchange wave")

The 2-wave ("sodium wave")


Proof

Notation for proof:

Let

superscripts + and - on the Riemann invariants and

subscripts + and - on the concentrations

denote the constant states + and - upstream and downstream of the 1- and 2-wave, respectively, i.e. let the superscripts on wj specify the state at which the Riemann invariant is evaluated:

w1- = [OAc]-, w1+ = [OAc]+,.................... (12)

or w1- = [OH]-, w1+ = [OH]+, .................... (12')

and w2- = [Na+]-, w2+ = [Na+]+..................... (13)

Similarly, let

α+OAcOH = [OAc]+[OH-]+/([OAc-]+[OH]+) .................... (14)

be the separation factor evaluated at c+.

Let a concentration without index be the concentration

We choose as Riemann invariants

w1 = [OAc], .................... (15)

or w1 = [OH], .................... (15')

and w2 = [Na+]..................... (16)

This choice is based on the following argument: According to (10), one solution, the 1-wave, has the properties

ξ ≠ vw/φ and [Na+]' = 0.................... (17)

i.e. the wave is retarded, and propagates the surface composition jump, as [OAc] and [OH] change ([OAc]' ≠ 0 ≠ [OH]'). The other solution of (10) is

ξ = vw/φ and [OAc]' = [OH]' = 0. ....................(18)

Fig. 1 gives the values of the Riemann invariants in the fan of waves. From this pattern of Riemann invariants we will derive the equations of the waves.


1-Wave

The invariance (12) of w2 across the 1-wave together with the electroneutrality in water, (5), tells us

[Na+]-= [OAc-] + [OH-] - [H+].................... (19)

We can use the H+/OH- equilibrium, (3), to eliminate [H+] from (19), solve for [OH-] and get the representations of the 1-wave:

[OH-] = a + √[a2 + Kw] ..........(1-wave).................... (20)

where a = 1/2 ([Na+]- - [OAc-]) and √ stands for the square root. At pH > 7 the 1-wave (19) is approximately a line with slope -1:

[OH-] = [Na+]- - [OAc-] (for [OH-] >> [H+])..................... (20')

Note that this slope reflects the 1/1 OAc-/OH- exchange process.


2-Wave

The other solution of (10) is

w2' = [Na+]' ≠ 0 and ξ = vw/φ ..................... (21)

The 1-Riemann invariant is constant:

w1 = w1+

or using the expression for w1

[OAc] = [OAc]+ and [OH] = [OH]+.................... (22)

Then due to the invariance of w1,

α+OAcOH = [OAc][OH-]+/([OAc-]+[OH-])..................... (14')

By assumption the separation factor is independent from the concentrations:

α+OAcOH = αOAcOH = 2.96.................... (23)

Solving this relationship for [OH-] gives us the equation of the 2-wave

[OH-] = [OH-]+ [OAc-]/[OAc-]+..........(2-wave)..................... (24)


Note that some other than the 1/1 exchange stoichiometry would show up in a separation factor involving higher powers of the concentration ratios. This non-linearity would then be passed on to the shape of the non-retarded wave.

Corollary

The neutral complex HOAc is without influence on the representation of the waves in concentration space (it does affect the retardation of the 1-wave).

Proof: [HOAc] did not enter into the deductions of the 1- and 2-wave equations.

3-Component System: Na - OAc - Cl - OH

When Cl- ions are added to the water of the previous system, we obtain 2 more variables, [Cl-] and [Cl], and one more constraint (25). Therefore the system has 3 components, now:


Definition 4:
Number of species and constraints


Definition 5:
Components and conservation laws


Theorem 3:
Waves

The addition of chloride ions introduces one more retarded wave. What has been said in Theorem 1 about the retarded wave can be generalized here to the two retarded waves. The 3-wave here is the analog of the 2-wave characterized in Theorem 1:


Fan of waves in 3-component system Na-OAc-Cl-OH
FIGURE 3 : Fan of waves in 3-component system Na-OAc-Cl-OH. c-= {[OAc-]-, [Cl-]-, [Na+]-}, c+= {[OAc-]+, [Cl-]+, [Na+]+}. Across the non-retarded (3-)wave concentration of sodium changes from [Na+]+ to [Na+]-. Concentration vectors on the 3- or 2-wave will be referred to as c = {[OAc-], [Cl-], [Na+]}.


Waves in concentration space of Na-OAc-Cl-OH system
FIGURE 4. Some 3-waves (eq. (33), heavy lines) and surface (36) for [Na+]- = 0.1 mol/L with 1- and 2-waves. Asymptotic equations (42), (44) and (51), (52) used for 1-, 2-waves, respectively. 2-waves is the starlike fan of lines.


Projection of 1- and 2-waves of Fig. 4 onto OAc-Cl plane
FIGURE 5. Same as Fig. 4 (with 3-waves omitted), view is perpendicular onto {[OAc-], [Cl-]} plane. 1-waves are more densely spaced lines with negative slope. Their upward bending is a consequence of the HOAc complex. Waves exist on the right side of the line intersecting the axes at 0.1 mol/L because of the water dissociation, the H+/OH- equilibrium.

The 1-and 2-rarefaction waves ("ion exchange waves")

The 3-wave ("sodium wave")
Click here to jump over proof of this theorem.
Proof

The sum of the set of conservation laws (9), (11), (28) is

(vw/φ - ξ)([OAc-] + [Cl-] + [OH-] - [H+])' - ξ/φ ([OAc] + [Cl] + [OH])'= 0
........(29)

Electroneutrality on the surface (27) requires that the sum of all adsorbed concentrations equaling the ion exchange capacity XT be constant. Thus ([OAc] + [Cl] + [OH])' = 0. Therefore (29) is simply

(vw/φ - ξ)([OAc-] + [Cl-] + [OH-] - [H+])' = 0 .................... (29')

which is the conservation law for Na+(see electroneutrality condition (26)). We will choose

w3 = [Na+]..................... (30)

Na+ is a component we could have chosen instead of one of the others, OAc and Cl.


3-wave

The proof of point 2 of the theorem on the 3-wave is analogous to the corresponding proof in the chloride free system. Let c = {[OAc-], [Cl-], [OH-]} lie on the 3-wave (see Fig. 3). Across a non-retarded wave, ξ = vw/φ, the concentrations on the surface cannot change (see conservation laws (9, 11, 28):

[OAc]' = 0, or [OAc] = [OAc]+.................... (31)
[Cl]' = 0, or [Cl] = [Cl]+
[OH]' = 0, or [OH] = [OH]+

The equation for the 3-wave in three-dimensional concentration space c = {[OAc-], [Cl-], [Na+]} is given by the invariance of the separation factors

α+OAcOH = αOAcOH.................... (32)
α+ClOH = αClOH

where the separation factors are evaluated at the compositions c+ and c, downstream and upstream of the 3-wave, respectively. Using the invariance relations (30) in (31) gives us the equations for the 3-wave:

[OH-] = [OH-]+/[OAc-]+ [OAc-].................... (33)

[Cl-] = [Cl-]+/[OAc-]+ [OAc-].

Note that like in the 2 component system separation factors that are non-linear in the ratios of the concentrations would pass this non-linearity on to the shape of the waves.


1- and 2-rarefaction wave

The 1- and 2-waves have to meet the condition (29') for non redarded speeds ξ: vw/φ ≠ ξ:

[Na+] = [OAc-] + [Cl-] + [OH-] - [H+] = [Na+]- ,.................... (34)

Again, like in the OAc - OH system above, the water dissociation equilibrium (3) is used to eliminate [H+],

[OAc-] + [OH-] + [Cl-] - Kw/[OH-] = [Na+]-, ....................(35)

and the resulting expression is solved for [OH-],

[OH-] = a + √[a2 + Kw].................... (36)

with

a = ([Na+]- - [OAc-] - [Cl-])/2

Fig. 4 shows the surface (36) in concentration space. It can be envisaged as a surface attached to and having a bend near the line [OAc-] + [Cl-] = [Na+]- in the {[OAc-],[ Cl-]} plane.

Since the exact Riemann invariants for the subsystem with constant [Na+] cannot be guessed from physico-chemical principles, I will calculate the 1- and 2-waves in two ways:



Eigenvectors of the Jacobian

The adsorption isotherms Ci = [OAc], [Cl] and [OH] can be constructed combining the definitions (4), (25) of the separation factors and of the ion exchange capacity (27):

Ci([OAc-], [Cl-], [OH-]) = XT ai OH ci / d.................... (37)

with

d = αOAcOH[OAc-] + αClOH[Cl-] + [OH-]

where cj = [OAc-], [Cl-] or [OH-].

The elements of the Jacobian matrix C'([OAc-], [Cl-], [OH-]) of the isotherms are the partial derivatives of the isotherms with respect to the concentrations [OAc-], [Cl-] and [OH-] [Gruber, 1990]. The k-waves ck are everywhere tangential to the k-eigenvectors rk of this matrix C'. I have constructed the k-waves from the k-eigenvectors with a suitable Euler or Runge-Kutta method [Gruber, 1995].

As long as the waves are (straight) lines (which is predominantly the case in the systems discussed here), the fastest and most convenient integration procedure is the Euler method. The k-wave is composed of segments

ck(n) = ck(n-1) + Delta rk(ck(n-1)).................... (38)

where n is the segment number and Delta is the length of a segment pointing in the direction of the k-eigenvector rk(ck(n-1)) evaluated at the beginning of the wave segment ck(n-1).


Asymptotic Solution 1: Low pH Approximation of 1- and 2-Waves:
(Negligible Adsorption of OH-)


If the adsorbed OH- ions can be neglected, i.e. if

[OAc] + [Cl] >> [OH].................... (39)

which is equivalent to

αOAcOH[OAc-] + αClOH[Cl-] >> [OH-],.................... (39')

the conservation law for [OH-] can be dropped from the set (9), (11), (28). The remaining set describes exchange between acetate and chloride, which can be solved similarly as the exchange between acetate and hydroxyl in the Na-OAc-OH system:

(vw/φ - ξ) ([OAc-] {1 + β [H+]})' - ξ/φ [OAc]'= 0 .................... (9)

(vw/φ - ξ) [Cl-]' - ξ/φ [Cl]'= 0 .................... (28)

The sum of these conservation laws (abbreviating [Ac] = [OAc] {1 + β [H+]})

(vw/φ - ξ) ([Ac] + [Cl-])' - ξ/φ ([OAc] + [Cl])'= 0

can be solved, because -for negligible [OH]- OAc and Cl cover the entire ion exchange resin surface ([OAc] + [Cl] = XT):

(vw/φ - ξ) ([Ac] + [Cl-])' = 0,.................... (40)

For ξ ≠ vw/φ, the 1-wave,

w1 = ([Ac] + [Cl-])' = 0.................... (41)

which can be integrated -as above:

[Ac] = [Ac]- - [Cl-] + [Cl-]-..................... (42)

Plugging electroneutrality in water, (34) with constant [Na+] = [Na+]- across the 1- and 2- wave, (see Fig. 3),

[OAc-] + [Cl-] - [H+] = [Na+]-,

(the latter equation includes the low pH approximation (39)).

into (42) and solving for [Cl-] we arrive at the low pH approximation for the 1-wave

[Cl] = (β[Na+]-/b - 1)[OAc-] + (b*-[OAc]- + [Cl]-) / b.................... (43)

with

b = 1 + β[OAc-], b-* = 1 + β Kw/[OH-]-

[OH-]- = a- + √[a2 + Kw].................... (36)

a- = ([Na+]- - [OAc-]- - [Cl-]-)/2


[OH-] = f([OAc-], [Cl-]) is the surface (36) with the 1- and 2-waves.

The fast (2-)wave in the system (9) and (28) is non-retarded (it is non-retarded only in the framwork of the low pH approximation). It is a line through the point {[OAc-]+, [Cl-]+} and the origin of concentration space {[OAc-], [Cl-]} = {0, 0}.

[OAc] = [OAc-]+/[Cl-]+ [Cl-].................... (44)

The proof is analogous to the one given for the 2-component system.

(43) and (44) are the waves in the low pH region (right part of Fig. 5).


Asymptotic Solution 2: High pH Approximation
(Negligible HOAc complex)


In the high pH approximation the presence of protons i.e. H+ and HOAc, is neglected.

[OH-] >> [H+] - [HOAc].................... (45)
[OAc-] >> [HOAc]..................... (46)

Thus, the set of conservation laws (9), (11) and (28) assumes the following simple form

(vw/φ - ξ)[OAc-]'- ξ/φ [OAc]' = 0.................... (9')

(vw/φ - ξ)[OH-]' - ξ/φ [OH]' = 0 ....................(11')

(vw/φ - ξ)[Cl-]' - ξ/φ [Cl]' = 0 .................... (28)


We get a Langmuir type adsorption isotherm due to

of the form

Ci([OAc-], [Cl-], [OH-]) = XT ai OH ci / dl.................... (47)

with

dl = 1 + (αOAcOH - 1)[OAc-] + (αClOH - 1)[Cl-],

ci = ci/[Na+]-.

The analytical solution of this type of equations (9'), (11'), (28) and (35') has been published in the literature [Rhee, Aris, Amundson, 1989] and been applied to the present set of equations by Gruber [1996].

The solution of this 2-component system are waves that lie on lines tangential to a parabola (see Fig. 6)

Cllp = k - [OAcp]/k + 2 √[-[OAcp]/k] (lower branch of parabola).................... (48)

Clup = k -[OAcp]/k - 2 √[-[OAcp]/k] (upper branch of parabola).................... (49)

where

k = - (1 - aOAcCl)/(1 - αOAcOH) [Na+]-..................... (50)

The equations of the 1-waves are

[Cl-] = Clup + dClup/d[OAcp] ([OAc-] - [OAcp]) .................... (51)

and of the 2-waves

[Cl-] = Cllp + dCllp/d[OAcp] ([OAc-] - [OAcp]).................... (52)

The waves (tangents to the parabola) are shown in Fig. 6.

Physical interpretation of the grid of 1- and 2-waves: Buffers

For the following discussion, note that the parabola touches the concentration axes at [OAc-] = - k and [OH-] = k .


Parabola and 1- and 2-waves in Na-OAc-Cl-OH system

FIGURE 6. Parabola (48), (49) and calculated asymptotic 1- and 2-waves (51), (52), in a 2-component system with conservation laws (9'), (11'), (28) and constant normality [Na+]- = 0.1 mol/L, (35'). These high pH asymptotic waves tangential to the parabola and the exact waves (i.e. the waves (38), derived from the Jacobian) coincide in the triangle bounded by the line [OAc-] + [Cl-] = [Na+]- = 0.1 mol/L. 1-waves are the lines with negative slope. k = 0.306 [Na+]- . aOAcCl = 0.4, αOAcOH = 2.95.

Comparison with Experiment

In 1984, Barbara J. Bennett performed experiments with the 2- and 3-component systems discussed here and interpreted the results on the basis of Helfferich's version (the original version) of multicomponent chromatography [Helfferich and Bennett, 1984a, b]. She used the Rohm & Haas Amberlite IRA-400 strong-base anion exchange resin. Chemical equilibrium between adsorber resin and the liquid phase was assumed.

Fig. 7 shows a portion of concentration space with the calculated grid of 1- and 2-waves and the experimentally determined 1-, 2- and 3-waves ("observed composition path"). Since measurement errors were not reported, I am not able to decide whether the observed 3-wave actually stops short of the surface [Na+]- in which the 1- and 2-waves run. If one allows an uncertainty of 0.2 pH units (or less in unison with uncertainties of [OAc-] and [Cl-]), the 3-wave runs into the required surface.

Fig. 8 allows a comparison of the observed composition path and the calculated grid of the 1- and 2-waves. Helfferich and Bennett comment that


Grid of 1- and 2-waves together with experimental composition path in Na-OAc-OH system

same as Fig. 7, different viewpoint

FIGURE 7. Calculated grid of 1- and 2-rarefaction waves (see Fig. 4) and experimentally observed waves (path through concentration space represented by heavy curve). Note that the calculated 3-wave (connecting the origin of concentration space {[OAc-], [Cl-], [OH-]} = {0, 0, 0} with the surface of the 1- and 2-wave grid) has been omiited to avoid clutter. Upper and lower part of the figure differ only in the point from which concentration space is viewed.


Grid of 1- and 2-waves together with experimental composition path projected onto OAc-Cl plane in Na-OAc-OH system

FIGURE 8. Experimentally determined 1- and 2-wave (observed composition path, heavy line) entered into Fig. 6. Main part of deviation between observation and multicomponent chromatography as presented in this paper is attributed to concentration dependent separation factors not considered in this analysis [Helfferich and Bennett, 1984a and b].

Conclusions

The grid of k-waves in concentration space is calculated The grid topology reflects the non-linearity introduced by the chemical interactions and reveals the relevant processes, i.e.

When the breakthrough curves c(x = l, t) of the 2- and 3-component ion exchange system are drawn in concentration space, the resulting composition path aligns to a degree with the grid of k-waves that seems satisfactory for environmental assessments.

Acknowledgement

I wish to thank Friedrich G. Helfferich of The Pennsylvania State University, University Park, PA, USA for having encouraged this work and helped me find these experimental data. The data have been generated by Barbara J. Bennett in the framework of her thesis for the Degree of Master of Science in the Department of Chemical Engineering at Penn State University.

Notation

[Ac].......... concentration of acetate (= [OAc-] + [HOAc]) in water (mol/L of water),
αClOH.......... constant (i.e. concentration independent) separation factor (25) (dimensionless),
αOAcOH.......... constant (i.e. concentration independent) separation factor (4) (dimensionless),
β.......... HOAc complex formation constant (= [HOAc]/[OAc-]) (dimensionless),
[Cl-].......... concentration of chloride in water (mol/L of water),
[Cl].......... concentration of chloride on adsorbing surfaces (mol/L of system ("IX bed")),
[Clp].......... Cl- concentration on parabola (mol/L of water),
c.......... concentration vector in water (= {[OAc-], [Cl-], [OH-]}) (mol/L of water),
cj ..........normalized concentration of component j (= cj/[Na+]- ) (mol/L of water),
concentration space ..........{[OAc-], [Cl-], [OH-]} space,
d.......... adsorption isotherm denominator (= αOAcOH[OAc-] + αClOH[Cl-] + [OH-], (37)) (mol/L of water),
dl.......... Langmuir adsorption isotherm denominator (= 1 + αOAcOH[OAc-] + αClOH[Cl-], (47)) (dimensionless),
φ (or f).......... fraction of system volume filled with water (0.4),
[H+].......... concentration of free protons in water (mol/L of water),
[HOAc].......... concentration of HOAc complex in water (mol/L of water),
Kw.......... dissociation constant of water (mol2/L2 of water),
k.......... size parameter of parabola (= - (1 - aOAcCl)/(1 - αOAcOH) [Na+]-) (mol/L of water),
k.......... index of k-rarefaction wave. An n0 component system has n0 k-rarefaction waves (k = 1, ... n0),
l.......... length of the ion exchange column (m),
l.......... superscript assigned to quantities referring to lower part of parabola,
[Na+].......... concentration of sodium in water (mol/L of water),
[Na+]-.......... concentration of sodium in column feed (mol/L of water),
[Na+]+.......... concentration of sodium in column at t = 0 (pre-equilibrant) (mol/L of water),
n0.......... number of chemical components and waves in system, it is also the dimension of concentration space,
[OAc-].......... concentration of OAc- in water (mol/L of water),
[OAcp].......... OAc- concentration on parabola (mol/L of water),
[OAc].......... concentration of OAc- on ion exchange resin (mol/L of system ("IX bed")),
[OH-].......... concentration of OH- in water (mol/L of water),
[OH].......... concentration of OH- on ion exchange resin (mol/L of system ("IX bed")),
p.......... index assigned to point cp = {[OAcp-], [Clp-] on parabola,
√[].......... square root of the expression in [],
t.......... time variable (s),
u.......... superscript assigned to quantities referring to upper part of parabola,
vw.......... flux of water (m3 of water per m2 of flow channel cross section. Note: channels have tiny cross sections),
vw/φ .......... speed of the water in the porous medium (m/s)
wk.......... k-Riemann invariant (the invariant that changes across the k-rarefaction wave),
XT.......... concentration of adsorbing sites on ion exchange resin (mol/L of system ("IX bed")),
x.......... spatial variable (m),
ξ, ξ(c).......... speed of concentration vector c (m/s),
-.......... subscript assigned to concentration at t = 0, x <= 0,
+.......... subscript assigned to concentration at t = 0, x > 0,
'.......... derivative with respect to ξ, meaning the change across the wave.

References

Gruber, J., Contaminant accumulation during transport through porous media, Water Resour. Res., 26, 99 - 107, 1990. (an enhanced version of my paper.)

Gruber, J., Waves in a two-component system: The oxide surface as a variable charge adsorbent, Ind. Eng. Chem. Res., 34, 2769-2781, 1995. (a modified version of my paper.)

Gruber, J. Advective transport of interacting solutes: the chromatographic model, in: U. Förstner and W. Calmano (eds.), "Sediments and Toxic Substances", Environmental Sciences Series, Springer, Heidelberg, 1996. (a modified version of my paper.)

Helfferich, F.G. and B.J. Bennett, Weak electrolytes, polybasic acids, and buffers in anion exchange columns: I. Sodium acetate and sodium carbonate systems, Reactive Polymers, 3, 51 - 66, 1984 a.

Helfferich, F.G. and B.J. Bennett, Weak electrolytes, polybasic acids, and buffers in anion exchange columns: II. Sodium acetate chloride system, Solvent Extraction and Ion Exchange, 2, 1151 - 1184, 1984 b.

Helfferich, F. and G. Klein, Multicomponent Chromatography - Theory of Interference, Marcel Dekker, New York, 1970. (out of print, available from University Microfilms, Ann Arbor, Michigan, USA)
Morel, F.M.M. and J.G. Hering, Principles and Applications of Aquatic Chemistry, John Wiley and Sons, New York, 1993.

Oelkers, E.H., Physical and chemical properties of rocks and fluids for chemical mass transport calculations, in: Reactive Transport in Porous Media, P.C. Lichtner, C.I. Steefel and E.H. Oelkers (eds.), Reviews in Mineralogy, 34, 131-191, 1996.

Rhee, H.-K., R. Aris and N.R. Amundson, First-Order Partial Differential Equations: Volume II, Theory and Application of Hyperbolic Systems of Quasilinear Equations, Prentice-Hall, Englewood Cliffs, New Jersey, NJ 07632, 1989.

Stumm, W. and J.J. Morgan, Aquatic Chemistry, 2nd ed., John Wiley, 1981.

Appendix: The Shape of Waves

How the adsorption stoichiometry determines the shapes of waves in 2-component systems.


Location of this page

last modified: December 15, 2013
Page compiled by Joachim Gruber