Aimen E. Gheribi,Patrice Chartrand
Abstract:In aluminium electrolysis cells, a ledge of frozen electrolyte is formed, attached to the sides of the cell. The control of the side ledge thickness is essential in ensuring a reasonable lifetime for the cells. Numerical modelling of the side ledge thickness requires an accurate knowledge of the thermal transport properties as a function of temperature. Unfortunately, there is a considerable lack of experimental data for the large majority of the phases constituting the side ledge. The aim of this work is to provide, for each phase possibly present in the side ledge, a formulation of the thermal conductivity as a function of both temperature and size. To achieve this, we consider reliable physical models linking the density of the lattice vibration energy and the phonon mean free path to key parameters: the high temperature limit of the Debye temperature and the Gruneisen constant. These model parameters can be obtained by simultaneous fitting of (i) the heat capacity, (ii) the thermal expansion tensor coefficient and (iii) the adiabatic elastic constants, on relevant physical models. Where data is missing, first principles (ab initio) calculations are used to determine directly the model parameters. For compounds for which data is available, the model's predictions are found to be in very good agreement with the reported experimental data.
What problem does this paper attempt to address?
This paper attempts to solve the problem of the thermal conductivity of the side ledge in aluminum electrolytic cells, especially the thermal conductivity of each phase at different temperatures and grain sizes. Specifically, the paper aims to provide a mathematical expression for the thermal conductivity of each phase that may exist in the side ledge, with temperature and size as variables.
### Research Background and Problem
During the aluminum electrolysis process, a layer of frozen electrolyte will form on the side wall, which is called the side ledge. Controlling the thickness of this side ledge is crucial for ensuring the reasonable service life of the electrolytic cell. Numerical simulation of the side ledge thickness requires an accurate understanding of its heat transfer characteristics as a function of temperature. However, experimental data for most of the phases that make up the side ledge are very scarce.
### Main Objectives
The main objective of the paper is to provide a formula for the thermal conductivity for each phase that may appear in the side ledge, which depends on temperature and size. To this end, the authors considered reliable physical models, linking the lattice vibration energy density and the phonon mean - free path to key parameters, such as the high - temperature limit of the Debye temperature and the Grünneisen constant. These model parameters can be obtained by simultaneously fitting the following physical quantities:
- Specific heat capacity
- Coefficient of thermal expansion
- Adiabatic elastic constant
For compounds lacking data, the model parameters are directly determined using ab initio calculations. For compounds with experimental data, the model prediction results are in very good agreement with the reported experimental data.
### Formula Summary
1. **Temperature and Size Dependence of Thermal Conductivity**:
\[
K_{\text{ms}}=\psi[\varphi_1(x, T), \varphi_2(x, T), \cdots, \varphi_n(x, T), K_1(T, d_1, p_1), K_2(T, d_2, p_2), \cdots, K_n(T, d_n, p_n), P]+\Psi_{\text{ms}}
\]
where:
- \(x\) is the composition vector
- \(T\) is the average thermodynamic temperature
- \(\varphi_i\) is the volume fraction of the \(i\) - th phase
- \(K_i(T, d_i, p_i)\) is the thermal conductivity of the \(i\) - th phase, which depends on temperature, average grain size and average porosity
- \(P\) is the total intergranular porosity
- \(\Psi_{\text{ms}}\) takes into account all other microstructural parameters
2. **Lattice Thermal Conductivity**:
\[
K_{\text{lat}}(T)=A\times\bar{m}\frac{\theta_D^3}{n^{2 / 3}}\times\frac{1}{\gamma_\infty^2}(\infty - 0.514\gamma_\infty + 0.228)\times\delta(T)e^{\phi/T}
\]
where:
- \(A\) is the constant \(6.51982\times 10^{-3}(\frac{k_B}{\hbar})^3\)
- \(\bar{m}\) is the average mass per atom
- \(\theta_D\) is the Debye temperature
- \(n\) is the number of ions in the unit cell
- \(\gamma_\infty\) is the high - temperature limit of the Grünneisen constant
- \(\delta(T)\) is the temperature - dependent average atomic volume
- \(\phi\) is a constant, assumed to be 0 for fluorides, oxides and carbides
3. **Grain Size Effect**:
\[
K_{\text{lat}}(T, d)=K_{\text{lat}}^\infty(T)\left[1-\sqrt{\frac{\sigma(T)}{d}}\right]