Time-dependent artificial viscosity

From: Leonard Romano <leonard.romano_at_tum.de>
Date: Mon, 26 Jul 2021 22:04:52 +0200

Dear gadget-list members,

I noticed some possibly major issues with the current implementation of
time-dependent artificial viscosity:

 1. The decay velocity is computed during the density loop, using Csnd
    even though Csnd depends on the density that is still to be computed
 2. Alpha is integrated during the density loop, which in particular
    means that it will be integrated in time as many times as loops are
    necessary to reach a good density estimate. I think that this should
    be moved somewhere outside the density loop to ensure correct

I am worried that with these issues, it might not work as intended. I
suggest for the decay velocity to compute it during the force loop and
integrating Alpha during the kick operations.


Leonard Romano, B.Sc.(レオナルド・ロマノ)
Physics Department
Technical University of Munich (TUM), Germany
Theoretical Astrophysics Group
Department of Earth and Space Science
Graduate School of Science, Osaka University, Japan
he / him / his
Received on 2021-07-26 22:04:57

This archive was generated by hypermail 2.3.0 : 2023-01-10 10:01:32 CET