Hydrostatic equilibrium boundary conditon

[kyriakosTapinou]

I've been trying to establish a hydrostatic equilibrium with constant gravity. I've used the Rayliegh-Taylor example as a starting point. Using the same style of pressure boundary condition where the ghost cell pressure is integrated from the domain pressure value is leading to pressure waves emitted on the boundary. Do you have any suggestions or resources you can point me to?

In the current case it is all hydrodynamic and there are no fields. My pressure and density profile are not linear throughout the domain.

[tomidakn]

I guess the pressure waves are emitted because the boundary conditions are not sufficiently accurate and the ghost cells are (slightly) out of equilibrium. You can quantify it by measuring the pressure gradient force between the first active cell and first ghost cell. I suggest you to improve the accuracy of the extrapolation in the calculation of the ghost cell values, as the current one implemented in the RT pgen does not look very accurate.

Alternatively, you may directly apply boundary conditions on the flux. However, you need to modify the code by yourself as we do not provide such an interface at this moment. The fluxes from the Riemann solver are stored in the x(1/2/3)flux arrays, and you can modify them before they are applied in Hydro::AddFluxDivergence.

[kyriakosTapinou]

Thanks for the response - I'm trying to implement the following remedy now but I think it is something you would deem to specific to be introduced into the main branch.

https://www.aanda.org/articles/aa/abs/2019/11/aa36387-19/aa36387-19.html