Implementation of contamination by irradiated materials

[prisma34]

Add to thermochemical properties packages the contamination modules.
The new modules will be named <heavymetal_name>_contamination.py and will be placed within this path: properties/<heavymetal_name>_thermochemical_properties.

The work will be organized according to the following tasks:

• Implementation of lead related properties.
• Implementation of bismuth related properties.
• Implementation of LBE related properties.
• Updating doc accordingly.
• Implementation of the corresponding tests.

[prisma34]

After closing #184, #185, #186, three actions are still to be made:

• Implementation of metallic compounds' vapour pressure correlations.
• Import of Iodine and Caesium lead-related correlations in lbe_contamination.py: this is to implement each correlations only once.
• Implementing Henry constant of Lead and Bismuth in lbe_thermochemical.py.

[ilariacasalbore]

After reviewing the implemented correlations, here are a few observations:

Polonium in Lead:
• Correlation Abakumov 1994 for PbPo vapour pressure should be Abakumov 1974a, 1982 (5.16, section 5.2.1.1.5.2.1);
• Eq. 5.23 and 5.26, section 5.2.1.1.5.2.3 for the Henry constant and activity coefficient could be added: defined for Po in LBE but valid for lead also (5.26 considered the best available estimate of γ for Po in Pb, see section 5.2.1.1.5.3 and addendum).

Cadmium and Thallium in Lead:
• Correlation for the activity coefficient provided at a single temperature, to decide whether implementing them or not (eq. 5.36, section 5.2.1.1.7.2 and eq. 5.40, section 5.2.1.1.8.2).

Iodine in Lead:
• To check whether log in eq. 5.45 (5.2.1.1.9.1) is log10 or natural log;
• Eq. 5.44 and 5.45 are complementary: to decide if both can be provided as default equations, for the relative temperature ranges.

Cesium in lead:
• To check whether log in eq. 5.54 (5.2.1.1.11.3) is log10 or natural log;
• Alternative equation for the thermodynamic activity coefficient could be added: eq. 5.52 (5.2.1.1.11.3), with temperature range 953-1104 K.

Polonium in Bismuth:
• Default correlation for the thermodynamic activity coefficient could be changed to eq. 5.21: considered the best available estimate of γ, see section 5.2.1.1.5.3;

Cadmium and Thallium in Bismuth:
• Correlation for the activity coefficient provided at a single temperature, to decide whether implementing them or not (eq. 5.36, section 5.2.1.1.7.2 and eq. 5.41, section 5.2.1.1.8.2).

Iodine in Bi:
• To check temperature range of vapour pressure correlation.

Polonium in LBE:
• Range of eq. 5.23 and 5.25 could be extended, see the addendum.

Mercury in LBE:
• To check whether log in eq. 5.34 (5.2.1.1.6.2) is log10 or natural log;
• To check range of correlations for Henry constant and thermodynamic activity coefficient (see also addendum).

Cadmium in LBE:
• To range of eq. 5.38 (5.2.1.1.7.2) and check whether log is log10 or natural log;
• Correlation for vapour pressure could be derived from correlations for henry constant and thermodynamic activity coefficient.

Thallium in LBE:
• To check range of Henry constant and thermodynamic activity coefficient;
• To check whether log in Henry constant correlation is log10 or natural log;
• Correlation for vapour pressure could be derived from correlations for henry constant and thermodynamic activity coefficient.

Iodine in LBE:
• To check whether correlations implemented for Henry constant are correct for Iodine in LBE (most correlations are defined for Iodine in Pb) and the given temperature ranges.

Cesium in LBE:
• Additional correlation for the thermodynamic activity coefficient could be added (it is not between the recommended correlations).