I am trying to develop magnetodynamics solver for 2D case (x-y plane), with assumed current density J and magnetic vector potential A in the x-y plane, ie. with null z-component (Jz=0 and Az=0)
In the Elmer ModelsManual document, there is sentence (Bold) from the 15.1. - Introduction paragraph witch concerns me:
Does this means that edge elements in some sense are not yet fully implemented for 2D elements of 2D mesh?"This module may be used to solve a version of the Maxwell equations in the A-V form. The approximation
of the associated vector potential variable ~A is here done by using edge element basis functions, while the
classic Lagrange interpolation is applied to compute the scalar potential V . The use of edge elements LIMITS
the applicability of this solver to 3D problems. In addition to performing the computations in the time
domain, the analogous version of the equations may also be solved in the frequency domain. Furthermore,
an additional solver may be called to produce nodal approximations of derived fields after the two potentials
have been obtained."
In EdgeElementInfo function (ElementDescription module) it seems to me that edge basis functions are calculated for 2D elements too.
The only concern is 5191 line in code (I am not entirely sure what traces of 2-D edge elements refers to exactly)
Is my attempt to develop this solver for any reason unfeasible in Elmer (or maybe redundant)?IF (cdim == 2 .AND. dim==1) THEN
CALL Warn('EdgeElementInfo', 'Traces of 2-D edge elements have not been implemented yet')
RETURN
END IF
Thanks!
Cheers,
Stje