Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-09-25 Thread Raymond Smith
ut with a >>>> thermal diffusivity is simply not the correct representation of your >>>> physical situation. The actual conservation equation for thermal energy is >>>> (neglecting convection and source terms) >>>> > > >>>> > >

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-09-25 Thread Daniel Wheeler
; from the correct conserved quantities. I would guess that if you run it >>> with any rho, c_p that are the same in the concrete and insulator, you will >>> get the correct result with the "alpha form" because, being both constant >>> and uniform, they can be swept int

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Raymond Smith
than unity. I have realised that the >> equation should be specified as >> > > >> > > >> > > eqn = TransientTerm(coeff=C) == DiffusionTerm(coeff=k) >> > > >> > > >> > > where the volumetric heat capacity is C = rh

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Raymond Smith
ho*c_p) has a value other than unity. I have realised that the > equation should be specified as > > > > > > > > > eqn = TransientTerm(coeff=C) == DiffusionTerm(coeff=k) > > > > > > > > > where the volumetric heat capacity is C = rho * c_p

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Guyer, Jonathan E. Dr.
thermal > > conductivity. For a steady case, this equation reduces to > > 'DiffusionTerm(coeff=k)==0' and gives a correct result. > > > > > > I have attached a figure with the updated comparison of FiPy and TEMP/W for > > an insulated concrete slab, show

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Conor Fleming
iffusionTerm respectively. > > > Many thanks, > > Conor > > > From: fipy-boun...@nist.gov<mailto:fipy-boun...@nist.gov> > [fipy-boun...@nist.gov<mailto:fipy-boun...@nist.gov>] on behalf of Conor > Fleming [conor.flem...@eng.ox.ac.uk<mailto:conor.flem

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Raymond Smith
to > add a note to the documentation page 'examples.diffusion.mesh1D', > explaining that for some applications, e.g. the heat equation, it is > appropriate to separate the (thermal) diffusivity into two portions, which > act on the TransientTerm and DiffusionTerm respectively. >

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Guyer, Jonathan E. Dr.
, explaining that > for some applications, e.g. the heat equation, it is appropriate to separate > the (thermal) diffusivity into two portions, which act on the TransientTerm > and DiffusionTerm respectively. > > > Many thanks, > > Conor > > > From: fipy-bou

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Raymond Smith
nd DiffusionTerm respectively. > > > > Many thanks, > > Conor > > > ------------------ > *From:* fipy-boun...@nist.gov [fipy-boun...@nist.gov] on behalf of Conor > Fleming [conor.flem...@eng.ox.ac.uk] > *Sent:* 08 August 2014 17:24 > *To:* fipy@nist.gov &

RE: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-18 Thread Conor Fleming
uk] Sent: 08 August 2014 17:24 To: fipy@nist.gov Subject: RE: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient Hi Kris, Thank you for the prompt response. You are right - altering the insulation conductivity in the FiPy m

RE: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-08 Thread Conor Fleming
@nist.gov Subject: Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient Conor, if you reduce the thermal conductivity in the insulation to about 0.1, the fipy solution looks about like the other model (the knee in T is about at

Re: Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-08 Thread Kris Kuhlman
Conor, if you reduce the thermal conductivity in the insulation to about 0.1, the fipy solution looks about like the other model (the knee in T is about at 400 degrees C). Is there an issue with how your compute or specify this in fipy or the other model? Kris On Fri, Aug 8, 2014 at 9:35 AM, C

Unexpected result, possibly wrong, result solving 1D unsteady heat equation with spatially-varying diffusion coefficient

2014-08-08 Thread Conor Fleming
Hi, I am using FiPy to determine the depth of heat penetration into concrete structures due to fire over a certain period of time. I am solving the unsteady heat equation on a 1D grid, and modelling various scenarios, e.g. time-dependent temperature boundary condition, temperature-dependent dif