Bidirectional exchange

Author: davidscn

channel-transport-particles/fluid-openfoam/constant/transportProperties

@@ -8,4 +8,7 @@ FoamFile
 
 transportModel  Newtonian;
 
-nu              nu [ 0 2 -1 0 0 0 0 ] 1;
+nu              [0 2 -1 0 0 0 0] 1.5e-05; // = air at 20°C m^2/s
+rho             1.2; // kg/m^3
+g               (0 -9.81 0);
+

channel-transport-particles/fluid-openfoam/system/controlDict

@@ -6,7 +6,7 @@ FoamFile
     object      controlDict;
 }
 
-application     pimpleFoam;
+application     AndersonJacksonFoam;
 
 startFrom       startTime;
 

channel-transport-particles/fluid-openfoam/system/fvSchemes

@@ -18,7 +18,7 @@ gradSchemes
 
 divSchemes
 {
-    default         none;
+    default         Gauss linear;
     div(phi,U)      bounded Gauss upwind;
     div((nuEff*dev2(T(grad(U))))) Gauss linear;
 }

channel-transport-particles/fluid-openfoam/system/preciceDict

@@ -12,6 +12,14 @@ participant Fluid;
 
 modules (FF);
 
+FF
+{
+  nameAlpha alphaSolid;
+  nameExplicitMomentum ExplicitMomentum;
+  nameImplicitMomentum ImplicitMomentum;
+};
+
+
 interfaces
 {
   Interface1
@@ -22,11 +30,15 @@ interfaces
 
     readData
     (
+        ImplicitMomentum
+        ExplicitMomentum
+        Alpha
     );
 
     writeData
     (
         Velocity
+        PressureGradientFull
     );
   };
-};
+};

channel-transport-particles/precice-config.xml

@@ -1,43 +1,84 @@
 <?xml version="1.0" encoding="UTF-8" ?>
-<precice-configuration experimental="true">
+<precice-configuration experimental="True">
   <log>
     <sink
-      filter="%Severity% > debug and %Rank% = 0"
+      type="stream"
+      output="stdout"
       format="---[precice] %ColorizedSeverity% %Message%"
-      enabled="true" />
+      filter="(%Severity% > debug) and not ((%Rank% != 0))" />
   </log>
 
-  <profiling mode="fundamental" synchronize="false" />
+  <profiling mode="off" />
 
-  <data:vector name="Velocity" />
+  <data:scalar name="Alpha" waveform-degree="2" />
+  <data:vector name="ExplicitMomentum" waveform-degree="2" />
+  <data:scalar name="ImplicitMomentum" waveform-degree="2" />
+  <data:vector name="Velocity" waveform-degree="0" />
+  <data:vector name="PressureGradientFull" waveform-degree="0" />
 
-  <mesh name="Fluid-Mesh" dimensions="2">
+  <mesh name="Fluid-Mesh" dimensions="3">
     <use-data name="Velocity" />
+    <use-data name="PressureGradientFull" />
+    <use-data name="Alpha" />
+    <use-data name="ExplicitMomentum" />
+    <use-data name="ImplicitMomentum" />
   </mesh>
 
   <participant name="Fluid">
     <provide-mesh name="Fluid-Mesh" />
     <write-data name="Velocity" mesh="Fluid-Mesh" />
+    <write-data name="PressureGradientFull" mesh="Fluid-Mesh" />
+    <read-data name="Alpha" mesh="Fluid-Mesh" />
+    <read-data name="ExplicitMomentum" mesh="Fluid-Mesh" />
+    <read-data name="ImplicitMomentum" mesh="Fluid-Mesh" />
   </participant>
 
   <participant name="Particles">
-    <receive-mesh name="Fluid-Mesh" from="Fluid" api-access="true" />
+    <receive-mesh name="Fluid-Mesh" from="Fluid" safety-factor="0.5" api-access="true" />
     <read-data name="Velocity" mesh="Fluid-Mesh" />
+    <read-data name="PressureGradientFull" mesh="Fluid-Mesh" />
+    <write-data name="Alpha" mesh="Fluid-Mesh" />
+    <write-data name="ExplicitMomentum" mesh="Fluid-Mesh" />
+    <write-data name="ImplicitMomentum" mesh="Fluid-Mesh" />
     <mapping:rbf-pum-direct
       direction="read"
       from="Fluid-Mesh"
       constraint="consistent"
       vertices-per-cluster="30">
-      <basis-function:compact-polynomial-c4 support-radius="1.5" />
+      <basis-function:compact-polynomial-c4 support-radius="0.5" />
     </mapping:rbf-pum-direct>
+    <mapping:coarse-graining
+      direction="write"
+      to="Fluid-Mesh"
+      constraint="conservative"
+      radius="0.05" />
   </participant>
 
   <m2n:sockets acceptor="Fluid" connector="Particles" exchange-directory=".." />
 
-  <coupling-scheme:parallel-explicit>
+  <coupling-scheme:serial-explicit>
+    <participants second="Fluid" first="Particles" />
     <time-window-size value="0.005" />
     <max-time value="0.25" />
-    <participants first="Fluid" second="Particles" />
-    <exchange data="Velocity" mesh="Fluid-Mesh" from="Fluid" to="Particles" />
-  </coupling-scheme:parallel-explicit>
+    <exchange data="Velocity" mesh="Fluid-Mesh" from="Fluid" to="Particles" substeps="false" />
+    <exchange
+      data="PressureGradientFull"
+      mesh="Fluid-Mesh"
+      from="Fluid"
+      to="Particles"
+      substeps="false" />
+    <exchange data="Alpha" mesh="Fluid-Mesh" from="Particles" to="Fluid" substeps="true" />
+    <exchange
+      data="ExplicitMomentum"
+      mesh="Fluid-Mesh"
+      from="Particles"
+      to="Fluid"
+      substeps="true" />
+    <exchange
+      data="ImplicitMomentum"
+      mesh="Fluid-Mesh"
+      from="Particles"
+      to="Fluid"
+      substeps="true" />
+  </coupling-scheme:serial-explicit>
 </precice-configuration>