Measuring the Effects of Artificial Viscosity in SPH Simulations of Rotating Fluid Flows

TL;DR

This paper introduces a practical, accurate local method to evaluate artificial viscosity effects in SPH simulations of rotating fluids, improving upon previous approaches by providing detailed, position-dependent analysis.

Contribution

The study develops a simple, broadly applicable technique based on entropy creation to assess artificial viscosity effects locally in SPH simulations, surpassing prior methods in accuracy and versatility.

Findings

The new method accurately characterizes viscous effects across the disc.

It reveals the limitations of previous approximate methods.

Explicit dependencies of effective viscosity on flow parameters are estimated.

Abstract

A commonly cited drawback of SPH is the introduction of spurious shear viscosity by the artificial viscosity term in situations involving rotation. Existing approaches for quantifying its effect include approximate analytic formulae and disc-averaged be- haviour in specific ring-spreading simulations, based on the kinematic effects produced by the artificial viscosity. These methods have disadvantages, in that they typically are applicable to a very small range of physical scenarios, have a large number of simplifying assumptions, and often are tied to specific SPH formulations which do not include corrective (e.g., Balsara) or time-dependent viscosity terms. In this study we have developed a simple, generally applicable and practical technique for evaluating the local effect of artificial viscosity directly from the creation of specific entropy for each SPH particle. This local approach is simple and quick to implement, and it al- lows a detailed characterization of viscous effects as a function of position. Several advantages of this method are discussed, including its ease in evaluation, its greater accuracy and its broad applicability. In order to compare this new method with ex- isting ones, simple disc flow examples are used. Even in these basic cases, the very roughly approximate nature of the previous methods is shown. Our local method pro- vides a detailed description of the effects of the artificial viscosity throughout the disc, even for extended examples which implement Balsara corrections. As a further use of this approach, explicit dependencies of the effective viscosity in terms of SPH and flow parameters are estimated from the example cases. In an appendix, a method for the initial placement of SPH particles is discussed which is very effective in reducing numerical fluctuations.