[Gmsh] Boundary Conditions.

Christophe Geuzaine geuzaine at acm.caltech.edu
Tue Oct 19 16:39:31 CEST 2004


Stott CH Mr (PG/R - Maths & Statistics) wrote:
> Hi,
> 
> I've searched through the manuals and the mailing list archives, but 
> can't find an answer to this.
> 
> Say I define a square domain...
> 
>  lc=0.1;
> 
>  Point(1) = {0, 0, 0, 0.1*lc}; Point(2) = {1, 0, 0, lc}; Point(3) = {0, 
> 1, 0, lc}; Point(4) = {1, 1, 0, lc};
>  Line(1) = {1,2}; Line(2) = {2,4}; Line(3) = {4,3}; Line(4) = {3,1};
> 
>  Line Loop(5) = {1,2,3,4};
>  Plane Surface(6) = {5};
> 
> I want to use this with a diffusion solver, with two types of boundary 
> conditions (Dirichlet and Neumann). So I set up some Physical Entities,
> 
>  Physical Line(10) = {1,2,3};
>  Physical Line(11) = {4};
> 
>  Physical Surface(100) = {6};
> 
> Say I want to have a zero Neumann condition on Physical Line 10, and a 
> non-zero Dirichlet condition on Physical Line 11. Do I have to hard code 
> the numerical values of these into my solver? or is there someway of 
> attaching extra data to a Physical Line?

No: the only purpose of the physical entities is to tag some of the 
elements in your mesh with specific integers. It's up to you to give a 
meaning to these tags in your solver.

Christophe

-- 
Christophe Geuzaine
Applied and Computational Mathematics, Caltech
geuzaine at acm.caltech.edu - http://geuz.org