On Jan 16, 2013, at 6:35 AM, John Assael wrote: > Thank you very much for your help! > Do you know how can I set only a point (a small circle) as the heat source > and not the whole side?
examples.diffusion.mesh1D, examples.diffusion.electrostatics, and examples.phase.anisotropy all show how to set a CellVariable to have a value in only one region. You can either use Q = CellVariable(..., value=0.) Q.setValue(Q0, where=(x - cx)**2 + (y - cy)**2 <= r) or Q = Q0 * ((x - cx)**2 + (y - cy)**2 <= r) Note: due to http://matforge.org/fipy/ticket/447, be sure to use x, y = mesh.cellCenters[0], mesh.cellCenters[1] _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
