In this work, three variant Time Domain Random Walk (TDRW) algorithms were developed for the problem of solute transport in a single fracture-matrix system, where an arbitrary inlet boundary condition can be applied. One approach performs an additional evaluation of integral in terms of injection boundary and the solution to a Dirac delta case. One method makes use of two functions, dependent on specific boundary conditions, to estimate the particle arrival time. The other additionally introduces the concept of solute injection time, resulting from the injection boundary, into the calculation ...