[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.
Thanks
-------------- 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>