Page 1 of 1

Periodic Boundaries problem

Posted: Wed Aug 17, 2016 11:00 am
by amir
Hello everyone,

I use a periodic boundary condition for my simulation with Code_Saturne 4.0.4. I have successfully used it for all boundaries.

I am now testing for a periodic solution in one direction and non-periodic in another.
For me the easiest solution would be if Saturne would allow for non-periodic in Y and periodic in X.

I tried to remove Y periodic faces without success (incomplete boundary condition error). I also put the translation vector equals to zero for Y, another error appeared.

So I am not sure if I can do it with Saturne and appreciate any help you can offer on this issue.

Many thanks
Kind regards
Amir

Re: Periodic Boundaries problem

Posted: Wed Aug 17, 2016 8:51 pm
by Yvan Fournier
Hello,

Yes, you can use periodicity in a single direction when needed. You must not define a periodicity with 0 distance.

When you have undefined boundary condition issues, did you check the recommended postprocessing (as recommended on the listing file) ?

First, did you define boundary conditions for faces which are not periodic ?

Regards,

Yvan

Re: Periodic Boundaries problem

Posted: Thu Aug 18, 2016 8:48 am
by amir
Hello Yvan,
Thanks for your reply,

Actually I tied both, periodicity distance equals zero for Y and using the boundary condition.

First, I put a zero distance for periodicity for Y without boundary condition; I got some error like that:

#############
Signal SIGSEGV (accès à une zone mémoire interdite) intercepté !

Pile d'appels :
1: 0x2b7ade0d8491 <cs_join_perio_split_update+0x411> (libsaturne.so.0)
2: 0x2b7ade0fe6cf <cs_join_update_mesh_after_split+0x5c8f> (libsaturne.so.0)
3: 0x2b7ade0a6852 <cs_join_all+0x20a2> (libsaturne.so.0)
4: 0x2b7added478f <cs_preprocess_mesh+0x12f> (libsaturne.so.0)
5: 0x2b7adde4373d <cs_run+0x11d> (libsaturne.so.0)
6: 0x2b7adde43512 <main+0x112> (libsaturne.so.0)
7: 0x2b7ae0ea0b15 <__libc_start_main+0xf5> (libc.so.6)
8: 0x403619 <> (cs_solver)
Fin de la pile
#############

Then I removed periodicity for Y and defined a boundary condition on it, the error was:

#############
Signal SIGSEGV (accès à une zone mémoire interdite) intercepté !

Pile d'appels :
1: 0x406308 <cs_user_boundary_conditions_+0x3f8> (cs_solver)
2: 0x2abd83cc49a0 <tridim_+0x3640> (libsaturne.so.0)
3: 0x2abd83ba3196 <caltri_+0x3006> (libsaturne.so.0)
4: 0x2abd83b8ba1e <cs_run+0x3fe> (libsaturne.so.0)
5: 0x2abd83b8b512 <main+0x112> (libsaturne.so.0)
6: 0x2abd86be8b15 <__libc_start_main+0xf5> (libc.so.6)
7: 0x403669 <> (cs_solver)
Fin de la pile
#############

At the end, I put the periodicity equals to zero for Y and set a boundary condition on it, the error was the same as the first case.

Enclosed please find 3 listing files.
Kind regards,
Amir

Re: Periodic Boundaries problem

Posted: Thu Aug 18, 2016 10:35 am
by Yvan Fournier
Hello,

In any case, you must not define a periodicity in y (a peruodicity with 0 distance).

For the othet case, I won't be able to check before a week.

Regards,

Yvan