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 "i omega" has been fogotten in the VectorHelmholtz implementation?rot mu^-1 rot E - omega^2 eps E = i omega J
If you assume a spatially constant current density over an infinite large room then it induces an electric field
If you set, for example, 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.E = - i / ( omega eps ) J
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