diff --git a/test/freeflow/navierstokes/donea/CMakeLists.txt b/test/freeflow/navierstokes/donea/CMakeLists.txt
index 394df47f735e2213b71c7304b38e532767bab4ee..71bdac37d33b9b31e9258210d4af12f93097d58a 100644
--- a/test/freeflow/navierstokes/donea/CMakeLists.txt
+++ b/test/freeflow/navierstokes/donea/CMakeLists.txt
@@ -1,3 +1,6 @@
+add_input_file_links()
+dune_symlink_to_source_files(FILES grids)
+
 dumux_add_test(NAME test_ff_stokes_donea_momentum
                LABELS freeflow navierstokes
                SOURCES main_momentum.cc
@@ -30,6 +33,31 @@ dumux_add_test(NAME test_ff_stokes_donea_nocaching
                              --command "${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea_nocaching params.input
                              -Problem.Name test_ff_stokes_donea_nocaching -Problem.UseNeumann false")
 
+dumux_add_test(NAME test_ff_stokes_donea_donut_nocaching
+              LABELS freeflow navierstokes
+              TARGET test_ff_stokes_donea_nocaching
+              CMAKE_GUARD "( dune-alugrid_FOUND AND HAVE_UMFPACK )"
+              COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+              CMD_ARGS       --script fuzzy
+                             --files ${CMAKE_SOURCE_DIR}/test/references/test_ff_stokes_donea_donut-reference.vtu
+                                     ${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea_donut_nocaching-00001.vtu
+                             --command "${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea_nocaching params.input
+                             -Grid.File grids/donut.dgf
+                             -Problem.Name test_ff_stokes_donea_donut_nocaching -Problem.UseNeumann false")
+
+dumux_add_test(NAME test_ff_stokes_donea_donut_twisted_nocaching
+              LABELS freeflow navierstokes
+              TARGET test_ff_stokes_donea_nocaching
+              CMAKE_GUARD "( dune-alugrid_FOUND AND HAVE_UMFPACK )"
+              COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+              CMD_ARGS       --script fuzzy
+                             --files ${CMAKE_SOURCE_DIR}/test/references/test_ff_stokes_donea_donut-reference.vtu
+                                     ${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea_donut_twisted_nocaching-00001.vtu
+                             --command "${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea_nocaching params.input
+                             -Grid.File grids/donut_twisted.dgf
+                             -Problem.Name test_ff_stokes_donea_donut_twisted_nocaching -Problem.UseNeumann false"
+                             --zeroThreshold {"connectivity":1e6})
+
 dumux_add_test(NAME test_ff_stokes_donea
               LABELS freeflow navierstokes
               SOURCES main.cc
@@ -42,4 +70,27 @@ dumux_add_test(NAME test_ff_stokes_donea
                              --command "${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea params.input
                              -Problem.Name test_ff_stokes_donea -Problem.UseNeumann false")
 
-dune_symlink_to_source_files(FILES "params.input")
+dumux_add_test(NAME test_ff_stokes_donea_donut
+               LABELS freeflow navierstokes
+               TARGET test_ff_stokes_donea
+               CMAKE_GUARD "( dune-alugrid_FOUND AND HAVE_UMFPACK )"
+               COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+               CMD_ARGS       --script fuzzy
+                               --files ${CMAKE_SOURCE_DIR}/test/references/test_ff_stokes_donea_donut-reference.vtu
+                                       ${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea_donut-00001.vtu
+                               --command "${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea params.input
+                               -Grid.File grids/donut.dgf
+                               -Problem.Name test_ff_stokes_donea_donut -Problem.UseNeumann false")
+
+dumux_add_test(NAME test_ff_stokes_donea_donut_twisted
+               LABELS freeflow navierstokes
+               TARGET test_ff_stokes_donea
+               CMAKE_GUARD "( dune-alugrid_FOUND AND HAVE_UMFPACK )"
+               COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+               CMD_ARGS       --script fuzzy
+                               --files ${CMAKE_SOURCE_DIR}/test/references/test_ff_stokes_donea_donut-reference.vtu
+                                       ${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea_donut_twisted-00001.vtu
+                               --command "${CMAKE_CURRENT_BINARY_DIR}/test_ff_stokes_donea params.input
+                               -Grid.File grids/donut_twisted.dgf
+                               -Problem.Name test_ff_stokes_donea_donut_twisted -Problem.UseNeumann false"
+                               --zeroThreshold {"connectivity":1e6})
diff --git a/test/freeflow/navierstokes/donea/grids/donut.dgf b/test/freeflow/navierstokes/donea/grids/donut.dgf
new file mode 100644
index 0000000000000000000000000000000000000000..716399c2c98929b3b4f93ab7126db31777ad47c1
--- /dev/null
+++ b/test/freeflow/navierstokes/donea/grids/donut.dgf
@@ -0,0 +1,205 @@
+DGF
+% Elements = 84  |  Vertices = 108
+
+VERTEX
+1.0000000000000001e-01 1.0000000000000001e-01
+2.0000000000000001e-01 1.0000000000000001e-01
+1.0000000000000001e-01 2.0000000000000001e-01
+2.0000000000000001e-01 2.0000000000000001e-01
+2.0000000000000001e-01 0.0000000000000000e+00
+3.0000000000000004e-01 0.0000000000000000e+00
+3.0000000000000004e-01 1.0000000000000001e-01
+4.0000000000000002e-01 0.0000000000000000e+00
+4.0000000000000002e-01 1.0000000000000001e-01
+5.0000000000000000e-01 0.0000000000000000e+00
+5.0000000000000000e-01 1.0000000000000001e-01
+3.0000000000000004e-01 2.0000000000000001e-01
+4.0000000000000002e-01 2.0000000000000001e-01
+5.0000000000000000e-01 2.0000000000000001e-01
+0.0000000000000000e+00 2.0000000000000001e-01
+0.0000000000000000e+00 3.0000000000000004e-01
+1.0000000000000001e-01 3.0000000000000004e-01
+0.0000000000000000e+00 4.0000000000000002e-01
+1.0000000000000001e-01 4.0000000000000002e-01
+2.0000000000000001e-01 3.0000000000000004e-01
+2.0000000000000001e-01 4.0000000000000002e-01
+0.0000000000000000e+00 5.0000000000000000e-01
+1.0000000000000001e-01 5.0000000000000000e-01
+2.0000000000000001e-01 5.0000000000000000e-01
+3.0000000000000004e-01 3.0000000000000004e-01
+4.0000000000000002e-01 3.0000000000000004e-01
+3.0000000000000004e-01 4.0000000000000002e-01
+4.0000000000000002e-01 4.0000000000000002e-01
+5.0000000000000000e-01 3.0000000000000004e-01
+5.0000000000000000e-01 4.0000000000000002e-01
+3.0000000000000004e-01 5.0000000000000000e-01
+4.0000000000000002e-01 5.0000000000000000e-01
+6.0000000000000009e-01 0.0000000000000000e+00
+6.0000000000000009e-01 1.0000000000000001e-01
+7.0000000000000007e-01 0.0000000000000000e+00
+7.0000000000000007e-01 1.0000000000000001e-01
+6.0000000000000009e-01 2.0000000000000001e-01
+7.0000000000000007e-01 2.0000000000000001e-01
+8.0000000000000004e-01 0.0000000000000000e+00
+8.0000000000000004e-01 1.0000000000000001e-01
+8.0000000000000004e-01 2.0000000000000001e-01
+9.0000000000000002e-01 1.0000000000000001e-01
+9.0000000000000002e-01 2.0000000000000001e-01
+6.0000000000000009e-01 3.0000000000000004e-01
+6.0000000000000009e-01 4.0000000000000002e-01
+7.0000000000000007e-01 3.0000000000000004e-01
+7.0000000000000007e-01 4.0000000000000002e-01
+6.0000000000000009e-01 5.0000000000000000e-01
+7.0000000000000007e-01 5.0000000000000000e-01
+8.0000000000000004e-01 3.0000000000000004e-01
+9.0000000000000002e-01 3.0000000000000004e-01
+8.0000000000000004e-01 4.0000000000000002e-01
+9.0000000000000002e-01 4.0000000000000002e-01
+1.0000000000000000e+00 2.0000000000000001e-01
+1.0000000000000000e+00 3.0000000000000004e-01
+1.0000000000000000e+00 4.0000000000000002e-01
+8.0000000000000004e-01 5.0000000000000000e-01
+9.0000000000000002e-01 5.0000000000000000e-01
+1.0000000000000000e+00 5.0000000000000000e-01
+0.0000000000000000e+00 6.0000000000000009e-01
+1.0000000000000001e-01 6.0000000000000009e-01
+2.0000000000000001e-01 6.0000000000000009e-01
+0.0000000000000000e+00 7.0000000000000007e-01
+1.0000000000000001e-01 7.0000000000000007e-01
+2.0000000000000001e-01 7.0000000000000007e-01
+3.0000000000000004e-01 6.0000000000000009e-01
+4.0000000000000002e-01 6.0000000000000009e-01
+3.0000000000000004e-01 7.0000000000000007e-01
+4.0000000000000002e-01 7.0000000000000007e-01
+5.0000000000000000e-01 6.0000000000000009e-01
+5.0000000000000000e-01 7.0000000000000007e-01
+0.0000000000000000e+00 8.0000000000000004e-01
+1.0000000000000001e-01 8.0000000000000004e-01
+2.0000000000000001e-01 8.0000000000000004e-01
+1.0000000000000001e-01 9.0000000000000002e-01
+2.0000000000000001e-01 9.0000000000000002e-01
+3.0000000000000004e-01 8.0000000000000004e-01
+4.0000000000000002e-01 8.0000000000000004e-01
+3.0000000000000004e-01 9.0000000000000002e-01
+4.0000000000000002e-01 9.0000000000000002e-01
+5.0000000000000000e-01 8.0000000000000004e-01
+5.0000000000000000e-01 9.0000000000000002e-01
+2.0000000000000001e-01 1.0000000000000000e+00
+3.0000000000000004e-01 1.0000000000000000e+00
+4.0000000000000002e-01 1.0000000000000000e+00
+5.0000000000000000e-01 1.0000000000000000e+00
+6.0000000000000009e-01 6.0000000000000009e-01
+7.0000000000000007e-01 6.0000000000000009e-01
+6.0000000000000009e-01 7.0000000000000007e-01
+7.0000000000000007e-01 7.0000000000000007e-01
+8.0000000000000004e-01 6.0000000000000009e-01
+9.0000000000000002e-01 6.0000000000000009e-01
+1.0000000000000000e+00 6.0000000000000009e-01
+8.0000000000000004e-01 7.0000000000000007e-01
+9.0000000000000002e-01 7.0000000000000007e-01
+1.0000000000000000e+00 7.0000000000000007e-01
+6.0000000000000009e-01 8.0000000000000004e-01
+6.0000000000000009e-01 9.0000000000000002e-01
+7.0000000000000007e-01 8.0000000000000004e-01
+7.0000000000000007e-01 9.0000000000000002e-01
+6.0000000000000009e-01 1.0000000000000000e+00
+7.0000000000000007e-01 1.0000000000000000e+00
+8.0000000000000004e-01 8.0000000000000004e-01
+9.0000000000000002e-01 8.0000000000000004e-01
+8.0000000000000004e-01 9.0000000000000002e-01
+9.0000000000000002e-01 9.0000000000000002e-01
+1.0000000000000000e+00 8.0000000000000004e-01
+8.0000000000000004e-01 1.0000000000000000e+00
+#
+
+CUBE
+0 1 2 3
+4 5 1 6
+5 7 6 8
+7 9 8 10
+1 6 3 11
+6 8 11 12
+8 10 12 13
+14 2 15 16
+15 16 17 18
+2 3 16 19
+16 19 18 20
+17 18 21 22
+18 20 22 23
+3 11 19 24
+11 12 24 25
+19 24 20 26
+24 25 26 27
+12 13 25 28
+25 28 27 29
+20 26 23 30
+26 27 30 31
+9 32 10 33
+32 34 33 35
+10 33 13 36
+33 35 36 37
+34 38 35 39
+35 39 37 40
+39 41 40 42
+13 36 28 43
+28 43 29 44
+36 37 43 45
+43 45 44 46
+44 46 47 48
+37 40 45 49
+40 42 49 50
+45 49 46 51
+49 50 51 52
+42 53 50 54
+50 54 52 55
+46 51 48 56
+51 52 56 57
+52 55 57 58
+21 22 59 60
+22 23 60 61
+59 60 62 63
+60 61 63 64
+23 30 61 65
+30 31 65 66
+61 65 64 67
+65 66 67 68
+66 69 68 70
+62 63 71 72
+63 64 72 73
+72 73 74 75
+64 67 73 76
+67 68 76 77
+73 76 75 78
+76 77 78 79
+68 70 77 80
+77 80 79 81
+75 78 82 83
+78 79 83 84
+79 81 84 85
+47 48 86 87
+69 86 70 88
+86 87 88 89
+48 56 87 90
+56 57 90 91
+57 58 91 92
+87 90 89 93
+90 91 93 94
+91 92 94 95
+70 88 80 96
+80 96 81 97
+88 89 96 98
+96 98 97 99
+81 97 85 100
+97 99 100 101
+89 93 98 102
+93 94 102 103
+98 102 99 104
+102 103 104 105
+94 95 103 106
+99 104 101 107
+#
+BOUNDARYDOMAIN
+default 1
+#
+
+#
diff --git a/test/freeflow/navierstokes/donea/grids/donut_twisted.dgf b/test/freeflow/navierstokes/donea/grids/donut_twisted.dgf
new file mode 100644
index 0000000000000000000000000000000000000000..2fa3f76c3db837284b704f6d9764b8395125cb3a
--- /dev/null
+++ b/test/freeflow/navierstokes/donea/grids/donut_twisted.dgf
@@ -0,0 +1,205 @@
+DGF
+% Elements = 84  |  Vertices = 108
+
+VERTEX
+1.0000000000000001e-01 1.0000000000000001e-01
+2.0000000000000001e-01 1.0000000000000001e-01
+1.0000000000000001e-01 2.0000000000000001e-01
+2.0000000000000001e-01 2.0000000000000001e-01
+2.0000000000000001e-01 0.0000000000000000e+00
+3.0000000000000004e-01 0.0000000000000000e+00
+3.0000000000000004e-01 1.0000000000000001e-01
+4.0000000000000002e-01 0.0000000000000000e+00
+4.0000000000000002e-01 1.0000000000000001e-01
+5.0000000000000000e-01 0.0000000000000000e+00
+5.0000000000000000e-01 1.0000000000000001e-01
+3.0000000000000004e-01 2.0000000000000001e-01
+4.0000000000000002e-01 2.0000000000000001e-01
+5.0000000000000000e-01 2.0000000000000001e-01
+0.0000000000000000e+00 2.0000000000000001e-01
+0.0000000000000000e+00 3.0000000000000004e-01
+1.0000000000000001e-01 3.0000000000000004e-01
+0.0000000000000000e+00 4.0000000000000002e-01
+1.0000000000000001e-01 4.0000000000000002e-01
+2.0000000000000001e-01 3.0000000000000004e-01
+2.0000000000000001e-01 4.0000000000000002e-01
+0.0000000000000000e+00 5.0000000000000000e-01
+1.0000000000000001e-01 5.0000000000000000e-01
+2.0000000000000001e-01 5.0000000000000000e-01
+3.0000000000000004e-01 3.0000000000000004e-01
+4.0000000000000002e-01 3.0000000000000004e-01
+3.0000000000000004e-01 4.0000000000000002e-01
+4.0000000000000002e-01 4.0000000000000002e-01
+5.0000000000000000e-01 3.0000000000000004e-01
+5.0000000000000000e-01 4.0000000000000002e-01
+3.0000000000000004e-01 5.0000000000000000e-01
+4.0000000000000002e-01 5.0000000000000000e-01
+6.0000000000000009e-01 0.0000000000000000e+00
+6.0000000000000009e-01 1.0000000000000001e-01
+7.0000000000000007e-01 0.0000000000000000e+00
+7.0000000000000007e-01 1.0000000000000001e-01
+6.0000000000000009e-01 2.0000000000000001e-01
+7.0000000000000007e-01 2.0000000000000001e-01
+8.0000000000000004e-01 0.0000000000000000e+00
+8.0000000000000004e-01 1.0000000000000001e-01
+8.0000000000000004e-01 2.0000000000000001e-01
+9.0000000000000002e-01 1.0000000000000001e-01
+9.0000000000000002e-01 2.0000000000000001e-01
+6.0000000000000009e-01 3.0000000000000004e-01
+6.0000000000000009e-01 4.0000000000000002e-01
+7.0000000000000007e-01 3.0000000000000004e-01
+7.0000000000000007e-01 4.0000000000000002e-01
+6.0000000000000009e-01 5.0000000000000000e-01
+7.0000000000000007e-01 5.0000000000000000e-01
+8.0000000000000004e-01 3.0000000000000004e-01
+9.0000000000000002e-01 3.0000000000000004e-01
+8.0000000000000004e-01 4.0000000000000002e-01
+9.0000000000000002e-01 4.0000000000000002e-01
+1.0000000000000000e+00 2.0000000000000001e-01
+1.0000000000000000e+00 3.0000000000000004e-01
+1.0000000000000000e+00 4.0000000000000002e-01
+8.0000000000000004e-01 5.0000000000000000e-01
+9.0000000000000002e-01 5.0000000000000000e-01
+1.0000000000000000e+00 5.0000000000000000e-01
+0.0000000000000000e+00 6.0000000000000009e-01
+1.0000000000000001e-01 6.0000000000000009e-01
+2.0000000000000001e-01 6.0000000000000009e-01
+0.0000000000000000e+00 7.0000000000000007e-01
+1.0000000000000001e-01 7.0000000000000007e-01
+2.0000000000000001e-01 7.0000000000000007e-01
+3.0000000000000004e-01 6.0000000000000009e-01
+4.0000000000000002e-01 6.0000000000000009e-01
+3.0000000000000004e-01 7.0000000000000007e-01
+4.0000000000000002e-01 7.0000000000000007e-01
+5.0000000000000000e-01 6.0000000000000009e-01
+5.0000000000000000e-01 7.0000000000000007e-01
+0.0000000000000000e+00 8.0000000000000004e-01
+1.0000000000000001e-01 8.0000000000000004e-01
+2.0000000000000001e-01 8.0000000000000004e-01
+1.0000000000000001e-01 9.0000000000000002e-01
+2.0000000000000001e-01 9.0000000000000002e-01
+3.0000000000000004e-01 8.0000000000000004e-01
+4.0000000000000002e-01 8.0000000000000004e-01
+3.0000000000000004e-01 9.0000000000000002e-01
+4.0000000000000002e-01 9.0000000000000002e-01
+5.0000000000000000e-01 8.0000000000000004e-01
+5.0000000000000000e-01 9.0000000000000002e-01
+2.0000000000000001e-01 1.0000000000000000e+00
+3.0000000000000004e-01 1.0000000000000000e+00
+4.0000000000000002e-01 1.0000000000000000e+00
+5.0000000000000000e-01 1.0000000000000000e+00
+6.0000000000000009e-01 6.0000000000000009e-01
+7.0000000000000007e-01 6.0000000000000009e-01
+6.0000000000000009e-01 7.0000000000000007e-01
+7.0000000000000007e-01 7.0000000000000007e-01
+8.0000000000000004e-01 6.0000000000000009e-01
+9.0000000000000002e-01 6.0000000000000009e-01
+1.0000000000000000e+00 6.0000000000000009e-01
+8.0000000000000004e-01 7.0000000000000007e-01
+9.0000000000000002e-01 7.0000000000000007e-01
+1.0000000000000000e+00 7.0000000000000007e-01
+6.0000000000000009e-01 8.0000000000000004e-01
+6.0000000000000009e-01 9.0000000000000002e-01
+7.0000000000000007e-01 8.0000000000000004e-01
+7.0000000000000007e-01 9.0000000000000002e-01
+6.0000000000000009e-01 1.0000000000000000e+00
+7.0000000000000007e-01 1.0000000000000000e+00
+8.0000000000000004e-01 8.0000000000000004e-01
+9.0000000000000002e-01 8.0000000000000004e-01
+8.0000000000000004e-01 9.0000000000000002e-01
+9.0000000000000002e-01 9.0000000000000002e-01
+1.0000000000000000e+00 8.0000000000000004e-01
+8.0000000000000004e-01 1.0000000000000000e+00
+#
+
+CUBE
+0 1 2 3
+4 5 1 6
+5 7 6 8
+7 9 8 10
+1 6 3 11
+11 12 6 8 % this element is twisted: original order is 6 8 11 12
+8 10 12 13
+14 2 15 16
+15 16 17 18
+2 3 16 19
+16 19 18 20
+17 18 21 22
+18 20 22 23
+3 11 19 24
+11 12 24 25
+19 24 20 26
+24 25 26 27
+12 13 25 28
+25 28 27 29
+20 26 23 30
+26 27 30 31
+9 32 10 33
+32 34 33 35
+10 33 13 36
+33 35 36 37
+34 38 35 39
+35 39 37 40
+39 41 40 42
+13 36 28 43
+28 43 29 44
+36 37 43 45
+43 45 44 46
+44 46 47 48
+37 40 45 49
+40 42 49 50
+45 49 46 51
+49 50 51 52
+42 53 50 54
+50 54 52 55
+46 51 48 56
+51 52 56 57
+52 55 57 58
+21 22 59 60
+22 23 60 61
+59 60 62 63
+60 61 63 64
+23 30 61 65
+30 31 65 66
+61 65 64 67
+65 66 67 68
+66 69 68 70
+62 63 71 72
+63 64 72 73
+72 73 74 75
+64 67 73 76
+67 68 76 77
+73 76 75 78
+76 77 78 79
+68 70 77 80
+77 80 79 81
+75 78 82 83
+78 79 83 84
+79 81 84 85
+47 48 86 87
+69 86 70 88
+86 87 88 89
+48 56 87 90
+56 57 90 91
+57 58 91 92
+87 90 89 93
+90 91 93 94
+91 92 94 95
+70 88 80 96
+80 96 81 97
+88 89 96 98
+96 98 97 99
+81 97 85 100
+97 99 100 101
+89 93 98 102
+93 94 102 103
+98 102 99 104
+102 103 104 105
+94 95 103 106
+99 104 101 107
+#
+BOUNDARYDOMAIN
+default 1
+#
+
+#
diff --git a/test/freeflow/navierstokes/donea/properties.hh b/test/freeflow/navierstokes/donea/properties.hh
index 5f9e89b2a996fac671e363f880904cfc4e1d115c..9247d3c52fe291dd155a654281f91653bce1b193 100644
--- a/test/freeflow/navierstokes/donea/properties.hh
+++ b/test/freeflow/navierstokes/donea/properties.hh
@@ -28,7 +28,11 @@
 #define ENABLECACHING 0
 #endif
 
+#if HAVE_DUNE_ALUGRID
+#include <dune/alugrid/grid.hh>
+#else
 #include <dune/grid/yaspgrid.hh>
+#endif
 
 #include <dumux/discretization/fcstaggered.hh>
 #include <dumux/discretization/cctpfa.hh>
@@ -63,7 +67,11 @@ struct FluidSystem<TypeTag, TTag::DoneaTest>
 
 template<class TypeTag>
 struct Grid<TypeTag, TTag::DoneaTest>
+#if HAVE_DUNE_ALUGRID
+{ using type = Dune::ALUGrid<2, 2, Dune::cube, Dune::nonconforming>; };
+#else
 { using type = Dune::YaspGrid<2>; };
+#endif
 
 template<class TypeTag>
 struct EnableGridGeometryCache<TypeTag, TTag::DoneaTest> { static constexpr bool value = ENABLECACHING; };
diff --git a/test/references/test_ff_stokes_donea_donut-reference.vtu b/test/references/test_ff_stokes_donea_donut-reference.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ef7df03702382f4b6bd129f52b7974de3065e066
--- /dev/null
+++ b/test/references/test_ff_stokes_donea_donut-reference.vtu
@@ -0,0 +1,172 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <UnstructuredGrid>
+    <Piece NumberOfCells="84" NumberOfPoints="108">
+      <CellData Scalars="p" Vectors="velocity_liq (m/s)">
+        <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
+          0.1275 0.1875 0.2275 0.2475 0.185876 0.226068 0.246941 0.0475 0.0475 0.129124 0.128932 0.0475
+          0.128059 0.1875 0.227043 0.187957 0.2275 0.247105 0.2475 0.187895 0.2275 0.2475 0.2275 0.248059
+          0.228932 0.1875 0.189124 0.1275 0.247895 0.2475 0.227957 0.2275 0.2275 0.1875 0.125876 0.187043
+          0.126068 0.0475 0.0475 0.187105 0.126941 0.0475 0.0475 0.126941 0.0475 0.126068 0.187105 0.2275
+          0.187043 0.2275 0.2475 0.0475 0.125876 0.1275 0.1875 0.227957 0.189124 0.228932 0.247895 0.248059
+          0.1875 0.2275 0.2475 0.2275 0.2475 0.2275 0.187895 0.128059 0.0475 0.187957 0.128932 0.0475
+          0.247105 0.246941 0.227043 0.226068 0.2475 0.2275 0.1875 0.129124 0.185876 0.1275 0.0475 0.1875
+        </DataArray>
+        <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+        </DataArray>
+        <DataArray type="Float32" Name="velocity_liq (m/s)" NumberOfComponents="3" format="ascii">
+          0.00313589 -0.00313589 0 0.00319571 -0.000834885 0 0.00488548 -0.000605641 0 0.00581017 -0.000217602 0
+          0.00653201 -0.00337585 0 0.00947902 -0.00245221 0 0.0111676 -0.000882814 0 0.000834885 -0.00319571 0
+          0.000605641 -0.00488548 0 0.00337585 -0.00653201 0 0.00245221 -0.00947902 0 0.000217602 -0.00581017 0
+          0.000882814 -0.0111676 0 0.00683621 -0.00683621 0 0.00986338 -0.00496603 0 0.00496603 -0.00986338 0
+          0.00716411 -0.00716411 0 0.0115909 -0.00178508 0 0.00841898 -0.00254547 0 0.00178508 -0.0115909 0
+          0.00254547 -0.00841898 0 0.00581017 0.000217602 0 0.00488548 0.000605641 0 0.0111676 0.000882814 0
+          0.00947902 0.00245221 0 0.00319571 0.000834885 0 0.00653201 0.00337585 0 0.00313589 0.00313589 0
+          0.0115909 0.00178508 0 0.00841898 0.00254547 0 0.00986338 0.00496603 0 0.00716411 0.00716411 0
+          0.00254547 0.00841898 0 0.00683621 0.00683621 0 0.00337585 0.00653201 0 0.00496603 0.00986338 0
+          0.00245221 0.00947902 0 0.000834885 0.00319571 0 0.000605641 0.00488548 0 0.00178508 0.0115909 0
+          0.000882814 0.0111676 0 0.000217602 0.00581017 0 -0.000217602 -0.00581017 0 -0.000882814 -0.0111676 0
+          -0.000605641 -0.00488548 0 -0.00245221 -0.00947902 0 -0.00178508 -0.0115909 0 -0.00254547 -0.00841898 0
+          -0.00496603 -0.00986338 0 -0.00716411 -0.00716411 0 -0.00841898 -0.00254547 0 -0.000834885 -0.00319571 0
+          -0.00337585 -0.00653201 0 -0.00313589 -0.00313589 0 -0.00683621 -0.00683621 0 -0.00986338 -0.00496603 0
+          -0.00653201 -0.00337585 0 -0.00947902 -0.00245221 0 -0.0115909 -0.00178508 0 -0.0111676 -0.000882814 0
+          -0.00319571 -0.000834885 0 -0.00488548 -0.000605641 0 -0.00581017 -0.000217602 0 -0.00254547 0.00841898 0
+          -0.00841898 0.00254547 0 -0.00716411 0.00716411 0 -0.00178508 0.0115909 0 -0.000882814 0.0111676 0
+          -0.000217602 0.00581017 0 -0.00496603 0.00986338 0 -0.00245221 0.00947902 0 -0.000605641 0.00488548 0
+          -0.0115909 0.00178508 0 -0.0111676 0.000882814 0 -0.00986338 0.00496603 0 -0.00947902 0.00245221 0
+          -0.00581017 0.000217602 0 -0.00488548 0.000605641 0 -0.00683621 0.00683621 0 -0.00337585 0.00653201 0
+          -0.00653201 0.00337585 0 -0.00313589 0.00313589 0 -0.000834885 0.00319571 0 -0.00319571 0.000834885 0
+        </DataArray>
+        <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+        </DataArray>
+        <DataArray type="Float32" Name="pressureExact" NumberOfComponents="1" format="ascii">
+          0.1275 0.1875 0.2275 0.2475 0.1875 0.2275 0.2475 0.0475 0.0475 0.1275 0.1275 0.0475
+          0.1275 0.1875 0.2275 0.1875 0.2275 0.2475 0.2475 0.1875 0.2275 0.2475 0.2275 0.2475
+          0.2275 0.1875 0.1875 0.1275 0.2475 0.2475 0.2275 0.2275 0.2275 0.1875 0.1275 0.1875
+          0.1275 0.0475 0.0475 0.1875 0.1275 0.0475 0.0475 0.1275 0.0475 0.1275 0.1875 0.2275
+          0.1875 0.2275 0.2475 0.0475 0.1275 0.1275 0.1875 0.2275 0.1875 0.2275 0.2475 0.2475
+          0.1875 0.2275 0.2475 0.2275 0.2475 0.2275 0.1875 0.1275 0.0475 0.1875 0.1275 0.0475
+          0.2475 0.2475 0.2275 0.2275 0.2475 0.2275 0.1875 0.1275 0.1875 0.1275 0.0475 0.1875
+        </DataArray>
+        <DataArray type="Float32" Name="velocityExact" NumberOfComponents="3" format="ascii">
+          0.00290174 -0.00290174 0 0.00300586 -0.000423047 0 0.00442516 -0.000307978 0 0.00523741 -0.000111684 0
+          0.00627539 -0.00304805 0 0.00923849 -0.00221898 0 0.0109342 -0.000804684 0 0.000423047 -0.00300586 0
+          0.000307978 -0.00442516 0 0.00304805 -0.00627539 0 0.00221898 -0.00923849 0 0.000111684 -0.00523741 0
+          0.000804684 -0.0109342 0 0.0065918 -0.0065918 0 0.0097043 -0.00479883 0 0.00479883 -0.0097043 0
+          0.00706473 -0.00706473 0 0.0114855 -0.00174023 0 0.00836148 -0.00256193 0 0.00174023 -0.0114855 0
+          0.00256193 -0.00836148 0 0.00523741 0.000111684 0 0.00442516 0.000307978 0 0.0109342 0.000804684 0
+          0.00923849 0.00221898 0 0.00300586 0.000423047 0 0.00627539 0.00304805 0 0.00290174 0.00290174 0
+          0.0114855 0.00174023 0 0.00836148 0.00256193 0 0.0097043 0.00479883 0 0.00706473 0.00706473 0
+          0.00256193 0.00836148 0 0.0065918 0.0065918 0 0.00304805 0.00627539 0 0.00479883 0.0097043 0
+          0.00221898 0.00923849 0 0.000423047 0.00300586 0 0.000307978 0.00442516 0 0.00174023 0.0114855 0
+          0.000804684 0.0109342 0 0.000111684 0.00523741 0 -0.000111684 -0.00523741 0 -0.000804684 -0.0109342 0
+          -0.000307978 -0.00442516 0 -0.00221898 -0.00923849 0 -0.00174023 -0.0114855 0 -0.00256193 -0.00836148 0
+          -0.00479883 -0.0097043 0 -0.00706473 -0.00706473 0 -0.00836148 -0.00256193 0 -0.000423047 -0.00300586 0
+          -0.00304805 -0.00627539 0 -0.00290174 -0.00290174 0 -0.0065918 -0.0065918 0 -0.0097043 -0.00479883 0
+          -0.00627539 -0.00304805 0 -0.00923849 -0.00221898 0 -0.0114855 -0.00174023 0 -0.0109342 -0.000804684 0
+          -0.00300586 -0.000423047 0 -0.00442516 -0.000307978 0 -0.00523741 -0.000111684 0 -0.00256193 0.00836148 0
+          -0.00836148 0.00256193 0 -0.00706473 0.00706473 0 -0.00174023 0.0114855 0 -0.000804684 0.0109342 0
+          -0.000111684 0.00523741 0 -0.00479883 0.0097043 0 -0.00221898 0.00923849 0 -0.000307978 0.00442516 0
+          -0.0114855 0.00174023 0 -0.0109342 0.000804684 0 -0.0097043 0.00479883 0 -0.00923849 0.00221898 0
+          -0.00523741 0.000111684 0 -0.00442516 0.000307978 0 -0.0065918 0.0065918 0 -0.00304805 0.00627539 0
+          -0.00627539 0.00304805 0 -0.00290174 0.00290174 0 -0.000423047 0.00300586 0 -0.00300586 0.000423047 0
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+          0.1 0.1 0 0.2 0.1 0 0.1 0.2 0 0.2 0.2 0
+          0.2 0 0 0.3 0 0 0.3 0.1 0 0.4 0 0
+          0.4 0.1 0 0.5 0 0 0.5 0.1 0 0.3 0.2 0
+          0.4 0.2 0 0.5 0.2 0 0 0.2 0 0 0.3 0
+          0.1 0.3 0 0 0.4 0 0.1 0.4 0 0.2 0.3 0
+          0.2 0.4 0 0 0.5 0 0.1 0.5 0 0.2 0.5 0
+          0.3 0.3 0 0.4 0.3 0 0.3 0.4 0 0.4 0.4 0
+          0.5 0.3 0 0.5 0.4 0 0.3 0.5 0 0.4 0.5 0
+          0.6 0 0 0.6 0.1 0 0.7 0 0 0.7 0.1 0
+          0.6 0.2 0 0.7 0.2 0 0.8 0 0 0.8 0.1 0
+          0.8 0.2 0 0.9 0.1 0 0.9 0.2 0 0.6 0.3 0
+          0.6 0.4 0 0.7 0.3 0 0.7 0.4 0 0.6 0.5 0
+          0.7 0.5 0 0.8 0.3 0 0.9 0.3 0 0.8 0.4 0
+          0.9 0.4 0 1 0.2 0 1 0.3 0 1 0.4 0
+          0.8 0.5 0 0.9 0.5 0 1 0.5 0 0 0.6 0
+          0.1 0.6 0 0.2 0.6 0 0 0.7 0 0.1 0.7 0
+          0.2 0.7 0 0.3 0.6 0 0.4 0.6 0 0.3 0.7 0
+          0.4 0.7 0 0.5 0.6 0 0.5 0.7 0 0 0.8 0
+          0.1 0.8 0 0.2 0.8 0 0.1 0.9 0 0.2 0.9 0
+          0.3 0.8 0 0.4 0.8 0 0.3 0.9 0 0.4 0.9 0
+          0.5 0.8 0 0.5 0.9 0 0.2 1 0 0.3 1 0
+          0.4 1 0 0.5 1 0 0.6 0.6 0 0.7 0.6 0
+          0.6 0.7 0 0.7 0.7 0 0.8 0.6 0 0.9 0.6 0
+          1 0.6 0 0.8 0.7 0 0.9 0.7 0 1 0.7 0
+          0.6 0.8 0 0.6 0.9 0 0.7 0.8 0 0.7 0.9 0
+          0.6 1 0 0.7 1 0 0.8 0.8 0 0.9 0.8 0
+          0.8 0.9 0 0.9 0.9 0 1 0.8 0 0.8 1 0
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+          0 1 3 2 4 5 6 1 5 7 8 6
+          7 9 10 8 1 6 11 3 6 8 12 11
+          8 10 13 12 14 2 16 15 15 16 18 17
+          2 3 19 16 16 19 20 18 17 18 22 21
+          18 20 23 22 3 11 24 19 11 12 25 24
+          19 24 26 20 24 25 27 26 12 13 28 25
+          25 28 29 27 20 26 30 23 26 27 31 30
+          9 32 33 10 32 34 35 33 10 33 36 13
+          33 35 37 36 34 38 39 35 35 39 40 37
+          39 41 42 40 13 36 43 28 28 43 44 29
+          36 37 45 43 43 45 46 44 44 46 48 47
+          37 40 49 45 40 42 50 49 45 49 51 46
+          49 50 52 51 42 53 54 50 50 54 55 52
+          46 51 56 48 51 52 57 56 52 55 58 57
+          21 22 60 59 22 23 61 60 59 60 63 62
+          60 61 64 63 23 30 65 61 30 31 66 65
+          61 65 67 64 65 66 68 67 66 69 70 68
+          62 63 72 71 63 64 73 72 72 73 75 74
+          64 67 76 73 67 68 77 76 73 76 78 75
+          76 77 79 78 68 70 80 77 77 80 81 79
+          75 78 83 82 78 79 84 83 79 81 85 84
+          47 48 87 86 69 86 88 70 86 87 89 88
+          48 56 90 87 56 57 91 90 57 58 92 91
+          87 90 93 89 90 91 94 93 91 92 95 94
+          70 88 96 80 80 96 97 81 88 89 98 96
+          96 98 99 97 81 97 100 85 97 99 101 100
+          89 93 102 98 93 94 103 102 98 102 104 99
+          102 103 105 104 94 95 106 103 99 104 107 101
+        </DataArray>
+        <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+          4 8 12 16 20 24 28 32 36 40 44 48
+          52 56 60 64 68 72 76 80 84 88 92 96
+          100 104 108 112 116 120 124 128 132 136 140 144
+          148 152 156 160 164 168 172 176 180 184 188 192
+          196 200 204 208 212 216 220 224 228 232 236 240
+          244 248 252 256 260 264 268 272 276 280 284 288
+          292 296 300 304 308 312 316 320 324 328 332 336
+        </DataArray>
+        <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>