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

Reply via email to