From 3829795e266ac6ecc4ab4a558215ddb6d7da3a03 Mon Sep 17 00:00:00 2001
From: Lars von Wolff <lars.von-wolff@ians.uni-stuttgart.de>
Date: Fri, 9 Apr 2021 13:28:14 +0200
Subject: [PATCH] Reaction for second upscaled example
---
src/paper_thinstrip_mixed/tmixed_upscaledc.hh | 8 +++++++-
1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/src/paper_thinstrip_mixed/tmixed_upscaledc.hh b/src/paper_thinstrip_mixed/tmixed_upscaledc.hh
index 83aa75a..bea9c72 100644
--- a/src/paper_thinstrip_mixed/tmixed_upscaledc.hh
+++ b/src/paper_thinstrip_mixed/tmixed_upscaledc.hh
@@ -78,8 +78,14 @@ template<typename Param, typename PAdaptor, typename CGFS, typename CContainer,
auto phicOld = pAdapt.evalPhiCOld(xg);
auto phiSolid = pAdapt.evalPhiSolid(xg);
auto phiSolidOld = pAdapt.evalPhiSolidOld(xg);
+ RF source = 0;
+ if(xg > 0.1 && xg < 0.3)
+ {
+ source = (xg-0.1)*(0.3-xg);
+ }
r.accumulate(lfsv,0, eg.geometry().volume()*
- ((phic*c - phicOld*cOld + param.phys.uAst*(phiSolid-phiSolidOld))/param.time.dt));
+ ((phic*c - phicOld*cOld + param.phys.uAst*(phiSolid-phiSolidOld))/param.time.dt
+ + source));
//std::cout << "bye alpha" << std::endl;
} //alpha volume
--
GitLab