Russian Journal of Nonlinear Dynamics, 2021, vol. 17, no. 1, pp. 49-57. Full-texts are available at http://nd.ics.org.ru DOI: 10.20537/nd210105
NONLINEAR PHYSICS AND MECHANICS
MSC 2010: 78-04, 82B30, 82D10
Modelling of the Thermodynamic Properties of the Plasma Mixture
A. S. Polyukhin
Numerical modelling of the thermodynamic properties of plasma mixture is performed using the Thomas-Fermi model with two different approaches. For this purpose, a numerical algorithm, as well as program realization, is developed to solve the Thomas-Fermi equations with quantum-exchange corrections. For the first time a comparison between different methods for taking account of the heterogeneous composition of plasma is made and an algorithm for estimating the corrections for mixtures is developed.
Keywords: computational model, magnetized plasma, numerical method, thermodynamic property
1. Introduction
The main purpose of this work is to develop an algorithm for modelling thermodynamic properties of plasma mixtures in a magnetic field using the Thomas-Fermi model [1] with quantum and exchange corrections and to compare different approaches to the given problem. All calculations are made by employing two different approaches: by solving the system of equations of the Thomas -Fermi model with the condition of thermodynamic equilibrium and by using average ion approximation.
Received January 25, 2021 Accepted March 04, 2021
Alexander S. Polyukhin polyukhinas@gmail.com
Joint Institute for High Temperatures of the Russian Academy of Sciences
ul. Izhorskaya 2, Moscow, 125412 Russia
Moscow Institute of Physics and Technology
Institutskiy per. 9, Dolgoprudny, Moscow Region, 141701 Russia
Correct and accurate modeling of plasma thermodynamic parameters is a complex and time-consuming task. The analysis of the Thomas-Fermi model with quantum-exchange corrections, comparison with other models and its extension to the region of multicomponent systems was presented earlier, but not in such a wide range of parameters.
In this paper, we study the thermodynamic parameters of one-component and multicomponent plasma in an external magnetic field based on the Thomas-Fermi model with quantumexchange corrections. All data are calculated in the region of temperatures of 10 eV < T < 105 since it is a preferable area for the problems of fusion energy [2].
2. Numerical methods
The main idea of the first approach is solving the system of Thomas-Fermi equations [3] for each component of mixture (N substances) with the additional condition that the sum of components' volumes is equal to the volume of the whole system:
m„-
Pi
(2.1)
where mi is a mass fraction of the component, pi is partial density of the component, and p is average density. Using the condition of thermodynamic equilibrium that sets equality of the chemical potential of each component of the system, one can obtain a system of equations which can be used to determine potential
^ - a rl
dx2 " -"i 2 I
#i
(i = 1, 2, ..., N),
$(1) =
dx
Pi
Vu
x=1
0i(O) =
A
Or,
0i
P
Vi = Vj = —Q (V'i, i), r 0i(x) _ Vi(r) + Pi
r0i x
9
A
r0l = 1,388 ( -i Pi
1/3
n
h(x) = J
Vk dy
1 + exp(y - x) '
(2.2)
r0i,
where Vi(r) is electrostatic potential, 0 is reduced electrostatic potential, ¡j,i is chemical potential, is reduced chemical potential, r0i is the boundary radius of an atomic cell, Zi is the charge of a nucleus, Ai is atomic weight, and Ik(x) is the Fermi-Dirac function. Here and in what follows, all values with index i correspond to the simple substance of the mixture. Radii of atomic cells should be chosen using (2.1) and the condition of thermodynamic equilibrium. Using linearization and additional transformations, the system (2.2) can be solved numerically using iteratively the sweep method [4].
9
x
After estimation of potential 0 all basic thermodynamic functions of selected species can be calculated according to the Thomas-Fermi model by the following relations [5]:
Se —
4v/203/2rg
n
4
7
x2 L
3/3/2(-i?) ~ z I 1 J3/2 (2d)5/2
3 J V x
0
dx
- Z$'(0),
Pe
^sM-V),
Ee — Epe. + Ek
Je Epe + Eke, 1
Eke = ^'2 jx%,2
0
—
pe n2
where rj = —0(1), v = ^kVq, and Se, Pe, Ee are thermodynamic functions of the electron system. The pressure, internal energy, and entropy of the nuclear component are determined by the classical formulas of an ideal gas. In order to get thermodynamic properties of the system, the sum of the components of the mixture should be carried out:
;05/2
1
/3/^-3 Jx213/2
E=
S=
Y,ElhlAl
EM* ' l
SihiAi EM* '
(2-3)
(2-4)
where S is the entropy of the system, E is the internal energy of the system, hi is the mole fraction of the i component, Si is the entropy of the i component, and Ei is the internal energy of the i component. The pressure of the system can be found using the relation of the classical Thomas-Fermi model because the pressure of all components is the same due to the condition of thermodynamic equilibrium. Thermodynamic consistency is achieved because the pressure depends on chemical potential, which is the same for all species due to condition (2.2).
Introduction of quantum-exchange corrections in the Thomas-Fermi model for mixtures allows expanding the limit of applicability in the region of lower temperatures. The method that provides corrections is based on the use of the quantum mechanical approach in the derivation of the Thomas-Fermi equation [6]. The theory provides the following equation for correction of the potential:
Xi
n
r2 r
"'Oí-'1/2
x
Xl
n
_r2 y I r 0lxY
ti (x)
Xl(0)— 0, Xl(1) —1,
where xi is the reduced correction of the Thomas-Fermi potential and Y = I1/2(x)I'1 /2 (x) +
X r -1 2
+ 6 f I'1 ^ (t) dt. In this approach, the correction for pressure, internal energy and entropy
1
x
can be written as
AP
192
AE; = —rSi
f_ 2 3tt2'
1 1
I .!■ \i ) /1 , f ^^ dx + 2 f x2Y (^ | At
Lo
+ -
6n
-X-(0)+0,2690017Z5/3,
AS, =
2i? 2
xXi(x)/V2
fai(x)
dx + 4
x2F
fa(x)
dx
+ -
=Xi(0),
.o o
where AP,, AE,, AS, are quantum-exchange corrections of thermodynamic functions. Since the total contribution of corrections for thermodynamic functions requires sums such as in relations (2.3) and (2.4), it will inevitably violate the condition of thermodynamic equilibrium. For equilibrium, the system should have equal values of pressure and chemical potential for each species, which cannot be achieved because correction for these thermodynamic properties is individual for every component of the mixture. It should be mentioned that the values of these corrections depend greatly on temperature and density, so in the region of high temperatures and densities the corrections become negligible and the condition of thermodynamic equilibrium will be performed automatically.
Another approach of mixtures' modelling is the average ion approximation [7]. This method was first applied to a similar problem of determining the thermodynamic properties using the Saha ionization model. The main idea is to replace the original system of various ions with a new system consisting of ions averaged over all components:
Z = hiZi, A = hA.
To calculate thermodynamic properties and their corrections using the Thomas-Fermi model, an average ion approximation can be used as well by replacing the system of equations (2.2) with only one equation of an averaged cell. In this approximation, the whole mixture can be described by the following system:
d2(f) dx2
= axl
1/2
fa'(l) =
dfa dx
x=1
= e=-ri>
x"
n
~roh/2
«°> = W„' my W5
x I n
x(0)=0, x'(l) = l,
where fa is the reduced Thomas-Fermi potential for the averaged ion and x stands for reduced corrections of the Thomas -Fermi potential for the averaged ion.
Thermodynamic functions can be calculated using relations for the classical model with the Thomas-Fermi potential of the averaged ion. This approach satisfies the condition of thermo-dynamic equilibrium since each species will have equal values of average chemical potential, thus all thermodynamic properties for all species of the mixture will be equal. The accuracy of this approach should be established by comparing the results with the existing data.
1
1
3. Results and discussion
To model the thermodynamic properties of plasma mixtures that could be used in problems of thermonuclear fusion, two mixtures were chosen: 60 % gold + 40 % silver and 90 % beryllium + 10 % tungsten. The first mixture of heavy cores can be implemented as a target component for problems of inertial and magnetic fusion [8-11]. The second one can be used in the creation of diverters and blankets of tokamaks. It should be mentioned that calculations performed in this work could be applied to a system with a strong magnetic field as it does not affect the present computations [12].
Figure 1 corresponds to the determination of pressure in the isochoric process (p = 1 g/cm3) as a function of temperature using different approaches. It is seen that all methods give similar results. However, considering quantum-exchange corrections reduce the pressure in the region of low temperature, since these corrections have a negative value. Differences between the method of taking account of the heterogeneity of plasma using the average ion approximation and the approach based on solving the system of equations with the condition of thermodynamic equilibrium are not significant. This fact is in direct agreement with the theory since a negative pressure is described by repulsive forces of exchange interactions. Such forces are completely ignored in the classical Thomas-Fermi model because their contribution is not significant at high temperatures.
io7-105-
10"1
101 102 103 104 105
T, eV
Fig. 1. Pressure of the mixture of gold and silver, based on the Thomas-Fermi model with and without quantum-exchange correction at a density of 1 g/cm3.
Considering the solution method of the system of Thomas-Fermi equations with the additional condition of thermodynamic equilibrium, it can be inferred that, as mentioned above, taking account of quantum-exchange correction violates the condition of equilibrium. However, modelling shows that this approach insignificantly changes the result for heavy cores with a temperature higher than 10 eV.
The second mixture has a completely different picture of the distribution of thermodynamic functions. The main reason is variation in charge between the components of the mixture. Since the difference between the core charge of tungsten and that of beryllium is significant, the difference in averaging is much more serious. This fact leads to the results shown in Fig. 2. The
60 % (Au) + 40 % (Ag)
Average ion Multi mixture
---Average ion with corrections - - Multi mixture with corrections i i
106 10s
«S 104 o
a," io3 102 io1 10°
IO1 IO2 103 IO4 105
T, eV
Fig. 2. Pressure of the mixture of beryllium and tungsten, based on the Thomas-Fermi model with and without quantum-exchange correction at a density of 1 g/cm3.
method of solving equations of the Thomas -Fermi model with the additional condition of ther-modynamic equilibrium usually gives lower values as compared to the average ion approximation. In the region of high temperatures, quantum-exchange corrections are negligible, all methods give the same values.
The modelling values of internal energy for the first mixture are almost equal regardless of the method of averaging species and taking into account the quantum-exchange corrections. The results presented in Fig. 3 can be explained by the fact that the contribution of the quantum exchange corrections for the internal energy at the chosen density begins to appear at much lower temperatures (T < 10 eV), as shown in Fig. 4. Since the plasma for the problems of
107
106
M 105 a ^ 104
103 102
101 102 103 104 105
T, eV
Fig. 3. Internal energy of the mixture of gold and silver, based on the Thomas-Fermi model with and without quantum-exchange correction at a density of 1 g/cm3.
90 % (Be) + 10 % (W)
— Average ion Multi mixture
---Average ion with corrections
- - Multi mixture with corrections
=F
=F
60 % (Au) + 40 % (Ag)
/ X / /
/ /
-Average ion Multi mixture - ---Average ion with corrections - - Multi mixture with corrections i i
10 106 M105 tf 104 103 102
60 % (Au) + 40 % (Ag
Average ion Multi mixture
Average ion with corrections Multi mixture with corrections
10"
10° 101
10J T, eV
10d
104
105
Fig. 4. Internal energy of the mixture of gold and silver, based on the Thomas-Fermi model with and without quantum-exchange correction at a density of 1 g/cm3, in a wider temperature range.
thermonuclear fusion is not used at such low temperatures, the calculation gives similar values, regardless of whether or not the quantum-exchange corrections are taken into account.
Considering the values of internal energy in the second mixture, one can draw similar conclusions, as in the case of the calculation of pressure. The values obtained using the same method do not differ when considering the quantum-exchange corrections for the reasons described above. But the values obtained by different methods are different and this is seen in Fig. 5. As in the case of pressure, the difference between the charges of the beryllium and tungsten nuclei leads to a strong variation of internal energy values, achieved by two approaches of averaging.
60 % (Au) + 40 % (Ag)
10'
10<
10;
60
10'
10:
10:
10
-----
-Average ion Multi mixture ---Average ion with corrections Multi mixture with corrections " i i
101
102
103 T, eV
104
105
Fig. 5. The internal energy of the mixture of beryllium and tungsten, based on the Thomas-Fermi model with and without quantum-exchange correction at a density of 1g/cm3.
All the arguments related to the internal energy can also be referred to the entropy of the system described in Fig. 6. In the region of relatively low temperatures, significant errors in the calculated values of entropy with corrections may appear due to the use of a simple approximation for the Fermi-Dirac functions, while there are more complex analogues [13, 14]. It is shown in [15, 16] that the calculation of internal energy and entropy according to the Thomas-Fermi model is less accurate than that of pressure, so in the temperature range (T < 100 eV), quantum exchange corrections can give nonphysical results. Such an approximation of the Fermi-Dirac functions was used since the plasma temperatures suitable for problems of thermonuclear fusion are much higher than the value of 100 eV. Moreover, areas where quantum-exchange corrections reach the value of a thermodynamic function should not give a physical result, since the correction is derived under the assumption of being small relative to the thermodynamic function it corrects. At high temperatures (T > 10 keV), all the methods give similar results for evaluation of the entropy.
105
~ io4 >
u
fe 103 Kf 102
101
101 102 103 104 105
T, eV
Fig. 6. Entropy of the mixture of gold and silver based on the Thomas-Fermi model with and without quantum-exchange correction at a density of 1 g/cm3.
Based on the above, it can be concluded that the calculations performed using the two approaches give minor differences, especially in the region where quantum-exchange corrections are neglegible, but the average ion approximation has several advantages. It provides the condition of thermodynamic equilibrium in the field of applicability of quantum exchange-corrections and also has a simpler implementation.
4. Conclusion
As a result of the study, two algorithms have been developed and implemented for modelling the thermodynamic parameters of plasma mixture based on the Thomas -Fermi model with quantum-exchange corrections. The approach to solving the system of equations for the Thomas -Fermi potential with corrections has been proposed for the first time.
As an example of modelling the thermodynamic properties of a substance, gold-silver and beryllium-tungsten mixtures were considered. For them, the distribution of thermodynamic functions and their quantum-exchange corrections were obtained using different approaches. More-
60 % (Au) + 40 % (Ag)
__
— ■--~
1 / / -Average ion Multi mixture ---Average ion with corrections - - - Multi mixture with corrections i i
I 1 1
over, the differences between the approaches in the description of the plasma mixture were discussed. For the first time, the thermodynamic parameters of gold-silver and beryllium-tungsten
mixtures were calculated at the selected concentrations.
References
[1] Zel'dovich, Ya. B. and Raizer, Yu. P., Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena, New York: Dover, 2002, p. 178.
[2] Ryzhkov, S.V. and Chirkov, A. Yu., Alternative Fusion Fuels and Systems, Boca Raton, Fla.: CRC Press, Taylor & Francis Group, 2019.
[3] Shemyakin, O.P., Levashov, P. R., and Krasnova, P. A., TFmix: A High-Precision Implementation of the Finite-Temperature Thomas-Fermi Model for a Mixture of Atoms, Comput. Phys. Commun, 2019, vol. 235, pp. 378-387.
[4] Kuzenov, V. V. and Shumaev, V. V., Description of the Thermodynamic Properties of Plasma in Saha and Thomas-Fermi Models, Prikl. Fiz, 2015, no. 2, pp. 32-36 (Russian).
[5] Nikiforov, A.F., Novikov, V. G., and Uvarov, V. B., Quantum-Statistical Models of Hot Dense Matter: Methods for Computation Opacity and Equation of State, Progr. Math. Phys., vol.37, Basel: Birkhauser, 2005, p. 34.
[6] Kirzhnits, D. A., Quantum Corrections to the Thomas-Fermi Equation, JETP, 1957, vol.5, no. 1, pp. 64-71; see also: Zh. Eksper. Teoret. Fiz., 1957, vol.32, no. 1, pp. 115-123.
[7] Belov, A. A., Golovanov, R.V., Kalitkin, N. N., Kozlitin, I. A., Koryakin, I. A., and Kuzmina, L.V., TEFIS Database. Thermodynamic Functions of Substances, Preprint No. 219, Moscow: KIAM, 2018 (Russian).
[8] Kuzenov, V. V. and Ryzhkov, S.V., Calculation of Plasma Dynamic Parameters of the Magneto-Inertial Fusion Target with Combined Exposure, Phys. Plasmas, 2019, vol. 26, no. 9, Art. No. 092704.
[9] Kuzenov, V.V., Ryzhkov, S.V., and Shumaev, V. V., Computer Analysis of Transport, Optical and Thermodynamic Properties of Plasma, High Temp. Mater. Process., 2014, vol. 18, nos. 1-2, pp. 99-109.
[10] Kuzenov, V. V. and Ryzhkov, S.V., Radiation-Hydrodynamic Modeling of the Contact Boundary of the Plasma Target Placed in an External Magnetic Field, Prikl. Fiz., 2014, no. 3, pp. 26-30 (Russian).
[11] Kuzenov, V.V. and Ryzhkov, S.V., Mathematical Modeling of Plasma Dynamics for Processes in Capillary Discharges, Russian J. Nonlinear Dyn., 2019, vol. 15, no. 4, pp. 543-550.
[12] Liberman, M. A. and Johansson, B., Properties of Matter in Ultrahigh Magnetic Fields and the Structure of the Surface of Neutron Stars, Physics-Uspekhi, 1995, vol.38, no. 2, pp. 117-136; see also: Uspekhi Fiz. Nauk, 1995, vol. 165, no. 2, pp. 121-142.
[13] Golovanov, R. V. and Lutskii, K.I., Computation of the Integral Fermi-Dirac Function, Math. Models Comput. Simul., 2012, vol. 4, pp. 464-470.
[14] Dyachkov, S. A. and Levashov, P. R., Region of Validity of the Finite-Temperature Thomas-Fermi Model with Respect to Quantum and Exchange Corrections, Phys. Plasmas, 2014, vol.21, no. 5, 052702, 6 pp.
[15] Kuzenov, V. V., Ryzhkov, S. V., and Shumaev, V. V., Numerical Thermodynamic Analysis of Alloys for Plasma Electronics and Advanced Technologies, Probl. Atom. Sci. Tech., 2015, vol. 98, no. 4, pp. 53-56.
[16] Kuzenov, V. V., Ryzhkov, S.V., and Shumaev, V. V., Thermodynamic Properties of Magnetized Plasma Evaluated by Thomas-Fermi Model, Prikl. Fiz, 2014, no. 3, pp. 22-25.