Hello, In general, defining a singular problem (such as a point source - dirac delta) may lead to unexpected numerical results. An alternative may be to define a circle geometry with a hole at the arbitrary location that you wish. The hole may be small but not infinetesimal. The boundary of the hole being set to a fixed value. You may use gmsh to generate the mesh.
Pierre On Fri, Mar 01, 2013 at 12:27:21PM -0500, Raymond Smith wrote: > Hi Adrian, > > I'm not sure if this is what you're going for, but perhaps you could use an > implicit source to apply the "internal boundary condition" as outlined in > http://www.ctcms.nist.gov/fipy/documentation/USAGE.html#applying-internal-boundary-conditions. > If you're just trying to keep a particular point at a fixed value, maybe > that will work? > > Cheers, > Ray > > > On Fri, Mar 1, 2013 at 12:11 PM, Adrian Jacobo <[email protected] > > wrote: > > > Hi All, > > > > I'm trying to solve the diffusion equation with a point source at an > > arbitrary location inside a circle. My idea is to adapt the example > > examples.diffusion.circle<http://www.ctcms.nist.gov/fipy/examples/diffusion/generated/examples.diffusion.circle.html#module-examples.diffusion.circle>but > > I'm not sure how to proceed. > > I think the idea should be create a variable which is zero everywhere > > except for where the point source is located, and add this variable as the > > source. But I don't know how to do this. > > > > Any ideas or help will be greatly appreciated. > > > > Best, > > Adrian. > > > > > > _______________________________________________ > > fipy mailing list > > [email protected] > > http://www.ctcms.nist.gov/fipy > > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > > > > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] -- ----------------------------------------------------------- Pierre de Buyl Université libre de Bruxelles Physique des Systèmes Complexes et Mécanique Statistique Membre du conseil d'administration - ULB EuroSciPy 2013 co-chair http://www.euroscipy.org/ tél: +32 (0)2 650 57 92 web: http://homepages.ulb.ac.be/~pdebuyl/ ----------------------------------------------------------- _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
