“…With this method it is possible to incorporate complex spatial structures in the simulation domain. The features of the method range from single phase, single component flow [16,17] to multi component [18,19,20,21] and/or multi phase configurations [18,22,23,24]. Sophisticated data structures and algorithms can be included [25,26] as well as extended boundary conditions [27,22,28].…”