Next Article in Journal
Degradation Behavior and Lifetime Prediction of Polyurea Anti-Seepage Coating for Concrete Lining in Water Conveyance Tunnels
Next Article in Special Issue
Improving the Electroluminescence Properties of New Chrysene Derivatives with High Color Purity for Deep-Blue OLEDs
Previous Article in Journal
Exploring the Potential Application of an Innovative Post-Weld Finishing Method in Butt-Welded Joints of Stainless Steels and Aluminum Alloys
Previous Article in Special Issue
Study on Influence Factors of H2O2 Generation Efficiency on Both Cathode and Anode in a Diaphragm-Free Bath
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Benchtop Zone Refinement of Simulated Future Spent Nuclear Fuel Pyroprocessing Waste

1
Materials and Engineering Research Institute, Sheffield Hallam University, Sheffield S1 1WB, UK
2
The Henry Royce Institute and Department of Materials Science and Engineering, The University of Sheffield, Sheffield S1 3JD, UK
3
Department of Chemistry, University of Sheffield, Western Bank, Sheffield S10 2TN, UK
4
School of Engineering, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
5
National Nuclear Laboratory, Sellafield, Seascale, Cumbria CA20 1PG, UK
*
Author to whom correspondence should be addressed.
Materials 2024, 17(8), 1781; https://doi.org/10.3390/ma17081781
Submission received: 16 March 2024 / Revised: 5 April 2024 / Accepted: 7 April 2024 / Published: 12 April 2024

Abstract

:
The UK’s adoption of pyroprocessing of spent nuclear fuel as an alternative to the current aqueous processing routes requires a robust scientific underpinning of all relevant processes. One key process is the clean-up of the contaminated salt from the electroreducing and electrorefining processes. A proposed method for this clean-up is zone refining, whereby the tendency of the contaminants to remain in the liquid phase during melting and freezing is exploited to ‘sweep’ the contaminants to one end of the sample. Experiments were performed, utilising off-the-shelf laboratory equipment, to demonstrate the feasibility of zone refining for clean-up of electroreducing and electrorefining wastes. This was successful for the electrorefining simulant samples, with effective segregation coefficient, keff, values, which provide a measure of the degree of separation in the sample, between 0 and 1. Lower values indicate greater separation, with values of as low as 0.542 achieved here, corresponding to a reduction in RECl3 content from 10.0 wt.% to 8.4 wt.% (for 80% salt reuse). Due to difficulties in obtaining a fully homogeneous electroreducing simulant waste, it was not possible to demonstrate the feasibility of zone refining using the current experimental setup. Further research is required to elucidate the correct preparation conditions for production of homogeneous electroreducing waste simulants.

1. Introduction

Nuclear power has remained one of the most significant modes of power generation since its inception in the 1950s [1]. However, one of the key challenges in its continued use is waste management. Countries generally operate either an open or closed fuel cycle [2]. In an open fuel cycle, the nuclear fuel is used within the power plant and then is disposed of once it is spent. Conversely, in a closed fuel cycle, once the fuel is spent, it is removed from the power plant and reprocessed, with some or all of the fuel being reused to generate power [3]. Most countries which reprocess their spent nuclear fuel (SNF) utilise aqueous processes, e.g., PUREX (plutonium and uranium recovery by extraction), where solvent extraction combines the use of acids and organic complexants to dissolve the SNF and separate the actinides from the fission products [4]. However, there are some disadvantages to aqueous processes, such as their vulnerability to proliferation and radiation damage of the solvents [5].
Molten-salt-based electrochemical reprocessing (also known as pyroprocessing) solves these issues, as it is inherently proliferation-resistant and the inorganic electrolytes are not vulnerable to radiation damage [6]. The US developed an electrochemical process for the reprocessing of spent metallic fuel from sodium-cooled fast reactors, such as EBR-II [7]. Its application is being expanded to process spent fuels from oxide fuel reactors [6]. South Korea also utilises a molten-salt-based process, with the Korea Atomic Energy Institute (KAERI) beginning development of the process in 1997 [8]. Generally, the primary stage of pyroprocessing is electrorefining, where a metallic feedstock is dissolved in LiCl-KCl eutectic (LKE; 44:56 wt%) salt, which melts at 355 °C, and the re-useable actinides (mainly U and Pu) are deposited at multiple cathodes [4]. Due to the requirement of a metallic feedstock for electrorefining, an additional stage, known as electroreduction, is required for processing of oxide fuels [8]. In electroreduction, the oxide fuels are dissolved into a LiCl molten salt, with the oxygen transported to the anode and released [8]. In both processes, highly radioactive fission products build up in the salt as the useful actinides are electrochemically extracted, leaving the fission products to accumulate, requiring further treatment. These salts become rich in alkalis (Cs, Rb), alkaline earths (Ba, Sr), and rare earths (La, Ce, Pr…) [9]. To operate efficient and successful pyroprocessing of spent nuclear fuel, the LiCl and LKE molten salts used in electroreducing and electrorefining need to be periodically treated post-use to remove fission products and other impurities [4]. There are two primary reasons for this: (1) removal of the fission products means that a significant volume of the salts can be reused, increasing the efficiency of the system; and (2) chloride wastes tend to have low solubility in most common waste forms—reducing the volume of chloride waste through treatment will drastically reduce final waste volumes [4].
Recent research efforts, primarily out of South Korea and the USA, have worked to establish the applicability of zone refining, and other melt crystallisation techniques, for the treatment of pyroprocessing salt wastes [3,4]. The melt crystallisation techniques investigated include zone freezing [10], layer crystallisation (or coldfinger crystallisation) [10,11,12], the Czochralski method [13], and zone refining [6,13]. In general, this research has focused on the separation of CsCl and SrCl2 from LiCl, i.e., refinement of electroreducing waste salt where Cs-137 and Sr-90 will be the main fission product contaminants. This is due to the difficulty of separating CsCl and SrCl2 from LiCl using common techniques such as precipitation [10]. Cho et al. demonstrated separation efficiencies for Cs and Sr of >90% for both the zone freezing and layer crystallisation methods, provided the process rate (crucible rising velocity for zone freezing, crystal growth rate for layer crystallisation) was suitably slow (1.7 mm/h and <5 g/min, respectively) [10]. In general, the separation efficiency was similar for Cs and Sr. Cho et al. confirmed the efficacy of the layer crystallisation method, observing >90% separation efficiencies for Cs and Sr up to a crystal growth rate of 5 g/min, although the results suggest a greater efficiency of Cs separation compared to Sr separation [11]. Versey et al. also observed effective separation of Cs from LiCl using layer crystallisation, and they demonstrated that the efficiency can be maximised through increasing the cooling gas flow rate and crystal growth time [12]. The Czochralski method was investigated by Lee et al., with effective separation of Cs and Sr from LiCl observed, particularly at lower extraction speeds (2.5 mm/h) [13]. However, the calculated segregation coefficient (k) values suggest that Cs was separated more efficiently than Sr using this method. Shim et al. investigated zone refining for separation of Cs and Sr from both LiCl and LKE [14,15]. The efficiency of separation was inversely proportional to the freezing rate, and proportional to the number of passes. Cs was separated more effectively than Sr, but both elements were more difficult to separate from the LKE than from the LiCl salt. It was also found that the elements were more difficult to separate in LiCl-KCl-CsCl-SrCl2 compared to LiCl-KCl-CsCl or LiCl-KCl-SrCl2.
One key issue for melt crystallisation techniques for processing of pyroprocessing wastes is their relatively slow rates. For a recovery yield of 60% LiCl from an LiCl-CsCl-SrCl2 waste, zone refining would take around 80 h, zone freezing would take around 102 h, and the Czochralski method around 52 h [14]. However, the advantage of the zone-refining technique is that multiple heating zones can be implemented, which increases the throughput by reducing the number of passes required to refine the waste salt [14]. Whilst there have been several investigations into melt crystallisation techniques for electroreducing waste salt [6,13], the same cannot be said for its electrorefining counterpart, where various rare-earth radionuclides will be the main impurities. This is likely due to the efficacy of other separation methods, such as reactive precipitation and reactive distillation [4]. However, the study of zone refinement for electrorefining waste treatment is of value for comparison to other methods. One proposed method of molten salt waste treatment is zone refining [8]. This technique was developed by W. G. Pfann at Bell Labs in the early 1950s [16], and can be used on any system with a non-zero segregation coefficient, k, as calculated by
k = C S C L
where CS is the concentration of the impurity/impurities in the solid salt, and CL is the concentration of the impurity/impurities in the molten salt. A schematic of this method can be seen in Figure 1. A moving heater generates a moving molten zone, travelling from one end of the sample to the other. When k < 1, as the salt behind the moving molten zone solidifies, the impurities remain within the molten zone at a higher concentration than in the recrystallised salt. Thus, as the molten zone travels along the sample it ‘sweeps’ the impurities to the end of the sample [14]. This method has most commonly been used in the preparation of high-purity materials for the production of semi-conductors [17,18,19,20] and scintillators [21,22].
The theoretical distribution of an impurity in a 1D zone-refining system can be calculated through several different equations based on the region within the sample. For a 1D system in equilibrium with sample length, L; position, x; normalised position X = x/L; molten zone length, z; and normalised molten zone length, Z = z/L, the sample can be split into three regions: 0 ≤ X ≤ 1 − Z, 1 − ZX < 1, and X = 1 [15]. For a single-pass system, the impurity distribution in the 0 ≤ X ≤ 1 − Z region can be calculated by:
C S X = C 0 1 1 k e x p k X Z
where CS(X) is the concentration of the impurity at position X and C0 is the initial impurity concentration [16]. In the 1 − ZX < 1 region, no new material is added to the molten zone and the impurity distribution is calculated by [15]
C S X = C 0 1 1 k e x p k 1 Z Z 1 X 1 Z Z k 1
In the region X = 1, for k < 1, CS(X) approaches infinity; for k = 1, CS(X) = C0; and for k > 1, CS(X) = 0.
However, realistic systems are not in equilibrium, and so the segregation coefficient, k, must be replaced with an effective segregation coefficient, keff. This can be calculated through the Burton–Prim–Slichter (BPS) theory [23,24]:
k e f f = k k + 1 k e x p V δ D
where V is the velocity of the molten zone; δ is the diffusion layer thickness at the solidification interface; and D is the impurity diffusion coefficient in the liquid. Thus, the separation efficiency of zone refining is positively correlated with increasing impurity diffusion coefficient and decreasing molten zone velocity and diffusion layer thickness.
For a multi-pass, 1D, non-equilibrium system, the sample is broken up into four regions: X = 0, 0 < X < 1 – Z, 1 − ZX < 1, and X = 1. The impurity concentration at X = 0 after n passes of the heater, CSn(X), is given by
C S n = k e f f d x Z q = 0 M 1 1 C S q   d x n 1
where dx is an individual element of length; M1 is the control volume, equal in length to Z; and q is the index of the dx elements. In the 0 < X < 1 − Z region, the impurity distribution can be calculated by [25]
C S X n = C S X d x n + k d x Z C S X + Z d x n 1 C S X d x n
In this region, the molten zone length is constant and re-solidification of the salt occurs on the left-hand side of the molten zone. In the last zone length section, 1 − ZX < 1, the heater begins to travel past the end of the sample and so no further material is added to the molten zone as re-solidification continues to occur on the left-hand side of the molten zone. In this region, the impurity distribution is given by [26]
C S X n = 1 + 1 k 1 X d X C S X d x n
where dX is the displacement between each calculated data point. The impurity concentration where the molten zone fully passes the sample (X = 1) is calculated by [26]
C S 1 n = 1 d X C 0 X = 0 X = 1 d X C S X n
In this work, benchtop apparatus was used to conduct zone-refinement procedures using simulated electrorefiner and electroreducer wastes. These experiments were performed to demonstrate the feasibility of this zone-refining technique for the clean-up of electroreducing and electrorefining wastes. Comparing the experimental data with the theoretical behaviours described through Equations (1)–(8) will illustrate the effective segregation coefficient, keff, as a measure for the degree of separation in the sample that was achieved. By demonstrating the suitability of laboratory-scale instrumentation, more progress can be made towards resolving the issues of appropriate waste management [27] for the nuclear sector, enabling more widespread use of nuclear power stations.

2. Materials and Methods

2.1. Materials and Sample Preparation

Raw materials of appropriate purity were obtained from Better Equipped, Staffordshire, UK (BE) and Alfa Aesar, Ward Hill, MA, USA (AA) (see Appendix A for details). All chemicals were kept in their original, sealed containers before use. Upon being opened, the containers were then stored within vacuum desiccators (with silica gel desiccant) to reduce hydration of the chlorides, which are generally hygroscopic. All weighing and mixing of chlorides took place within an Inert Corp, Amesbury, MA, USA, PureLab HE glovebox. To produce a salt sample for testing, the appropriate chemicals were placed inside the glovebox and the correct amounts were weighed out into an alumina BS111 crucible boat from Almath Crucibles Ltd., Newmarket, UK. The crucible boat was removed from the glovebox and placed inside the tube furnace (see Section 3.2) at a predetermined position. The furnace was heated at 2 °C/min to a temperature above the melting temperatures of the waste salt simulants, held for several hours (depending on sample), and then cooled back to room temperature at 2 °C/min, all while under a steady, low flow of N2.
Table 1 shows the three simulant salt waste compositions derived through consultation with NNL, based on their published example compositions [9]. Electroreducer waste is generated during the reduction of spent oxide fuels (UO2, MOx) to produce a metallic feed for electrorefining. The ‘No-Cs’ and ‘Full-Cs’ variants are based on whether a high-temperature decladding process is or is not used, respectively, which may result in the volatilisation of Cs [28]. The electrorefiner waste is an estimate of the salt composition after the residual actinides have been removed. Zone-refining experiments were performed on the electrorefiner waste simulant and the Full-Cs electroreducer waste simulant.

2.2. Experimental Setup

Zone-refining experiments were performed in an Elite Furnaces, Market Harborough, UK, 12/100/750 tube furnace with a fused quartz tube and flowing nitrogen. Due to the layout of the heating elements within the furnace (one controllable heating zone in the centre of the furnace), zone-refining experiments were performed in accordance with the schematics shown in Figure 2. The sample boat is placed off-centre in the furnace and during the heating stage, the end of the sample closest to the centre begins to melt, and this melting front continues across the sample until it is entirely molten (molten stage). Once the sample is fully molten, the temperature of the furnace is decreased at a controlled rate, such that a freezing front passes from the end of the sample furthest from the centre of the furnace to the end of the sample closest to the centre of the furnace (freezing stage). This is equivalent to zone refining with a molten zone which covers the entire sample, i.e., Z = 1. The theoretical impurity distributions for Z = 1 experiments are discussed in Section 2.4.
Prior to beginning zone-refining experiments, the temperature profile of the furnace as a function of the controller set temperature and position must be measured. To achieve this, the temperature of the furnace was set to a temperature of 400 °C–700 °C in 50 °C increments, at a heating rate of 2 °C min−1 and allowed to equilibrate for 30 min. The temperature from the centre point to the end of the tube was then measured in 0.5 cm increments using a steel thermocouple (TC Direct, Uxbridge, UK, 310 Stainless Steel) and handheld reader. The measured temperature profiles are shown in Figure 3, and temperatures measured at 0.0 cm are reported in Appendix A, Table A3. The general trend in the profiles is for the measured temperature to decrease as the thermocouple is moved further from the centre of the furnace, as expected. There are some shifts in the data, e.g., at 9.5 cm in the 700 °C profile, which are likely due to a slight change in the position of the thermocouple within the cross-section of the tube. To obtain more consistent data, these shifts were corrected by a constant factor, i.e., where a profile increased instead of decreased, every data point after that was offset by the same value. The profiles were then fit with a second-order polynomial fit of the form a x 2 + b x + c . The parameters for each profile are provided in Table 2. The positions of 10 °C isotherms (640 °C, 630 °C, 620 °C, etc.) for each profile were noted, and the change in position of these isotherms between profiles was calculated. These isotherm position changes were then converted to rates of isotherm position change as a function of the change in furnace set temperature. These rates can then be used to estimate freezing-front velocities through Equation (9):
V = r I × r T
where rI is the rate of isotherm position change in cm/°C; and rT is the set cooling rate of the furnace in °C/s. The minimum value of rT for the furnace is 0.0017 °C/s (0.1 °C/min) and the highest recommended value is 0.0333 °C/s (2.0 °C/min) in order to avoid thermal shock of the alumina sample boats. The calculated average rates of isotherm position change vary between 0.103 cm/°C and 0.168 cm/°C. The estimated range of possible V values is 0.00018 cm/s–0.00559 cm/s (0.108–3.354 mm/min), which compares well with the values used previously by Shim et al. (0.25–3 mm/min) [14,15].
For the zone-refining experiments, the samples (produced according to Section 2.1) were placed in an alumina boat, and then, into the tube, with the axial position of the boat along the furnace being recorded. The end caps were then fixed to the fused quartz tube, and the system was purged with N2 for 30 min. A low flow of N2 was continued for the rest of each experiment. Once the system was purged, the set temperature of the furnace was increased at 2 °C min−1 up to the target temperature. Once the set temperature had been reached, it was held for a period of time to allow the temperature to equilibrate. The set temperature was then reduced at a constant rate back to room temperature, such that the sample froze from one end to the other. From the rate of temperature decrease and the measured temperature profiles at different set temperatures, an estimate of the velocity of the freezing front was obtained. Once the furnace reached room temperature, the sample was removed and stored inside a desiccator prior to preparation for analysis.

2.3. Initial Testing

Prior to performing zone refining on samples containing the electrorefining simulant, initial tests were performed to assess the experimental setup. A number of these initial tests were carried out, using LKE with Ba or Cs as a tracer of refinement. Zone refining was performed as described in Section 2.2. The details for each test are shown in Table A2. Once each sample had undergone the zone refining, 1–2 cm sections were taken from each end of the sample using a Top Tech Preciso low-speed saw with diamond-coated abrasive cutting wheel. These sections were then sectioned longitudinally using the saw. The longitudinal sections were mounted in Buehler EpoThin resin, ground to a P1200 SiC paper finish, and then carbon-coated and painted with Silver DAG to ensure sample conductivity. The resin-mounted samples were then analysed by SEM-EDX (FEI, Hillsboro, OR, USA, Quanta 650 SEM and Oxford Instruments, Abingdon, UK, X-Max 80 EDX detector) to determine sample homogeneity and ascertain the degree of refinement achieved.

2.4. Full Simulant Zone-Refining Experiments

2.4.1. Electrorefining Simulant

Full electrorefining simulant samples were prepared as described in Section 2.1, and zone refining was performed in accordance with Section 2.2. For zone refining, each sample was placed with its inner end 5.0 cm from the centre of the furnace. Table 3 shows the target temperature for each sample, the estimated temperatures at the centre and ends of the boat, and the rate of temperature decrease and estimated freezing-front velocity for each sample. Samples were held under vacuum in desiccators upon removal from the tube furnace. Note that the nomenclature “FS” in the sample names refers to the compositions being those of their respective full waste simulants, as given in Table 1, where “LKE-” denotes electrorefining samples and “L-” denotes electroreducing samples.
After removal from the desiccator, the samples were cross-sectioned every 0.5 cm along their length using the low-speed saw. Each section of sample was weighed, and then, dissolved in 5 mL of 1% HNO3 solution made with UHQ Type II water. The dissolved samples were sent to the MIDAS facility at the University of Sheffield for ICP-OES analysis to determine the concentrations of Li, K, Y, La, Ce, Pr, Nd, and Sm as a function of position along the sample. The mg/L concentration data received were normalised to the mass of each section and converted to an equivalent mg per kg of sample value. One section of each sample was retained for SEM-EDX analysis to assess sample homogeneity.

2.4.2. Electroreducing Simulant

The two electroreducer samples were prepared in the same manner as the electrorefining samples. The samples were melted for 5 h at 620 °C before undergoing zone refining (see parameters in Table 4). As with the electrorefining samples, the samples were positioned 5.0 cm from the centre of the tube furnace, held under vacuum upon removal from the furnace, and analysed using ICP-OES (for Li, K, Cs, Rb, Ba, and Sr) and SEM-EDX.

2.5. Impurity Distribution Modelling

As noted in Section 2.2, due to the arrangement of the heating elements within the furnace, the zone-refining experiments were performed with a normalised molten zone length, Z, of 1. This means that the entire molten sample was subject to normal freezing and can therefore be modelled using the normal-freezing equations, as defined by Shim [14,15]. The impurity distribution for a region of normal freezing is given by Equation (3). Substituting Z = 1 into Equation (3) and using keff for a non-equilibrium system gives
C S X = C 0 k e f f 1 X k e f f 1
Thus, the impurity distribution is a function of X, which is independent of the experimental conditions; C0, which is user-controlled and defined by the batching of the salt mixture; and keff. In turn, keff is a function of k, which is dependent on the salt system; V, which is user-controlled through the furnace setup; δ, which can be reduced through improved mixing; and D, which is a fundamental property of each element that is being separated from the base salt (LiCl or LKE). The quantitative effect of each of these parameters can be demonstrated by varying each one in turn and plotting the values returned by Equation (10). Figure A9, Figure A10, Figure A11 and Figure A12 in Appendix A show the effects of k, V, δ, and D, respectively, on the impurity distribution profile of a theoretical 1D system.
Figure A9, Figure A10, Figure A11 and Figure A12 show that CS(X) is most sensitive to variation in k, with a higher value of k increasing segregation, as expected, and less dependent on factors pertaining to a particular salt system or element being separated. This should enable the estimation of k by comparison of experimental and modelled data. The next most influential parameter on CS(X) is V, where a lower value of V is preferable for greater segregation of impurities; however, there must be a compromise between the efficiency of the process and the duration of the process. The sensitivity of CS(X) to δ is lower than for k and V, and as a measure of the degree of mixing within the molten zone it is likely that it will remain constant between different tests within the same experimental framework. The sensitivity of CS(X) to the diffusion coefficient of the impurity was assessed through values for five different rare-earth chlorides (LaCl3, CeCl3, NdCl3, DyCl3, and SmCl3) in LKE taken from the literature [6]. Although the values of D span 1–2 orders of magnitude, there is little change in CS(X), with the relative standard deviation between the five profiles being between 1 and 2%. This suggests that the efficiency of the segregation can be maintained despite varying salt waste compositions.

3. Results

3.1. Initial Testing

The first four test samples utilised BaCl2 as the tracer for refinement. However, there appeared to be issues with the BaCl2 dissolving within the LKE. This was confirmed by SEM-EDX (Figure 4), which showed small, highly concentrated regions of BaCl2 within a matrix of KCl (N.B., Li is not detectable through EDX due to a low Z number). CsCl was therefore used as the tracer of refinement for the last two initial tests. SEM-EDX on the final sample showed some evidence of refinement of the LKE, with the Cs concentration significantly greater closer to the centre of the furnace (Figure 5). The data are taken from spectra collected for 10 min at ×150 magnification from the whole area of the image. Based on this result, it was deemed reasonable to move forward to zone-refining experiments on full simulant samples.

3.2. Electrorefining Simulant Zone Refining

Of the six full electrorefining simulant samples, two suffered from significant volatilisation (LKE-FSb and LKE-FSd, as shown in Figure 6). These samples underwent zone refining with the slowest cooling rates, and therefore, remained at higher temperatures for longer periods of time, which is the likely cause of the increased volatilisation. An SEM-EDX analysis of the residue remaining on the alumina boat from the LKE-FSb sample showed that it consisted of rare-earth (RE) chlorides, with little K present (Figure 7). In addition, O is observed in the data, which might suggest the formation of RE oxychlorides. This could explain why this residue remains; RE oxychlorides are prone to form in the presence of moist air [29], which would occur to a greater extent when cooling over longer periods.
The post-refinement LKE-FSa sample appeared to have a consistent eutectic composition across the full length, as shown by the average K/Li ratio measured by ICP-OES across the sample (0.46 ± 0.01). This suggests that the batched sample was well melted and mixed during the melting and zone-refining runs. The experimentally determined distribution profiles for the REs are shown in Figure 8 and Figure 9. Each experimental profile is fitted with Equation (10) (with keff provided by Equation (4)). To achieve this, the squared residual, R2, for each data point was calculated by
R 2 = E X C X C X 2
where EX is the experimentally determined concentration at normalised position X; and CX is the calculated concentration at X. The Solver function in Microsoft Excel was then utilised to fit Equation (10) to the data, with a focus on minimising the average R2 value. C0 and k were allowed to refine freely, whilst V, D, and δ were kept constant. V was calculated from the temperature gradient and cooling rate (see Section 2.2); the D values for Ce, Nd, Sm, and La were taken from Sridharan et al. [6], whilst those for Y and Pr were estimated to be of the same order of magnitude; and δ was chosen to be 0.001 cm as this represents a well-mixed molten phase, as would be expected from a sample which is fully molten before refinement [14].
The average R2 values for all samples can be found in Table A4 and the fitting parameters for the REs are shown in Table A5. This provides an estimate of keff for each RE. The separation behaviour of the REs varied significantly, with Nd exhibiting the greatest separation (keff = 0.542). Pr, Sm, and La all had keff values between 0.714 and 0.801, but Ce and Y showed only minor separation, with keff values greater than 0.900.
The LKE-FSc sample also showed a consistent K/Li ratio across the sample (0.50 ± 0.01), again suggesting that it was reasonably homogeneous. As with the LKE-Fsa sample, Nd showed the greatest separation (Figure A1). The behaviour of the other elements was more similar, with their keff values falling between 0.737 and 0.852 (Table A6).
The average K/Li ratio for LKE-FSd was 0.56 ± 0.02, suggesting reasonable homogeneity. The data for Nd, Ce, and La are shown in Figure A3, and those for Y, Sm, and Pr are shown in Figure A4. Due to the volatilisation of the sample during zone refining, the mass of sample available per section for ICP-OES was limited. This resulted in many of the Ce data being below the limit of detection for the instrument, and thus, no information can be drawn from the remaining data. The model parameters are shown in Table A8. Nd and Sm had the lowest keff values at 0.647 and 0.600, respectively. Pr and Y had similar keff values (0.724 and 0.789), whilst La had a significantly higher value than the other elements (0.888).
The LKE-FSe sample had an average K/Li ratio of 0.51 ± 0.03, suggesting reasonable consistency of composition across the sample. The experimental data and modelled profiles are shown in Figure A5 and Figure A6, with model parameters displayed in Table A9. Nd, Sm, and Pr exhibited very similar profiles, with calculated keff values between 0.875 and 0.896. La and Y demonstrated flatter profiles with a lower degree of separation; the keff values were 0.943 and 0.961, respectively. It was not possible to obtain a fit of the Ce data, suggesting that little to no refinement took place.
The LKE-FSf sample had an average K/Li ratio of 0.49 ± 0.04, suggesting reasonable consistency of composition across the sample, although the variation in this value with position was greater than for the other samples. The experimental data and modelled profiles for the REs are shown in Figure A7 and Figure A8, with the parameters for the modelled profiles shown in Table A10. La and Ce had the lowest keff values, at 0.691 and 0.723, respectively. Nd, Sm, and Pr exhibited similar separation behaviour, with keff values between 0.825 and 0.868, but Y showed a much lower separation level, with a keff of 0.930.

3.3. Electroreducing Simulant Zone Refining

The appearance of both of the electroreducing samples is consistent post-refinement, shown in Figure 10, with volatilisation at lower levels than for some of the electrorefining samples. The same elemental concentration profile fitting method was applied to the samples as was used for the electrorefining samples. The residual data for each fit can be found in Table A5 in Appendix A.
The experimental data and modelled profiles for the L-FSa and L-FSb samples are shown in Figure 11 and Figure 12, respectively. The residuals for each fit were significantly higher for these samples compared to those for the electrorefining samples (see Table A4 and Table A5). There are also systematic deviations of the modelled profiles from the experimental data. In particular, for the regions X ≤ 0.40 and X > 0.80, the calculated values are consistently greater than the experimental values, whilst for 0.55 < X < 0.80 the calculated values are consistently underestimates when compared to the experimental values. Due to these factors, the keff values generated by fitting with Equation (10) are suspect and cannot reasonably be used for comparison with the calculation of DFi values for refined and waste salt compositions.

4. Discussion

The approach taken in this work regarding the keff values generated by fitting with Equation (10) was found to be unsuitable when considering the electroreducing wastes, as shown in Figure 11 and Figure 12. As such, further discussions will focus on the electrorefining simulant zone-refining efforts.
In the electrorefining simulant zone-refining data, the elemental concentration profiles for the RE elements show that some degree of zone refinement was achieved in all five analysed samples (due to volatilisation, ICP-OES analysis of the LKE-FSb sample was not possible). Using the calculated elemental concentration profiles for each element in each sample, decontamination factors and waste compositions can be calculated. The decontamination factor, DFi, for a given element, i, can be calculated by
D F i = a v e r a g e   c o n c e n t r a t i o n   o f   i   f o r   X R X < 1 a v e r a g e   c o n c e n t r a t i o n   o f   i   f o r   0 X X R
where XR is a predetermined ‘refinement boundary’. For 0 ≤ X < XR, the base salt has been refined to remove some of the impurities, and for XRX < 1, the impurities have been concentrated to produce a waste material. The average concentrations are calculated by evaluating the definite integral across the two intervals (0 ≤ X < XR and XRX < 1) and dividing by the interval width:
a v e r a g e   c o n c e n t r a t i o n   o f   i = 0 X R C S X X R   o r   X R 0.99999 C S X 0.99999 X R
X = 0.99999 is taken to be a reasonable estimate of X = 1, as there is a singularity in CS(X) at X = 1. DFi values for XR = 0.8 and 0.9 are shown in Table 5 and Table 6, respectively. In general, the minimum enrichment of an RE in the ‘waste’ portion of the sample for XR = 0.8 is 8% (Y, LKE-FSe) and the maximum is 286% (Nd, LKE-FSa). The equivalent values for XR = 0.9 are 11% (Y, LKE-FSe) and 360% (Nd, LKE-FSa).
Based on the average elemental concentrations calculated through Equation (13), the composition of the refined salt and the waste salt can be estimated. Li, K, and Cl are assumed to remain in the same proportions throughout each sample. The estimated compositions of the refined and waste salts for XR = 0.8 and 0.9 are shown in Table 7, Table 8, Table 9 and Table 10. Due to the variation in DFi between the REs, their relative proportions vary between the samples, but an estimate of the overall compositional change can be made using the total wt.% of RECl3. The total RECl3 content of the waste salt was between 11.4 and 16.2 wt.% for XR = 0.8, and between 12.1 and 20.3 wt.% for XR = 0.9. These contaminant-enriched waste salts could be sent for immobilisation, or for further processing, e.g., phosphate precipitation of the RE ions from the chloride salt.
The primary parameter varied between the six electrorefining simulant tests was the cooling rate of the furnace, and hence, the velocity of the freezing front, V. Based on Equation (4), we would expect keff to decrease with decreasing V, assuming k, δ, and D remain constant. For the example system described in Figure A9 (k = 0.5, L = 100 cm, Z = 1, C0 = 1000 ppm, δ = 0.01 cm, D = 1.01 × 10−5 cm2/s), the keff value would be expected to decrease from 0.800 to 0.688 as V decreases from 0.0014 to 0.0008 cm/s (the velocity range of the successful experiments). However, when the average keff values for the six elements for each test are plotted against the estimated V, the relationship is unclear (Figure 13). Whilst there may generally be a slight decrease in keff with V, the range of the three keff values at 0.0014 cm/s (0.144) is greater than the difference between keff at 0.0008 cm/s (0.730) and the average of the three values at 0.0014 cm/s (0.833). Further experiments at a range of V values, with duplicate or triplicate data, are required to determine whether the relationship is present.
Local maxima in the distribution profiles around the normalised position 0.6 to 0.8, presented in Figure A7 and Figure A8, are indicative of incomplete refinement due to the single-pass experimental setup used. Further zone-refinement processing and the use of multi-pass modes could improve the final distribution profiles and result in profiles more consistent with the modelled profiles. However, the variation in the refinement behaviour of the REs is of interest as this will affect the efficacy of zone refining and the composition of the waste generated. The REs are chemically similar, with all six elements present in the salt as 3+ ions. However, there is a slight difference in their coordination within their respective chlorides, with Nd, Ce, La, Sm, and Pr nine-coordinated, whilst Y is six-coordinated. Plotting the average keff for each ion against their atomic number divided by ionic radius (Z/IR) appears to provide a general linear trend of increasing separation with increasing Z/IR, although this should be caveated with an R2 value of 0.64 (Figure 14). Whether there is a physical explanation for such a relationship remains unclear and would require further work to elucidate.

5. Conclusions

Zone refinement of electrorefining simulant salt was partially achieved using an experimental setup consisting of off-the-shelf equipment. As far as the authors are aware, data on successful zone refining of electrorefining simulant salts using such equipment have not been published previously in the literature. The zone-refining experiments on electroreducing simulant salt samples were not successful, likely due to the contaminant salts (particularly BaCl2 and SrCl2) not dissolving into the base LiCl salt. This work sets out a framework for future work, whereby experimental data can be assessed through the appropriate theoretical models. Further work is required to identify the optimum processing parameters for the electroreducing simulant salt. The experiments were performed in a ‘normal freezing’ arrangement, which allowed for single-pass experiments only. It is anticipated that by improving on the methods herein, and expanding the methods into multi-pass models, as described by Spim [25], molten salt zone refinement research can be accelerated significantly.
Clean-up of salt wastes is a key pillar for successful implementation of pyroprocessing for used nuclear fuel. It determines salt reuse rates and waste volumes and compositions, and thus, the work in this report is a vital part of understanding a potential process within salt clean-up. Further work is required to elucidate the optimal processing conditions for preparation of homogeneous, anhydrous, oxygen-free salt samples for testing. Of particular interest is the identification of temperatures/times which achieve homogeneity whilst minimising salt volatilisation. This may also necessitate further investment in laboratory-scale salt-handling facilities, such as dry gloveboxes containing milling, sectioning, grinding, and polishing equipment, and atmosphere-controlled furnaces.
The concentration of REEs at the bottom of the electrorefining samples and observation of REE-rich residues left after volatilisation are of interest for further investigation. The former may suggest that density of other effects could potentially be exploited to increase separation by utilising a vertical, rather than horizontal, zone-refining technique. The latter observation could indicate preferential volatilisation of the LKE over the REE chlorides, potentially opening up a new route for processing of this waste. Further work is required to investigate these potential opportunities further.

Author Contributions

Conceptualization, P.A.B., M.D.O., C.A.S., M.T.H. and D.M.; methodology, D.J.B. and P.A.B.; formal analysis, D.J.B. and A.S.; investigation, D.J.B. and A.S.; data curation, D.J.B., W.D. and C.M.; writing—original draft preparation, A.S. and D.J.B.; writing—review and editing, A.S., D.J.B. and P.A.B.; supervision, P.A.B.; project administration, P.A.B.; funding acquisition, P.A.B., M.D.O. and C.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded under the £46m Advanced Fuel Cycle Programme as part of the Department for Business, Energy and Industry Strategy’s (BEIS) £505m Energy Innovation Programme. This research utilised the HADES/MIDAS National Nuclear User Facilities at the University of Sheffield, established with financial support from EPSRC and BEIS, under grant numbers EP/T011424/1. This work was supported by the Henry Royce Institute for Advanced Materials, funded through EPSRC grants EP/R00661X/1, EP/S019367/1, EP/P02470X/1 and EP/P025285/1. Authors AS, MDO, CAS and PAB acknowledge funding from EPSRC under grant EP/X011607/1.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Details of reagents purchased and used in this work.
Table A1. Details of reagents purchased and used in this work.
ChemicalSupplierPurity (%)Notes
LiClBE99Anhydrous
KClBE99.5Anhydrous
SrCl2·6H2OBE98Dehydrated at 350 °C for 60 h before use
BaCl2·2H2OBE99Dehydrated at 350 °C for 60 h before use
YCl3AA99.9
LaCl3AA99.9
CeCl3AA99.5
NdCl3AA99.9
PrCl3AA99.99
SmCl3AA99.5
RbClAA99.8
CsClAA99.9
Table A2. Experimental details for initial tests of zone-refining apparatus.
Table A2. Experimental details for initial tests of zone-refining apparatus.
SampleTracerMelting Temp (°C)Melting Duration (h)Z-R Max Temp (°C)Z-R Max Temp Hold (h)Z-R Cooling Rate (°C/min)Sample End Position (cm)
LKE-Ba22 wt.% BaCl2N/AN/A10000.51.05.0
LKE-Ba2_b2 wt.% BaCl2400250021.05.0
LKE-Ba2_c2 wt.% BaCl2500450010.55.0
LKE-Ba55 wt.% BaCl2500450020.45.0
LKE-Cs55 wt.% CsCl500850020.25.0
LKE-Cs5_b5 wt.% CsCl500850080.25.0
Table A3. Measured furnace temperatures at horizontal position 0.0 cm.
Table A3. Measured furnace temperatures at horizontal position 0.0 cm.
Set Temperature (°C)Measured Temperature (°C)
400353
450385
500446
550495
600545
650593
700642
Table A4. Average R2 values for fitting of Equation (10) to electrorefining simulant data.
Table A4. Average R2 values for fitting of Equation (10) to electrorefining simulant data.
SampleR2 (Nd)R2 (Ce)R2 (La)R2 (Y)R2 (Sm)R2 (Pr)
LKE-FSa0.0040.0110.0020.0200.0360.011
LKE-FSc0.0510.1400.0400.0390.0560.046
LKE-FSd0.036-0.0120.0240.0430.032
LKE-FSe0.029-0.0230.0120.0270.023
LKE-FSf0.0190.0720.0390.0100.0330.018
Table A5. Average R2 values for fitting of Equation (10) to electroreducing simulant data.
Table A5. Average R2 values for fitting of Equation (10) to electroreducing simulant data.
SampleR2 (Ba)R2 (Sr)R2 (Cs)R2 (Rb)
L-FSa0.0770.0830.1150.110
L-FSb0.0570.0570.0710.065
Table A6. Model parameters of LKE-FSa experimental profiles.
Table A6. Model parameters of LKE-FSa experimental profiles.
ElementkV (cm/s)δ (cm)D (cm2/s)keff
Nd0.5270.00140.0012.40 × 10−50.542
La0.7780.00140.0011.01 × 10−50.801
Ce0.9320.00140.0011.28 × 10−50.938
Y0.8970.00140.0013.00 × 10−50.901
Sm0.6900.00140.0011.23 × 10−50.714
Pr0.7090.00140.0012.00 × 10−50.724
Figure A1. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSc sample.
Figure A1. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSc sample.
Materials 17 01781 g0a1
Figure A2. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSc sample compared to modelled profiles.
Figure A2. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSc sample compared to modelled profiles.
Materials 17 01781 g0a2
Table A7. Model parameters of LKE-FSc experimental profiles.
Table A7. Model parameters of LKE-FSc experimental profiles.
ElementkV (cm/s)δ (cm)D (cm2/s)keff
Nd0.6560.00110.0012.40 × 10−50.666
La0.7700.00110.0011.01 × 10−50.789
Ce0.8290.00110.0011.28 × 10−50.841
Y0.730.00110.0013.00 × 10−50.737
Sm0.840.00110.0011.23 × 10−50.852
Pr0.770.00110.0012.00 × 10−50.780
Figure A3. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSd sample compared to modelled profiles.
Figure A3. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSd sample compared to modelled profiles.
Materials 17 01781 g0a3
Figure A4. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSd sample compared to modelled profiles.
Figure A4. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSd sample compared to modelled profiles.
Materials 17 01781 g0a4
Table A8. Model parameters of LKE-FSd experimental profiles.
Table A8. Model parameters of LKE-FSd experimental profiles.
ElementkV (cm/s)δ (cm)D (cm2/s)keff
Nd0.6390.00080.0012.40 × 10−50.647
La0.8800.00080.0011.01 × 10−50.888
Ce-----
Y0.7760.00080.0013.00 × 10−50.789
Sm0.5800.00080.0011.23 × 10−50.600
Pr0.7080.00080.0012.00 × 10−50.724
Figure A5. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSe sample compared to modelled profiles.
Figure A5. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSe sample compared to modelled profiles.
Materials 17 01781 g0a5
Figure A6. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSe sample compared to modelled profiles.
Figure A6. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSe sample compared to modelled profiles.
Materials 17 01781 g0a6
Table A9. Model parameters of LKE-FSe experimental profiles.
Table A9. Model parameters of LKE-FSe experimental profiles.
ElementkV (cm/s)δ (cm)D (cm2/s)keff
Nd0.8680.00140.0012.40 × 10−50.875
La0.9350.00140.0011.01 × 10−50.943
Ce-----
Y0.9590.00140.0013.00 × 10−50.961
Sm0.8790.00140.0011.23 × 10−50.893
Pr0.8900.00140.0012.00 × 10−50.896
Figure A7. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSf sample compared to modelled profiles.
Figure A7. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSf sample compared to modelled profiles.
Materials 17 01781 g0a7
Figure A8. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSf sample compared to modelled profiles.
Figure A8. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSf sample compared to modelled profiles.
Materials 17 01781 g0a8
Figure A9. Dependence of CS(X) on k in a single-pass, normal-freezing experiment. Parameters used: L = 100 cm, Z = 1, C0 = 1000 ppm, V = 0.0001 cm/s, δ = 0.01 cm, D = 1.01 × 10−5 cm2/s.
Figure A9. Dependence of CS(X) on k in a single-pass, normal-freezing experiment. Parameters used: L = 100 cm, Z = 1, C0 = 1000 ppm, V = 0.0001 cm/s, δ = 0.01 cm, D = 1.01 × 10−5 cm2/s.
Materials 17 01781 g0a9
Figure A10. Dependence of CS(X) on V in a single-pass, normal-freezing experiment. Parameters used: k = 0.5, L = 100 cm, Z = 1, C0 = 1000 ppm, δ = 0.01 cm, D = 1.01 × 10−5 cm2/s.
Figure A10. Dependence of CS(X) on V in a single-pass, normal-freezing experiment. Parameters used: k = 0.5, L = 100 cm, Z = 1, C0 = 1000 ppm, δ = 0.01 cm, D = 1.01 × 10−5 cm2/s.
Materials 17 01781 g0a10
Figure A11. Dependence of CS(X) on δ in a single-pass, normal-freezing experiment. Parameters used: k = 0.5, L = 100 cm, Z = 1, C0 = 1000 ppm, V = 0.0001 cm/s, D = 1.01 × 10−5 cm2/s.
Figure A11. Dependence of CS(X) on δ in a single-pass, normal-freezing experiment. Parameters used: k = 0.5, L = 100 cm, Z = 1, C0 = 1000 ppm, V = 0.0001 cm/s, D = 1.01 × 10−5 cm2/s.
Materials 17 01781 g0a11
Figure A12. Dependence of CS(X) on D in a single-pass, normal-freezing experiment. The D values are those noted by Sridharan et al. for rare-earth chlorides [6]. Parameters used: k = 0.5, L = 100 cm, Z = 1, C0 = 1000 ppm, V = 0.0001 cm/s, δ = 0.01 cm.
Figure A12. Dependence of CS(X) on D in a single-pass, normal-freezing experiment. The D values are those noted by Sridharan et al. for rare-earth chlorides [6]. Parameters used: k = 0.5, L = 100 cm, Z = 1, C0 = 1000 ppm, V = 0.0001 cm/s, δ = 0.01 cm.
Materials 17 01781 g0a12
Table A10. Model parameters of LKE-FSf experimental profiles.
Table A10. Model parameters of LKE-FSf experimental profiles.
ElementkV (cm/s)δ (cm)D (cm2/s)keff
Nd0.8160.00140.0012.40 × 10−50.825
La0.6610.00140.0011.01 × 10−50.691
Ce0.7010.00140.0011.28 × 10−50.723
Y0.9200.00140.0013.00 × 10−50.930
Sm0.8510.00140.0011.23 × 10−50.868
Pr0.8290.00140.0012.00 × 10−50.848

References

  1. International Energy Agency—Data and Statistics. Available online: https://www.iea.org/data-and-statistics?country=WORLD&fuel=Energy (accessed on 7 January 2024).
  2. Netter, P. Reprocessing of Spent Oxide Fuel from Nuclear Power Reactors; Woodhead Publishing Limited: Cambridge, UK, 2012; pp. 459–500. [Google Scholar]
  3. Dey, P.K.; Bansal, N.K. Spent fuel reprocessing: A vital link in Indian nuclear power program. Nucl. Eng. Des. 2006, 236, 723–729. [Google Scholar] [CrossRef]
  4. Riley, B.J. Electrochemical Salt Wasteform Development: A Review of Salt Treatment and Immobilization Options. Ind. Eng. Chem. Res. 2020, 59, 9760–9774. [Google Scholar] [CrossRef]
  5. Brambilla, G.; Sartorelli, A. Electrochemical Reprocessing of Nuclear. Fuels. Patent GB 2045511A, 29 October 1980. [Google Scholar]
  6. Sridharan, K.; Martin, S.; Mohammadian, M.; Sager, J.; Allen, T.; Simpson, M. Thermal properties of LiCl-KCl molten salt for nuclear waste separation. Trans. Am. Nucl. Soc. 2012, 106, 1240–1241. [Google Scholar]
  7. Gutknecht, T.Y.; Fredrickson, G.L. Fundamental Properties of Salts; Idaho National Lab.: Idaho Falls, ID, USA, 2012. [Google Scholar]
  8. Park, G.I.; Jeon, M.K.; Choi, J.H.; Lee, K.R.; Han, S.Y.; Kim, I.T.; Cho, Y.Z.; Park, H.S. Recent progress in waste treatment technology for pyroprocessing at KAERI. J. Nucl. Fuel Cycle Waste Technol. 2019, 17, 279–298. [Google Scholar] [CrossRef]
  9. Harrison, M. Waste from Pyro-processing: A Salty Challenge for Immobilisation Science. Nucl. Futur. 2021, 17, 39–48. [Google Scholar]
  10. Cho, Y.Z.; Park, G.H.O.; Lee, H.S.U.; Kim, I.N.T.; Han, D.S. Concentration of cesium and strontium elements involved in a LiCl waste salt by a melt crystallization process. Nucl. Technol. 2010, 171, 325–334. [Google Scholar] [CrossRef]
  11. Cho, Y.Z.; Ahn, B.G.; Eun, H.C.; Jung, J.S.; Lee, H.S. Melt crystallization process treatment of LiCl salt waste generated from electrolytic reduction process of spent oxide fuel. Energy Procedia 2011, 7, 525–528. [Google Scholar] [CrossRef]
  12. Versey, J.R.; Phongikaroon, S.; Simpson, M.F. Separation of CsCl from LiCl-CsCl molten salt by cold finger melt crystallization. Nucl. Eng. Technol. 2014, 46, 395–406. [Google Scholar] [CrossRef]
  13. Lee, H.S.; Oh, G.H.; Lee, Y.S.; Kim, I.T.; Kim, E.H.; Lee, J.H. Concentrations of CsCl and SrCl2 from a Simulated LiCl Salt Waste Generated by pyroprocessing by Using Czochralski Method. J. Nucl. Sci. Technol. 2009, 46, 392–397. [Google Scholar] [CrossRef]
  14. Shim, M.; Choi, H.G.; Yi, K.W.; Hwang, I.S.; Lee, J.H. Separation of CsCl and SrCl2 from a ternary CsCl-SrCl2-LiCl via a zone refining process for waste salt minimization of pyroprocessing. J. Nucl. Mater. 2016, 480, 403–410. [Google Scholar] [CrossRef]
  15. Shim, M.; Choi, H.G.; Choi, J.H.; Yi, K.W.; Lee, J.H. Separation of Cs and Sr from LiCl-KCl eutectic salt via a zone-refining process for pyroprocessing waste salt minimization. J. Nucl. Mater. 2017, 491, 9–17. [Google Scholar] [CrossRef]
  16. Pfann, W.G. Zone Melting. Science 1962, 135, 1101–1109. Available online: http://www.jstor.org/stable/1709243 (accessed on 2 February 2024). [CrossRef] [PubMed]
  17. Ezheiyan, M.; Sadeghi, H.; Tavakoli, M.H. Thermal Analysis Simulation of Germanium Zone Refining Process Assuming a Constant Radio-Frequency Heating Source. Chin. Phys. Lett. 2016, 33, 058102. [Google Scholar] [CrossRef]
  18. Ezheiyan, M.; Sadeghi, H. Simulation for purification process of high pure germanium by zone refining method. J. Cryst. Growth 2017, 462, 1–5. [Google Scholar] [CrossRef]
  19. Ito, T.; Ozaki, Y.; Tomioka, Y. Purification of β-Ga2O3 crystals by the zone refining method. Jpn. J. Appl. Phys. 2019, 58, 11. [Google Scholar] [CrossRef]
  20. Zhang, X.; Friedrich, S.; Friedrich, B. Separation behavior of arsenic and lead from antimony during vacuum distillation and zone refining. J. Mater. Res. Technol. 2020, 9, 4386–4398. [Google Scholar] [CrossRef]
  21. Mcalexander, T.W. Zone Refining of Raw Materials for Metal Halide Scintillator Crystal Growth. Master’s Thesis, University of Tennessee, Knoxville, TN, USA, 2014. [Google Scholar]
  22. Swider, S.; Lam, S.; Motakef, S.; Donohoe, E.; Coers, L.; Taylor, S.; Spencer, S. Impurity segregation in zone-refined precursors for crystalline halide scintillators. Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2015, 784, 5–8. [Google Scholar] [CrossRef]
  23. Pfann, W.G. Zone melting. Metall. Rev. 1957, 2, 29–76. [Google Scholar] [CrossRef]
  24. Kirk-Othmer. Encyclopedia of Chemical Technology, 3rd ed.; Wiley: New York, NY, USA, 1978. [Google Scholar]
  25. Spim, J.A.; Bernadou, M.J.S.; Garcia, A. Numerical modeling and optimization of zone refining. J. Alloys Compd. 2000, 298, 299–305. [Google Scholar] [CrossRef]
  26. Zhang, X.; Friedrich, S.; Liu, B.; Huang, T.; Friedrich, B. Computation-assisted analyzing and forecasting on impurities removal behavior during zone refining of antimony. J. Mater. Res. Technol. 2020, 9, 1221–1230. [Google Scholar] [CrossRef]
  27. IPCC. Climate Change 2014: Mitigation of Climate Change; IPCC: Geneva, Switzerland, 2014. [Google Scholar]
  28. MacDonald-Taylor, J.; Wright, J.; Edmondson, M.; Harrison, M.; Taylor, R. Pyro-processing: A hot topic in a sustainable fuel cycle. Nucl. Futur. 2021, 17, 46–54. [Google Scholar]
  29. Gaede, D. Chlorination and Selective Vaporization of Rare Earth Elements. Master’s Thesis, University of Montana, Missoula, MT, USA, 2016. [Google Scholar]
Figure 1. Schematic of zone-refinement apparatus.
Figure 1. Schematic of zone-refinement apparatus.
Materials 17 01781 g001
Figure 2. Schematics showing normal-freezing process of zone refining.
Figure 2. Schematics showing normal-freezing process of zone refining.
Materials 17 01781 g002
Figure 3. Measured temperature profiles within the Elite tube furnace as a function of set temperature and position.
Figure 3. Measured temperature profiles within the Elite tube furnace as a function of set temperature and position.
Materials 17 01781 g003
Figure 4. Cl, K, and Ba EDX maps of LKE-Ba5 (100× magnification).
Figure 4. Cl, K, and Ba EDX maps of LKE-Ba5 (100× magnification).
Materials 17 01781 g004
Figure 5. Cs/K, K/Cl, and Cs/Cl ratios of the LKE-Cs5_b sample, as measured by SEM-EDX, as a function of position along the sample.
Figure 5. Cs/K, K/Cl, and Cs/Cl ratios of the LKE-Cs5_b sample, as measured by SEM-EDX, as a function of position along the sample.
Materials 17 01781 g005
Figure 6. Images of electrorefining simulant samples following zone refinement.
Figure 6. Images of electrorefining simulant samples following zone refinement.
Materials 17 01781 g006
Figure 7. EDX maps of unvolatilised residue from LKE-FSb zone-refined sample.
Figure 7. EDX maps of unvolatilised residue from LKE-FSb zone-refined sample.
Materials 17 01781 g007
Figure 8. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSa sample compared to modelled profiles.
Figure 8. Experimental distribution profiles for Nd, Ce, and La for the LKE-FSa sample compared to modelled profiles.
Materials 17 01781 g008
Figure 9. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSa sample compared to modelled profiles.
Figure 9. Experimental distribution profiles for Y, Sm, and Pr for the LKE-FSa sample compared to modelled profiles.
Materials 17 01781 g009
Figure 10. Images of electroreducing simulant samples following zone refinement.
Figure 10. Images of electroreducing simulant samples following zone refinement.
Materials 17 01781 g010
Figure 11. Experimental distribution profiles for Ba, Sr, Cs, and Rb for the L-FSa sample compared to modelled profiles.
Figure 11. Experimental distribution profiles for Ba, Sr, Cs, and Rb for the L-FSa sample compared to modelled profiles.
Materials 17 01781 g011
Figure 12. Experimental distribution profiles for Ba, Sr, Cs, and Rb for the L-FSb sample compared to modelled profiles.
Figure 12. Experimental distribution profiles for Ba, Sr, Cs, and Rb for the L-FSb sample compared to modelled profiles.
Materials 17 01781 g012
Figure 13. Average keff for each electrorefining sample vs. estimated freezing-front velocity, V.
Figure 13. Average keff for each electrorefining sample vs. estimated freezing-front velocity, V.
Materials 17 01781 g013
Figure 14. Average keff for each element against atomic number/ionic radius (Z/IR).
Figure 14. Average keff for each element against atomic number/ionic radius (Z/IR).
Materials 17 01781 g014
Table 1. Simulant pyroprocessing salt waste compositions.
Table 1. Simulant pyroprocessing salt waste compositions.
Component (wt.%)Electroreducer Waste
(No-Cs)
Electroreducer Waste
(Full-Cs)
Electrorefiner Waste
LiCl84.1584.1539.60
KCl--50.40
Li2O0.850.85-
BaCl29.724.92-
CsCl-7.33-
RbCl1.400.73-
SrCl23.882.02-
YCl3--0.50
LaCl3--1.20
CeCl3--2.30
PrCl3--1.10
NdCl3--3.90
SmCl3--1.00
Table 2. Fit parameters for tube furnace temperature profiles.
Table 2. Fit parameters for tube furnace temperature profiles.
400 °C450 °C500 °C550 °C600 °C650 °C700 °C
a−0.107−0.055−0.094−0.077−0.007−0.176−0.104
b−2.283−4.4049−3.862−4.409−6.396−3.738−4.235
c353.79398.33449.10498.46552.29595.50643.80
Table 3. Zone-refining parameters for electrorefining simulant samples.
Table 3. Zone-refining parameters for electrorefining simulant samples.
SampleTarget Set Temp. (°C)Est. Temp. at
Inner End (°C)
Est. Temp. at Outer End (°C)Set Cooling Rate (°C/min)Est. Freezing-Front Vel., V (cm/s)
LKE-FSa5504744100.50.0014
LKE-FSb5504744100.10.0003
LKE-FSc5504744100.40.0011
LKE-FSd5504744100.30.0008
LKE-FSe5504744100.50.0014
LKE-FSf5504744100.50.0014
Table 4. Zone-refining parameters for electroreducing simulant samples.
Table 4. Zone-refining parameters for electroreducing simulant samples.
SampleTarget Set Temp. (°C)Est. Temp. at Inner End (°C)Est. Temp. at Outer End (°C)Set Cooling Rate (°C/min)Est. Freezing-Front Vel., V (cm/s)
L-FSa7706906210.50.0012
L-FSb7706906210.50.0012
Table 5. Calculated DFi values for zone refining of electrorefining simulant samples, XR = 0.8.
Table 5. Calculated DFi values for zone refining of electrorefining simulant samples, XR = 0.8.
LKE-FSaLKE-FScLKE-FSdLKE-FSeLKE-FSf
DFNd2.862.082.181.291.47
DFLa1.521.561.261.121.96
DFCe1.131.39--1.82
DFY1.231.731.561.081.15
DFSm1.851.652.451.251.31
DFPr1.811.721.811.241.37
Table 6. Calculated DFi values for zone refining of electrorefining simulant samples, XR = 0.9.
Table 6. Calculated DFi values for zone refining of electrorefining simulant samples, XR = 0.9.
LKE-FSaLKE-FScLKE-FSdLKE-FSeLKE-FSf
DFNd3.602.472.611.381.62
DFLa1.691.751.341.162.30
DFCe1.171.52--2.10
DFY1.291.981.751.111.20
DFSm2.151.873.011.321.41
DFPr2.091.972.091.311.49
Table 7. Estimated electrorefining refined salt composition, XR = 0.8.
Table 7. Estimated electrorefining refined salt composition, XR = 0.8.
wt.%LKE-FSaLKE-FScLKE-FSdLKE-FSeLKE-FSf
LiCl40.2940.1740.1439.7540.04
KCl51.2751.1351.0950.6050.96
YCl30.480.440.450.490.49
LaCl31.091.081.141.171.01
CeCl32.242.132.302.301.98
PrCl30.950.960.951.051.02
NdCl32.843.213.153.683.56
SmCl30.850.890.770.950.94
Total100.00100.00100.00100.00100.00
Total RECl38.448.708.769.659.00
Table 8. Estimated electrorefining waste salt composition, XR = 0.8.
Table 8. Estimated electrorefining waste salt composition, XR = 0.8.
wt.%LKE-FSaLKE-FScLKE-FSdLKE-FSeLKE-FSf
LiCl36.8837.3237.4338.9837.84
KCl46.9347.4947.6449.6148.16
YCl30.590.750.700.530.56
LaCl31.651.681.441.321.97
CeCl32.542.972.302.303.59
PrCl31.711.651.711.301.40
NdCl38.116.676.874.775.24
SmCl31.581.461.901.191.24
Total100.00100.00100.00100.00100.00
Total RECl316.1915.1914.9211.4014.00
Table 9. Estimated electrorefining refined salt composition, XR = 0.9.
Table 9. Estimated electrorefining refined salt composition, XR = 0.9.
wt.%LKE-FSaLKE-FScLKE-FSdLKE-FSeLKE-FSf
LiCl40.1140.0039.9939.7039.91
KCl51.0550.9150.9050.5350.79
YCl30.490.460.470.490.49
LaCl31.121.121.161.181.06
CeCl32.262.192.302.302.07
PrCl30.991.000.991.071.05
NdCl33.093.403.363.763.67
SmCl30.900.920.830.970.96
Total100.00100.00100.00100.00100.00
Total RECl38.859.089.119.779.30
Table 10. Estimated electrorefining waste salt composition, XR = 0.9.
Table 10. Estimated electrorefining waste salt composition, XR = 0.9.
wt.%LKE-FSaLKE-FScLKE-FSdLKE-FSeLKE-FSf
LiCl35.0735.9736.0838.6836.85
KCl44.6345.7845.9149.2346.90
YCl30.630.900.810.550.59
LaCl31.901.951.551.372.44
CeCl32.653.322.302.304.35
PrCl32.071.972.071.401.56
NdCl311.128.408.775.205.95
SmCl31.931.722.501.281.35
Total100.00100.00100.00100.00100.00
Total RECl320.3018.2518.0112.0916.25
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Scrimshire, A.; Backhouse, D.J.; Deng, W.; Mann, C.; Ogden, M.D.; Sharrad, C.A.; Harrison, M.T.; McKendrick, D.; Bingham, P.A. Benchtop Zone Refinement of Simulated Future Spent Nuclear Fuel Pyroprocessing Waste. Materials 2024, 17, 1781. https://doi.org/10.3390/ma17081781

AMA Style

Scrimshire A, Backhouse DJ, Deng W, Mann C, Ogden MD, Sharrad CA, Harrison MT, McKendrick D, Bingham PA. Benchtop Zone Refinement of Simulated Future Spent Nuclear Fuel Pyroprocessing Waste. Materials. 2024; 17(8):1781. https://doi.org/10.3390/ma17081781

Chicago/Turabian Style

Scrimshire, Alex, Daniel J. Backhouse, Wei Deng, Colleen Mann, Mark D. Ogden, Clint A. Sharrad, Mike T. Harrison, Donna McKendrick, and Paul A. Bingham. 2024. "Benchtop Zone Refinement of Simulated Future Spent Nuclear Fuel Pyroprocessing Waste" Materials 17, no. 8: 1781. https://doi.org/10.3390/ma17081781

APA Style

Scrimshire, A., Backhouse, D. J., Deng, W., Mann, C., Ogden, M. D., Sharrad, C. A., Harrison, M. T., McKendrick, D., & Bingham, P. A. (2024). Benchtop Zone Refinement of Simulated Future Spent Nuclear Fuel Pyroprocessing Waste. Materials, 17(8), 1781. https://doi.org/10.3390/ma17081781

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop