[Getdp] Postprozessing:Curl of Curl
Kubicek Bernhard
Bernhard.Kubicek at arsenal.ac.at
Wed Aug 16 09:06:21 CEST 2006
Thank you very much for your reply.
My initial intention was this:
Having a current distribution, one calculates the magnetic vector potential.
Afterwards, in the post-processing, I would have liked to calculate the difference between the original current and the one obtained by rot rot A, to obtain a local error estimation, e.g. for mesh refinement.
Clearly, this was quite stupid, since a second derivative of a the linear finit element function space is mostly undefined.
Thank you very much again,
Bernhard Kubicek
-----Ursprüngliche Nachricht-----
Von: Andre Nicolet [mailto:andre.nicolet at fresnel.fr]
Gesendet: Freitag, 11. August 2006 12:15
An: Kubicek Bernhard
Betreff: Re: [Getdp] Postprozessing:Curl of Curl
It does not make sense to compute the Curl Curl of a quantity obtained
via the weak formulation of a second order operator. A is in the edge
element space, Curl A is in a subspace of the face element space but a
further derivation does not make sense in this framework. Your options are:
1) you do not really need Curl Curl A...
2) you need Curl Curl A "in a weak sense" i.e. in an integral the
integrand of which is something like w.Curl Curl A where w is some
(vector) weight factor so that you can integrate by part to obtain
something like Curl w.Curl A for the integrand...
3) completely change your formulation...
Think about this!
Andre
Kubicek Bernhard a écrit :
> Hi!
> Does anyone know how to postprocess a current that is defined by a
> formulation-calculated magnetic vector potential? I tried
> 1/mu0*{Curl {Curl A}}
> , with A being a BF_Edge, Form 1,
> but getDP complained the line in the preprocessing.
>
> Thank you, and a nice greetings,
> Bernhard
>
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp
>
--
---------------------------------------------------------------------------------
Prof. André Nicolet
Université Paul Cézanne
Institut Fresnel (UMR CNRS 6133)
Domaine Universitaire de Saint-Jérôme
F13397 Marseille cedex 20
France
tel: +33 4 91 28 87 73
fax: +33 4 91 67 44 28
E-mail: andre.nicolet at fresnel.fr
FOUNDATIONS OF PHOTONIC CRYSTAL FIBRES http://www.fresnel.fr/perso/renversez/foundations-photonic-crystal-fibres
http://www.icpress.co.uk/physics/p367.html