From 50f60f798aec36d5854829f7ad41588b9d8564ee Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <schaller@strw.leidenuniv.nl>
Date: Sat, 18 May 2019 17:48:38 +0100
Subject: [PATCH] When doing stellar or AGN feedback, set the particles'
viscosity to the maximal value.
---
src/black_holes/EAGLE/black_holes_iact.h | 3 +++
src/feedback/EAGLE/feedback_iact.h | 3 +++
2 files changed, 6 insertions(+)
diff --git a/src/black_holes/EAGLE/black_holes_iact.h b/src/black_holes/EAGLE/black_holes_iact.h
index 3174dea57a..0db85b3924 100644
--- a/src/black_holes/EAGLE/black_holes_iact.h
+++ b/src/black_holes/EAGLE/black_holes_iact.h
@@ -134,6 +134,9 @@ runner_iact_nonsym_bh_feedback(const float r2, const float *dx, const float hi,
hydro_set_physical_internal_energy(pj, xpj, cosmo, u_new);
hydro_set_drifted_physical_internal_energy(pj, cosmo, u_new);
+
+ /* Impose maximal viscosity */
+ hydro_set_viscosity_alpha_max_feedback(pj);
}
}
diff --git a/src/feedback/EAGLE/feedback_iact.h b/src/feedback/EAGLE/feedback_iact.h
index f1096beec2..4ad145e77d 100644
--- a/src/feedback/EAGLE/feedback_iact.h
+++ b/src/feedback/EAGLE/feedback_iact.h
@@ -300,6 +300,9 @@ runner_iact_nonsym_feedback_apply(const float r2, const float *dx,
hydro_set_physical_internal_energy(pj, xpj, cosmo, u_new);
hydro_set_drifted_physical_internal_energy(pj, cosmo, u_new);
+ /* Impose maximal viscosity */
+ hydro_set_viscosity_alpha_max_feedback(pj);
+
/* message( */
/* "We did some heating! id %llu star id %llu probability %.5e " */
/* "random_num %.5e du %.5e du/ini %.5e", */
--
GitLab