High Level Waste in Salt Repositories: Experiments and Simulations of Evaporation in the Underground - 17167
Advances in simulation of coupled heat and mass transfer in granular salt using the Finite Element Heat and Mass Transfer Code (FEHM) is presented in this paper. The modeling focuses on reproducing the results of an in-situ subsurface evaporation experiment beginning during the summer of 2015 carried out in the underground at the Waste Isolation Pilot Plant (WIPP). Run-of-mine (granular) salt was placed in a metal pan, and sample mass, as well as ambient temperature and humidity, were measured continuously on a data logger through May 2016. Changes in moisture content of the salt are interpreted from the change in mass, and these data provide insight into evaporation and water retention characteristics from salt in the WIPP underground. To allow simulation of the experimental data, FEHM was modified to include a variable humidity boundary condition. The new humidity boundary condition can adjust streams of inflowing air to reflect fractions of air and water at prescribed temperature and pressure conditions that in turn represent relative humidity. Attributes for the flow such as temperature, air mass fraction, and humidity can be specified for inflowing air. The in-situ experiment provided data for validation of FEHM code development to accurately compute water-vapor pressure lowering due to dissolved salt and capillary pressure for pore water stored granular WIPP salt. Using this new boundary condition capabilities, a numerical representation of the pan experiment was created. The numerical domain is a 2-dimensional square grid of 30 by 30 cm axes. Salt and air properties are assigned based on measured values. A low flow rate of air is specified across the pile of salt to represent drift ventilation airflow. Humidity and temperature of the flowing air are fixed to measured values, and the simulated pile of salt is allowed to take up or release water based only on gradients in humidity and water-vapor pressure. Mass changes of the modeled salt pan were compared to measured mass changes in the bench-scale experiment. Directionality of the modeled changes tracks well with measured values, with good correlation between peak and valley inflection points. This shows that in situations with airflow where the relative humidity boundary can apply, FEHM can produce accurate values if humidity, air flow, and initial conditions are well established. Finally, this work supports phased-field-thermal testing in salt planned for the underground at WIPP that is expected to begin in early 2017. (authors)