Dear all,
I developed a coupled HM model for porous media with many ODE interfaces. I found the model ran very slow. I am wondering whether these ODE interfaces contribute to the Jacobian matrix. Is it possible to neglect these contributions while keeping the accuracy of the coupled model. Do you have any ideas to deal with such a situation?
Thank you so much for your help.