I am trying to implement a 3D Cartesian exponentially varying viscosity benchmark on ASPECT. I did get the code running. However, my velocity and pressure L2 norm stay almost the same with increasing refinements (only uniform applied). I am wondering if anyone has any suggestion. Here I also attach some of my data.
refinement Velocity L2 norm Pressure L2 norm
3 2.832740e-02 1.420663e+04
4 2.271762e-02 1.420674e+04
5 2.261151e-02 1.420675e+04
6 2.260930e-02 1.420675e+04
I have a feeling that if my implementations of benchmark and boundary condition are correct, then the velocity L2 norm could be due to the solver while the pressure L2 norm could be due to the pressure normalization.
Thanks in advance,