WIP: aste Replay Mode tutorial on elastic flap case

CHT/heat_exchanger/buoyantSimpleFoam-CalculiX/precice-config.xml

@@ -1,4 +1,4 @@
-<?xml version='1.0' encoding='UTF-8'?>
+<?xml version="1.0"?>
 <precice-configuration>
   <log>
       <sink filter="%Severity% > debug and %Rank% = 0" format="---[precice] %ColorizedSeverity% %Message%" enabled="true"/>

FSI/3D_Tube/OpenFOAM-CalculiX/precice-config.xml

@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding='UTF-8'?>
+<?xml version="1.0"?>
 
 <precice-configuration>
 

FSI/cylinderFlap/OpenFOAM-CalculiX/Fluid/system/controlDict

@@ -10,8 +10,8 @@ FoamFile
 
 // Note: With OpenFOAM v1806 and OpenFOAM 6, the DyM solvers
 // were marked deprecated and merged into their respective standard solvers.
-application         pimpleFoam;    // OpenFOAM v1806, OpenFOAM 6, or newer
-// application         pimpleDyMFoam; // OpenFOAM v1712, OpenFOAM 5.x, or older
+// application         pimpleFoam;    // OpenFOAM v1806, OpenFOAM 6, or newer
+application         pimpleDyMFoam; // OpenFOAM v1712, OpenFOAM 5.x, or older
 
 startFrom           startTime;
 

FSI/cylinderFlap/OpenFOAM-FEniCS-bigdt/OpenFOAM-FEniCS/.gitignore

@@ -0,0 +1,13 @@
+Fluid
+!Fluid/0
+!Fluid/0.orig
+!Fluid/constant
+!Fluid/system
+!Fluid/Fluid.foam
+!Fluid/precice-adapter-config.yml
+Solid
+!Solid/cyl-flap.py
+!Solid/precice-adapter-config-fsi-s.json
+*.log
+precice-fenics-events.json
+precice-Fluid-events.json

FSI/cylinderFlap/OpenFOAM-FEniCS/Postprocessing/plot_tip.py

@@ -10,9 +10,10 @@
 import matplotlib.pyplot as plt
 
 fenics_watchpoint = np.loadtxt('../precice-fenics-watchpoint-point1.log', skiprows = 1)
-ccx_watchpoint = np.loadtxt('/home/richyrich/Bachelorarbeit/OpenFOAM-CalculiX/precice-Fluid-watchpoint-point1.log', skiprows = 1)
+#ccx_watchpoint = np.loadtxt('/home/richyrich/Bachelorarbeit/OpenFOAM-CalculiX/precice-Fluid-watchpoint-point1.log', skiprows = 1)
 
 
-plt.plot(fenics_watchpoint[:,0],fenics_watchpoint[:,8])
+plt.plot(fenics_watchpoint[:,0],fenics_watchpoint[:,7])
 
+plt.show()
 #plt.plot(ccx_watchpoint[:,0],ccx_watchpoint[:,8])

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

@@ -249,4 +249,3 @@ def avg(x_old, x_new, alpha):
 plt.ylabel("Tip displacement")
 plt.show()
 
-

FSI/flap_perp/ASTE/.gitignore

@@ -0,0 +1,5 @@
+preCICE-output/
+Solid-fenics*
+*.log
+vtk*
+Solid/FSI-S

FSI/flap_perp/ASTE/Allclean

@@ -0,0 +1,64 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+echo "Cleaning..."
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+# Participant 1: Fluid
+Participant1="Fluid"
+cd ${Participant1}
+    # Clean the case
+    cleanCase
+    rm -rfv 0
+    # Create an empty .foam file for ParaView
+    # Note: ".foam" triggers the native OpenFOAM reader of ParaView.
+    # Change to ".OpenFOAM" to use the OpenFOAM reader provided with OpenFOAM.
+    touch ${Participant1}.foam
+cd ..
+# Remove the log files
+rm -fv ${Participant1}_blockMesh.log
+rm -fv ${Participant1}_checkMesh.log
+rm -fv ${Participant1}_potentialFoam.log
+rm -fv ${Participant1}_decomposePar.log
+rm -fv ${Participant1}.log
+rm -fv ${Participant1}_reconstructPar.log
+
+# Participant 2: Solid
+Participant2="Solid"
+cd ${Participant2}
+cd FSI-S
+    # Clean the case
+	rm -fv *.vtu
+	rm -fv *.pvd
+    rm -fv *.log
+cd ..
+cd ..
+# Remove the log files
+rm -fv spooles.out
+rm -fv ${Participant2}.log
+
+# Remove the precice-events-files
+rm -fv precice-fenics-events.json
+rm -fv precice-Fluid-events.json
+
+# Remove the preCICE-related log files
+echo "Deleting the preCICE log files..."
+rm -fv \
+    precice-*.log \
+    precice-postProcessingInfo.log
+
+# Output files for preCICE versions before 1.2:
+rm -fv \
+    iterations-${Participant1}.txt iterations-${Participant2}.txt \
+    convergence-${Participant1}.txt convergence-${Participant2}.txt \
+    Events-${Participant1}.log Events-${Participant2}.log \
+    EventTimings-${Participant1}.log EventTimings-${Participant2}.log
+
+# Remove the preCICE address files
+rm -rfv precice-run
+rm -fv .*.address
+
+echo "Cleaning complete!"
+#------------------------------------------------------------------------------

FSI/flap_perp/ASTE/README.md

@@ -0,0 +1,95 @@
+## Overview
+
+This tutorial is an example how the Artificial Solver Testing Environment (aste) can be used. The idea is that you run a coupled simulation with two regular solvers and save the coupling data in every timestep. Then, this data is converted to aste-format. Finally, aste substitutes one of the solvers and the simulation can be "replayed" using only one solver and aste with the prior results.
+
+The base simulation for this tutorial is the OpenFOAM-FEniCS perpendicular flap tutorial. aste then substitutes OpenFOAM.
+
+## Requirements
+
+To run this tutorial you need to install the following components:
+- [preCICE](https://github.com/precice/precice/wiki/Get-preCICE)
+- [aste](https://github.com/precice/aste/tree/develop)
+- [FEniCS](https://fenicsproject.org/)
+- [FEniCS-Adapter](https://github.com/precice/fenics-adapter)
+- OpenFOAM, e.g. [OpenFOAM 7](https://openfoam.org/version/7/)
+- [OpenFOAM Adapter](https://github.com/precice/openfoam-adapter/wiki/Building) matching the OpenFOAM version.
+
+Make sure to add aste/build to the ```PATH``` such that the python scripts and ```preciceMap``` can be found from anywhere on your system.
+
+## Step-by-Step explanations
+
+### Generating vtk output during a simulation
+
+The base case for this tutorial is the OpenFOAM-FEniCS perpendicular flap tutorial. So let's start in the the root-directory of that simulation.
+To generate a vtk output in preCICE, add the following line to ```precice-config.xml``` after [line 41](https://github.com/precice/tutorials/blob/develop/FSI/flap_perp/OpenFOAM-FEniCS/precice-config.xml#L41) in the Solid participant:
+```<export:vtk directory="preCICE-output" />```
+
+Then, the simulation can be run in two terminals with ```bash ./runFluid``` and ```python3 Solid/perp_flap.py```.
+The exports are in the ```preCICE-output``` directory.
+
+### Converting the output to aste format.
+
+Copy the ```preCICE-output``` folder to the root directory of the aste tutorial. 
+To convert the files to the correct format, open the ```preCICE-output``` folder and run
+
+```precice_to_aste.py -n 500 -f Forces0 --datadim 3```
+
+### Replay of the simulation with aste
+
+#### The quick way to run 
+Copy perp_flap.py from OpenFOAM-FEniCS/Solid to aste/Solid
+
+run in two terminals ```python3 Solid/perp-flap.py``` and ```preciceMap -v -c precice-config.xml -p A --mesh preCICE-output/Solid-fenics --vectordata```
+
+Read on if you want to know what to change in the configuration files starting from the OpenFOAM-FEniCS tutorial.
+
+#### The way to do it yourself
+
+The aste-FEniCS simulation has its own ```precice-config.xml``` and ```precice-adapter-config-fsi-s.json```. 
+However, these can be generated by modifiying the files from the OpenFOAM-FEniCS case. 
+
+First, copy the files from the OpenFOAM-FEniCS tutorial into the respective directory of the aste-tutorial. 
+
+For ```precice-adapter-config-fsi-s.json``` we remove the line ```"write_data_name": "Displacements0",``` because we don't write data to aste.
+Plus, we change to the ```read_data_name``` from ```Forces0``` to ```Data``` since aste works with the general data name ```Data```.
+
+In ```precice-config.xml``` we need more changes:
+
+- Again, we rename ```Forces0``` to ```Data``` throughout the whole file.
+- We delete all lines where ```Displacements0``` occurs, since we are only coupling in the other direction.
+
+
+Then, you substitute the Fluid meshes and the Fluid participant with aste:
+
+- Rename ```Fluid-Mesh-Faces``` to ```MeshA``` throughout the whole file.
+- Delete the ```mesh``` ```"Fluid-Mesh-Nodes"``` and all lines where ```"Fluid-Mesh-Nodes"``` occurs.
+- Rename the participant ```Fluid``` to ```A```
+
+Last, you need to change the coupling scheme to an explicit scheme:
+
+- Change ```serial-implicit``` to ```serial-explicit```.
+- Explicit schemes don't iterate such that they don't need convergence measures, extrapolation or acceleration. Delete everything in the coulping scheme except for 
+```
+      <participants first="A" second="fenics" />
+      <max-time value="5.0" />
+      <time-window-size value="1.e-2" />
+      <exchange data="Data" mesh="Solid" from="A" to="fenics" />
+```
+
+### Run
+
+Finally, run the tutorial in two shells with ```python3 Solid/perp-flap.py``` and ```preciceMap -v -c precice-config.xml -p A --mesh preCICE-output/Solid-fenics --vectordata```.
+
+### Play Around
+
+If you want to explore more possibilities of the Replay-Mode here are some ideas:
+
+- Export the output at the Fluid participant. Instead of ```preCICE-output/Solid-fenics.*``` the files containing the forces will be ```preCICE-output/Fluid-Mesh-Faces-Fluid.*```
+
+- Export the forces in the [Calculix-Openfoam-tutorial](https://github.com/precice/tutorials/tree/master/FSI/flap_perp/OpenFOAM-CalculiX) with the same geometry and feed them to FEniCS with aste. 
+
+- If you are familiar with the OpenFOAM-Adapter you can also substitute FEniCS with aste such that aste writes previously exported displacements to OpenFOAM.
+
+
+
+