I constructed a simple 2-D model with an elastic crust, an elastic loading body and a viscoelastic mantle as the picture below shows.
My target is to see the subsidence rate of the surface under the locale loading body after a half million years. However, when I tried to insert gravity following section 7.17.4 of the manual, the velocity output
of the domain turns out to be asymmetrical and strange under a symmetrical loading after 20 000 years. The figures below show the velocity results after 15 000 year’s simulation and 500 000 year’s simulation, respectively.
Do you know what’s wrong with the simulation or how can I fulfill my target reasonably?