Thanks for your fast reply!
About your first point: Yes, I currently use a FeSystem composed of two 3D
Nedelec fes.
The second point I agree with you. Setting the additional values should not
be a problem since I know the analytical solution of the problem so this
should not be a problem concerning sulvability of the system. However you
are right, it should also not be necessary to set these values in the first
place so I can just skip this part.
As a kind of an extension to my first question: I have a function which
computes the analytical solution for a given position. Is there a function
to compute the best approximation for a given element type (like nedelec)?
Again, thank you for your time,
Kind regards,
Pascal
Am Dienstag, 16. Januar 2018 18:45:15 UTC+1 schrieb Wolfgang Bangerth:
>
> On 01/16/2018 01:41 AM, Pascal Kraft wrote:
> > I am currently using a FeSystem composed of two 3D Fields (real and
> > imaginary E-field) and I want to impose Dirichlet conditions of the sort
> > E(x,y,z) = E_{in}(x,y,z) on the input interface (in an xy-plane).
> > Earlier the z-component had been 0 so I did not run into real problems
> > and could use the project_boundary_values_curl_conforming_l2 function
> > (for the x and y components) and my own code for the z-components
> > (setting them 0 in the first cell layer). Now however I would like to
> > impose different values on the z-component. Since my cells have variable
> > thickness in z-direction (i.e. hanging node constraints appear) I would
> > like to use some library function for the interpolation.
> > Are there any feasible solutions to this?
> > Ideally I would like to keep using
> > project_boundary_values_curl_conforming_l2 for the tangential components
> > and just add some code for the z-components.
>
> The question doesn't quite have enough information:
> * Are you using Nedelec elements for the two fields?
> * If the input face is in the xy-plane, and you are using a finite
> element that is only curl-conforming (such as the Nedelec element),
> isn't it correct that you can only prescribe *tangential* components for
> boundary values? In other words, you cannot prescribe z-components anyway.
>
> Best
> W.
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: [email protected]
> <javascript:>
> www: http://www.math.colostate.edu/~bangerth/
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
