Re: thermal contact between two regions

2018-07-23 Thread Daniel Wheeler
On Mon, Jul 23, 2018 at 1:06 PM, Drew Davidson 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 >

Re: thermal contact between two regions

2018-07-23 Thread Drew Davidson
Hello Dr. Wheeler, '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

Re: thermal contact between two regions

2018-07-23 Thread Daniel Wheeler
On Tue, Jul 17, 2018 at 4:43 PM, Drew Davidson 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

Re: thermal contact between two regions

2018-07-17 Thread Drew Davidson
Hello, I continued to stare at https://www.mail-archive.com/fipy@nist.gov/msg02641.html. I realized Keff must be effective thermal conductivity as described in: Salazar, Agust n. “On Thermal Diffusivity.” *European Journal of Physics* 24, no. 4 (July 1, 2003): 351–58.

Re: thermal contact between two regions

2018-07-12 Thread Daniel Wheeler
I think you can do this by aligning the mesh with the contact region and specifying the diffusion coefficient at the contact faces. On Thu, Jul 12, 2018 at 12:28 PM, Drew Davidson wrote: > Hello, > > > I’m sorry for not being more specific. I was trying to incorporate >

Re: thermal contact between two regions

2018-07-12 Thread Drew Davidson
Hello, I’m sorry for not being more specific. I was trying to incorporate https://www.mail-archive.com/fipy@nist.gov/msg02626.html by reference. Sketch of problem: (1) Region 1 --* thermal contact resistance * -- Region 2 (2) Temperature at (1)

Re: thermal contact between two regions

2018-07-12 Thread Daniel Wheeler
On Wed, Jul 11, 2018 at 7:00 PM, Drew Davidson wrote: > > 1. Transient heat conduction rather than steady state. Same boundary > conditions but initial condition can be zero temperature eveywhere. Add a TransientTerm to the equation and step through the solution with time steps and sweeps if

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

2013-05-17 Thread Jonathan Guyer
On May 16, 2013, at 1:23 PM, Daniel Wheeler daniel.wheel...@gmail.com wrote: I believe the correct diffusion coefficient across the face that approximates a thin layer is D = Kc * (1 / (epsilon + (1 - epsilon) * D_ratio)) where D_ratio = Kc / 2 * (1 / K1 + 1 / K2) and

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

2013-05-16 Thread Frederic Durodie
Hi, Thanks for your proposed 1D solution. It bears part of the correct solution which should look like : from X=-L/2 the temperature increases to 2.5 deg C at X = 0; then there is a temperature jump to 3.5 deg C and the temperature increases linearly to 6 deg C at the X=L/2 boundary. So I'm a

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

2013-05-16 Thread Jonathan Guyer
On May 16, 2013, at 6:39 AM, Frederic Durodie frederic.duro...@googlemail.com wrote: So I'm a bit puzzled as to why for X from -L/2 to 0 the temperature is 0. : it should increase linearly to 2.5 deg C. I suspect the problem is that divergence of the flux at the cell just left of the

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

2013-05-16 Thread Daniel Wheeler
On Fri, May 10, 2013 at 2:11 AM, Frederic Durodie frederic.duro...@googlemail.com wrote: 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

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

2013-05-13 Thread Jonathan Guyer
On May 10, 2013, at 2:11 AM, Frederic Durodie frederic.duro...@googlemail.com wrote: 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. I don't know whether the

how to implement a thermal contact between two regions ?

2013-05-10 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