dolfin team mailing list archive

[Evan Lezar <question156430@xxxxxxxxxxxxxxxxxxxxx>]

Question #156430 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/156430

    Status: Open => Answered

Evan Lezar proposed the following answer:
Pablo,

I have been using FEniCS for solving waveguide cutoff and dispersion
problems, and can send you a few code samples if you like. In my problems I
use homogeneous Dirichlet BCs, but this should not be too much of a problem
to extend to the non-homogeneous case.

As a start, I have attached a simple example for calculating the TE cutoff
modes of a hollow rectangular waveguide. To extend this to the
non-homogeneous case, you would have to use a non-zero expression when
instantiating the DirichletBC class.

Feel free to contact me if you have any questions.

Regards
--
Evan Lezar

Computational Electromagnetics Group
Department of Electrical and Electronic Engineering
University of Stellenbosch
Stellenbosch
South Africa

www.evanlezar.com

GoogleTalk: evanlezar
Skype: evanlezar


On Sun, May 8, 2011 at 12:11 AM, Pablo Suarez <
question156430@xxxxxxxxxxxxxxxxxxxxx> wrote:

> Question #156430 on DOLFIN changed:
> https://answers.launchpad.net/dolfin/+question/156430
>
> Description changed to:
> Hello,
>
>      I was wondering if anyone can answer this question.
> I need to solve for the follow PDE:
>
>        \nabla \times \nabla E - k^2 E = 0
>       (curl(curl(E)) - k^2 E = 0
>
> with BC   n \times E = F
> E is a vector field
> F is a vector field
> k^2 is a constant (k^2 is not an eigenvalue)
> this is in the interior of a 3D surface (presumably smooth)
> n is the outward normal  to the surface.
>
> This is the cavity problem for Maxwell Equation's  with non-homgenous
> Dirichlet BC.
> How can I incorporate this condition into my code?
> I have followed the example in undocumented demos curl curl. If anyone can
> help it will be greatly appreciated.
>
> Best wishes,
>
> Pablo
>
> --
> You received this question notification because you are a member of
> DOLFIN Team, which is an answer contact for DOLFIN.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin
> Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help   : https://help.launchpad.net/ListHelp
>

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.