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.

Energy Loss Function of Solids Assessed by Ion Beam Energy-Loss Measurements: Practical Application to Ta2O5

We present a study where the energy loss function of Ta2O5, initially derived in the optical limit for a limited region of excitation energies from reflection electron energy loss spectroscopy (REELS) measurements, was improved and extended to the whole momentum and energy excitation region through a suitable theoretical analysis using the Mermin dielectric function and requiring the fulfillment of physically motivated restrictions, such as the f- and KK-sum rules. The material stopping cross section (SCS) and energy-loss straggling measured for 300–2000 keV proton and 200–6000 keV helium ion beams by means of Rutherford backscattering spectrometry (RBS) were compared to the same quantities calculated in the dielectric framework, showing an excellent agreement, which is used to judge the reliability of the Ta2O5 energy loss function. Based on this assessment, we have also predicted the inelastic mean free path and the SCS of energetic electrons in Ta2O5.