On Tue, Jul 17, 2018 at 4:43 PM, Drew Davidson <davidson...@gmail.com> wrote:
>
> I still need FiPy code for:
>
>
> dAP1/(dAP1+dAP2)
>
>
> and
>
>
> dAP2/(dAP1+dAP2)
>
>
> where dAP1 and dAP2 were distances from cell center to cell face for cells
> on either side of the interface. If these FiPy expressions are unavailable,
> I would think assuming .5 is going to be OK...
I think it's outlined in the link that you gave.
Kcontact = hc * mesh._cellDistances
Kavg = Kcell.harmonicFaceValue
K.setValue(Kcontact * Kavg / (Kavg + Kcontact * (1 - dx /
mesh._cellDistances)), where=mesh.physicalFaces['thermal contact'])
m._cellDistances is the distance between each cell and `dx` is the
volume of the cell for a 1D problem.
--
Daniel Wheeler
_______________________________________________
fipy mailing list
fipy@nist.gov
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
