I am interested in modeling a heat equation problem with a Gaussian heat source, so I need a way to define a source as a function of X and Y. A search turned up a previous 1D case from the older version in which something like this was done:
X = mesh.getCellCenters()[0] eq = (TransientTerm() == DiffusionTerm() + X) However, when I try to do this I get the following error: "TypeError: The coefficient must be rank 0 for a rank 0 solution variable." Did 3.0 change the way this works, or is the 2D case different, or there a problem somewhere else that I'm missing? Thanks, Adam _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
