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. In addition to performing the computations in the time

the applicability of this solver to 3D problems

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