“…α f,s = k f,s /(ρc) f,s , t * is time, µ is the dynamic viscosity, K is permeability, ρ is density, g is the modulus of the gravitational acceleration g, β is the thermal expansion coecient, T 0 is the reference temperature, e z is the unit vector for the direction z, ϕ is porosity, and c is the heat capacity per unit mass. The set of boundary conditions can be expressed as [11,12,13] …”