2.2.2 ONSAGER's Theorem

A rigorous description of thermal influences on the electrical current and vice versa has been presented by ONSAGER2.19 in 1931 [83,84]. His theory discusses the relations of reciprocity of reversible and irreversible processes, where the coupling of the electrical and the thermal subsystems are investigated. For instance, if the electrical driving force is denoted as $ {\mathbf{{X}}}_1{=}{\mathbf{{E}}}$ and the thermodynamic driving force $ {\mathbf{{X}}}_2$ is expressed as,

$\displaystyle {\mathbf{{X}}}_2 {=}- {\frac{1}{T}}{{\boldmath {\nabla}}}{T},$ (2.61)

where $ T$ has been identified as the absolute temperature by CARNOT2.20 [61,85], the corresponding equation system can be formulated with independent equations as
$\displaystyle {\mathbf{{X}}}_1 {=}R_1 {\mathbf{{J}}}_1$     (2.62)
$\displaystyle {\mathbf{{X}}}_2 {=}R_2 {\mathbf{{J}}}_2,$     (2.63)

where $ R_1$ and $ R_2$ are the electrical resistivity and the thermal ``heat resistance'', respectively. The heat resistance is also called thermal resistance $ R_{\mathrm{th}}$ in this thesis. The quantities $ {\mathbf{{J}}}_1$ and $ {\mathbf{{J}}}_2$ are the electrical and the thermal current, respectively. The thermal current density $ {\mathbf{{J}}}_2$ is also called heat flow density $ {{\mathbf{{q}}}_{\mathrm{th}}}$ . Several thermodynamic experiments over the last 150 years have shown that the electrical current is not independent of the temperature. Therefore, equations (2.62) and (2.63) are coupled. Introducing the standard notation, these equations can be adapted by cross coefficients $ R_{12}$ and $ R_{21}$ and represent the ONSAGER relations
$\displaystyle {\mathbf{{X}}}_1 {=}R_{11} {\mathbf{{J}}}_1 + R_{12} {\mathbf{{J}}}_2,$     (2.64)
$\displaystyle {\mathbf{{X}}}_2 {=}R_{21} {\mathbf{{J}}}_1 + R_{22} {\mathbf{{J}}}_2.$     (2.65)

For this equation system, THOMSON proposed the relation

$\displaystyle R_{12} {=}R_{21},$ (2.66)

which is also called ``reciprocity theorem'' of the ONSAGER relations. However, (2.66) implies that this relation follows from symmetric principles of thermodynamic theory. Hence, the reciprocity theorem neglects the loss during heat conduction and energy conversion and relation (2.66) assumes a balanced energy flow between the two subsystems. Thus, a steady stage is assumed with the request of (2.66) [61], where equilibrium conditions are applicable only within short range.

The principle of microscopic reversibility in (2.66) is less general than the second fundamental law of thermodynamics [83]. For further investigated coupled systems, the currents $ R_{ij}$ may have different signs due to the different directions of the energy flows. Therefore, (2.66) is not sufficient enough to fulfill the second law of thermodynamics. Hence, the necessary condition for the equation system consisting of (2.64) and (2.65) to guarantee the second law with $ {{\partial}_{t}}{S} \geq 0$ yields

$\displaystyle R_{12} + R_{21} \leq 2 \sqrt{R_{11}R_{22}}.$ (2.67)

This necessary condition has been originally proposed by BOLTZMANN in 1887 [86].

Writing the ONSAGER relations (2.64) and (2.65) as functions of driving forces $ {\mathbf{{X}}}_i$

$\displaystyle {\mathbf{{J}}}_1 {=}L_{11} {\mathbf{{X}}}_1 + L_{12} {\mathbf{{X}}}_2$     (2.68)
$\displaystyle {\mathbf{{J}}}_2 {=}L_{21} {\mathbf{{X}}}_1 + L_{22} {\mathbf{{X}}}_2,$     (2.69)

where the necessary condition of type (2.67) remains valid accordingly for $ L_{ij}$ as

$\displaystyle L_{12} + L_{21} \leq 2 \sqrt{L_{11}L_{22}}.$ (2.70)

To consider the ONSAGER relations in terms of energy, (2.68) and (2.69) can be multiplied by $ {\mathbf{{X}}}_1$ and $ {\mathbf{{X}}}_2$ , respectively, leading to

$\displaystyle {\mathbf{{J}}}_1 \cdot {\mathbf{{X}}}_1 {=}L_{11} {\mathbf{{X}}}_1 \cdot {\mathbf{{X}}}_1 + L_{12} {\mathbf{{X}}}_2 \cdot {\mathbf{{X}}}_1,$     (2.71)
$\displaystyle {\mathbf{{J}}}_2 \cdot {\mathbf{{X}}}_2 {=}L_{21} {\mathbf{{X}}}_1 \cdot {\mathbf{{X}}}_2 + L_{22} {\mathbf{{X}}}_2 \cdot {\mathbf{{X}}}_2.$     (2.72)

These equations represent the products of the driving forces $ {\mathbf{{X}}}_i$ and displacements of types of flow $ {\mathbf{{J}}}_i$ . The result of (2.71) and (2.72) can be described as the dissipated energy per volume and per time and reads

$\displaystyle T\theta {=}\sum_{i} {\mathbf{{J}}}_{i}\cdot{\mathbf{{X}}}_{i},$ (2.73)

where $ \theta$ is the entropy generation rate per unit volume and follows from the second law of thermodynamics (2.35)
$\displaystyle {{\partial}_{t}}{S} \geq 0 \quad\Longrightarrow\quad T {{\partial}_{t}}{S}$   $\displaystyle \geq 0,$ (2.74)
$\displaystyle T\theta {=}{\frac{\partial(T {{\partial}_{t}}{S})}{\partial V}} {=}\sum_{i} {\mathbf{{J}}}_{i}\cdot{\mathbf{{X}}}_{i}$   $\displaystyle \geq 0,$ (2.75)
$\displaystyle \Longrightarrow \quad \theta {=}\frac{\partial({{\partial}_{t}}S)}{\partial V} {=}{\frac{1}{T}}\sum_{i} {\mathbf{{J}}}_{i}\cdot{\mathbf{{X}}}_{i}$   $\displaystyle \geq 0,$ (2.76)

where the entropy generation rate $ {{\partial}_{t}}{S}$ can be determined by the sum of the power densities of all contributing subsystems. Hereby, the parts of the sums can be identified as the power densities of the participating systems which are determined by chemical reactions, the power loss due to heat transfer and JOULE's self-heating, and the power loss due to diffusion processes. The power density of chemical reactions can be expressed by scalar-valued quantities as

$\displaystyle T\theta^{\mathrm{chem}} {=}\sum_{i} J^{\mathrm{chem}}_{i} X^{\mat...
...frac{1}{V}} {{\partial}_{t}}{\xi}_j \right) \left(- \sum_j \nu_j \mu_j \right),$ (2.77)

where $ J^{\mathrm{chem}}_{i}$ is determined by the chemical reaction rate $ {{\partial}_{t}}{\xi}_j$ per unit volume $ V$ . The chemical driving force is represented by $ X^{\mathrm{chem}}_{i}=-\sum_j\nu_j\mu_j$ , where $ \mu_j$ denotes the chemical potential and $ \nu_j$ the stoichiometric coefficient of the participating atom. The electrical power density can be identified by
$\displaystyle T\theta^{\mathrm{el}} {=}\sum_{i}
{\mathbf{{J}}}^{\mathrm{el}}_{i...
...bf{{X}}}^{\mathrm{el}}_{i} {=}\sum_{i} {\mathbf{{J}}}_i\cdot{\mathbf{{E}}}_i{=}$      
$\displaystyle {=}- \sum_{i} {\mathbf{{J}}}_i \cdot {{\boldmath {\nabla}}}{\varphi_i},$     (2.78)

where $ {\mathbf{{J}}}$ and $ {\mathbf{{E}}}$ are the electrical current density and the electrical field, respectively. The electric field can be expressed by the spatial gradient of an electrical potential $ \varphi$ . For electro-magnetical subsystems, the power density has to be appropriately adapted as discussed in Section 2.2.3. Another important contribution to the global entropy increase is the power loss due to thermal heat flow, which can be expressed in terms of (2.73) by
$\displaystyle T\theta^{\mathrm{heat}} {=}\sum_{i}
{\mathbf{{J}}}^{\mathrm{heat}...
...cdot
\left(T_i {{\boldmath {\nabla}}}{\left({\frac{1}{T_i}} \right)} \right){=}$      
$\displaystyle {=}\sum_{i} {{{\mathbf{{q}}}_{\mathrm{th}}}}_{,i}\cdot \left( -
{\frac{1}{T_i}} {{\boldmath {\nabla}}}{T_i} \right),$     (2.79)

where $ {{\mathbf{{q}}}_{\mathrm{th}}}$ represents the local heat flux density. The second term in (2.79) depicts the thermal driving force according to FOURIER's empirical law. For diffusion processes, the power density can be identified as
$\displaystyle T\theta^{\mathrm{diff}}{=}
\sum_i{\mathbf{{J}}}_i^{\mathrm{diff}}...
...eft( - T_i {{\boldmath {\nabla}}}{\left({\frac{\mu_i}{T_i}}\right)}
\right) {=}$      
$\displaystyle {=}
\sum_i{\mathbf{{J}}}_i^{\mathrm{mol}}\cdot \left( - {{\boldmath {\nabla}}}{\mu_i} +
{\frac{\mu_i}{T_i}}{{\boldmath {\nabla}}}{T_i} \right),$     (2.80)

where $ {\mathbf{{J}}}^{\mathrm{mol}}_i$ is the mole number per unit area and time of the contributing species $ i$ . The driving force of diffusion processes is determined by the gradient of the chemical potential and by the gradient of the temperature.

Since the power density of diffusion processes has been determined for no external forces, an extension for the applied electrical field has to be made by introducing an additional term that depicts the force acting on charged particles inside the simulation domain. Hence, equation (2.80) has to be modified as

$\displaystyle T\theta {=}\sum_i{\mathbf{{J}}}_i^{\mathrm{mol}} \cdot \left( - {...
...\boldmath {\nabla}}}{T_i} - Z^*_i n_i {{\boldmath {\nabla}}}{\varphi_i}\right),$ (2.81)

where $ Z^*_i$ is the effective valence charge of the species $ i$ , $ n_i$ is the species concentration per mole, and $ \varphi_i$ is the corresponding electrical potential.

To conclude ONSAGER's thermodynamical treatment, the overall power density is thus given by the sum of the power densities of all contributing subsystems as

$\displaystyle T\theta$   $\displaystyle {=}T \sum_i\theta_i {=}T \left(\theta^{\mathrm{chem}} + \theta^{\mathrm{el}} + \theta^{\mathrm{heat}} + \theta^{\mathrm{diff}} \right)$ (2.82)
$\displaystyle T\theta$   $\displaystyle {=}- {\frac{1}{V}} {{\partial}_{t}}{\xi}_j \sum_j \nu_j \mu_j$  
    $\displaystyle - \sum_{i} {\mathbf{{J}}}_i \cdot {{\boldmath {\nabla}}}{\varphi_i}$  
    $\displaystyle + \sum_{i} {\mathbf{{q}}}_i\cdot \left( - {\frac{1}{T_i}}{{\boldmath {\nabla}}}{T_i} \right)$  
    $\displaystyle + \sum_i{\mathbf{{J}}}_i^{\mathrm{mol}} \cdot \left( - {{\boldmat...
...\boldmath {\nabla}}}{T_i} - Z^*_i n_i {{\boldmath {\nabla}}}{\varphi_i}\right),$ (2.83)

where the thermodynamic power density $ T\theta$ is determined by the contributing chemical reactions, the electrical burden, heat flows, and molar diffusion processes.


Stefan Holzer 2007-11-19