Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

多孔介质的干燥是由多孔结构内液相的蒸发和运动共同控制的。表面粗糙度引起的接触角滞后会影响多相流，例如液滴的接触线运动、排水过程中的相分布以及恒定接触半径模式下液滴干燥后形成的咖啡环。然而，接触角滞后对多孔介质中液体干燥的影响仍然是一个悬而未决的问题。格子玻尔兹曼模型 (LBM) 是一种先进的数值方法，越来越多地用于研究包括干燥在内的相变问题。在本文中，基于规定接触角的几何公式​​方案，我们在两相赝势LBM框架内实现了接触角磁滞模型。通过模拟位于平面和曲面上的液滴，测试和验证了在大范围内指定和自动测量接触角的能力和准确性。随后，通过模拟平面和曲面上的液滴干燥来验证所提出的接触角滞后模型。然后，研究了两个连接毛细管的干燥，考虑了不同接触角滞后范围对干燥动力学的影响。最后，该模型应用于研究双孔隙多孔介质的干燥，其中比较了有和没有接触角滞后的相分布和干燥速率。所提出的模型被证明能够准确处理不同的接触角滞后范围，并捕获不同多孔介质（包括平面和弯曲几何形状）干燥过程中的物理机制。