MARFA user's manual: Migration analysis of radionuclides in the far field

The computer code Migration Analysis of Radionuclides in the Far Field (MARFA) uses a particle-based Monte Carlo method to simulate the transport of radionuclides in a sparsely fractured geological medium. The algorithm uses non-interacting particles to represent packets of radionuclide mass. These particles are moved through the system according to rules that mimic the underlying physical transport and retention processes. The physical processes represented in MARFA include advection, longitudinal dispersion, Fickian diffusion into an infinite or finite rock matrix, equilibrium sorption, decay, and in-growth. Because the algorithm uses non-interacting particles, the transport and retention processes are limited to those that depend linearly on radionuclide concentration. Multiple non-branching decay chains of arbitrary length are supported, as is full heterogeneity in the transport and retention properties. Two variants of the code are provided. These two versions differ in how particles are routed through the computational domain. In MARFA 3.2.3, transport is assumed to occur along a set of trajectories or pathways that originate at radionuclide source locations. The trajectories are intended to represent the movement of hypothetical, advectively transported groundwater tracers and are typically calculated by pathline tracing in a discrete fracture network flow code. The groundwater speed and retention properties along each pathway may change in time, but the pathway trajectories are fixed. MARFA 3.3.1 allows the transport effects of changing flow directions to be represented by abandoning the fixed pathways and performing node routing within MARFA. (orig.)