Two slabs representing the fluid and solid domain are connected via an interface (see Figure , top left). The dimensions of both slabs are identical and set to 8.0 mm × 2.0 mm × 10.5 mm (Height x Width x Length). Energy and the two model species A and B are only transported via diffusion. Dirichlet boundary conditions for the temperature (T surf,s = 500 K, T surf,f = 350 K) and the mass fractions (w surf,s,A = 0.00, w surf,f,A = 0.01) are applied at the opposing surfaces of the two slabs. The remaining walls are set as adiabatic and impermeable, respectively. The thermal conductivity and the diffusion coefficients differ between the solid (λ s = 0.6 W m -1 K -1 , D s,A,B = 5 × 10 -6 m 2 s -1 ) and the fluid (λ f = 0.3 W m -1 K -1 , D f,A,B = 10 × 10 -6 m 2 s -1 ) phase, but are constant within the respective domain. Moreover, the density is considered as invariable and chemical reaction is excluded, so that the simulation can be compared with the analytical solution of the problem described in the supplementary information (SI Section S6). A structured mesh (total cell count: 4.28×10 4 ) containing only orthogonal cells and an unstructured mesh (total cell count: 2.66 × 10 4 ) with polyhedral cells is used for the numerical calculation.