diff --git a/src/paper_thinstrip_mixed/chns.ini b/src/paper_thinstrip_mixed/chns.ini index 7723e3c1b3c675901a00b3e7c303ffca63b912db..18c2e381670816c3ea30006b552e004a536826d3 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 5acf69a20caffa1d4869ab7994fbdc39bd624a80..7723e3c1b3c675901a00b3e7c303ffca63b912db 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 bc7c4efb3ebed18440163473194b2b2bb0976765..340ffb3ccd346b493e1e0e51892e263eb3ca0940 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 a7620a26fa635df401d22467e2f813900d210e40..3d611bd33291a1399662d987849425e3135ba52f 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