Sorption

While electrochemical problems predominantly require the modelling of charge converting heterogeneous electrochemical reactions, the concentrations of reactants are also significantly changed by non charge-transferring processes such as sorption.

Sorption is a general term for the mechanisms of adsorption and absorption. Adsorption describes the process of molecules being bound to a surface by means of adhesion. Absorption refers to molecules which enter a bulk volume—in electrochemical applications, this is most commonly an electrolyte that is modelled as a solid in Simcenter STAR-CCM+.

The sorption process can be written in the form of a reaction. For the example of water and steam/vapour at an interface, this sorption process can be written as:

1. EQUATION_DISPLAY
H 2 O ( fluid ) H 2 O ( sorbed )
(4164)
The surface sorption rate r s u r f of the process is modeled depending on the imbalance from equilibrium:
2. EQUATION_DISPLAY
r s u r f = r s u r f , 0 ( K eq c fluid a sorbed )
(4165)

where r s u r f , 0 is the sorption rate constant, c f l u i d is the concentration of the reactant, a s o r b e d is the activity of the sorbed species, and K eq is an equilibrium constant, which allows for discontinuous concentrations at interfaces—as in Henry's law.

The surface rates determine the surface flux of chemical species i given as:
3. EQUATION_DISPLAY
N n , s , i = v i r s u r f
(4166)

where v i is the stoichiometric coefficient of species i .

Different methods can be used to parameterize Eqn. (4165).

Springer

To allow for the most general representation, Simcenter STAR-CCM+ models the reaction rate in terms of chemical activities rather than concentrations to account for deviations from ideal solution behavior. As a consequence, Eqn. (4165) becomes:

4. EQUATION_DISPLAY
r s u r f = r s u r f , 0 (   C eq a s o r b e d )
(4167)

where C eq is an equilibrium-related term.

The equilibrium concentration term for the fluid species is computed as:

5. EQUATION_DISPLAY
C eq = ξ λ eq
(4168)

where ξ is the equilibrium concentration constant of the reaction.

λ eq is the dimensionless solute content at equilibrium and is usually computed in terms of chemical activities a :
6. EQUATION_DISPLAY
a = { ( c R T P sat ) for gas phase c  for liquid/intercalated phase
(4169)

where P s a t is the saturation pressure.

The Springer reaction formulation [836] is tailored to model water transport in Polymer Electrolyte Membrane (PEM) fuel cells. The dimensionless solute content at equilibrium λ eq becomes a dimensionless water content at equilibrium, and is expressed as:
7. EQUATION_DISPLAY
λ eq = A a 3 γ + B a 2 γ + C a γ + D a 1
(4170)

where γ is the rate exponent. For a Multi-Component Gas continuum in Simcenter STAR-CCM+, the dimensionless water content at equilibrium is cropped at a = 1 , which corresponds to a saturated gas (for example, water vapour). A , B , C , and D are coefficients that you specify. The default coefficient values in Simcenter STAR-CCM+ are from Springer [836] at 30°C for water.

Wu Li Berg
Similar to the Springer formulation, the Wu Li Berg reaction formulation [838] is tailored to model water transport in PEM fuel cells. The dimensionless water content at equilibrium is defined as:
8. EQUATION_DISPLAY
λ eq = 0.3 + 6 a γ ( 1 tanh ( a γ 0.5 ) ) + 3.9 a γ ( 1 + tanh ( a γ 0.89 0.23 ) )
(4171)
If the reactant activity a exceeds 1.0, the activity is limited to a value of 1.0.
Henry
Unlike the above fuel cell related expressions, Henry's law is not formulated in terms of reference concentrations and dimensionless solute contents—the surface sorption rate is expressed directly as:
9. EQUATION_DISPLAY
r s u r f = r s u r f , 0 ( H e a r γ a p γ )
(4172)
where H e is Henry's constant, a r is the activity of the reactant species, and a p is the activity of the product species, calculated as:
10. EQUATION_DISPLAY
a = { ( c R T P r e f ) for gas phase c  for liquid/intercalated phase
(4173)
where P r e f is the reference pressure ([1.0 atm]).

Sorption: Electrochemical Reaction Heating

The Electrochemical Reaction Heating model accounts for heat contributions which are due to the sorption process (for heat contributions due to reversible and irreversible electrochemical processes, see Eqn. (4143)).

Simcenter STAR-CCM+ calculates the heat contributions that are released during the sorption process according to:
11. EQUATION_DISPLAY
q˙A=ΣiNn,s,iHiWi
(4174)
where Nn,s,i represents the molar concentration flux of species i from the solid to the liquid, Hi is the enthalpy of species i, and Wi is the molar mass of species i.