[Chama-developers] gamma

[John A. T. <tu...@la...>]

so I happened to be talking with Doug tonight as he was leaving,
and I mentioned the fact that we were finding that the method for
computing the face conductivities in T was pretty bad, and he
responded that it was surprising because it was just a harmonic
average, right?

well, no - so I showed him where we were looking in the code to
find that it was using what we're calling the "gamma" scheme:

    k1*k2 (V1 + V2)                / V2   V1 \ -1
    ---------------  =  (V1 + V2) |  -- + --  |
     k1*V1 + k2*V2                 \ D1   D2 /

and he said "well, looks like a bug - and it's probably mine"

interesting!  that's been in there for a while!

he pointed out that we have some two-material conduction stuff in
the test suite that they thought was working, so I'll look into
that...

-JT