Syrthes coupled simulation crashes

[AdamLarat]

Hi,
Now that I have been successfully able to run the Syrthes' coupling tutorial 3_DISKS, I have tried to report the procedure to my own test case. Both simulation (Syrthes and Saturne) work independently.
However, once I try to run the coupled simulation, the Syrthes process crashes with the following error:

Code: Select all

Error during SYRTHES execution
==============================

ple_locator.c:3439: Fatal error.

Locator trying to use distant space dimension 1684426904
with local space dimension 0

syrthes: syr_cfd_coupling.c:219: _syr_error_handler: Assertion `0' failed.
[po21210:50962] *** Process received signal ***
[po21210:50962] Signal: Aborted (6)
[po21210:50962] Signal code:  (-6)
[po21210:50962] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x13140)[0x7fb2fcb01140]
[po21210:50962] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x141)[0x7fb2fc952d51]
[po21210:50962] [ 2] /lib/x86_64-linux-gnu/libc.so.6(abort+0x123)[0x7fb2fc93c537]
[po21210:50962] [ 3] /lib/x86_64-linux-gnu/libc.so.6(+0x2240f)[0x7fb2fc93c40f]
[po21210:50962] [ 4] /lib/x86_64-linux-gnu/libc.so.6(+0x316d2)[0x7fb2fc94b6d2]
[po21210:50962] [ 5] ./syrthes(+0x3c965)[0x55cff6239965]
[po21210:50962] [ 6] /home/alarat/Codes/Carambarre/libs/code_saturne-8.0.4//lib/libple.so.2(ple_error+0x7e)[0x7fb2fcc3d62e]
[po21210:50962] [ 7] /home/alarat/Codes/Carambarre/libs/code_saturne-8.0.4//lib/libple.so.2(ple_locator_extend_search+0x222)[0x7fb2fcc42062]
[po21210:50962] [ 8] /home/alarat/Codes/Carambarre/libs/code_saturne-8.0.4//lib/libple.so.2(ple_locator_set_mesh+0x278)[0x7fb2fcc48b28]
[po21210:50962] [ 9] ./syrthes(+0x3e080)[0x55cff623b080]
[po21210:50962] [10] ./syrthes(+0x3f023)[0x55cff623c023]
[po21210:50962] [11] ./syrthes(+0x3b1ad)[0x55cff62381ad]
[po21210:50962] [12] ./syrthes(+0xddd6)[0x55cff620add6]
[po21210:50962] [13] ./syrthes(+0x5fba)[0x55cff6202fba]
[po21210:50962] [14] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xea)[0x7fb2fc93dd7a]
[po21210:50962] [15] ./syrthes(+0x628a)[0x55cff620328a]
[po21210:50962] *** End of error message ***

Do you have any idea of what is going wrong? Especially concerning these "space dimensions" !
You can find the whole test case attached.
Thanks a lot!

[Yvan Fournier]

Hello,

The logging of the space dimension seems to be a bug, but I suspect your issue is that you have a 2D domain in Syrthes, and a 3D one in code_saturne. In the 3D disks example (whose reference setup should be updated), the "2D projection" option in the conjugate heat transfer section is set, which allows aligning the 3D computation with the 2D one.

Best regards,

Yvan

[AdamLarat]

Hi Yvan,
Thank you for your answer.
In my setup.xml Saturne's config file, I have the following section : 

Code: Select all

 <conjugate_heat_transfer>
      <external_coupling>
        <syrthes>
          <projection_axis>Z</projection_axis>
          <selection_criteria>1</selection_criteria>
          <syrthes_name>SYRTHES</syrthes_name>
          <verbosity>0</verbosity>
          <visualization>1</visualization>
        </syrthes>
        <syrthes_instances>
          <instance name="SOLID">
            <coupled_boundary label="Barre1_P1"/>
          </instance>
        </syrthes_instances>
      </external_coupling>
    </conjugate_heat_transfer>

With the projection line "<projection_axis>Z</projection_axis>".
Is this what you are talking about or is it something else ?

[AdamLarat]

All right! I have been able to run the simulation by rearranging the <conjugate_heat_transfer> as:

Code: Select all

    <conjugate_heat_transfer>
      <external_coupling>
        <projection_axis>Z</projection_axis>
        <syrthes_instances>
          <instance name="SOLID">
            <coupled_boundary label="Barre1_P1"/>
          </instance>
        </syrthes_instances>
        <selection_criteria>1</selection_criteria>
        <syrthes_name>SYRTHES</syrthes_name>
        <tolerance>0.1</tolerance>
        <verbosity>0</verbosity>
        <visualization>1</visualization>
      </external_coupling>
    </conjugate_heat_transfer>

It is still not clear for me why, but at least it is working!
Thanks!