[Getdp] Multiple natural boundary conditions

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


Hi,

Thanks a lot, it works too.

I've an other similar problem, but getdp create a NAN temperature field, I 
can't explain it. The geometry is different but still with 3 BCs, the main 
different is that now i define all the dimension of my geometry im meter, 
so the coordinates are small.


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/2ad19f9d/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: case.msh
Type: model/mesh
Size: 814459 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20141104/2ad19f9d/attachment-0001.msh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: case.pro
Type: application/octet-stream
Size: 2487 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20141104/2ad19f9d/attachment-0001.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: case.geo
Type: application/octet-stream
Size: 1327 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20141104/2ad19f9d/attachment-0001.geo>