Analytical modelling of bacterial migration and accumulation with rate-limited sorption in discrete fractures

An analytical model for bacterial accumulation in a discrete fractllre has been developed. The transport and accumlllation processes incorporate into the model include advection, dispersion, rate-limited adsorption, rate-limited desorption, irreversible adsorption, attachment, detachment, growth and first order decay botl1 in sorbed and aqueous phases. An analytical solution in Laplace space is derived and nlln1erically inverted. The model is implemented in the code BIOFRAC vvhich is written in Fortran 99. The model is derived for two phases, Phase I, where adsorption-desorption are dominant, and Phase II, where attachment-detachment are dominant. Phase I ends yvhen enollgh bacteria to fully cover the substratllm have accllillulated. The model for Phase I vvas verified by comparing to the Ogata-Banks solution and the model for Phase II was verified by comparing to a nonHomogenous version of the Ogata-Banks solution. After verification, a sensitiv"ity analysis on the inpllt parameters was performed. The sensitivity analysis was condllcted by varying one inpllt parameter vvhile all others were fixed and observing the impact on the shape of the clirve describing bacterial concentration verSllS time. Increasing fracture apertllre allovvs more transport and thus more accllffilliation, "Vvhich diminishes the dllration of Phase I. The larger the bacteria size, the faster the sllbstratum will be covered. Increasing adsorption rate, was observed to increase the dllration of Phase I. Contrary to the aSSllmption ofllniform biofilm thickness, the accllffilliation starts frOll1 the inlet, and the bacterial concentration in aqlleous phase moving towards the olitiet declines, sloyving the accumulation at the outlet. Increasing the desorption rate, redllces the dliration of Phase I, speeding IIp the accllmlilation. It was also observed that Phase II is of longer duration than Phase I. Increasing the attachment rate lengthens the accliffililation period. High rates of detachment speeds up the transport. The grovvth and decay rates have no significant effect on transport, althollgh increases the concentrations in both aqueous and sorbed phases are observed. Irreversible adsorption can stop accllillulation completely if the vallIes are high.

Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.