If you want to couple more than two participants, there are two options: You can combine multiple normal coupling schemes (composition) or you can use a fully-implicit multi-coupling scheme. On this page, we explain both options.

Composition of bi-coupling schemes

preCICE allows us to combine multiple coupling schemes by only specifying them one after the other in the configuration:

<coupling-scheme:parallel-explicit>
  <participants first="A" second="B"/>
  ...
</coupling-scheme:parallel-explicit>
<coupling-scheme:parallel-explicit>
  <participants first="C" second="A"/>
  ...
</coupling-scheme:parallel-explicit>

For this example, all three participants are executed in parallel to one another, whereas A exchanges data with B and C, but B and C do not exchange data with each other.

Numerically, it only makes sense to either only have explicit schemes or to combine one implicit scheme with several explicit ones. Combining two implicit schemes is forbidden. To find out more, you can read this paper. If you want to resolve more than one strong interaction, you need a fully-implicit multi-coupling.

Still, even with only explicit schemes, you can do very strange combinations. Some of them will not work, but run into a deadlock. To better understand this, let us modify the example slightly: We add an interaction between B and C and switch to serial coupling schemes everywhere.

<coupling-scheme:serial-explicit>
    <participants first="A" second="B" />
    ...
    <exchange data="AB-Data" mesh="..." from="A" to="B" />
  </coupling-scheme:serial-explicit>

  <coupling-scheme:serial-explicit>
    <participants first="B" second="C" />
    ...
    <exchange data="BC-Data" mesh="..." from="B" to="C" />
  </coupling-scheme:serial-explicit>

  <coupling-scheme:serial-explicit>
    <participants first="C" second="A" />
    ...
    <exchange data="CA-Data" mesh="..." from="C" to="A" />
  </coupling-scheme:serial-explicit>

The mesh over which the data is communicated plays no role. We get a circular dependency:

A -- first to second --> B;
B -- first to second --> C;
C -- first to second --> A;

All three participants are a second participant in a serial coupling scheme, in which they receive data. Thus, they all wait for data in initialize, which is, however, only sent in the advance call of another participant, as explained on the coupling flow pages.

There are two easy fixes here: Either switch to parallel coupling scheme somewhere (better everywhere, otherwise you could still screw up the situation) or swap first and second in any of the coupling schemes. For example in the last coupling scheme:

- <participants first="C" second="A" />
+ <participants first="A" second="C" />

Now, A is the first participant in its two coupling schemes and does not receive any data in initialize. It can move to advance and send the data the other two participants are waiting for.

This exact setup is tested in the circular integration test.

Want to read more? Have a look at a blog post on compositional coupling.

Fully-implicit multi-coupling

In a fully-implicit multi-coupling, an arbitrary number of solvers are executed in parallel to each other in an implicit fashion.

<coupling-scheme:multi>
  <participant name="A" control="yes"/>
  <participant name="B" />
  <participant name="C" />
  ...
</coupling-scheme:multi>

Exactly one participant needs to take control of the coupling. preCICE computes the convergence measures and the acceleration on this participant.

Data being exchanged between multiple solvers must be configured indentically in terms of data initialization and substeps. All other tags are similar to a normal implicit coupling.

To find out more about multi coupling, you can also read Section 3.8 in Benjamin’s thesis or have a look at the tutorial of multiple perpendicular flaps in a channel.