From cd6274807b0d744b7657c201bf06d4c4bbb00916 Mon Sep 17 00:00:00 2001
From: Sven Marcus <s.marcus@outlook.de>
Date: Thu, 11 Feb 2021 15:25:35 +0100
Subject: [PATCH] Finish Poiseuille Test
---
Python/poiseuille/test_poiseuille_l2.py | 14 +++++++++-----
1 file changed, 9 insertions(+), 5 deletions(-)
diff --git a/Python/poiseuille/test_poiseuille_l2.py b/Python/poiseuille/test_poiseuille_l2.py
index a3a731853..e3dd226ff 100644
--- a/Python/poiseuille/test_poiseuille_l2.py
+++ b/Python/poiseuille/test_poiseuille_l2.py
@@ -35,9 +35,14 @@ class TestPoiseuilleFlow(unittest.TestCase):
self.viscosities)
def test_poiseuille_flow(self):
- self.skipTest("This test is not implemented correctly yet")
plt.ion()
+ channel_height = 10
+ number_of_nodes = [16, 32, 64]
+ number_of_timesteps = [40_000, 80_000, 160_000]
+ viscosities = [5e-3, 1e-2, 2e-2]
+ l2_norm_results = []
+
physical_params = PhysicalParameters()
runtime_params = RuntimeParameters()
@@ -53,9 +58,8 @@ class TestPoiseuilleFlow(unittest.TestCase):
runtime_params.number_of_timesteps = timesteps
kernel.forcing_in_x1 = forcing
- grid_params = create_grid_params_with_nodes_in_column(nodes, delta_x)
- l2_error = get_l2_error_for_simulation(grid_params, physical_params, runtime_params, kernel)
- normalized_l2_errors.append(l2_error)
+ self.assertLessEqual(l2_norm_results[1], l2_norm_results[0] / 2)
+ self.assertLessEqual(l2_norm_results[2], l2_norm_results[1] / 2)
nodes_as_log = [np.log10(node) for node in self.number_of_nodes]
l2_norms_as_log = [np.log10(l2) for l2 in normalized_l2_errors]
@@ -77,7 +81,7 @@ def get_l2_error_for_simulation(grid_params, physical_params, runtime_params, ke
heights = get_heights(output_folder, runtime_params)
numerical_results = get_numerical_results(runtime_params, output_folder)
- analytical_results = get_analytical_results(grid_params, physical_params, kernel, heights)
+ analytical_results = get_analytical_results(physical_params, heights, grid_params.number_of_nodes_per_direction[2])
plt.plot(heights, numerical_results)
plt.plot(heights, analytical_results)
--
GitLab