Skip computation of collision sources of an ion species with itself

trixi-framework · Dec 18, 2024 · 5c772f6 · 5c772f6
1 parent 5dabad1
commit 5c772f6
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/equations/ideal_glm_mhd_multiion.jl b/src/equations/ideal_glm_mhd_multiion.jl
@@ -334,6 +334,9 @@ function source_terms_collision_ion_ion(u, x, t,
         S_q3 = zero(eltype(u))
         S_E = zero(eltype(u))
         for l in eachcomponent(equations)
+            # Do not compute collisions of an ion species with itself
+            k == l && continue
+
             rho_l, v1_l, v2_l, v3_l, p_l = get_component(l, prim, equations)
             T_l = p_l / (rho_l * gas_constants[l])