From cb4b857c26c42d1efdc9c6d040989ebf5e0f5404 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Maximilian=20H=C3=B6rl?=
<maximilian.hoerl@mathematik.uni-stuttgart.de>
Date: Wed, 20 May 2020 12:18:34 +0200
Subject: [PATCH] modify problem 1
---
dune/mmdg/problems/mmdgproblem1.hh | 40 +++++++++++++++++++++++++++++-
1 file changed, 39 insertions(+), 1 deletion(-)
diff --git a/dune/mmdg/problems/mmdgproblem1.hh b/dune/mmdg/problems/mmdgproblem1.hh
index f35a61f..e1caa9b 100644
--- a/dune/mmdg/problems/mmdgproblem1.hh
+++ b/dune/mmdg/problems/mmdgproblem1.hh
@@ -37,14 +37,52 @@ public:
//aperture d of the fracture at position pos
Scalar aperture (const Coordinate& pos) const
{
+ return d_ * (pos[1] + 0.01);
return d_;
}
+ //interface source term at position pos
+ Scalar qInterface (const Coordinate& pos) const
+ {
+ const Scalar de = aperture(pos);
+ const Scalar dPrime2 = gradAperture(pos) * gradAperture(pos);
+ const Scalar dPrimePrime = 0.0;
+
+ return (-de * dPrimePrime - dPrime2) * de;
+
+ return Scalar(0.0);
+ }
+
+ //tangential gradient of the aperture d of the fracture at position pos
+ Coordinate gradAperture (const Coordinate& pos) const
+ {
+ Coordinate gradD(0.0);
+ gradD[1] = d_;
+ return gradD;
+
+ return Coordinate(0.0);
+ }
+
+ //tangential permeability per aperture tensor of the interface at position pos
+ Matrix Kparallel (const Coordinate& pos) const
+ {
+ Matrix permeability(0.0);
+
+ for (int i = 0; i < dim; i++)
+ {
+ permeability[i][i] = 1.0 * aperture(pos);
+ }
+
+ //return identity matrix
+ return permeability;
+ }
+
+
//returns the recommended quadrature order to compute an integral
//over x * boundary(x)
int quadratureOrderBoundary () const
{
- return 2;
+ return 10;//2;
}
private:
--
GitLab