“…(ii) The evolution of liquid surface, y = h ( x , z , t ), is described by the kinematic boundary condition [25] ∂th+bold-italicu⋅bold∇h=bold-italicu⋅bold-italicjboldfalse^. (iii) Following Beavers & Joseph [26], an empirical boundary condition for shear stress is imposed at the liquid–porous interface, y = − d , ∂yuH=(αBJκ)false(uH−upHfalse),1emv=vp, p=pp, where α BJ is an empirical dimensionless coefficient, which, in fact, depends on the local structure of the porous material close to the liquid–porous interface [26,27]. In fact, this empirical boundary condition (2.8) is considered to take into account the jump of velocity in liquid and porous layers due to the boundary layer developed by the viscous friction effect in the vicinity of the liquid–porous interface.…”