diff --git a/tutorials/fluidSolidInteraction/3dTube/Allclean b/tutorials/fluidSolidInteraction/3dTube/Allclean index 4cd967994..0e8d3b37a 100755 --- a/tutorials/fluidSolidInteraction/3dTube/Allclean +++ b/tutorials/fluidSolidInteraction/3dTube/Allclean @@ -20,6 +20,7 @@ else fi rm -rf case.foam postProcessing residuals +rm -rf radialDisplacement.pdf fsiConvergence.pdf axialDisplacement.pdf # Convert case version to FOAM EXTEND solids4Foam::convertCaseFormatFoamExtend . @@ -37,6 +38,6 @@ fi # Link files back to default with pimpleFluid for file in $(find ./0 ./constant ./system -name '*.pimpleFluid') -do +do ln -vnsf ${file##*/} ${file%.*} done diff --git a/tutorials/fluidSolidInteraction/3dTube/Allrun b/tutorials/fluidSolidInteraction/3dTube/Allrun index 80fd3a57a..d4c576535 100755 --- a/tutorials/fluidSolidInteraction/3dTube/Allrun +++ b/tutorials/fluidSolidInteraction/3dTube/Allrun @@ -56,3 +56,13 @@ solids4Foam::runApplication -s fluid blockMesh -region fluid # Run solver solids4Foam::runApplication solids4Foam + +# Create plots +if command -v gnuplot &> /dev/null +then + echo "Generating plots" + gnuplot -p displacement.gnuplot &> /dev/null + gnuplot -p fsiConvergence.gnuplot &> /dev/null +else + echo "Please install gnuplot if you would like to generate the plots" +fi diff --git a/tutorials/fluidSolidInteraction/3dTube/displacement.gnuplot b/tutorials/fluidSolidInteraction/3dTube/displacement.gnuplot new file mode 100644 index 000000000..7027fcf04 --- /dev/null +++ b/tutorials/fluidSolidInteraction/3dTube/displacement.gnuplot @@ -0,0 +1,27 @@ +set terminal pdfcairo enhanced color solid + +# Axial displacement + +set output "axialDisplacement.pdf" +set title "3dTube: Axial Displacement at Point A" +set xlabel "Time, t [s]" +set ylabel "Axial Displacement [mm]" +set grid +set key top right + +set yrange [-0.1:0.1] + +plot \ + "postProcessing/0/solidPointDisplacement_displacement.dat" u 1:(1000*$4)\ + w l lc "red" lw 0.5 lt 1 t "Axial displacement" + +# Radial displacement + +set output "radialDisplacement.pdf" +set title "3dTube: Radial Displacement at Point A" +set ylabel "Radial Displacement [mm]" +set yrange [-0.1:0.15] + +plot \ + "postProcessing/0/solidPointDisplacement_displacement.dat" u 1:(1000*$3)\ + w l lc "red" lw 0.5 lt 1 t "Radial displacement" diff --git a/tutorials/fluidSolidInteraction/3dTube/fsiConvergence.gnuplot b/tutorials/fluidSolidInteraction/3dTube/fsiConvergence.gnuplot new file mode 100644 index 000000000..63a9328e7 --- /dev/null +++ b/tutorials/fluidSolidInteraction/3dTube/fsiConvergence.gnuplot @@ -0,0 +1,14 @@ +set terminal pdfcairo enhanced color solid + +set output "fsiConvergence.pdf" + +set title "3dTube: Number of FSI iterations" +set xlabel "Time, t [s]" +set ylabel "Number of FSI iterations" +set grid +set key top right + +set yrange [0:80] + +plot \ + "postProcessing/0/fsiConvergenceData.dat" u 1:2 w lp ps 0.3 pt 5 lc "red" lw 0.5 lt 1 t "Number of iterations" diff --git a/tutorials/fluidSolidInteraction/3dTube/system/controlDict b/tutorials/fluidSolidInteraction/3dTube/system/controlDict index 5f9224b3a..1137c6800 100644 --- a/tutorials/fluidSolidInteraction/3dTube/system/controlDict +++ b/tutorials/fluidSolidInteraction/3dTube/system/controlDict @@ -76,6 +76,12 @@ functions CofR (0.5 0.1 0); region fluid; } + fsiConvData + { + type fsiConvergenceData; + region fluid; + } + } DebugSwitches