“…[]. Among the methods reviewed, three are most relevant to this work: a weighted residual method based on mortar spaces [ Balhoff et al ., ; Mehmani et al ., ; Sun et al ., ], an iterative approach with a “handshake” domain [ Battiato et al ., ], and a semianalytical coupling technique that involves numerical solution of a global system of discretized equations that includes both domains [ Roubinet and Tartakovsky , ]. The mortar method utilizes a finite element space of reduced dimension and a Lagrange multiplier approach to compute matching boundary conditions for the adjacent pore‐scale and continuum‐scale regions in which computations are simultaneously performed [ Arbogast et al ., ; Peszynska et al ., ].…”