From 60a16f1b5dc279bd86973fb0b05118b4b92a8870 Mon Sep 17 00:00:00 2001
From: Lars von Wolff <lars.von-wolff@ians.uni-stuttgart.de>
Date: Thu, 27 May 2021 10:05:05 +0200
Subject: [PATCH] Other testcase
---
src/paper_thinstrip_mixed/chns.ini | 8 ++++----
.../{chns_reactive.ini => chns_hyperbolic.ini} | 16 ++++++++--------
src/paper_thinstrip_mixed/tmixed_boundary.hh | 2 +-
src/paper_thinstrip_mixed/tmixed_upscaledc.hh | 2 +-
4 files changed, 14 insertions(+), 14 deletions(-)
rename src/paper_thinstrip_mixed/{chns_reactive.ini => chns_hyperbolic.ini} (85%)
diff --git a/src/paper_thinstrip_mixed/chns.ini b/src/paper_thinstrip_mixed/chns.ini
index 7723e3c..18c2e38 100644
--- a/src/paper_thinstrip_mixed/chns.ini
+++ b/src/paper_thinstrip_mixed/chns.ini
@@ -4,7 +4,7 @@ D = 0.02
rho1 = 1.0
rho2 = 1.0
mu=0.005
-ReactionRate=0.5
+ReactionRate=2#0.5
#ReactionLimiter = 1000
SurfaceTensionScaling=0.1
VelDissipationRate = 1000
@@ -20,7 +20,7 @@ YResolution = 200
XResolution = 400
[output]
-filename = comparison_hyperbolic_bugfix
+filename = comparison_diffusion
[Phasefield]
eps=0.03
@@ -32,7 +32,7 @@ Sigma2=1
Sigma3=1
[Time]
-tMax = 1
+tMax = 4
dt = 0.0001
dtMax = 0.0001
dtMin = 0.000015
@@ -45,7 +45,7 @@ saveintervall = 0.005
pConst=0.05
uConst=0.5
uInflow=0.3
-ionSource=0.0
+ionSource=1.0
[Refinement]
etaRefine = 0.2
diff --git a/src/paper_thinstrip_mixed/chns_reactive.ini b/src/paper_thinstrip_mixed/chns_hyperbolic.ini
similarity index 85%
rename from src/paper_thinstrip_mixed/chns_reactive.ini
rename to src/paper_thinstrip_mixed/chns_hyperbolic.ini
index 5acf69a..7723e3c 100644
--- a/src/paper_thinstrip_mixed/chns_reactive.ini
+++ b/src/paper_thinstrip_mixed/chns_hyperbolic.ini
@@ -4,7 +4,7 @@ D = 0.02
rho1 = 1.0
rho2 = 1.0
mu=0.005
-ReactionRate=2#0.5
+ReactionRate=0.5
#ReactionLimiter = 1000
SurfaceTensionScaling=0.1
VelDissipationRate = 1000
@@ -16,11 +16,11 @@ level = 0
filename = grids/rectangle_periodic_coarse.msh
#YScalingFactor = 0.5
YScaling = 1
-YResolution = 100#200
-XResolution = 100#400
+YResolution = 200
+XResolution = 400
[output]
-filename = comparison_diffusion
+filename = comparison_hyperbolic_bugfix
[Phasefield]
eps=0.03
@@ -32,9 +32,9 @@ Sigma2=1
Sigma3=1
[Time]
-tMax = 4
-dt = 0.002#0.0001
-dtMax = 0.002#0.0001
+tMax = 1
+dt = 0.0001
+dtMax = 0.0001
dtMin = 0.000015
saveintervall = 0.005
#dtIncreaseFactor=1.2
@@ -45,7 +45,7 @@ saveintervall = 0.005
pConst=0.05
uConst=0.5
uInflow=0.3
-ionSource=1.0
+ionSource=0.0
[Refinement]
etaRefine = 0.2
diff --git a/src/paper_thinstrip_mixed/tmixed_boundary.hh b/src/paper_thinstrip_mixed/tmixed_boundary.hh
index bc7c4ef..340ffb3 100644
--- a/src/paper_thinstrip_mixed/tmixed_boundary.hh
+++ b/src/paper_thinstrip_mixed/tmixed_boundary.hh
@@ -176,7 +176,7 @@ public:
//bool inMiddle = (pref.TempXpos > 0.4) && (pref.TempXpos < 0.6);
//RF mod = 0.1*std::sin(2*3.141592*pref.TempXpos);//inMiddle?0.3:0;
- RF d_sine = /*0;*/0.15*sin(pref.TempXpos*2*M_PI);
+ RF d_sine = 0;//0.15*sin(pref.TempXpos*2*M_PI);
//std::cout << "TempXPos" << pref.TempXpos;
RF rf = pref.tw.dw.eval_shape(1/eps*(x[0] - dSolid));
diff --git a/src/paper_thinstrip_mixed/tmixed_upscaledc.hh b/src/paper_thinstrip_mixed/tmixed_upscaledc.hh
index a7620a2..3d611bd 100644
--- a/src/paper_thinstrip_mixed/tmixed_upscaledc.hh
+++ b/src/paper_thinstrip_mixed/tmixed_upscaledc.hh
@@ -81,7 +81,7 @@ template<typename Param, typename PAdaptor, typename CGFS, typename CContainer,
RF source = 0;
if(xg > 0.1 && xg < 0.3)
{
- source = (xg-0.1)*(0.3-xg)/0.04;
+ source = (xg-0.1)*(0.3-xg)/0.016;
}
r.accumulate(lfsv,0, eg.geometry().volume()*
((phic*c - phicOld*cOld + param.phys.uAst*(phiSolid-phiSolidOld))/param.time.dt
--
GitLab