Accurate and efficient calculation of response times for groundwater flow

We study measures of the amount of time required for transient flow in heterogeneous porous media to effectively reach steady state, also known as the response time. Here, we develop a new approach that extends the concept of mean action time. Previous applications of the theory of mean action time to estimate the response time use the first two central moments of the probability density function associated with the transition from the initial condition, at $t=0$, to the steady state condition that arises in the long time limit, as $t \to \infty$. This previous approach leads to a computationally convenient estimation of the response time, but the accuracy can be poor. Here, we outline a powerful extension using the first $k$ raw moments, showing how to produce an extremely accurate estimate by making use of asymptotic properties of the cumulative distribution function. Results are validated using an existing laboratory-scale data set describing flow in a homogeneous porous medium. In addition, we demonstrate how the results also apply to flow in heterogeneous porous media. Overall, the new method is: (i) extremely accurate; and (ii) computationally inexpensive. In fact, the computational cost of the new method is orders of magnitude less than the computational effort required to study the response time by solving the transient flow equation. Furthermore, the approach provides a rigorous mathematical connection with the heuristic argument that the response time for flow in a homogeneous porous medium is proportional to $L^2/D$, where $L$ is a relevant length scale, and $D$ is the aquifer diffusivity. Here, we extend such heuristic arguments by providing a clear mathematical definition of the proportionality constant.Comment: 22 pages, 3 figures, accepted version of paper published in Journal of Hydrolog

New homogenization approaches for stochastic transport through heterogeneous media

The diffusion of molecules in complex intracellular environments can be strongly influenced by spatial heterogeneity and stochasticity. A key challenge when modelling such processes using stochastic random walk frameworks is that negative jump coefficients can arise when transport operators are discretized on heterogeneous domains. Often this is dealt with through homogenization approximations by replacing the heterogeneous medium with an $\textit{effective}$ homogeneous medium. In this work, we present a new class of homogenization approximations by considering a stochastic diffusive transport model on a one-dimensional domain containing an arbitrary number of layers with different jump rates. We derive closed form solutions for the $k$th moment of particle lifetime, carefully explaining how to deal with the internal interfaces between layers. These general tools allow us to derive simple formulae for the effective transport coefficients, leading to significant generalisations of previous homogenization approaches. Here, we find that different jump rates in the layers gives rise to a net bias, leading to a non-zero advection, for the entire homogenized system. Example calculations show that our generalized approach can lead to very different outcomes than traditional approaches, thereby having the potential to significantly affect simulation studies that use homogenization approximations.Comment: 9 pages, 2 figures, accepted version of paper published in The Journal of Chemical Physic

On the non-vanishing of the Collins mechanism for single spin asymmetries

The Collins mechanism provides a non-perturbative explanation for the large single spin asymmetries found in hard semi-inclusive reactions involving a transversely polarized nucleon. However, there are seemingly convincing reasons to suspect that the mechanism vanishes, and indeed it does vanish in the naive parton model where a quark is regarded as an essentially 'free' particle. We give an intuitive analysis which highlights the difference between the naive picture and the realistic one, and shows how the Collins mechanism arises when the quark is described as an off-shell particle by a field in interaction. A typographical error is corrected in this version.Comment: 15 pages, 2 figure

Raman scattering in a Heisenberg {\boldmath S=1/2S=1/2} antiferromagnet on the triangular lattice

We investigate two-magnon Raman scattering from the $S=1/2$ Heisenberg antiferromagnet on the triangular lattice, considering both the effect of renormalization of the one-magnon spectrum by 1/S corrections and final-state magnon-magnon interactions. The bare Raman intensity displays two peaks related to one-magnon van-Hove singularities. We find that 1/S self-energy corrections to the one-magnon spectrum strongly modify this intensity profile. The central Raman-peak is significantly enhanced due to plateaus in the magnon dispersion, the high frequency peak is suppressed due to magnon damping, and the overall spectral support narrows considerably. Additionally we investigate final-state interactions by solving the Bethe-Salpeter equation to $O(1/S)$. In contrast to collinear antiferromagnets, the non-collinear nature of the magnetic ground state leads to an irreducible magnon scattering which is retarded and non-separable already to lowest order. We show that final-state interactions lead to a rather broad Raman-continuum centered around approximately twice the 'roton'-energy. We also discuss the dependence on the scattering geometry.Comment: 7 pages, 5 figure

Calculating how long it takes for a diffusion process to effectively reach steady state without computing the transient solution

Mathematically, it takes an infinite amount of time for the transient solution of a diffusion equation to transition from initial to steady state. Calculating a \textit{finite} transition time, defined as the time required for the transient solution to transition to within a small prescribed tolerance of the steady state solution, is much more useful in practice. In this paper, we study estimates of finite transition times that avoid explicit calculation of the transient solution by using the property that the transition to steady state defines a cumulative distribution function when time is treated as a random variable. In total, three approaches are studied: (i) mean action time (ii) mean plus one standard deviation of action time and (iii) a new approach derived by approximating the large time asymptotic behaviour of the cumulative distribution function. The new approach leads to a simple formula for calculating the finite transition time that depends on the prescribed tolerance $\delta$ and the $(k-1)$th and $k$th moments ($k \geq 1$) of the distribution. Results comparing exact and approximate finite transition times lead to two key findings. Firstly, while the first two approaches are useful at characterising the time scale of the transition, they do not provide accurate estimates for diffusion processes. Secondly, the new approach allows one to calculate finite transition times accurate to effectively any number of significant digits, using only the moments, with the accuracy increasing as the index $k$ is increased.Comment: 17 pages, 2 figures, accepted version of paper published in Physical Review

Fast computation of effective diffusivities using a semi-analytical solution of the homogenization boundary value problem for block locally-isotropic heterogeneous media

Direct numerical simulation of diffusion through heterogeneous media can be difficult due to the computational cost of resolving fine-scale heterogeneities. One method to overcome this difficulty is to homogenize the model by replacing the spatially-varying fine-scale diffusivity with an effective diffusivity calculated from the solution of an appropriate boundary value problem. In this paper, we present a new semi-analytical method for solving this boundary value problem and computing the effective diffusivity for pixellated, locally-isotropic, heterogeneous media. We compare our new solution method to a standard finite volume method and show that equivalent accuracy can be achieved in less computational time for several standard test cases. We also demonstrate how the new solution method can be applied to complex heterogeneous geometries represented by a grid of blocks. These results indicate that our new semi-analytical method has the potential to significantly speed up simulations of diffusion in heterogeneous media.Comment: 29 pages, 4 figures, 5 table