From c1aa0ef443ea9eb6b5b78a4f61fdff665e9ab320 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Fri, 5 Feb 2021 20:27:28 +0100
Subject: [PATCH] Updating multi drop meshes to better accuracy.
---
Meshes/multi_drop_gauss.geo | 10 +++++-----
Meshes/multi_drop_gauss_far.geo | 10 +++++-----
Meshes/multi_drop_gauss_med.geo | 10 +++++-----
3 files changed, 15 insertions(+), 15 deletions(-)
diff --git a/Meshes/multi_drop_gauss.geo b/Meshes/multi_drop_gauss.geo
index 80ddb2a..a0c98d9 100644
--- a/Meshes/multi_drop_gauss.geo
+++ b/Meshes/multi_drop_gauss.geo
@@ -5,19 +5,19 @@ Physical Surface("surfacedomain") = {1};
Physical Volume("volumedomain") = {1};
Field[3] = MathEval;
-Field[3].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[3].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[4] = MathEval;
-Field[4].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[4].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[5] = MathEval;
-Field[5].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[5].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[6] = MathEval;
-Field[6].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[6].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.25)*(z-0.25))-0.25)))";
Field[7] = MathEval;
-Field[7].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[7].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.25)*(z-0.25))-0.25)))";
Field[8] = Min;
Field[8].FieldsList = {3,4,5,6,7};
diff --git a/Meshes/multi_drop_gauss_far.geo b/Meshes/multi_drop_gauss_far.geo
index e9a25ee..d9e4f86 100644
--- a/Meshes/multi_drop_gauss_far.geo
+++ b/Meshes/multi_drop_gauss_far.geo
@@ -5,19 +5,19 @@ Physical Surface("surfacedomain") = {1};
Physical Volume("volumedomain") = {1};
Field[3] = MathEval;
-Field[3].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[3].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[4] = MathEval;
-Field[4].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[4].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[5] = MathEval;
-Field[5].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[5].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[6] = MathEval;
-Field[6].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[6].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.25)*(z-0.25))-0.25)))";
Field[7] = MathEval;
-Field[7].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[7].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.25)*(z-0.25))-0.25)))";
Field[8] = Min;
Field[8].FieldsList = {3,4,5,6,7};
diff --git a/Meshes/multi_drop_gauss_med.geo b/Meshes/multi_drop_gauss_med.geo
index f189c78..5a28df7 100644
--- a/Meshes/multi_drop_gauss_med.geo
+++ b/Meshes/multi_drop_gauss_med.geo
@@ -5,19 +5,19 @@ Physical Surface("surfacedomain") = {1};
Physical Volume("volumedomain") = {1};
Field[3] = MathEval;
-Field[3].F = ".5 - (.5-0.01)*exp(-(1./0.015) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[3].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[4] = MathEval;
-Field[4].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[4].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[5] = MathEval;
-Field[5].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[5].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
Field[6] = MathEval;
-Field[6].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[6].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.25)*(z-0.25))-0.25)))";
Field[7] = MathEval;
-Field[7].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.5)*(z-0.5))-0.25)))";
+Field[7].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.25)*(z-0.25))-0.25)))";
Field[8] = Min;
Field[8].FieldsList = {3,4,5,6,7};
--
GitLab