On Mon, Jul 23, 2018 at 1:06 PM, Drew Davidson <davidson...@gmail.com> wrote: > > 'dx' is actually the thickness of the thermal contact region. It is not the > volume of the cell. If it were a volume, the ratio dx/mesh._cellDistances > would not be dimensionless. I was seeking a thermal contact region of zero > thickness (dx=0), giving a true thermal contact resistance.
Sorry, my mistake on that one. > I let this thread drop because the Salazar paper I referenced gave me > expressions for effective thermal conductivity and effective thermal > diffusivity, and this allowed me to arrive at a reasonable-looking FiPy > result for transient thermal behavior. My result still deviates a bit from > the exact solution at steady state, and this deviation persists despite mesh > refinement, but I am busy with other things and have called it good enough > for now. Thanks for following up and sharing your code. -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]