From 34573d7d63593932cd7933755f9ee94790c470f4 Mon Sep 17 00:00:00 2001 From: JackHsu Date: Tue, 30 Jul 2019 13:32:44 +0800 Subject: [PATCH] Corrected acciden removal. --- acc/mhd_solver/stencil_process.sps | 3 +++ 1 file changed, 3 insertions(+) diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps index 4b8bae9..a9e9747 100644 --- a/acc/mhd_solver/stencil_process.sps +++ b/acc/mhd_solver/stencil_process.sps @@ -181,6 +181,9 @@ heat_conduction( in Scalar ss, in Scalar lnrho) { (gamma - AcReal(1.)) * laplace(lnrho); const Vector second_term = gamma * inv_cp_sound * gradient(ss) + (gamma - AcReal(1.)) * gradient(lnrho); + const Vector third_term = gamma * (inv_cp_sound * gradient(ss) + + gradient(lnrho)) + grad_ln_chi; + const Scalar chi = AC_THERMAL_CONDUCTIVITY / (exp(value(lnrho)) * cp_sound); return cp_sound * chi * (first_term + dot(second_term, third_term));