[Getdp] Multiple natural boundary conditions

Hamid Badi h.badi at bbs-slama.com
Wed Nov 5 01:02:34 CET 2014


Thanks a lot

Le mardi 4 novembre 2014 17:12:43 UTC+1, michae... at infineon.com a écrit :
>
>  Hi,
>
>  
>
> the point is that in GetDP it’s not allowed to subtract or add anything 
> from/to
>
> a Dof in a single Galerkin statement. Therefore you have to split your 2nd 
>
> Galerkin into two:
>
>  
>
> Then it works J
>
>  
>
> Cheers,
>
> Michael
>
>  
>
>  
>  
> *From:* getdp [mailto:getdp-... at ace20.montefiore.ulg.ac.be <javascript:>] *On 
> Behalf Of *Hamid Badi
> *Sent:* Tuesday, November 04, 2014 1:14 PM
> *To:* getdp-... at googlegroups.com <javascript:>
> *Subject:* [Getdp] Multiple natural boundary conditions
>  
>  
>  
> Hello,
>  
>  
>  
> My problem is common, i'd like to solve the heat equation with, for 
> example, 2 natural Robin BC and one adibatique BC : -lambda.dnT = h.(T-Ti) 
> on Gamma_i with i={1,2} and lambda.dnT = 0 on Gamma_3. 
>  
> I write the weak form with u in Form0 of the heat equation leads to the 
> classical form :
>   
>  
>   
> Equation{
>       Galerkin{ [lambda[] * Dof{Grad u}, {Grad u}];
>         In Omega; Jacobian JVol; Integration I1;}
>    
>   Galerkin{ [ h[]* (Dof{u} - T[]), {u}];
>        In Ambiance; Jacobian JSur; Integration I1;}
> }
>   
>  
>  
>  
>  
> In my problem there is no essential BC so i do not have Constraint block.
>  
> The computation give a 0 temperature field, I do not understand what's 
> wrong.
>  
>  
>  
> Moreover : 
>  
> 1 - I wan to store in a file the temperature field respects to the nodes 
> but PostOperation block computes over elements.
>  
> 2 - How can I compute the integral of the flux over Gamma_i, to do so, I 
> need for each boundary edge : the flux and the normal vector
>  
>  
>  
> Best regards.
>   
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20141104/6decdbc4/attachment-0001.html>