2D waveguide using Nedelec elements

[Rolf Wester]

Hi,

I would like to compute eigen-solutions of optical waveguides with z-dependence ~ exp(-i k_z z) using the deal.ii Nedelec elements. I found some resources on this topic but no one that worked for me. I would be very appreciative for a link to a deal.ii example that solves the optical waveguide eigen-problem or any other help.

Thank you in anticipation
Regards
Rolf

[Wolfgang Bangerth]

On 08/17/2017 08:59 AM, Rolf Wester wrote:
>
> I would like to compute eigen-solutions of optical waveguides with
> z-dependence ~ exp(-i k_z z) using the deal.ii Nedelec elements. I found
> some resources on this topic but no one that worked for me. I would be
> very appreciative for a link to a deal.ii example that solves the
> optical waveguide eigen-problem or any other help.

Rolf -- I don't know whether anyone has solved this problem already (I
suspect that there are people), but in any case, there is at least no
code that's available to me. Absent such a code, you'll have to write
one yourself. There are of course example programs that solve
eigenproblems that you can use as a basis, and you only have to change
the underlying equation you want to solve. Have you tried? What concrete
problems did you find?

Best
W.

--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

[Rolf Wester]

Hi Wolfgang,

I still try to implement it. There is one point that I could not figure out until now. I need PETScWrappers::BlockSparseMatrix for the eigensolver and I don't know how to initialize these with the block information. I have two function spaces, FE_Nedelec<dim>(order-1) and FE_DGQ<dim>(order).

Kind regards
Rolf

[Jean-Paul Pelteret]

Hi Rolf,

Step-55 has PETSc compatibility and constructs a block system. You can look there for information as to how to initialise the block system.

I hope that this helps.

Regards,

Jean-Paul