New question #152565 on DOLFIN: https://answers.launchpad.net/dolfin/+question/152565
I'm trying to solve a standard electromagnetic problem to test my code. The problem is an infintesimal electric dipole. The linear form involving current would look something like: V = FunctionSpace(mesh, "Nedelec 1st kind H(curl)", order) u = TrialFunction(V) v = TestFunction(V) L = dot(j*k0*Z0*J_imp, v)*dx where (ignoring the difficulties of the imaginary part) j*k0*Z0 is essentially a constant, J_imp is a function that represents the impressed current and v is the testing function. This is fine if J_imp is a volumetric current distribution that can easily be represented as a function in V. An infintesimal electric dipole is modelled as a current point source. In other words J_imp is a spatial dirac delta function. This kind of source is also used to excite some antenna models, which is what I'm working towards. How could one get dolfin to integrate such a dirac-delta function? Thanks Neilen -- You received this question notification because you are a member of DOLFIN Team, which is an answer contact for DOLFIN. _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp

