First thing I would look at is convergence on each time step. Is the model iterating and converging? It sounds like it is not. This could be because your convergence settings are too loose or too tight, or because your discretization is too coarse.

Edited January 31, 2017 by woodward

Share this post

Link to post

Share on other sites

I have checked and the numerical process is converging! I tried to change the Run option and the iteration parameters, I also converted the heterogeneous materials in homogeneous ones: the heads remain as in the initial condition with the exception of the side where the transient boundary conditions are applied.

Yes the unsaturated curves can cause problems. I also recommend NOT using the spline interpolation, as sometimes this results in a non-monotonic curve. Be careful about the curve spacing and how far away from the water table you specify it (as the last value will be used for all higher tensions as well).