I am still quite new to Elmer and wondering if this is actually a bug or just because I have overseen the factor by myself. If it is not a bug, I am of course happy to learn where my mistake is.

According to Elmer Modules manual, the VectorHelmholtz module should solve the time harmonic Maxwell equations, which are basically (17.3):

Is it possible that the prefactor infront of the current density "rot mu^-1 rotE- omega^2 epsE= i omegaJ

*i omega*" has been fogotten in the VectorHelmholtz implementation?

If you assume a spatially constant current density over an infinite large room then it induces an electric field

If you set, for example,E= - i / ( omega eps )J

*J = 1 A/m^2*and insert

*eps*, you would expect for a frequency of

*10 GHz*an induced field of

*~ -i 1.8 V/m*. In all my simulations I however get an induced field of

*~ -2.9e-11 V/m*, which is off by 11 orders of magnitude, to be precise by a factor of

*i * 2 * pi * 10GHz*.

Furthermore the induced field seems to have the wrong phase. Despite it beeing phase shifted by -Pi/2 due to the "-i" factor, it is phase shifted by Pi in my simulation output.

Assuming that the prefactor "i omega" has been just forgotten in the VectorHelmholtz implementation, all phases and field magnitudes come out correctly.

I have actually also checked

*VectorHelmholtz.F90*and not found a factor of i*omega infront of the current density, though I might have overseen it. I am actually a bit confused by this fluid dynamics formulation...

Elmer Version: 8.3

OS: Ubuntu 16.04 LTS

How to reproduce: The simplest way to reproduce this, is by running the attached .sif file, on a domain of size 0.09x0.03x0.4m³ meshed with maxh=0.003m (~ 0.1*lambda)

Of course this is not an infinite large room and J is not constant. Though, for the rot terms to become (such) important, we would expect to see wave propagation or a standing wave in the simulation output. You can see this by for example changing the sine in z-direction to a cosine. Then J has a sharp discontinuity at the simulation domain boundaries and you should see a standing wave.

Thank you,

Regards,

Intex