I'm modeling ice flow using a user defined creep strain rate tensor. I'm getting close to the expected velocities but the stresses are much too high and have an odd distribution that I think indicates some numerical error (see attached screenshot of the xy component of the deviatoric stress tensor). I've tried a few different solvers (MUMPS and PARADISO are the ones I've tried that solve in a reasonable time) and I've tried decreasing the mesh size and have the same problem. I'd really appreciate any ideas on what to try to fix this. I've attached the mph file.
↧