Maybe that is not the best place for the following question, if such, let me know and sorry.

I would like to model the flow of expanding bentonite in a fracture following the given reference. In order to take into account the swelling pressure, they modified the Stokes equation (momentum equation) by adding the swelling pressure explicitly (equation 4-6 of reference): \begin{equation} \frac{\partial \bf{u}}{\partial t} -\nabla \cdot \bf{\tau} + \nabla (p+p_{sw}) = 0 \\ \end{equation}

After that they add the static friction at the wall (equation 4-7 in reference): \begin{equation} \frac{\partial \bf{u}}{\partial t} -\nabla \cdot \bf{\tau} + \nabla (p+p_{sw}(1-2\mu)) = 0 \\ \end{equation}

the friction coefficient is \begin{equation} \mu \end{equation} I wonder if such assumption in a 3D model would be correct. I assume that if I solve the system with the last equation, the nodes that are not close to the wall are also going to take into account the friction coefficient. I think it has more sense to add the Friction force into the wall as a boundary condition. I thought to use a navier slip condition.

Does it make sense what I am saying, or I am wrong.


Best Regards



Pont, A., Coene, E., & Idiart, A. (2020). Bentonite erosion project. Svensk kärnbränslehantering AB [Swedish Nuclear Fuel and Waste Management Company]}



Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.