From 9e8bc2c39956b9bef7441585b302a852fd79564e Mon Sep 17 00:00:00 2001
From: Martin Schneider <martin.schneider@iws.uni-stuttgart.de>
Date: Thu, 12 Dec 2019 12:00:01 +0100
Subject: [PATCH] [box][tests] Add box forchheimer test
---
.../1p/isothermal/CMakeLists.txt | 11 ++++
.../test_1p_forchheimer_box-reference.vtu | 63 +++++++++++++++++++
2 files changed, 74 insertions(+)
create mode 100644 test/references/test_1p_forchheimer_box-reference.vtu
diff --git a/test/porousmediumflow/1p/isothermal/CMakeLists.txt b/test/porousmediumflow/1p/isothermal/CMakeLists.txt
index 090267217b..c5c3982128 100644
--- a/test/porousmediumflow/1p/isothermal/CMakeLists.txt
+++ b/test/porousmediumflow/1p/isothermal/CMakeLists.txt
@@ -44,6 +44,17 @@ dumux_add_test(NAME test_1p_forchheimer_tpfa
--command "${CMAKE_CURRENT_BINARY_DIR}/test_1p_forchheimer_tpfa params_forchheimer.input -Problem.Name test_1p_forchheimer_tpfa"
--zeroThreshold {"velocity_liq \(m/s\)":1e-12})
+dumux_add_test(NAME test_1p_forchheimer_box
+ LABELS porousmediumflow 1p
+ SOURCES main.cc
+ COMPILE_DEFINITIONS TYPETAG=OnePTestBox FORCHHEIMER=1
+ COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ CMD_ARGS --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/test_1p_forchheimer_box-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/test_1p_forchheimer_box-00001.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p_forchheimer_box params_forchheimer.input -Problem.Name test_1p_forchheimer_box"
+ --zeroThreshold {"velocity_liq \(m/s\)":1e-12})
+
# a gstat test (becaue it's a random permeability field we can't test against a reference solution)
dumux_add_test(NAME test_1p_gstat
LABELS porousmediumflow 1p
diff --git a/test/references/test_1p_forchheimer_box-reference.vtu b/test/references/test_1p_forchheimer_box-reference.vtu
new file mode 100644
index 0000000000..cd8c417b90
--- /dev/null
+++ b/test/references/test_1p_forchheimer_box-reference.vtu
@@ -0,0 +1,63 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="20" NumberOfPoints="33">
+ <PointData Scalars="p" Vectors="velocity_liq (m/s)">
+ <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
+ 200000 200000 190000 190000 200000 190000 180000 180000 180000 170000 170000 170000
+ 160000 160000 160000 150000 150000 150000 140000 140000 140000 130000 130000 130000
+ 120000 120000 120000 110000 110000 110000 100000 100000 100000
+ </DataArray>
+ <DataArray type="Float32" Name="velocity_liq (m/s)" NumberOfComponents="3" format="ascii">
+ 0 0.236926 0 0 0.236926 0 2.54118e-17 0.236926 0 2.54118e-17 0.236926 0
+ 0 0.236926 0 2.54118e-17 0.236926 0 2.57874e-17 0.236926 0 3.44067e-17 0.236926 0
+ 4.3026e-17 0.236926 0 -9.49765e-17 0.236926 0 -6.04993e-17 0.236926 0 -2.6022e-17 0.236926 0
+ 1.31051e-16 0.236926 0 1.03806e-17 0.236926 0 -1.1029e-16 0.236926 0 1.18368e-16 0.236926 0
+ -1.95407e-17 0.236926 0 -1.5745e-16 0.236926 0 -7.91435e-18 0.236926 0 -1.65337e-17 0.236926 0
+ -2.5153e-17 0.236926 0 -5.39001e-17 0.236926 0 -6.49394e-18 0.236926 0 4.09122e-17 0.236926 0
+ -2.07381e-17 0.236926 0 1.80488e-17 0.236926 0 5.68356e-17 0.236926 0 2.10432e-17 0.236926 0
+ 2.10432e-17 0.236926 0 2.10432e-17 0.236926 0 1.99153e-17 0.236926 0 1.99153e-17 0.236926 0
+ 1.99153e-17 0.236926 0
+ </DataArray>
+ </PointData>
+ <CellData Scalars="process rank">
+ <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
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 0.5 0 0 0 0.1 0 0.5 0.1 0
+ 1 0 0 1 0.1 0 0 0.2 0 0.5 0.2 0
+ 1 0.2 0 0 0.3 0 0.5 0.3 0 1 0.3 0
+ 0 0.4 0 0.5 0.4 0 1 0.4 0 0 0.5 0
+ 0.5 0.5 0 1 0.5 0 0 0.6 0 0.5 0.6 0
+ 1 0.6 0 0 0.7 0 0.5 0.7 0 1 0.7 0
+ 0 0.8 0 0.5 0.8 0 1 0.8 0 0 0.9 0
+ 0.5 0.9 0 1 0.9 0 0 1 0 0.5 1 0
+ 1 1 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 2 3 7 6
+ 3 5 8 7 6 7 10 9 7 8 11 10
+ 9 10 13 12 10 11 14 13 12 13 16 15
+ 13 14 17 16 15 16 19 18 16 17 20 19
+ 18 19 22 21 19 20 23 22 21 22 25 24
+ 22 23 26 25 24 25 28 27 25 26 29 28
+ 27 28 31 30 28 29 32 31
+ </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
+ </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
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
--
GitLab