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 ]

Reply via email to