Spherical benchmarks: mismatch between analytical and computed pressure values

[gthyagi]

Hi @knepley,

I am running spherical benchmarks from the following papers:

1. Paper1
2. Paper2

I successfully reproduced the benchmark results from the first paper, but I'm encountering difficulties with the second. The numerical pressure values deviate significantly from the analytical ones. The primary difference between the benchmarks in these two papers lies in the initial conditions, specifically the magnitude of the velocity on the outer surface and the magnitude of the density difference.

For both benchmarks, I am using identical solver settings.

Analytical (left) and numerical (right) pressure solution:

Analytical (left) and numerical (right) velocity solution:

Here are the scripts I am using:

Paper 1: Benchmark script Paper 2: Benchmark script