Solving Maxwell's Equations in time-harmonic form. More...
|
Modules | |
| Two Dimensions | |
| Three Dimensions | |
Detailed Description
Solving Maxwell's Equations in time-harmonic form.
The equations solved are the time-harmonic Maxwell's equations:
![\[ - \varepsilon i \omega \vec E + \text{curl} \vec H = \sigma \vec E + \vec J \]](form_167.png)
![\[ \mu i \omega \vec H + \text{curl} \vec E = 0 \]](form_168.png)
with the divergence conditions (assuming no charge density)
![\[ \text{div} \mu \vec H = 0 \qquad \text{and} \qquad \text{div} \varepsilon \vec E = 0. \]](form_169.png)
Extracting 
from the second equation and inserting in into the first formally yields
![\[ \text{curl}( \mu^{-1} \text{curl} \vec E) - \omega^2 (\varepsilon + \frac{\sigma}{i\omega}) \vec E = -i \omega \vec J, \]](form_171.png)
the so-called electric source problem. The appropriate space is
![\[ H(\text{curl}; D) := \{\vec u \in L^2(D)^3 : \text{curl} \vec u \in L^2(D)^3 \}. \]](form_172.png)
However, the curl-curl form above is not defined for electric fields in 
. This is resolved in the variational formulation of the electric source problem:
Find 
such that
![\[ \int_D (\mu^{-1} \text{curl} \vec E \cdot \text{curl} \vec v - \omega^2 \tilde\varepsilon \vec E \cdot \vec v) \, d\vec x = \int_D \vec f \cdot \vec v \, d\vec x \quad \forall \vec v \in H_0(\text{curl}; D), \]](form_175.png)
where 
is with the prefect conductor boundary conditions 
, 
and 
. The associated bilinear operator of this variational formulation is not elliptic and the divergence condition is an independent constraint. Normally, Maxwell's equations are discretised using Nedelec's elements. However, there are good reasons why one would like to use standard H1-conforming FEM. This is possible using weighted regularization.