An incident flux expansion transport theory method suitable for coupling to diffusion theory methods in hexagonal geometry

► A new hybrid stochastic-deterministic transport theory method to couple with diffusion theory. ► The method is implemented in 2D hexagonal geometry. ► The new method produces excellent results when compared with Monte Carlo reference solutions. ► The method is fast, solving all test cases in less than 12 s.