“…Indeed, a global γ = 1 is often a default in artificial compressibility studies (e.g., [30,56,122,123]), although higher values are also found (e.g., [62,78,124]). A wide range of studies fail to document the values used for γ; likely because it formally controls convergence rates rather than accuracy (e.g., [73,77,81,83,120])-although it can also affect the simulation result for upwind schemes [124]. Methods have been developed to select an optimum global γ based on flow characteristics and/or model time step (e.g., [63,64,[125][126][127]), or using mesh spacing [72], or as a preconditioning parameter for a solution matrix [128].…”