Diffuse interface modeling of non‐isothermal Stokes‐Darcy flow with immersed transmissibility conditions

The coupling between free and porous medium flows has received significant attention since it plays an important role in a wide range of problems from fluid‐soil interactions to biofluid dynamics. However, modeling this coupled process remains a difficult task as it often involves a domain decomposition algorithm in conjunction with a special treatment at the interface. The problem can become more challenging under non‐isothermal conditions because it requires the iterative procedure at every time step to simultaneously meet the transient mass continuity, force equilibrium, and energy balance for the entire system. This article presents a diffuse interface framework for modeling non‐isothermal Stokes‐Darcy flow and the corresponding finite element formulation that bypasses the need for explicitly splitting the domain into two, which enables the unified treatment for distinct regions with different hydrothermal flow regimes. To achieve this goal, we employ the Allen‐Cahn type phase field model to generate the diffuse geometry, where the solution field can be seen as a regularized approximation of the Heaviside indicator function, allowing us to transfer the interface conditions into a set of immersed boundary conditions. Our formulation suggests that the isothermal operator splitting strategy can be adopted without compromising accuracy if the heat and mass transfer processes are decoupled by assuming that the density and viscosity of the phase constituents are independent to the temperature. Numerical examples are also introduced to verify the implementation and to demonstrate the model capacity.