[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] assembling equations containing curls
From: |
Yves Renard |
Subject: |
Re: [Getfem-users] assembling equations containing curls |
Date: |
Mon, 10 Sep 2007 09:49:15 +0200 |
User-agent: |
KMail/1.9.5 |
Hi alex,
No, we did not add a special term to take the curl of a finite element method
directly into account in the generic assembly description. It should be of
course possible to add it quite easily but for the moment, you can use a term
like
't=comp(vGrad(#1)); V(#1)+=t(:,3,2)-t(:,2,3); W(#1)+=t(:,1,3)-t(:,3,1);
Z(#1)+=t(:,2,1)-t(:,1,2,)';
Does this answer your question ?
Yves.
> hi everybody!
>
> thanks for the answer, but sorry for my bad questioning. i think the
> vector identity mentioned will work.
> unfortunatly i forgot to say that i also have terms like integral
> curl(A) to deal with. So my main question is, if there exists a function
> like vRot(#i), like the other comp commands, or way to deal directly
> with curls.
>
> regards and thanks in advance, alex
>
> Richard George wrote:
> > I guess you get an authorative answer from someone like Yves Renard
> > who wrote getfem, but
> >
> > could you express curl(curl(A)) in terms of available functions using
> > an identity such as
> >
> > curl(curl(A)) = grad(div(A)) - vector_laplacian(A)
> >
> > would that make sense / allow enough numerical accuracy?
> >
> > Richard George
> >
> > =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
> >
> > On Sun, 9 Sep 2007, Alex Magnes wrote:
> >> Hi everybody!
> >>
> >>
> >> Is there any possibility to assemble a term like integral curl(curl A)
> >> where A is a vector field? I just found examples dealing with gradients
> >> and divergences.
> >>
> >> regards and thanks in advance, alex
> >>
> >>
> >>
> >>
> >> _______________________________________________
> >> Getfem-users mailing list
> >> address@hidden
> >> https://mail.gna.org/listinfo/getfem-users
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users
--
Yves.
-------------------------------------------------------------------------
Yves Renard (address@hidden) tel : (33) 04.72.43.80.11
Pole de Mathematiques, fax : (33) 04.72.43.85.29
Institut Camille Jordan - CNRS UMR 5208
INSA de Lyon, Universite de Lyon
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard
-------------------------------------------------------------------------