Skip to content
Snippets Groups Projects
Commit e7b0f927 authored by Lars von Wolff's avatar Lars von Wolff
Browse files

State of 1st example for 2f1s_thinstrip

parent 13875593
No related branches found
No related tags found
No related merge requests found
...@@ -17,7 +17,7 @@ filename = grids/rectangle_periodic_coarse.msh ...@@ -17,7 +17,7 @@ filename = grids/rectangle_periodic_coarse.msh
#YScalingFactor = 0.5 #YScalingFactor = 0.5
YScaling = 1 YScaling = 1
YResolution = 200 YResolution = 200
XResolution = 100 XResolution = 400
[output] [output]
filename = comparison_hyperbolic filename = comparison_hyperbolic
...@@ -33,8 +33,8 @@ Sigma3=1 ...@@ -33,8 +33,8 @@ Sigma3=1
[Time] [Time]
tMax = 1 tMax = 1
dt = 0.001 dt = 0.0001
dtMax = 0.001 dtMax = 0.0001
dtMin = 0.000015 dtMin = 0.000015
saveintervall = 0.005 saveintervall = 0.005
#dtIncreaseFactor=1.2 #dtIncreaseFactor=1.2
......
...@@ -121,8 +121,8 @@ public: ...@@ -121,8 +121,8 @@ public:
nsGo = std::make_unique<NS_GO>(gfs.ns,gfs.nsCC,gfs.ns,gfs.nsCC,*nsLop,mbe); nsGo = std::make_unique<NS_GO>(gfs.ns,gfs.nsCC,gfs.ns,gfs.nsCC,*nsLop,mbe);
chGo = std::make_unique<CH_GO>(gfs.ch,gfs.chCC,gfs.ch,gfs.chCC,*chLop,mbe); chGo = std::make_unique<CH_GO>(gfs.ch,gfs.chCC,gfs.ch,gfs.chCC,*chLop,mbe);
// Linear solver // Linear solver
chLs = std::make_unique<CH_LS>(*chGo,500,100,3,1); chLs = std::make_unique<CH_LS>(*chGo,500,100,3,0);
nsLs = std::make_unique<NS_LS>(*nsGo,1000,200,3,1); nsLs = std::make_unique<NS_LS>(*nsGo,1000,200,3,0);
// Nonlinear solver // Nonlinear solver
nsNewton = std::make_unique<NS_Newton>(*nsGo,nsV,*nsLs); nsNewton = std::make_unique<NS_Newton>(*nsGo,nsV,*nsLs);
...@@ -131,8 +131,8 @@ public: ...@@ -131,8 +131,8 @@ public:
chNewton = std::make_unique<CH_Newton>(*chGo,chV,*chLs); chNewton = std::make_unique<CH_Newton>(*chGo,chV,*chLs);
chNewton->setParameters(param.chNewtonTree); chNewton->setParameters(param.chNewtonTree);
printVector(chV,10,"chV"); //printVector(chV,10,"chV");
printVector(nsV,10,"nsV"); //printVector(nsV,10,"nsV");
updatePhiC(); updatePhiC();
updatePhiSolid(); updatePhiSolid();
...@@ -174,7 +174,7 @@ public: ...@@ -174,7 +174,7 @@ public:
Dune::FieldVector<double,1> integral; Dune::FieldVector<double,1> integral;
Dune::PDELab::integrateGridFunction(productFunction,integral,2); Dune::PDELab::integrateGridFunction(productFunction,integral,2);
kf = integral[0]; kf = integral[0];
std::cout << "Kf = " << kf << std::endl; //std::cout << "Kf = " << kf << std::endl;
} }
void updatePhiC() void updatePhiC()
......
...@@ -146,7 +146,7 @@ template<typename Param, typename PAdaptor, typename CGFS, typename CContainer, ...@@ -146,7 +146,7 @@ template<typename Param, typename PAdaptor, typename CGFS, typename CContainer,
// contribution to residual on inside and outside elements // contribution to residual on inside and outside elements
auto dudn = (outsideC-insideC)/distance; auto dudn = (outsideC-insideC)/distance;
RF peclet = 5; RF peclet = 5; //TODO THIS HAS TO BE FIXED IN OPTIONS!!!!!!!!!!
auto DiffusionCoeff = 1/peclet * 0.5*(pAdapt.evalPhiC(inside_global) + pAdapt.evalPhiC(outside_global)); auto DiffusionCoeff = 1/peclet * 0.5*(pAdapt.evalPhiC(inside_global) + pAdapt.evalPhiC(outside_global));
r_i.accumulate(lfsv_i,0,-dudn*DiffusionCoeff + flux); r_i.accumulate(lfsv_i,0,-dudn*DiffusionCoeff + flux);
r_o.accumulate(lfsv_o,0, dudn*DiffusionCoeff - flux); r_o.accumulate(lfsv_o,0, dudn*DiffusionCoeff - flux);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment