WSRC-MS-2001-00603

Prevention of Solids Formation During Canyon Processing
of Legacy Nuclear Materials at the Savannah River Site

William D. Rhodes and William. J. Crooks, III
Westinghouse Savannah River Company
Aiken, SC 29808

Jerry D. Christian, Consultant
Electrode Specialties Company.
Idaho Falls, ID 83404

This document was prepared in conjunction with work accomplished under Contract No. DE-AC09-96SR18500 with the U.S. Department of Energy.

DISCLAIMER

This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.

This report has been reproduced directly from the best available copy.

Available for sale to the public, in paper, from:  U.S. Department of Commerce, National Technical Information Service, 5285 Port Royal Road, Springfield, VA 22161,  phone: (800) 553-6847,  fax: (703) 605-6900,  email:  orders@ntis.fedworld.gov   online ordering:  http://www.ntis.gov/support/index.html

Available electronically at  http://www.osti.gov/bridge/

Available for a processing fee to U.S. Department of Energy and its contractors, in paper, from: U.S. Department of Energy, Office of Scientific and Technical Information, P.O. Box 62, Oak Ridge, TN 37831-0062,  phone: (865 ) 576-8401,  fax: (865) 576-5728,  email:  reports@adonis.osti.gov

Abstract

The Savannah River Site (SRS) Environmental Management (EM) nuclear material stabilization program includes the dissolution of legacy materials from various DOE sites. The SRS canyon facilities were designed to dissolve and process spent nuclear fuel and targets. As the processing of typical materials is completed, unusual and exotic nuclear materials are being targeted for stabilization. These unusual materials are often difficult to dissolve using historical flowsheet conditions and require more aggressive dissolver solutions. Solids must be prevented in the dissolver to avoid expensive delays associated with the build-up of insoluble material in downstream process equipment. Moreover, it is vital to prevent precipitation of all solids, especially plutonium-bearing solids, since their presence in dissolver solutions raises criticality safety issues.

To prevent precipitation of undesirable solids in aqueous process solutions, the accuracy of computer models to predict the formation of precipitate formation requires incorporation of plant specific fundamental data. Incorporation of these data into a previously developed thermodynamic computer program, that applies the Pitzer correlation to derive activity coefficient parameters, has enhanced the ability to minimize the potential of unwanted precipitation under a wide variety of operating conditions at any DOE site working with EM materials in aqueous solutions.

Introduction

Research and development focused on SRS canyon dissolver precipitation issues was important during the Sand, Slag, and Crucible (SS&C) campaign of 1997. During this campaign, white solids were identified as potassium tetrafluoroborate (KBF4), indicating a reduction in soluble boron, a neutron adsorbing poison. The conditions that shift the equilibria towards precipitation are qualitatively understood in terms of Le Chatelier’s principle by considering the following equation:

K+ (aq) + H3BO3 (aq) + 4 F- (aq) + 3 H+ (aq) = KBF4 (s) + 3 H2O (l) (1)

The Idaho National Engineering and Environmental Laboratory (INEEL) developed expertise in aqueous fluoride chemistry as a result of processing naval nuclear fuels at the Idaho Nuclear Technology and Engineering Center (INTEC, formerly the Idaho Chemical Processing Plant, ICPP). This process included dissolution in hydrofluoric and nitric acids that incorporated boron as a soluble neutron poison for criticality control. This processing need required development of a thermodynamic speciation program for predicting multiple fluoride species equilibrium concentrations in representative plant solutions. As a result of the SS&C campaign issues, the INEEL model was used to predict residue dissolution of calcium fluoride in the presence of boric acid and to control corrosion of the stainless steel dissolver vessel. However, the INEEL speciation program thermodynamic data are applicable at ionic strength conditions for the INEEL process solutions, i.e., do not have activity coefficient data. Therefore, application to SRS solutions with high ionic strength requires specific chemical species information. To this end, the INEEL speciation computer program is being updated with new basic chemical data in order to better predict and avoid solids production precipitation in aqueous process solutions at SRS.

The objective of the project is to incorporate activity coefficients into the speciation program that has been developed to calculate individual component concentrations in acidic aqueous fluoride systems. This will enable accurate predictions of solubilities of potentially precipitating species in plant solutions and provide the ability to calculate solution adjustments to assure stability. In order to do this, solubility and activity coefficient data must be fitted to a suitable activity coefficient model and its parameters determined. Then, the fitted model can be used to calculate the activity coefficients for process solution compositions. The computer program has potential applications at DOE sites working with EM materials in aqueous solutions.

Modeling to Address Precipitation in the Canyon Dissolver

In laboratory tests to support the Sand, Slag, and Crucible (SS&C) campaign and the Mark 42 Fuel Tube campaign, the presence of high concentration of fluoride ions in boric acid/nitric acid solutions led to the formation of a white solid (see Table 1). The white solids were collected from flowsheet simulations, and were identified as KBF4.

Table 1.
Identification of KBF4 Precipitate in SRS Dissolver Simulation Tests.

Date

Test

[HNO3]o

(M)

[F-]o

(M)

[B]o

(g/L)

Observation

Dec. 1997

SS&C - simulation

9.3

0.30

2.5

Unidentified white solid

May 1998

SS&C - simulation

1.0

0.23

1.7

White solid, KBF4

Nov. 1998

SS&C – test 1
SS & C – test 2

8.8

0.32

1.6

No solids

8.6

0.40

2.2

KBF4 (s)

Feb. 1999

Mark 42 – simulation
[Al] = 0.44 M

1.0

0.40

2.5

No solids

1.0

0.50

2.5

KBF4 (s), minor

1.0

0.60

2.5

KBF4 (s), more

 

Without known KBF4 activity coefficients at the conditions evaluated, the INEEL program under predicted the saturation of KBF4, as shown in Table 2.

Table 2.
Modeling Results for KBF4 Experiments, 20oC

[HNO3]0

[KF]0

[B]0

[Al]o

Mark-42: simulation

Using INEEL program, Calculated:

M

M

(g/L)

(M)

Observation

Calculated

1.0

0.50

2.5

0.44

KBF4 (s), few

[BF4-] = 41.2% of saturation
(i.e. no precipitation is predicted).

1.0

0.60

2.5

0.44

KBF4 (s), some

Predicts saturated KBF4.
Calculated KSP= 1.2724 x 10-3
[vs. literature: 1.27 x 10-3]
Precipitate composition:
1.9% of K+
7.6% of F (4.9% of F as KBF4)


In recent years, the INEEL modeling capability has been expanded with the incorporation of complexation equilibrium calculations into a free energy minimization program with a database for over 15,000 compounds. To apply the model to new applications, the user incorporates data for the performance of phase equilibrium calculations. The INEEL program applies the Pitzer model,, a widely used model for which parameters have been extensively tabulated for various salts and acids. For applications to multielectrolyte solutions, data from both single and binary salt solutions are required to obtain ion interaction parameters for all ions in solution. Figure 1 shows the prediction capability of Pitzer single-salt equation parameters for NaNO3 activity coefficients. Without the plant specific activity coefficients, the The Pitzer equation is suitable to about 6 molal, but not must be evaluated on a case-by case basis at higher ionic strengths.

The INEEL model incorporates multiple fluoride complexation constants and solubilities of fluoride species (e.g. aluminum fluoride and zirconium fluoride) that are involved in multiple complexation equilibria. This enables predictions of conditions (e.g. reagent concentrations and temperatures) that assure solution stability. The INEEL model will be applied to evaluate SRS dissolver solution compositions and predict equilibrium concentrations and the possible formation of undesirable solids. However, at ionic strengths pertinent to SRS plant solutions and specifically to the KBF4 solubility product and activity coefficient determinations, the model needs improvement via incorporation of experimental data.

At SRS, the following interactions are important: KBF4 – NaNO3 (no common ion), KBF4 – NaBF4 (common anion), and KBF4 – KNO3 (common cation). Based on solubility measurements at plant solution ionic strengths, binary and ternary KBF4 activity coefficient parameters were determined. These data enable solubility extrapolation to zero ionic strength and determination of Pitzer parameters.

Once the salt solubilities have been determined as a function of ionic strength, the activity coefficients are calculated as follows. For the general salt dissolution, Eq. (2), the concentration equilibrium constant and thermodynamic equilibrium constants are obtained by Eqs. (3) and (4).

Figure 1 Pitzer Coefficients for NaNO3 Single Salt Equation Fitted to Hamer and Wu Data
Figure 1
Pitzer Coefficients for NaNO3 Single Salt Equation Fitted to Hamer & Wu Data

AxBy = xAz+ + yBz- (2)

Equation (3)

Equation (4)

Here, C is molar concentration, a is activity, g is mean activity coefficient [(Equation)1/n ], and n is x + y. Let S0 and g 0 be the solubility and mean activity coefficient, respectively, of the salt in pure H2O and S and g be the corresponding values in a solution with added electrolyte (with no common ion) that increases the ionic strength :

[m = ½Equation] (5)

KTh = xxyyS0g 0 = xxyySg (6)

so that Sg = S0g 0. Taking logarithms, we have

log S = log S0g 0 – log g (7)

Once S0g 0 is known, the activity coefficient at a given ionic strength is obtained from the measured solubility. To obtain S0g 0, log S is plotted against (m )1/2. The plot is extrapolated to (m )1/2 = 0. The intercept gives S0 at zero ionic strength, where g 0 = 1 and, thus, we have the value of S0g 0. Eq. (7) is solved for g from the measured solubility at each ionic strength.

A commercial free energy minimization program, HSC Chemistry® for Windows, provides the capability of inputting enthalpy of formation, entropy, and heat capacity terms for individual species. Simple activity coefficient expressions or the values can also be inputted. In the case of experimental solubility constants, thermodynamic data are expressed for the reaction; individual species values are not provided. The INEEL model possesses general equations and methodology to convert equilibrium constants into a consistent set of thermodynamic parameters for use in the HSC database and program. Based on the experimental solubility data, the activity coefficients are obtained from the INEEL model. The plant solution stability is evaluated with the application of the HSC program. Solution compositions can be varied to determine the concentration limit at which precipitation will begin.

Determination of Binary and Ternary Activity Coefficients

Various well-established thermodynamic methods are known for determining the activity coefficients of electrolytic solutions. These methods include vapor pressure, freezing point depression, boiling point elevation, osmotic pressure, solubility, and electromotive force measurements. Activity coefficients of KBF4 as a function of ionic strength will be determined by simple solubility measurements at various ionic strengths. Specifically, the determination of KBF4 binary and ternary activity coefficient parameters was based on KBF4 solubility measurements as a function of the ionic strength of an adjuster (NaNO3, NaBF4, and KNO3). The fluoborate ion hydrolyzes slightly to give yield H3BO3 and HF. Therefore, chemical additions (small additions of HF and H3BO3 at levels that will not contribute to ion interactions) were made to the test solutions, preventing hydrolysis of BFEquation that would otherwise occur to about 3.7%. These data, along with literature values of Pitzer parameters for interactions of Na+-Equation, K+-Equation, Na+-Equation, and K+-Na+ enable evaluation of all pertinent two-salt interaction parameters yielding KBF4 activity coefficients as a function of ionic strength. The KBF4 solution was analyzed for B and K concentration by inductively couple plasma-atomic emission spectroscopy (ICP-AES).

The solubility of the boron in the NaNO3 and KNO3 ionic strength solutions is shown in Figure 32; whereas, the solubility of potassium in NaNO3 and NaBF4 ionic solutions in shown in Figure 43. As expected, both cases the show a drastic dramatic drop in solubility with the interfering ionic species. (The boron concentration at zero ionic strength addition is slightly greater than the potassium concentration because of the H3BO3 added to prevent Equationhydrolysis; this will be accounted for in deriving solubility products.)

Summary

With the objective of preventing precipitation of undesirable solids during aggressive SRS dissolution processes of EM materials, the INEEL computer program is being updated with new basic chemical data resulting in a better ability to predict and avoid solids production in aqueous process solutions at SRS. The basic chemical data includes solubility, activity coefficients, and solubility products of potassium tetrafluoroborate (KBF4) at ionic strengths expected in process solutions. This program will calculate the equilibrium position for a given starting dissolver solution composition and the solution stability is determined. Solution compositions can be varied to determine the concentration limit at which precipitation will begin in a dissolver solution.

Figure 2 Anion Interaction Effects Reduce B Solubility
Figure 2.
Anion Interaction Effects Reduce B Solubility.

 

Figure 3 Cation Interaction Effects Reduce K Solubility
Figure 3.
Cation Interaction Effects Reduce K Solubility.

This effort will enable the avoidance of aqueous solution concentrations that may cause solids formation. Processing of off-normal material will have less chance to produce unwanted solids that stop work, and potentially delay processing canyon facilities campaigns. Delays would in turn result in higher overall life cycle operating costs. Moreover, less time will be needed in the

This effort will enable the avoidance of aqueous solution concentrations that may cause solids formation. Processing of off-normal material will have less chance to produce unwanted solids that stop work, and potentially delay processing canyon facilities campaigns. Delays would in turn result in higher overall life cycle operating costs. Moreover, less time will be needed in the future to determine the concentration of dissolution solutions since the region where no solids form is better defined. Results should immediately impact dissolution of Rocky Flats material by shortening the time it takes to determine dissolving solutions. In the long run, schedules to dissolve off-normal material or process aqueous solutions that are stored throughout the complex shall be favorably impacted.

Future Work

Future efforts shall continue to address the INEEL speciation model with the incorporation of experimentally determined mercury fluoride (HgF2) and plutonium fluoride (PuF4) solubilities. The solubility and activity coefficients of PuF4 will be determined by measuring the solubility of ThF4 as a function of ionic strength, and correcting the determined complexation constants and activity coefficients to PuF4 using the Born equation. Other specific solids of interest for future work include calcium fluoride, boric acid, aluminum nitrate, plutonium salts, and insoluble oxides from stainless steel corrosion products.

Hydrofluoric acid (HF) is an important species in modeling complexation equilibria and solubilities of fluoride salts in process solutions. In the current INEEL speciation program model, its activity coefficient is assumed to be unity. This has been adequate for INEEL process solutions in which the free HF and HNO3 concentrations have been less than 0.1 and1.8 molar, respectively. However, at greater concentrations of each, as occur in SRS process solutions, the activity coefficient of HF increases dramatically. Accordingly, the activity coefficients of HF as a function of ionic strength will be determined by measuring the partial pressure above a solution by infrared spectroscopy or by an alternate alternative transpiration technique.

References:

  1. K. S. Pitzer, J. Phys. Chem., 77, 268 (1973).
  2. K. S. Pitzer, "Ion Interaction Approach: Theory and Data Correlation," Chapter 3 in Activity Coefficients in Electrolyte Solutions, 2nd Edition, K. S. Pitzer, editor, CRC Press, Boca Raton, 1991.
  3. W. J. Hamer and Y. Wu, "Osmotic Coefficients and Mean Activity Coefficients of Uni-univalent Electrolytes in Water at 25°C," J. Phys. Chem. Ref. Data, Vol. 1, No. 4, 1047-1099 (1972).
  4. A. Roine, Outokumpo HSC Chemistry® for Windows, Chemical Reaction and Equilibrium Software with Extensive Thermochemical Database, Version 4.0, Outokumpo Research Oy Information Service, P O Box 60, FIN-28101 Pori, Finland; available from ESM Software, Hamilton, Ohio.
  5. C.F. Prutton and A.H. Maron, Fundamentals of Physical Chemistry (The MacMillan Company, 5th Edition, New York, 1957), 510.
  6. R. F. Platford, "Osmotic and activity coefficients of simple borates in aqueous solution at 25°", Canadian Journal of Chemistry, 47, 2271-2273 (1969).