Other testcase

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

src/paper_thinstrip_mixed/chns_reactive.ini → 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

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));

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