Hallo Markus,
Orientation vector is defined based on experiment you want to reproduce.
For instance, modeling any finite length dipole of specified power can
be done as:
/*f*(*p*) = a * *n* * \delta (*p* - p_s);/
*n* = {n_x n_y n_z} -- orientation vector;
a = I * L; I -- current, L -- length
Same can be done for magnetic dipoles (antennas).
So this is quite convenient and many real sources can be well
approximated by this vector function.
However, I'm not sure now that the dot product of orientation vector and
shape function is correct representation.
The dot product gives us scalar, whereas we need vector with each
component assigned to individual field components of FE_Nedelec. Is this
still correct?
Thanks.
On 25.04.2012 18:32, Markus Bürg wrote:
Hello Alexander,
then you are right and inserting the dot product of orientation vector
and shape function should work. Where do you get the orientation
vector from? Should this be another argument of the function?
Best Regards,
Markus
On 25.04.2012 15:21, Alexander Grayver wrote:
Markus,
My previous email wasn't clear.
For EM, we can have some vector that defines orientation of the unit
dipole, then source is defined as following
/f(p) = n * \delta (p - p_s);/ n = {n_x n_y n_z} -- orientation vector
Elements of this vector should then be assigned to corresponding
cell's components (i.e. E{x,y,z}).
For instance, x-directed unit dipole would give us 3 component vector
with only 1st element == 1.
--
Regards,
Alexander
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
