Hi Mike Thanks for your answer and sorry about my long response time. I think this source idea is actually what I need.
What do you mean by rank-0 variable? I have seen the term (and rank-1) in the fipy manual but not any explanation as to its meaning. Is there a way (during the solve loop) of observing the error? Thanks David On Thu, Jun 4, 2009 at 4:23 PM, Mike Atkins <[email protected]> wrote: > I suppose I am looking for a spatial and time dependent source, which >> can represent the inflow of heat energy over my specific time interval >> I did notice the reference to sources on pages 52,53 of the manual but I >> do not know how >> to extrapolate this information to my specific problem. >> > > In order to apply this method, you should only need to make a few changes > to your second attempt. You can try something like this: > > x, y = m2D.getCellCenters() > cv = CellVariable(mesh=m2D, value = initialTemp) > heater = CellVariable(mesh=m2D) > heater.setValue(4000/(1+(x-0.5)**2+(y-0.5)**2)**2,where=(x-0.5)**2 +\ > (y-0.5)**2 < 0.04) > > You'll Also need to change the equation to incorporate the heater term: > > eq = TransientTerm() == ImplicitDiffusionTerm(coeff=D)+(t<.03)*heater > > The way I outlined above is just one of many ways you can incorporate a > source term with a rank-0 CellVariable to describe the heating. > > -Mike Atkins > > -- David Wende home +972-8-9353488 work +972-2-5886116 mobile +972-54-234-6479
