[Getdp] Periodic Boundary Conditions

Lars Rindorf lhr at com.dtu.dk
Wed Feb 6 09:36:56 CET 2008

Dear Raymond
I don't know what you exactly are trying to model, but it seems to me that your Constraint{ ...}'s were wrong. I have changed your constraint and functionspace, and now the script seems to work according to what you have described.
I have attached the .pro file
Hope this helps you.
Lars Rindorf


From: getdp-bounces at geuz.org on behalf of Raymond Michaux
Sent: Tue 2/5/2008 9:39 PM
To: getdp at geuz.org
Subject: [Getdp] Periodic Boundary Conditions

Dear all,
I'am trying to use the "Periodic Boundary Conditions" in a simplified model of permanent magnet 
linear motor, according to an information given by Christophe Geuzaine on the same topic (see mailing 
list June 2005).
In order to check the method, I used a model with magnet centered on the teeth, I used the vector 
potential formulation and I ran the same problem with two sets of  boundary conditions :
    1. A Dirichlet a = 0 on three sides
    2. A Dirichlet a = 0 on upper border and periodic conditions between left and right side
For both cases, there was no error reported.
As the geometry is a symmetrical one, both methods should be equivalent but the calculation results are
obviously different.
Could you help me. The files are attached for tests. The file ".don" contains the main dimensions 
and is included in both ".geo" and ".pro". For the second case, please modify ".don" and ".pro" as 
indicated inside the files.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MoteurULg.pro
Type: application/octet-stream
Size: 3450 bytes
Desc: MoteurULg.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20080206/a52f8e10/attachment.pro>