support WR bindings and fenics adapter

Author: BenjaminRodenberg

FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py

@@ -90,10 +90,11 @@ def remaining_boundary(x, on_boundary):
 # get the adapter ready
 
 # read fenics-adapter json-config-file)
-adapter_config_filename = "../tools/precice-adapter-config-fsi-s.json"
+adapter_config_filename = "../tools/precice-adapter-config-fsi-s_wr.json"
+other_adapter_config_filename = "../tools/precice-adapter-config-excite.json"
 
 # create Adapter
-precice = Adapter(adapter_config_filename)
+precice = Adapter(adapter_config_filename, other_adapter_config_filename)
 
 # create subdomains used by the adapter
 clamped_boundary_domain = AutoSubDomain(left_boundary)
@@ -107,6 +108,7 @@ def remaining_boundary(x, on_boundary):
                                 dimension=dim,
                                 dirichlet_boundary=clamped_boundary_domain)
 
+precice_dt = float(precice_dt)
 fenics_dt = precice_dt  # if fenics_dt == precice_dt, no subcycling is applied
 #fenics_dt = 0.02  # if fenics_dt < precice_dt, subcycling is applied
 dt = Constant(np.min([precice_dt, fenics_dt]))
@@ -238,6 +240,7 @@ def avg(x_old, x_new, alpha):
             displacement_out << (u_n, t)
 
         u_tip.append(u_n(0.6, 0.2)[1])
+        print(u_n(0.6, 0.2))
         time.append(t)
 
 
@@ -251,3 +254,4 @@ def avg(x_old, x_new, alpha):
 
 output_file = open("subiteration_out.txt", "a")
 output_file.write("{};{}\n".format(u_tip[-1], dt(0)))
+precice.finalize()

FSI/cylinderFlap/OpenFOAM-FEniCS/tools/excitement.py

@@ -4,7 +4,6 @@
 import numpy as np
 import precice_future as precice
 
-force_excitement = np.array([0, 10, 0])
 omega_excitement = 1
 
 # beam geometry
@@ -18,7 +17,7 @@
 y_attack = .5 * (y_bottom + y_top)
 
 
-def compute_force(time):
+def compute_force(time, force_excitement):
     omega = omega_excitement
     force_0 = force_excitement
     return force_0 * np.sin(omega * time)
@@ -40,7 +39,7 @@ def compute_force(time):
 write_data_name = "Forces0"
 read_data_name = "Displacements0"
 
-n = 1
+n_vertices = 1
 
 solver_process_index = 0
 solver_process_size = 1
@@ -51,7 +50,13 @@ def compute_force(time):
 mesh_id = interface.get_mesh_id(mesh_name)
 
 dimensions = interface.get_dimensions()
-vertices = np.array([x_attack, y_attack, 0])
+force_excitement = np.zeros((n_vertices, dimensions))
+force_excitement[:, 0] = 0
+force_excitement[:, 1] = 10
+
+vertices = np.zeros((n_vertices, dimensions))
+vertices[:, 0] = x_attack
+vertices[:, 1] = y_attack
 
 vertex_ids = interface.set_mesh_vertices(mesh_id, vertices)
 write_data_id = interface.get_data_id(write_data_name, mesh_id)
@@ -61,7 +66,7 @@ def compute_force(time):
 t = 0
 
 if interface.is_action_required(precice.action_write_initial_data()):
-    forces = compute_force(t + dt)
+    forces = compute_force(t + dt, force_excitement)
     interface.write_block_vector_data(write_data_id, vertex_ids, forces)
     interface.fulfilled_action(precice.action_write_initial_data())
 
@@ -73,7 +78,7 @@ def compute_force(time):
         print("DUMMY: Writing iteration checkpoint")
         interface.fulfilled_action(precice.action_write_iteration_checkpoint())
 
-    forces = compute_force(t+dt)
+    forces = compute_force(t+dt, force_excitement)
     print("sending forces: {}".format(forces))
     interface.write_block_vector_data(write_data_id, vertex_ids, forces)
     dt = interface.advance(dt)

FSI/cylinderFlap/OpenFOAM-FEniCS/tools/precice-adapter-config-excite.json

@@ -0,0 +1,12 @@
+{
+  "solver_name": "fenics",
+  "config_file_name": "../tools/precice-config.xml",  
+  "interface": {
+      "coupling_mesh_name": "Excitement",
+      "write_data_name": "Forces0",
+      "read_data_name": "Displacements0"
+    },
+  "waveform": {
+      "n_substeps": 1
+    }
+}

FSI/cylinderFlap/OpenFOAM-FEniCS/tools/precice-adapter-config-fsi-s_wr.json

@@ -0,0 +1,12 @@
+{
+  "solver_name": "fenics",
+  "config_file_name": "../tools/precice-config_wr.xml",  
+  "interface": {
+      "coupling_mesh_name": "Solid",
+      "write_data_name": "Displacements",
+      "read_data_name": "Forces"
+    },
+  "waveform": {
+      "n_substeps": 1
+    }
+}

FSI/cylinderFlap/OpenFOAM-FEniCS/tools/precice-config.xml

@@ -6,7 +6,7 @@
       <sink type="stream" output="stdout"  filter= "%Severity% > debug" format="preCICE:%ColorizedSeverity% %Message%" enabled="true" />
     </log>
 
-    <solver-interface dimensions="3" >
+    <solver-interface dimensions="2" >
 
         <data:vector name="Forces0"/>
         <data:vector name="Displacements0"/>
@@ -26,8 +26,8 @@
             <use-mesh name="Solid" from="fenics"/>
             <write-data name="Forces0" mesh="Excitement-Mesh"/>
             <read-data name="Displacements0" mesh="Excitement-Mesh"/>
-            <mapping:petrbf-thin-plate-splines direction="write" from="Excitement-Mesh" to="Solid" constraint="conservative" z-dead="true"/>
-            <mapping:petrbf-thin-plate-splines direction="read" from="Solid" to="Excitement-Mesh" constraint="consistent" z-dead="true"/>
+            <mapping:petrbf-thin-plate-splines direction="write" from="Excitement-Mesh" to="Solid" constraint="conservative" />
+            <mapping:petrbf-thin-plate-splines direction="read" from="Solid" to="Excitement-Mesh" constraint="consistent" />
         </participant>
 
         <participant name="fenics">

FSI/cylinderFlap/OpenFOAM-FEniCS/tools/precice-config_wr.xml

@@ -0,0 +1,51 @@
+<?xml version="1.0"?>
+
+<precice-configuration>
+
+    <log>
+      <sink type="stream" output="stdout"  filter= "%Severity% > debug" format="preCICE:%ColorizedSeverity% %Message%" enabled="true" />
+    </log>
+
+    <solver-interface dimensions="2" >
+
+        <data:vector name="Forces1"/>
+        <data:vector name="Displacements1"/>
+
+        <mesh name="Excitement-Mesh">
+            <use-data name="Forces1"/>
+            <use-data name="Displacements1"/>
+        </mesh>
+
+        <mesh name="Solid">
+            <use-data name="Forces1"/>
+            <use-data name="Displacements1"/>
+        </mesh>
+
+        <participant name="Excitement">
+            <use-mesh name="Excitement-Mesh" provide="yes"/>
+            <use-mesh name="Solid" from="fenics"/>
+            <write-data name="Forces1" mesh="Excitement-Mesh"/>
+            <read-data name="Displacements1" mesh="Excitement-Mesh"/>
+            <mapping:petrbf-thin-plate-splines direction="write" from="Excitement-Mesh" to="Solid" constraint="conservative" />
+            <mapping:petrbf-thin-plate-splines direction="read" from="Solid" to="Excitement-Mesh" constraint="consistent" />
+        </participant>
+
+        <participant name="fenics">
+            <use-mesh name="Solid" provide="yes"/>
+            <read-data name="Forces1" mesh="Solid"/>
+            <write-data name="Displacements1" mesh="Solid"/>
+            <watch-point mesh="Solid" name="point1" coordinate="0.6;0.2;0." />
+        </participant>
+
+        <m2n:sockets from="Excitement" to="fenics" distribution-type="gather-scatter"/>
+
+        <coupling-scheme:serial-explicit>
+            <participants first="Excitement" second="fenics" />
+            <max-time value="1"/>
+            <timestep-length value="0.1" />
+            <exchange data="Forces1" mesh="Solid" from="Excitement" to="fenics"/>
+            <exchange data="Displacements1" mesh="Solid" from="fenics" to="Excitement" initialize="0"/>
+        </coupling-scheme:serial-explicit>
+    </solver-interface>
+
+</precice-configuration>