From e414d20142e6e10b8fcafb099028643add3c3480 Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Fri, 29 Jun 2018 17:15:22 +0200
Subject: [PATCH] [test][kepsilon] Implement test for wallfunctions

---
 test/freeflow/ransnc/CMakeLists.txt           |  18 +-
 .../test_flatplate2cni_wallfunction.input     |  33 ++
 ...st_flatplate_kepsilon2cni_wallfunction.vtu | 292 ++++++++++++++++++
 3 files changed, 340 insertions(+), 3 deletions(-)
 create mode 100644 test/freeflow/ransnc/test_flatplate2cni_wallfunction.input
 create mode 100644 test/references/test_flatplate_kepsilon2cni_wallfunction.vtu

diff --git a/test/freeflow/ransnc/CMakeLists.txt b/test/freeflow/ransnc/CMakeLists.txt
index f2db9e649d..7db4c65c00 100644
--- a/test/freeflow/ransnc/CMakeLists.txt
+++ b/test/freeflow/ransnc/CMakeLists.txt
@@ -5,8 +5,11 @@ dune_add_test(NAME test_flatplate_kepsilon2c
               COMPILE_ONLY CMAKE_GUARD HAVE_UMFPACK)
 target_compile_definitions(test_flatplate_kepsilon2c PUBLIC "KEPSILON=1")
 
-dune_add_test(NAME test_flatplate_kepsilon2cni
-              SOURCES test_flatplate.cc
+add_executable(test_flatplate_kepsilon2cni EXCLUDE_FROM_ALL test_flatplate.cc)
+target_compile_definitions(test_flatplate_kepsilon2cni PUBLIC "KEPSILON=1" "NONISOTHERMAL=1")
+
+dune_add_test(NAME test_flatplate_kepsilon2cni_twolayer
+              TARGET test_flatplate_kepsilon2cni
               CMAKE_GUARD HAVE_UMFPACK
               COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
               CMD_ARGS       --script fuzzy
@@ -14,7 +17,16 @@ dune_add_test(NAME test_flatplate_kepsilon2cni
                                      ${CMAKE_CURRENT_BINARY_DIR}/test_flatplate2cni-00036.vtu
                              --command "${CMAKE_CURRENT_BINARY_DIR}/test_flatplate_kepsilon2cni
                                        test_flatplate2cni.input")
-target_compile_definitions(test_flatplate_kepsilon2cni PUBLIC "KEPSILON=1" "NONISOTHERMAL=1")
+
+dune_add_test(NAME test_flatplate_kepsilon2cni_wallfunction
+              TARGET test_flatplate_kepsilon2cni
+              CMAKE_GUARD HAVE_UMFPACK
+              COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+              CMD_ARGS       --script fuzzy
+                             --files ${CMAKE_SOURCE_DIR}/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu
+                                     ${CMAKE_CURRENT_BINARY_DIR}/test_flatplate2cni_wallfunction-00025.vtu
+                             --command "${CMAKE_CURRENT_BINARY_DIR}/test_flatplate_kepsilon2cni
+                                       test_flatplate2cni_wallfunction.input")
 
 dune_add_test(NAME test_flatplate_komega2c
               SOURCES test_flatplate.cc
diff --git a/test/freeflow/ransnc/test_flatplate2cni_wallfunction.input b/test/freeflow/ransnc/test_flatplate2cni_wallfunction.input
new file mode 100644
index 0000000000..53b88eeca2
--- /dev/null
+++ b/test/freeflow/ransnc/test_flatplate2cni_wallfunction.input
@@ -0,0 +1,33 @@
+[TimeLoop]
+DtInitial = 1e-2 # [s]
+TEnd = 1000 # [s]
+
+[Grid]
+Verbosity = true
+Positions0 = 0.0 6.0
+Positions1 = 0.0 0.5
+Cells0 = 10
+Cells1 = 5
+
+[Problem]
+Name = test_flatplate2cni_wallfunction # name passed to the output routines
+InletVelocity = 0.1 # [m/s]
+EnableGravity = false
+
+[RANS]
+TurbulentSchmidtNumber = 0.7
+TurbulentPrandtlNumber = 0.85
+
+[KEpsilon]
+YPlusThreshold = 10. # should be small (10-30) for coarse grids
+
+[Assembly]
+NumericDifferenceMethod = 0
+
+[Newton]
+MaxSteps = 10
+TargetSteps = 8
+MaxRelativeShift = 1e-5
+
+[Vtk]
+AddVelocity = true
diff --git a/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu b/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu
new file mode 100644
index 0000000000..27b0b3c411
--- /dev/null
+++ b/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu
@@ -0,0 +1,292 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <UnstructuredGrid>
+    <Piece NumberOfCells="50" NumberOfPoints="66">
+      <CellData Scalars="p" Vectors="dv_x/dx_">
+        <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000
+        </DataArray>
+        <DataArray type="Float32" Name="rhoMolar" NumberOfComponents="1" format="ascii">
+          42.3689 42.27 42.179 42.0948 42.0167 41.9441 41.8765 41.8134 41.7544 41.6995 42.4677 42.4525
+          42.4329 42.4103 42.3854 42.3591 42.3319 42.304 42.2758 42.249 42.4757 42.4736 42.4702 42.4654
+          42.4594 42.4523 42.4442 42.4352 42.4255 42.4156 42.4764 42.4762 42.4756 42.4748 42.4736 42.472
+          42.4699 42.4675 42.4646 42.4615 42.4765 42.4765 42.4764 42.4762 42.476 42.4756 42.4751 42.4745
+          42.4736 42.4726
+        </DataArray>
+        <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
+          1.227 1.22413 1.22149 1.21905 1.21678 1.21467 1.21271 1.21087 1.20916 1.20757 1.22983 1.22936
+          1.22878 1.22812 1.22739 1.22663 1.22583 1.22503 1.22421 1.22343 1.22972 1.22972 1.22967 1.22956
+          1.22942 1.22923 1.22901 1.22877 1.22849 1.22821 1.23008 1.23004 1.23001 1.22998 1.22994 1.22989
+          1.22983 1.22975 1.22967 1.22958 1.23011 1.23011 1.2301 1.23008 1.23007 1.23005 1.23003 1.23001
+          1.22998 1.22994
+        </DataArray>
+        <DataArray type="Float32" Name="v_x/v_x,max" NumberOfComponents="1" format="ascii">
+          0.990969 0.973296 0.956439 0.940778 0.926124 0.912446 0.899661 0.887712 0.876656 0.868418 0.999584 0.998119
+          0.995507 0.992304 0.988612 0.98463 0.980453 0.976158 0.971933 0.96861 1 0.999957 0.999637 0.999166
+          0.99846 0.9976 0.996604 0.995491 0.994336 0.993356 0.999961 1 1 0.999998 0.999893 0.999752
+          0.999576 0.999366 0.999141 0.998933 0.999907 0.999898 0.999931 1 1 1 1 1
+          1 1
+        </DataArray>
+        <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
+          0.00053604 0.000470416 0.000412942 0.000356245 0.000300092 0.000244574 0.000189766 0.000135716 8.23746e-05 2.86046e-05 0.000530008 0.000470289
+          0.000413008 0.000356306 0.000300144 0.000244621 0.000189815 0.000135772 8.24722e-05 2.96917e-05 0.000525566 0.000470018 0.000413027 0.00035635
+          0.000300193 0.000244673 0.000189869 0.000135835 8.26262e-05 3.09723e-05 0.000522663 0.0004698 0.000413057 0.000356402 0.000300248 0.00024473
+          0.000189926 0.000135896 8.27676e-05 3.19309e-05 0.000521235 0.000469677 0.000413073 0.000356432 0.000300281 0.000244764 0.00018996 0.000135934
+          8.28511e-05 3.24247e-05
+        </DataArray>
+        <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
+          1.44385e-05 1.44992e-05 1.45555e-05 1.46078e-05 1.46565e-05 1.47021e-05 1.47447e-05 1.47846e-05 1.48221e-05 1.48571e-05 1.4378e-05 1.43872e-05
+          1.43991e-05 1.44128e-05 1.4428e-05 1.44441e-05 1.44608e-05 1.44779e-05 1.44953e-05 1.45118e-05 1.43717e-05 1.43733e-05 1.43756e-05 1.43786e-05
+          1.43823e-05 1.43867e-05 1.43917e-05 1.43973e-05 1.44033e-05 1.44094e-05 1.43727e-05 1.43728e-05 1.4373e-05 1.43735e-05 1.43742e-05 1.43752e-05
+          1.43764e-05 1.43779e-05 1.43796e-05 1.43815e-05 1.43728e-05 1.43728e-05 1.43728e-05 1.43729e-05 1.4373e-05 1.43732e-05 1.43735e-05 1.43738e-05
+          1.43743e-05 1.43749e-05
+        </DataArray>
+        <DataArray type="Float32" Name="nu_t" NumberOfComponents="1" format="ascii">
+          7.89256e-05 7.82221e-05 7.7372e-05 7.64328e-05 7.54439e-05 7.44331e-05 7.34201e-05 7.24193e-05 7.14416e-05 7.04928e-05 0.000126096 0.000117829
+          0.000110689 0.000104601 9.95182e-05 9.53974e-05 9.2188e-05 8.9834e-05 8.82758e-05 8.75691e-05 0.000130936 0.000126833 0.000123141 0.000119767
+          0.000116652 0.000113763 0.000111082 0.000108603 0.000106328 0.00010437 0.000131202 0.000127576 0.000124514 0.000121874 0.000119557 0.000117495
+          0.000115638 0.000113951 0.000112411 0.000111021 0.000131217 0.000127628 0.000124632 0.000122087 0.000119892 0.000117973 0.000116276 0.00011476
+          0.000113395 0.000112167
+        </DataArray>
+        <DataArray type="Float32" Name="l_w" NumberOfComponents="1" format="ascii">
+          0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.15 0.15
+          0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.25 0.25 0.25 0.25
+          0.25 0.25 0.25 0.25 0.25 0.25 0.35 0.35 0.35 0.35 0.35 0.35
+          0.35 0.35 0.35 0.35 0.45 0.45 0.45 0.45 0.45 0.45 0.45 0.45
+          0.45 0.45
+        </DataArray>
+        <DataArray type="Float32" Name="y^+" NumberOfComponents="1" format="ascii">
+          10.7535 10.7466 10.7359 10.723 10.7086 10.6933 10.6774 10.6612 10.6452 10.6307 32.396 32.4908
+          32.5576 32.604 32.6347 32.6529 32.6611 32.6611 32.6556 32.6509 54.0171 54.2037 54.3513 54.4694
+          54.5638 54.6384 54.6963 54.7401 54.7736 54.805 75.6188 75.888 76.1053 76.2844 76.4326 76.5554
+          76.6566 76.7395 76.8092 76.8755 97.2235 97.5699 97.8509 98.084 98.2788 98.4419 98.5787 98.6929
+          98.7912 98.8854
+        </DataArray>
+        <DataArray type="Float32" Name="u^+" NumberOfComponents="1" format="ascii">
+          31.9824 31.438 30.9437 30.4986 30.0953 29.7281 29.3927 29.0856 28.805 28.5932 32.2604 32.2398
+          32.2077 32.169 32.1259 32.08 32.0323 31.9836 31.9356 31.8921 32.2738 32.2991 32.3413 32.3914
+          32.4459 32.5025 32.5599 32.617 32.6717 32.7069 32.2726 32.3005 32.3531 32.4184 32.4925 32.5727
+          32.657 32.744 32.8296 32.8905 32.2708 32.2972 32.3508 32.4185 32.496 32.5807 32.6709 32.7647
+          32.8578 32.9256
+        </DataArray>
+        <DataArray type="Float32" Name="k" NumberOfComponents="1" format="ascii">
+          4.94091e-05 4.85322e-05 4.74831e-05 4.63372e-05 4.5146e-05 4.39444e-05 4.27564e-05 4.15987e-05 4.04831e-05 3.9415e-05 4.27605e-05 3.73789e-05
+          3.33351e-05 3.02767e-05 2.79635e-05 2.6228e-05 2.49508e-05 2.40457e-05 2.34453e-05 2.30279e-05 4.2569e-05 3.68613e-05 3.23891e-05 2.88229e-05
+          2.59407e-05 2.35887e-05 2.16574e-05 2.00672e-05 1.87565e-05 1.76601e-05 4.25671e-05 3.68498e-05 3.23561e-05 2.87511e-05 2.58088e-05 2.33726e-05
+          2.13309e-05 1.96027e-05 1.81276e-05 1.6859e-05 4.25675e-05 3.68507e-05 3.23575e-05 2.87522e-05 2.58076e-05 2.33658e-05 2.13139e-05 1.957e-05
+          1.80729e-05 1.6778e-05
+        </DataArray>
+        <DataArray type="Float32" Name="epsilon" NumberOfComponents="1" format="ascii">
+          2.7838e-06 2.71002e-06 2.62262e-06 2.52827e-06 2.43141e-06 2.33498e-06 2.24093e-06 2.15054e-06 2.06461e-06 1.98345e-06 1.30505e-06 1.0672e-06
+          9.03532e-07 7.88723e-07 7.07169e-07 6.48985e-07 6.07765e-07 5.79265e-07 5.6042e-07 5.45006e-07 1.24558e-06 9.64165e-07 7.6672e-07 6.24283e-07
+          5.19177e-07 4.40201e-07 3.80026e-07 3.33715e-07 2.97782e-07 2.68938e-07 1.24294e-06 9.57956e-07 7.56724e-07 6.10436e-07 5.01421e-07 4.18444e-07
+          3.54128e-07 3.03498e-07 2.63097e-07 2.3041e-07 1.24282e-06 9.57609e-07 7.56071e-07 6.09416e-07 4.99974e-07 4.16503e-07 3.51624e-07 3.00354e-07
+          2.59242e-07 2.2587e-07
+        </DataArray>
+        <DataArray type="Float32" Name="y^+_nom" NumberOfComponents="1" format="ascii">
+          13.3325 13.1583 12.965 12.7618 12.5548 12.3482 12.1449 11.947 11.756 11.5725 40.1657 39.7824
+          39.3176 38.8032 38.2609 37.7063 37.1501 36.6003 36.063 35.5436 66.9722 66.3682 65.6364 64.826
+          63.9706 63.0942 62.2139 61.3423 60.4889 59.6604 93.7548 92.9188 91.9074 90.7887 89.6097 88.4031
+          87.1924 85.995 84.8238 83.6862 120.541 119.467 118.168 116.733 115.222 113.677 112.128 110.596
+          109.099 107.646
+        </DataArray>
+        <DataArray type="Float32" Name="u^+_nom" NumberOfComponents="1" format="ascii">
+          25.7957 25.6758 25.6234 25.6261 25.6698 25.744 25.8411 25.9552 26.0834 26.2662 26.02 26.3306
+          26.6701 27.0297 27.4018 27.7806 28.1617 28.5412 28.9182 29.2966 26.0308 26.3791 26.7807 27.2166
+          27.6748 28.1466 28.6256 29.1065 29.5848 30.045 26.0298 26.3802 26.7904 27.2393 27.7145 28.2073
+          28.7109 29.2198 29.7277 30.2137 26.0283 26.3776 26.7886 27.2393 27.7175 28.2143 28.7231 29.2383
+          29.7533 30.246
+        </DataArray>
+        <DataArray type="Float32" Name="inNearWallRegion" 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
+        </DataArray>
+        <DataArray type="Float32" Name="isMatchingPoint" NumberOfComponents="1" format="ascii">
+          1 1 1 1 1 1 1 1 1 1 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="temperature" NumberOfComponents="1" format="ascii">
+          283.869 284.533 285.147 285.718 286.248 286.744 287.207 287.641 288.047 288.426 283.209 283.31
+          283.441 283.592 283.758 283.935 284.118 284.305 284.494 284.675 283.156 283.17 283.192 283.224
+          283.264 283.311 283.365 283.425 283.491 283.557 283.151 283.152 283.156 283.161 283.17 283.18
+          283.194 283.21 283.229 283.25 283.15 283.15 283.151 283.152 283.153 283.156 283.159 283.164
+          283.169 283.176
+        </DataArray>
+        <DataArray type="Float32" Name="lambda" NumberOfComponents="1" format="ascii">
+          0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535
+          0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535
+          0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535
+          0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535 0.0255535
+          0.0255535 0.0255535
+        </DataArray>
+        <DataArray type="Float32" Name="lambda_t" NumberOfComponents="1" format="ascii">
+          0.114414 0.113133 0.111665 0.110093 0.10847 0.106835 0.105213 0.103625 0.102084 0.100598 0.183221 0.171153
+          0.160711 0.151793 0.144335 0.138274 0.133537 0.130043 0.127703 0.126602 0.190361 0.184374 0.178984 0.174054
+          0.169498 0.165269 0.161342 0.157706 0.154366 0.151487 0.190681 0.185414 0.180965 0.177127 0.173757 0.170754
+          0.168047 0.165587 0.163337 0.161307 0.190695 0.18548 0.181128 0.177431 0.174241 0.171452 0.168985 0.16678
+          0.164794 0.163006
+        </DataArray>
+        <DataArray type="Float32" Name="X^H2O_gas" NumberOfComponents="1" format="ascii">
+          4.4882e-06 1.13316e-05 1.9234e-05 2.74364e-05 3.55055e-05 4.32067e-05 5.04251e-05 5.71163e-05 6.32768e-05 6.89298e-05 5.06258e-05 8.25272e-05
+          0.000102577 0.000115096 0.000122791 0.000127381 0.000129968 0.000131265 0.000131748 0.000131886 0.000504155 0.000419579 0.000358106 0.000312499
+          0.000277991 0.000251396 0.000230544 0.000213938 0.000200533 0.000189736 5.41871e-05 8.86178e-05 0.000110188 0.000123524 0.000131569 0.000136212
+          0.000138673 0.000139742 0.000139929 0.000139498 6.23616e-06 1.55258e-05 2.59594e-05 3.64931e-05 4.65666e-05 5.59013e-05 6.43821e-05 7.19872e-05
+          7.87394e-05 8.46231e-05
+        </DataArray>
+        <DataArray type="Float32" Name="x^H2O_gas" NumberOfComponents="1" format="ascii">
+          7.21491e-06 1.82158e-05 3.09189e-05 4.41042e-05 5.7075e-05 6.94544e-05 8.10575e-05 9.18131e-05 0.000101716 0.000110802 8.13802e-05 0.000132659
+          0.000164885 0.000185007 0.000197377 0.000204754 0.000208911 0.000210997 0.000211772 0.000211994 0.000810198 0.000674315 0.000575543 0.000502257
+          0.000446804 0.000404065 0.000370556 0.000343867 0.000322324 0.000304971 8.71048e-05 0.000142448 0.000177118 0.000198553 0.000211485 0.000218948
+          0.000222902 0.000224622 0.000224922 0.000224229 1.00248e-05 2.49581e-05 4.173e-05 5.86626e-05 7.48552e-05 8.98601e-05 0.000103492 0.000115717
+          0.00012657 0.000136027
+        </DataArray>
+        <DataArray type="Float32" Name="D^H2O_gas" NumberOfComponents="1" format="ascii">
+          2.28281e-05 2.29243e-05 2.30135e-05 2.30964e-05 2.31737e-05 2.32459e-05 2.33136e-05 2.3377e-05 2.34364e-05 2.34919e-05 2.27326e-05 2.27473e-05
+          2.27661e-05 2.27881e-05 2.28121e-05 2.28376e-05 2.28641e-05 2.28912e-05 2.29187e-05 2.29449e-05 2.27249e-05 2.27269e-05 2.27303e-05 2.27348e-05
+          2.27406e-05 2.27474e-05 2.27552e-05 2.27639e-05 2.27734e-05 2.27829e-05 2.27242e-05 2.27245e-05 2.2725e-05 2.27258e-05 2.27269e-05 2.27285e-05
+          2.27305e-05 2.27328e-05 2.27356e-05 2.27386e-05 2.27241e-05 2.27242e-05 2.27242e-05 2.27244e-05 2.27246e-05 2.2725e-05 2.27254e-05 2.27261e-05
+          2.27269e-05 2.27279e-05
+        </DataArray>
+        <DataArray type="Float32" Name="X^Air_gas" NumberOfComponents="1" format="ascii">
+          0.999996 0.999989 0.999981 0.999973 0.999964 0.999957 0.99995 0.999943 0.999937 0.999931 0.999949 0.999917
+          0.999897 0.999885 0.999877 0.999873 0.99987 0.999869 0.999868 0.999868 0.999496 0.99958 0.999642 0.999687
+          0.999722 0.999749 0.999769 0.999786 0.999799 0.99981 0.999946 0.999911 0.99989 0.999876 0.999868 0.999864
+          0.999861 0.99986 0.99986 0.999861 0.999994 0.999985 0.999974 0.999964 0.999953 0.999944 0.999936 0.999928
+          0.999921 0.999915
+        </DataArray>
+        <DataArray type="Float32" Name="x^Air_gas" NumberOfComponents="1" format="ascii">
+          0.999993 0.999982 0.999969 0.999956 0.999943 0.999931 0.999919 0.999908 0.999898 0.999889 0.999919 0.999867
+          0.999835 0.999815 0.999803 0.999795 0.999791 0.999789 0.999788 0.999788 0.99919 0.999326 0.999424 0.999498
+          0.999553 0.999596 0.999629 0.999656 0.999678 0.999695 0.999913 0.999858 0.999823 0.999801 0.999789 0.999781
+          0.999777 0.999775 0.999775 0.999776 0.99999 0.999975 0.999958 0.999941 0.999925 0.99991 0.999897 0.999884
+          0.999873 0.999864
+        </DataArray>
+        <DataArray type="Float32" Name="D_t" NumberOfComponents="1" format="ascii">
+          0.000112751 0.000111746 0.000110531 0.00010919 0.000107777 0.000106333 0.000104886 0.000103456 0.000102059 0.000100704 0.000180137 0.000168326
+          0.000158127 0.00014943 0.000142169 0.000136282 0.000131697 0.000128334 0.000126108 0.000125099 0.000187052 0.00018119 0.000175916 0.000171096
+          0.000166645 0.000162518 0.000158688 0.000155147 0.000151898 0.0001491 0.000187432 0.000182251 0.000177877 0.000174106 0.000170796 0.00016785
+          0.000165197 0.000162788 0.000160587 0.000158602 0.000187452 0.000182325 0.000178046 0.00017441 0.000171275 0.000168533 0.000166109 0.000163943
+          0.000161993 0.000160239
+        </DataArray>
+        <DataArray type="Float32" Name="dv_x/dx_" NumberOfComponents="3" format="ascii">
+          -0.00225394 0.667851 0 -0.00217095 0.669801 0 -0.00202187 0.671063 0 -0.0018661 0.671855 0
+          -0.00172642 0.672293 0 -0.00160067 0.672455 0 -0.00148574 0.672401 0 -0.00137459 0.672174 0
+          -0.0011415 0.671857 0 -0.000964413 0.67161 0 0.000522323 0.00452515 0 0.00040148 0.0134182 0
+          0.000256776 0.0218391 0 0.000153755 0.0296495 0 7.50583e-05 0.0368933 0 1.34847e-05 0.0436174 0
+          -3.51535e-05 0.0498627 0 -6.79046e-05 0.0556619 0 -7.05355e-05 0.0610108 0 -6.18771e-05 0.0649717 0
+          0.000727883 0.000188815 0 0.000714713 0.000946757 0 0.000683355 0.00227145 0 0.00064292 0.0039072 0
+          0.000601554 0.00575383 0 0.0005609 0.00774579 0 0.000521684 0.00983577 0 0.000483428 0.0119854 0
+          0.000410134 0.0141058 0 0.000355904 0.0157689 0 0.000732714 -4.67119e-05 0 0.000748586 -2.95359e-05 0
+          0.000750182 0.000149062 0 0.000734153 0.000423564 0 0.000714823 0.000785621 0 0.000693808 0.00122919 0
+          0.000671449 0.00174691 0 0.000643769 0.00232855 0 0.000560045 0.00293621 0 0.000492479 0.00345515 0
+          0.000721319 -5.44696e-05 0 0.000747342 -0.00010255 0 0.000758871 -6.9395e-05 0 0.000749013 1.72196e-06 0
+          0.000735825 0.000108935 0 0.000721115 0.000253745 0 0.000704911 0.000436624 0 0.000681644 0.000655289 0
+          0.0005979 0.000891114 0 0.000528885 0.00110955 0
+        </DataArray>
+        <DataArray type="Float32" Name="dv_y/dx_" NumberOfComponents="3" format="ascii">
+          0.000143034 0.00161815 0 -1.42142e-05 0.00160501 0 -1.63532e-05 0.00151698 0 -1.46173e-05 0.00142993 0
+          -1.30493e-05 0.00134807 0 -1.17878e-05 0.00127119 0 -1.0827e-05 0.00119864 0 -1.08167e-05 0.0011286 0
+          -3.08974e-05 0.00104199 0 -5.05871e-05 0.000623728 0 0.000267502 0.000213415 0 -1.26457e-05 0.000297876 0
+          -2.18851e-05 0.00033137 0 -2.11142e-05 0.000351569 0 -1.98432e-05 0.000364008 0 -1.86789e-05 0.000370969 0
+          -1.7823e-05 0.000373888 0 -1.95805e-05 0.00037283 0 -6.31092e-05 0.000354208 0 -0.000104566 0.000227268 0
+          0.000209229 -0.000684334 0 5.44495e-06 -0.000631257 0 -7.40444e-06 -0.000570022 0 -9.17765e-06 -0.000515007 0
+          -9.81598e-06 -0.000465721 0 -1.01411e-05 -0.000421332 0 -1.05167e-05 -0.000381064 0 -1.40967e-05 -0.000344354 0
+          -5.51578e-05 -0.000310641 0 -9.29006e-05 -0.000172744 0 0.000127223 -0.00071575 0 6.40631e-06 -0.000750451 0
+          -2.50999e-06 -0.000733241 0 -3.73076e-06 -0.000709849 0 -4.23069e-06 -0.000684475 0 -4.5695e-06 -0.000658205 0
+          -4.99325e-06 -0.000631546 0 -7.84335e-06 -0.000603939 0 -3.45076e-05 -0.000562029 0 -5.86001e-05 -0.000342362 0
+          4.24618e-05 -0.000708686 0 2.52984e-06 -0.000762847 0 -6.37462e-07 -0.000755203 0 -1.05003e-06 -0.000740376 0
+          -1.20856e-06 -0.000723035 0 -1.31952e-06 -0.000704111 0 -1.47248e-06 -0.000684035 0 -2.51047e-06 -0.000661862 0
+          -1.15616e-05 -0.00062004 0 -1.96784e-05 -0.00038651 0
+        </DataArray>
+        <DataArray type="Float32" Name="velocity_gas (m/s)" NumberOfComponents="3" format="ascii">
+          0.0993142 0.000135454 0 0.0979715 0.000128731 0 0.0967091 0.000118397 0 0.0955452 0.000109107 0
+          0.0944698 0.000100856 0 0.0934735 9.34475e-05 0 0.092549 8.67105e-05 0 0.0916906 8.04549e-05 0
+          0.0908995 7.373e-05 0 0.0903208 4.33776e-05 0 0.100178 0.000242722 0 0.10047 0.000240752 0
+          0.100659 0.000227547 0 0.100778 0.00021449 0 0.100844 0.00020221 0 0.100868 0.000190678 0
+          0.10086 0.000179795 0 0.100826 0.00016929 0 0.100779 0.000156298 0 0.100741 9.35581e-05 0
+          0.100219 0.000178137 0 0.100655 0.000188306 0 0.101077 0.000184671 0 0.101475 0.00017942 0
+          0.101848 0.000173658 0 0.102197 0.000167641 0 0.102522 0.000161488 0 0.102823 0.000155021 0
+          0.103102 0.000144572 0 0.103315 8.88313e-05 0 0.100215 0.000105855 0 0.100659 0.0001145 0
+          0.101114 0.000113543 0 0.10156 0.000111488 0 0.101995 0.000109066 0 0.102417 0.000106411 0
+          0.102827 0.000103582 0 0.103223 0.000100419 0 0.1036 9.41702e-05 0 0.103895 5.90102e-05 0
+          0.10021 3.49867e-05 0 0.100649 3.82156e-05 0 0.101107 3.80225e-05 0 0.10156 3.74506e-05 0
+          0.102006 3.67625e-05 0 0.102443 3.60003e-05 0 0.102871 3.5179e-05 0 0.103289 3.42333e-05 0
+          0.103689 3.21664e-05 0 0.104006 2.03594e-05 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
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+          0 0 0 0.6 0 0 0 0.1 0 0.6 0.1 0
+          1.2 0 0 1.2 0.1 0 1.8 0 0 1.8 0.1 0
+          2.4 0 0 2.4 0.1 0 3 0 0 3 0.1 0
+          3.6 0 0 3.6 0.1 0 4.2 0 0 4.2 0.1 0
+          4.8 0 0 4.8 0.1 0 5.4 0 0 5.4 0.1 0
+          6 0 0 6 0.1 0 0 0.2 0 0.6 0.2 0
+          1.2 0.2 0 1.8 0.2 0 2.4 0.2 0 3 0.2 0
+          3.6 0.2 0 4.2 0.2 0 4.8 0.2 0 5.4 0.2 0
+          6 0.2 0 0 0.3 0 0.6 0.3 0 1.2 0.3 0
+          1.8 0.3 0 2.4 0.3 0 3 0.3 0 3.6 0.3 0
+          4.2 0.3 0 4.8 0.3 0 5.4 0.3 0 6 0.3 0
+          0 0.4 0 0.6 0.4 0 1.2 0.4 0 1.8 0.4 0
+          2.4 0.4 0 3 0.4 0 3.6 0.4 0 4.2 0.4 0
+          4.8 0.4 0 5.4 0.4 0 6 0.4 0 0 0.5 0
+          0.6 0.5 0 1.2 0.5 0 1.8 0.5 0 2.4 0.5 0
+          3 0.5 0 3.6 0.5 0 4.2 0.5 0 4.8 0.5 0
+          5.4 0.5 0 6 0.5 0
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+          0 1 3 2 1 4 5 3 4 6 7 5
+          6 8 9 7 8 10 11 9 10 12 13 11
+          12 14 15 13 14 16 17 15 16 18 19 17
+          18 20 21 19 2 3 23 22 3 5 24 23
+          5 7 25 24 7 9 26 25 9 11 27 26
+          11 13 28 27 13 15 29 28 15 17 30 29
+          17 19 31 30 19 21 32 31 22 23 34 33
+          23 24 35 34 24 25 36 35 25 26 37 36
+          26 27 38 37 27 28 39 38 28 29 40 39
+          29 30 41 40 30 31 42 41 31 32 43 42
+          33 34 45 44 34 35 46 45 35 36 47 46
+          36 37 48 47 37 38 49 48 38 39 50 49
+          39 40 51 50 40 41 52 51 41 42 53 52
+          42 43 54 53 44 45 56 55 45 46 57 56
+          46 47 58 57 47 48 59 58 48 49 60 59
+          49 50 61 60 50 51 62 61 51 52 63 62
+          52 53 64 63 53 54 65 64
+        </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
+        </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
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>
-- 
GitLab