On Jan 17, 2013, at 12:58 PM, John Assael wrote: > from fipy import * > import time > > dx = dy = 10 > nx = ny = 100 > L = dx * nx > mesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny) > > # Fluid Properties > tc0=0.141 > den0=867. > cp0=1717 > > # Wire Properties > tcw=57.5 > denw=16690. > cpw=135. > Qw=1.5*10**9
> However the change of Temperature only takes place within the radius of the > wire and do not seem to affect the outer fluid. Do you know how can I fix > this? The scales of your parameters are completely different. You are supplying 1.5e9 W/m**3 in the wire, but the diffusion term only dissipates 5.75 W/m**3 for a 10 K temperature difference (k dT / dx**2) on the wire side; it's even smaller in the fluid. You either need terms of the same order or you need to accept that diffusion is not important at the rate that Q is pumping energy into your wire. _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
