Procedure for determining the number of thermal diffusion columns in square cascade for separation of Ne stable isotopes

NUCLEAR ENERGY SCIENCE AND ENGINEERING

Procedure for determining the number of thermal diffusion columns in square cascade for separation of Ne stable isotopes

Fatemeh Mansourzadeh ，

Mohammad Mahdi Shadman，

Javad Karimi Sabet，

Valiyollah Ghazanfari

Nuclear Science and Techniques

Vol.34, No.5

Article number 68

Published in print May 2023

Available online 13 May 2023

DOI：10.1007/s41365-023-01217-z

35100

The thermal diffusion column represents one method of separating stable isotopes. This method is advantageous for small-scale operations because of the simplicity of the apparatus and small inventory, especially in gas-phase operations. Consequently, it has attracted attention for its applicability in tritium and noble gas separation systems. In this study, the R cascade was used to design and determine the number of columns. A square cascade was adopted for the final design because of its flexibility, and calculations were performed to separate ²⁰Ne and ²²Ne isotopes. All the R cascades that enriched the Ne isotopes by more than 99% were investigated, the number of columns was determined, and the square cascade parameters were optimized using the specified columns. Additionally, a calculation code ''RSQ_CASCADE'' was developed. A unit separation factor of three was considered, and the number of studied stages ranged from 10 to 20. The results showed that the column separation power, relative total flow rate, and required number of columns were linearly related to the number of stages. The separation power and relative total flow decreased and the number of columns increased as the stage number increased. Therefore, a cascade of 85 columns is recommended to separate the stable Ne isotopes. These calculations yielded a 17-stage square cascade with five columns in each stage. By changing the stage cut, feed point, and cascade feed flow rate, the best parameters for the square cascade were determined according to the cascade and column separation powers. As the column separation power had a maximum value in cascade feed 50, it was selected for separating Ne isotopes.

Stable IsotopeThermal diffusionCascadeSeparationNeon

Introduction

Stable isotopes are widely applied as precursors for radioisotopes and used in fields such as medicine, electronics, mineral exploration, metallurgy, industry, and physics research [1, 2]. Understanding separation techniques would facilitate the production of a substantial of stable isotopes. Thermal diffusion is an important process because of its positive attributes, such as a high separation coefficient, simplicity, and cascade operation. Moreover, the apparatus involved is advantageous for small-scale operation because of its simplicity and small inventory, particularly for gas-phase operations. Consequently, it has attracted attention for its applicability in tritium and noble gas separation systems [3, 4]. This method produces expensive isotopes with high purity on a small scale (1–2 g of the target isotope per day at a concentration of over 99%). The temperature difference is the most crucial factor affecting mass transfer in a column. Separation using a temperature gradient was first conducted in 1938 by Classius and Dickel, who showed that high-efficiency component separation could be achieved by creating a convective flow using a temperature gradient [5]. Over an extended period, Clausius et al. separated small quantities of the stable isotopes of chlorine, nitrogen, oxygen, carbon, and successfully separated many isotopes, including ¹²⁶Xe (33.6%), ¹³²Xe (17.5%), ²¹Ne (99.6%), and ²⁰Ne (99.95%) [3, 5]. Other researchers have reported successful stable isotope separation using thermal diffusion columns. In 1960, Saxena and Watson optimized the physical parameters of a thermal diffusion column [6]. In 1968, Rutherford et al. designed a column sample to compare the theoretical and experimental results of Ar, CH₄, and Xe isotope separation and utilized a six-stage cascade to separate ¹²⁴Xe using 24 columns [2]. Subsequently, in 1969, they experimentally and theoretically studied the isotopic separation of O₂ and N₂ in a thermal diffusion column [7], and in 1970, they examined the separation of ²¹Ne column and cascade [8]. In 1972, Roger and Rutherford designed a cascade to separate ⁸⁵Kr using a thermal diffusion column [9], and in 1979, Charles designed and investigated the ideal cascade of thermal diffusion [10]. In 1981, Rutherford et al. analyzed the behavior of a cascade by changing the time required to separate multicomponent isotopes [11]. In 1982, Zieger et al. optimized intermediate isotope enrichment in multicomponent systems [12]. In 2003, Vasaru studied the separation of ¹³C in a seven-stage cascade using 19 thermal diffusion columns, and investigated the enrichment of ²⁰Ne and ²²Ne in an eight-stage cascade comprising eight columns and ⁸⁶Kr and ⁸⁷Kr in a 15-stage cascade of 35 columns [13]. From 1987 to 1997, Yamamoto et al. optimized parameters such as feed and pressure in a thermal diffusion column using analytical and experimental solutions for various noble gases [4, 14-17]. In 2020, Saiki et al. investigated hydrogen isotope separation via thermal diffusion in a two-dimensional cascade of numerous small cavities [18]. Furthermore, from 2014 to 2021, many researchers have studied numerical theories for multicomponent systems [19-22].

The thermal diffusion cascade was investigated in this study because of its application in noble gas isotope separation. The most notable approach to this problem involves determining the separating elements in a square cascade. Because of the importance of square cascades in separating multicomponent isotopes, determining the number of columns is necessary. To the best of our knowledge, although a considerable amount of research has been conducted, it has not been published. In the aforementioned studies, cascades with specific structures were used for isotope separation, and the number of columns for specific isotope separation was not specified.

The minimum total flow rate is the main criterion for determining the cascade stages and separation units, particularly when using analytical approaches [23]. This is a typical and efficient approach. In situations in which the separation factor is high, determining the number of separation units using a numerical method is more accurate than using an analytical method. When the difference between the analytical and numerical results is large, previous studies have suggested using the numerical method for determining the number of separation units. Notably, the error also increases when the number of isotopes increases. Therefore, numerical techniques are essential in cascade design [20].

This study investigated the main concepts involved in designing cascades for separating stable isotopes using a thermal diffusion column. This procedure was used as an example for Ne isotope separation, and the match abundance ratio cascade (R cascade) was used to estimate the number of columns required. We modified the numerical algorithm presented in Ref. [20] and compared the analytical and numerical results, after which the required parameters were determined to evaluate the number of columns in the R and square cascades. Owing to its flexibility, a square cascade was then designed using the calculated columns. By changing the stage cut, feed point, and cascade feed rate, the best parameters for the square cascade were determined according to the cascade and column separation powers. To this end, cascade modeling and design methods for thermal diffusion columns were identified, and calculation codes for the cascade design were prepared.

Research Theories

The common R and square cascades are shown in Fig. 1. The feed mixture (F) with composition z_i,F enters a cascade. The cascade then delivers the heavy product (W) and light product (P) with compositions x_i,w and y_i,p, respectively. Stage s receives the feed flow rate Z_s and delivers the up and down flow rates M_s and N_s, respectively [19-21].

Fig. 1

(a) Square, and (b) R cascades

The governing equations are listed in Table 1. The total number of stages is S, and the number of components in the feed mixture is N_c. In square cascades, the number of separation units and feed flow rate Z_s are constant in each stage (Eq. (1)). In contrast to the square cascade, the rate of return flows ε and ε' in the R cascade is zero [21].

Governing equations in the square and R cascades

Equation	Note	No.
$Z_{1} = Z_{2} = \dots = Z_{S} = Z$	Only in square cascade	(1)
$Z_{s} = M_{s - 1} + N_{s + 1}, s \neq s_{f}, s \neq S, s \neq 1$		(2)
$Z_{s} = M_{s - 1} + N_{s + 1} + F, s = s_{f}$		(3)
$Z_{s}_{} = M_{s - 1} + ε^{'}, s = S$	$ε^{'} = 0, in R cascade$	(4)
$Z_{s} = N_{2} + ε, s = 1$	$ε = 0, in R cascade$	(5)
$M_{s}^{} = Z θ_{s}$		(6)
$N_{s}^{} = Z (1 - θ_{s})$		(7)
$Z_{s} z_{i, s} = M_{s} y_{i, s} + N_{s} x_{i, s}$		(8)
$Z_{s} z_{i, s} = M_{s - 1} y_{i, s - 1} + N_{s + 1} x_{i, s + 1} + F z_{i, F}, s = s_{f}$		(9)
$Z_{s} z_{i, s} = M_{s - 1} y_{i, s - 1} + N_{s + 1} x_{i, s + 1}, s \neq s_{f}$		(10)
$Z_{s} z_{i, s} - M_{s - 1} y_{i, s - 1} - ε^{'} y_{i, P} = 0, s = S$	$ε^{'} = 0 in R cascade$	(11)
$Z_{s} z_{i, s} - N_{s + 1} x_{i, s + 1} - ε x_{i, w} = 0, s = 1$	$ε^{} = 0 in R cascade$	(12)
$α_{i j, s} = \frac{(y_{i, s} / y_{j, s})}{(x_{i, s} / x_{j, s})} = α_{0, s}^{(M_{j} - M_{i})}, (i = j - 1, j = 2, \dots, N_{c})$		(13)
$\sum_{i = 1}^{N_{c}} z_{i, s} = \sum_{i = 1}^{N_{c}} y_{i, s} = \sum_{i = 1}^{N_{c}} x_{i, s} = 1$		(14)
${R^{″}}_{(k_{1}, k_{2}), s + 1} = R_{(k_{1}, k_{2}), s} = {R^{'}}_{(k_{1}, k_{2}), s - 1}$	Only in R cascade	(15)
$R_{(k_{1}, k_{2}), s} = \frac{z_{k 1, n}}{z_{k 2, n}}, {R^{'}}_{(k_{1}, k_{2}), s} = \frac{y_{k 1, n}}{y_{k 2, n}}, {R^{″}}_{(k_{1}, k_{2}), s + 1} = \frac{x_{k 1, n}}{x_{k 2, n}}$		(16)

Eqs. 1–7 show the gas flow rates in cascade pipes [20-22]. Eqs. 2–5 show the mass flow balance at the mixing points from 1 to S. In these equations, the cut in each stage is θ_s. Eqs. 6 and 7 can be obtained using the mass flow balance in each stage and the definition of cut in stages ( $θ_{s} = \frac{M_{s}}{Z_{s}}$ ).

Equations 8–15 show the concentrations of isotopes in the cascade stages. Equation 8 shows the mass conservation of the isotope i^th in the stage s^th, and Eqs. 9–12 shows the mass balance of the component i^th at the mixing points. $α_{ij, s}$ is the separation factor between the isotopes i^th and j^th at the s^th stage (Eq. 13), and $α_{0}$ is the overall separation factor for the unit mass difference. M_i and M_j refer to the molar weights of the i^th and j^th isotopes, respectively. Finally, the concentrations must satisfy the constraints in Eq. 14 [20-22], that shows the sum of the concentrations of each component at the mixing points should be equal to 1.

In the R cascade, the match abundance ratio for the two isotopes k₁ and k₂ is defined using Eq. 15 [20, 21]. The parameters R, R $",$ and R $""$ are the abundance ratios in the feed, light, and heavy flows, respectively (Eq. 16). Previous studies have shown that the R cascade is suitable for thermal diffusion cascade design, such as the gas centrifuge cascade [10, 24].

The analytical approach to the R cascade is as follows [21, 23, 27]: $y_{i, P}^{} = \frac{ξ_{i}^{s_{f}} - 1}{ξ_{i}^{N_{F}} - ξ_{i}^{- (S - s_{f} + 1)}} z_{i, F}^{} / \sum_{j = 1}^{N_{c}} \frac{ξ_{j}^{s_{f}} - 1}{ξ_{j}^{N_{F}} - ξ_{j}^{- (S - s_{f} + 1)}} z_{j, F}^{}$ (17) $x_{i, W}^{} = \frac{1 - ξ_{i}^{- (S - s_{f} + 1)}}{ξ_{i}^{s_{f}} - ξ_{i}^{- (S - s_{f} + 1)}} z_{i, F}^{} / \sum_{j = 1}^{N_{C}} \frac{1 - ξ_{j}^{- (S - s_{f} + 1)}}{ξ_{j}^{s_{f}} - ξ_{j}^{- (S - s_{f} + 1)}} z_{j, F}^{}$ (18) such that $ξ_{i} = γ_{0}^{M^{*} - M_{i}}, and M^{*} = \frac{M_{k 1} + M_{k 2}}{2} .$ (19)

The number of stages (S) and feed stage ( $s_{f}$ ) with the given external isotope concentrations ( $y_{k, P}^{}, x_{k, W}^{}$ ) can be obtained by solving Eqs. (17) and (18). The feed flow rate at each stage (Z_s) can be obtained using Eq. 20. $Z_{s}^{} = {\begin{matrix} W \sum_{i = 1}^{N_{c}} x_{i, W}^{} \frac{ξ_{i} + 1}{ξ_{i} - 1} (ξ_{i}^{s} - 1), (1 \leq s \leq s_{f}) \\ P \sum_{i = 1}^{N_{c}} y_{i, P}^{} \frac{ξ_{i} + 1}{ξ_{i} - 1} (1 - ξ_{i}^{s - S - 1}), (1 \leq s \leq s_{f}) \end{matrix}$ (20)

Based on Eq. 20, the parameter $\sum_{s = 1}^{S} Z_{s}$ depends on $ξ_{i}$ , and according to Eq. 19 $ξ_{i}$ is a function of M^*. Hence, the cascade parameters for isotopes k₁ and k₂ can be determined. If M^* is used as an optimization variable and has a value between the actual components ( $M_{1} \leq M * \leq M_{N c}$ ), there is a continuous set of cascades.

2.1

Separation power

The separative power was calculated using Eq. (21). $Δ U = P V_{(y_{1, P}, y_{2, P}, \dots y_{N c, P})} + W V_{(x_{1, W}, x_{2, W}, \dots x_{N c, W})} - F V_{(z_{1, F}, z_{2, F}, \dots z_{N c, F})}$ (21) where F, P, and W are the flows of the feed, light, and heavy products, respectively, and $y_{i, P}$ , $x_{i, W}$ , and $z_{i, F}$ , are the corresponding isotopic concentrations. Moreover, $V (c_{1,} c_{2, \dots} c_{N c})$ is the separation potential of the multicomponent system.

There is no standard interpretation of the separation potential of multicomponent systems despite the many approaches devised for its evaluation. Thus far, relationships have been introduced for separation potentials based on different assumptions. If a small number of impurity isotopes are present in the feed and the two components are the main isotopes, the most suitable potential relation is [25]. $\begin{matrix} V (c_{1,} c_{2, \dots} c_{N c}) = {(\frac{Δ M_{N c 1}}{Δ M_{k_{1} k_{2}}})}^{2} (c_{k_{2}} - c_{k_{1}}) ln \frac{c_{k_{2}}}{c_{k_{1}}} \\ + \sum_{\begin{array}{l} k = 1 \\ k \neq k_{1}, k_{2} \end{array}}^{N_{c}} {(\frac{Δ M_{N c 1}}{Δ M_{k_{1} k}})}^{2} c_{k} ln \frac{c_{k}}{c_{k_{1}}} \end{matrix}$ (22) where $c_{i}$ is the isotopic concentration of the product, waste, or feed streams, and k₁ and k₂ are the indices of the main isotopes.

Calculation Procedure

First, different R cascades and their effects on the separation of ²⁰Ne and ²²Ne were studied, and a cascade design was performed. Consequently, the required number of columns was estimated using the R cascade. Subsequently, a square cascade with specific columns was determined and optimized for the separation of Ne isotopes.

In this regard, the computational code "RSQ_CASCADE" was developed, and the required number of columns was determined to be $y_{{}^{20}N e, P} > 99 %$ , and $x_{{}^{22}N e, W} > 99 %$ by changing the number of R-cascade stages from 10 to 20. The concentrations of Ne stable isotopes in the feed stream were $z_{1, F ({}^{20}N e)}$ =0.9048, $z_{2, F ({}^{21}N e)}$ =0.0027, and $z_{3, F ({}^{22}N e)}$ =0.0925. The components $k_{1}$ and $k_{2}$ in the R cascade were ${}^{20}N e$ and ${}^{22}N e$ , respectively.

The parameters that determine the number of stages were separation power, number of columns, and relative total flow (product flow ratio to total feed flow). In this study, the unit separation factor and feed flow rate of the thermal column were 50–60 cm³/s and 3, respectively. These parameters were obtained from experimental data [15].

The numerical procedure for the RSQ code is shown in Fig. 2. This code consists of two separate sections for modeling the R cascade (Fig. 2a) and the square cascade (Fig. 2b).

Fig. 2

Numerical procedure for (a) R cascade and (b) square cascade

3.1

R cascade

The differences between the numerical and analytical total flow values for various alphas are listed in Table 2. When the separation factor increases, the difference between analytical and numerical methods for determining the value of $\sum Z$ increases. Therefore, we used numerical methods to calculate the separation elements. Notably, these results are for Ne isotopes, and additional changes occur as the separation factor increases for other isotopes. For example, $\sum Z / P$ for the separation of ⁷⁸Kr using an analytical solution was 2791.2 [27], which is equal to 2863.63 for the numerical calculation and leads to a 2.5% error. This result corresponds to a unit separation factor of 1.1. Therefore, the number of separation units was numerically determined in this study.

Calculated total flow rate values for the numerical and analytical methods in different alphas

α₀	1.05	1.09	1.5	2	2.5	3
(∑Z), Analytical	30678.25	22530.74	10963.15	7044.151	5206.569	4361.049
(∑Z), Numerical	30684.34	22545.11	11112.93	7323.95	5565.853	4791.102
error, %	0.02	0.06	1.37	3.97	6.9	9.86

In the numerical approach, proportional to the initial values of the θs, the flow rates and the concentration of isotopes were determined based on Eqs. (2)–(7) [20]. Subsequently, nonlinear Eqs. (8)–(15) were solved using iterative methods. Therefore, the flow rates and isotope concentrations at each stage were determined. The trust-region method was used to model the R cascade and calculate the nonlinear equations. The solution method is very sensitive to preliminary guesses and often fails to find a physical solution even when the procedure converges. The continuation technique can obtain a solution with an arbitrary initial guess [20]. The R-cascade parameters were determined according to the proposed procedure using a continuation technique. For the mixing points, the value of $α_{0}$ and ${R^{'}}_{k_{1} k_{2}, s - 1}^{(i t r - 1)}$ in the algorithm of continuation technique was applied according to Eq. 23. The itr counts the number of iterations, and $0 < τ < 1$ controls the speed of this transition. This allows the value of ${R^{″}}_{k_{1} k_{2}, s + 1}^{(i t r - 1)}$ to gradually approach that of ${R^{'}}_{k_{1} k_{2}, s - 1}^{(i t r - 1)}$ . The error was calculated using Eq. 24. ${R^{'}}_{k_{1} k_{2}, s - 1}^{i t r} - {R^{″}}_{k_{1} k_{2}, s + 1}^{i t r} = [{R^{'}}_{k_{1} k_{2}, s - 1}^{(i t r - 1)} - {R^{″}}_{k_{1} k_{2}, s + 1}^{(i t r - 1)}] (1 - e^{- τ . i t r})$ (23) $e r r o r = m a x_{s ϵ S} (\frac{| z_{s}^{(i t r)} - z_{s}^{(itr - 1)} |}{z_{s}^{(n_{i})}}, \frac{| M_{s}^{(n_{i})} - M_{s}^{(n_{i} - 1)} |}{M_{s}^{(n_{i})}}, \frac{| N_{s}^{(n_{i})} - N_{s}^{(n_{i} - 1)} |}{N_{s}^{(n_{i})}})$ (24)

The R cascade was simulated as described in Ref. [20]; however, only the solving of nonlinear equations, and the equation for the continuation technique in it is based on the Newton method and the definition of enrichment factor. A comparison of the two algorithms showed that the convergence time and number of loops required were reduced. Replacing the Newton method with a trust region and changing the homotopy equation reduced the number of required iterations from 1322 to 47. If only the trust region was used without changing the homotopy equation, the iteration convergence was reduced from 1322 to 687.

3.2

Square cascade

As mentioned before, in the simulation of the square cascade, Z_s is known and provided as the input. Moreover, the parameters M_s, N_s, ε, and ε' are unknown. Therefore, there are 2S + 2 unknown parameters. The independent equations, which include the equations of flow at the stages and mixing points, are 2S. In solving these equations, two parameters must be determined. In a square cascade, the cut is the most significant operational parameter. Having the cut in the first stage and the cascade cut, the flow rates are obtained using Eqs. (1)–(7) [21]. Additionally, using nonlinear Eqs. (8)–(13), the concentration of isotopes is determined using the Q-iteration method explained briefly in this study.

By extracting the parameter $z_{i, s}$ from Eq. (8), and replacing it with Eqs. (9) and (10), Eqs. (25)–(28) are obtained in terms of $x_{i, s} .$ $(M_{1} q + N_{1} - ε) x_{i, 1} - N_{2} x_{i, 2} = δ (s, S_{F}) F z_{i, F}$ (25) $(N_{s} + M_{s} q) x_{i, s} - M_{s - 1} q x_{i, s - 1} - N_{s + 1} x_{i, s - 1} = δ (s, S_{F}) F z_{i, F}$ (26) $(M_{P} q + N_{P} - ε' q) x_{i, P} - M_{P - 1} q x_{i, P - 1} = δ (s, S_{F}) F z_{i, F}$ (27) $δ (s, S_{F}) = {\begin{matrix} 0, (s \neq S_{F}) \\ 1, (s = S_{F}) \end{matrix}$ (28)

In these relations, $(M_{1} q + N_{1} - ε) x_{i, 1} - N_{2} x_{i, 2} = δ (s, S_{F}) F z_{i, F}$ , and considering this definition for Eq. (13), Eq. (29) is obtained as follows: $q_{i, s} = q_{j, s} α_{0, s}^{(M_{j} - M_{i})} .$ (29)

First, the value of $q_{j, s}$ is estimated for all stages; subsequently, it is calculated in all stages using Eq. (29). Therefore, the nonlinear Eqs. (25)–(27) are linearized, and the value of $x_{i, s}$ is calculated for all stages. Subsequently, the value of $y_{i, s}$ is calculated from the definition of $q_{i, s}$ . Additionally, $z_{i, s}$ is determined by considering the conservation equations for the isotope i^th in the stage s^th. By applying Eq. 14, the error value is checked at each calculation step, $q_{i, s}$ is modified, and the calculations are repeated. More detailed explanations are provided in Ref. [26]. Although the principles remain the same, the equations differ.

Results and Discussion

4.1

Results for R cascade and column number estimation

The number of stages in the R cascade was changed from 10 to 20. Figure 3a shows that an increase in the number of cascade stages decreased the column separation power. This parameter varied from 27.48 to 26.95 (sccm). The column separation power for the 10-stage cascade was slightly higher than other R cascades.

Fig. 3

(a) column separation power, (b) column number, (c) light relative total flow, and (d) heavy relative total flow according to the number of stages

As shown in Fig. 3b, the number of columns increased from 85 in the 10-stage cascade to 105 in the 20-stage cascade. Despite the increase in the number of columns and stages, the column separation power did not change considerably. Moreover, the flow rates of light and heavy products largely remained constant. However, as shown in Fig. 3c and 3d, the relative total flow rate decreased as the number of stages increased. The light relative total flow is $P r o d u c t / \sum_{s = 1}^{S} Z_{s}$ , and the heavy relative total flow is $W a s t e / \sum_{s = 1}^{S} Z_{s}$ . Notably, the waste in this study is related to the flow rate of the ²²Ne product stream that withdraws from Stage one. This phenomenon indicates that the total flow rate increases as the number of stages increased. Therefore, the best choice is a cascade with 10 stages, fewer columns, and better separation power than other R cascades. Moreover, reducing the number of columns decreases operating costs while improving separation efficiency. Therefore, because of the same light and heavy products, the same enrichment, and greater column separation power the R cascade with 85 columns was selected.

4.2

Results for square cascade

A square cascade with 85 columns was investigated based on R cascade studies. The square cascade comprised 17 stages, with five columns in each stage. By changing the $θ_{s}$ , $s_{f}$ , and F, the square cascade parameters were determined according to the cascade separation power and column separation power. Owing to their flexibility, these parameters in the square cascade can change each time. Table 3 presents the main parameters of the optimal squared cascade.

Main square cascade parameters

P/F	θN	θ₂	θ₁	Feed (sccm)	s_f	Z/F	Opt. Cascade
0.9076	0.5907	0.4908	0.5000	25	11	10	1
0.9074	0.6567	0.4946	0.4900	41.67	16	6	2
0.9075	0.6900	0.4915	0.4900	50	16	5	3
0.9075	0.7400	0.4869	0.4900	62.5	10	4	4
0.9079	0.8233	0.4792	0.4900	83.33	10	3	5

As mentioned above, the optimal feed rate value was 50–60 sccm. Therefore, the input flow rate for all stages was 250–300 sccm. Thus, changes could only be observed regarding $s_{f}$ , $s_{f}$ , and $F_{}$ increasing the enrichment of both Ne-20 and Ne-22 to 99%. A suitable value for the stage cut was considered to be close to 0.5.

As the feed rate increased, the amounts of light and heavy products in the cascade increased. Moreover, with an increase in the cascade feed flow rate, the column separation power reached a maximum of 27.8 (cm³/min) at 50 sccm and then decreased (Fig. 4a). In contrast, as shown in Fig. 4b, the cascade separation power increased as cascade feed increased. Therefore, the higher the feed rate, the better the yield and separation capacity.

Fig. 4

Variation　in (a) column separation power and (b) cascade separation power according to the square cascade feed rate

In the cascade design, a pressure-regulating valve must be used at each stage to separate the desired isotopes according to the operating conditions of the column. Thus, adjusting the cut to a certain value is possible, and the separation can be performed. If the design constraints can be overcome, cascade 5 is more suitable because of its greater cascade separation power and higher product (with an increasing cascade feed rate, the cut deviation rate increases from 0.5). Otherwise, cascade 3 seems more suitable than other R cascades because of its higher column separation power.

Nevertheless, the correct and optimal use of the cascade in isotope separation depends on the preparation of the operational documentation of the column and the effects of the parameters on the separation factor.

Conclusion

In this study, a thermal diffusion column cascade was investigated to separate Ne stable isotopes. Owing to the difference between the numerical and analytical total flow values for various alphas, numerical method to calculate the separation elements in the R cascade. Therefore, two calculation codes were developed in MATLAB. Subsequently, the thermal column separation power and the light and heavy relative total flow rates were studied in the R cascade from stages 10–20 to determine the required column numbers. The results suggest a 10-stage R cascade with 85 columns to separate stable Ne isotopes. The R cascade is better than the square cascade for separating stable Ne isotopes because of its proximity to ideal conditions. Because different cuts are required in all stages, single-column tests should be performed according to these requirements.

In the long term, in separating other isotopes of noble gases, a square cascade should be used because of its flexibility. Subsequently, a square cascade was designed, and the calculations resulted in a square cascade with 17 stages and five columns for each stage.

References

[1]

L. Chen, R. Yan, X. Z. Kang, et al.,

Study on the production characteristics of 131I and 90Sr isotopes in a molten salt reactor

. Nucl. Sci. Tech. 32, 33 (2021). doi: 10.1007/s41365-021-00867-1