Re: how to implement a thermal contact between two regions ?

2013-05-16 Thread Frederic Durodie
'rgbcm'[sweep % 5]+['-','-.',':','--'][(sweep/5)%4], label = 'sweep %d' % (sweep+1)) # plot the theoretical solution pl.plot([-L/2,0.,0.,L/2], [0.,Ps/K1*L/2,Ps/K1*L/2+Ps/hc,Ps/K1*L/2+Ps/hc+Ps/K1*L/2], 'k

how to implement a thermal contact between two regions ?

2013-05-09 Thread Frederic Durodie
Dear FiPy users and developers, could you help me with how to implement a thermal contact, hc [W/(m^2.K)], between two regions : so there is like a discontinuity in the temperature between the two regions. In some cases I could implement it as a thin layer, dx [m], with a given thermal diffusion

Re: unexpected result using var.getFacGrad().constrain(|gradvar|, where=boundary)

2013-04-30 Thread Frederic Durodie
Thank you very much ! This indeed resolved the issues. Kind regards, Frederic ___ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]

unexpected result using var.getFacGrad().constrain(|gradvar|,where=boundary)

2013-04-29 Thread Frederic Durodie
In a simple 2D heat conduction problem I try to impose a power flux on part of the outer boundary and a fixed temperature on the remaining part. However the resulting temperature distribution is not what I expect and furthermore depends on how the problem geometry is rotated (alphadeg : I did this