Applying reactive models to column experiments to assess the hydrogeochemistry of seawater intrusion: optimising ACUAINTRUSION and selecting cation exchange coefficients with PHREEQC

Three sets of laboratory column experimental results concerning the hydrogeochemistry of seawater intrusion have been modelled using two codes: ACUAINTRUSION (Chemical Engineering Department, University of Alicante) and PHREEQC (U.S.G.S.). These reactive models utilise the hydrodynamic parameters determined using the ACUAINTRUSION TRANSPORT software and fit the chloride breakthrough curves perfectly. The ACUAINTRUSION code was improved, and the instabilities were studied relative to the discretisation. The relative square errors were obtained using different combinations of the spatial and temporal steps: the global error for the total experimental data and the partial error for each element. Good simulations for the three experiments were obtained using the ACUAINTRUSION software with slight variations in the selectivity coefficients for both sediments determined in batch experiments with fresh water. The cation exchange parameters included in ACUAINTRUSION are those reported by the Gapon convention with modified exponents for the Ca/Mg exchange. PHREEQC simulations performed using the Gains-Thomas convention were unsatisfactory, with the exchange coefficients from the database of PHREEQC (or range), but those determined with fresh water – natural sediment allowed only an approximation to be obtained. For the treated sediment, the adjusted exchange coefficients were determined to improve the simulation and are vastly different from those from the database of PHREEQC or batch experiment values; however, these values fall in an order similar to the others determined under dynamic conditions. Different cation concentrations were simulated using two different software packages; this disparity could be attributed to the defined selectivity coefficients that affect the gypsum equilibrium. Consequently, different calculated sulphate concentrations are obtained using each type of software; a smaller mismatch was predicted using ACUAINTRUSION. In general, the presented simulations by ACUAINTRUSION and PHREEQC produced similar results, making predictions consistent with the experimental data. However, the simulated results are not identical to the experimental data; sulphate (total S) is overpredicted by both models, most likely due to such factors as the kinetics of gypsum, the possible variations in the exchange coefficients due to salinity and the neglect of other processes.

Mineralogical analysis of ceramic tiles by FTIR: A quantitative attempt

A method for quantitative mineralogical analysis by ATR-FTIR has been developed. The method relies on the use of the main band of calcite as a reference for the normalization of the IR spectrum of a mineral sample. In this way, the molar absorptivity coefficient in the Lambert–Beer law and the components of a mixture in mole percentage can be calculated. The GAMS equation modeling environment and the NLP solver CONOPT (©ARKI Consulting and Development) were used to correlate the experimental data in the samples considered. Mixtures of different minerals and gypsum were used in order to measure the minimum band intensity that must be considered for calculations and the detection limit. Accordingly, bands of intensity lower than 0.01 were discarded. The detection limit for gypsum was about 7% (mol/total mole). Good agreement was obtained when this FTIR method was applied to ceramic tiles previously analyzed by X-ray diffraction (XRD) or mineral mixtures prepared in the lab.