You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: agglomeration_poisson/README.md
+1-15Lines changed: 1 addition & 15 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -1,4 +1,4 @@
1
-
# A posteriori error estimator for first order hyperbolic problems
1
+
# An agglomeration-based polytopic Discontinuous Galerkin solver for the Poisson equation
2
2
3
3
4
4
## Running the code:
@@ -13,13 +13,6 @@ on the command line to configure the program. After that you can compile with `m
13
13
14
14
on the command line.
15
15
16
-
### Parameter file:
17
-
18
-
If you run `./DG_advection_reaction parameters.prm`, an error message will tell you that a parameter file has been created for you. You can open it and change some useful parameters like the number of refinement cycles, the advection coefficient, and others. If you don't specify anything, then the default values used for the test case (see paragraph below) will be used.
19
-
20
-
21
-
22
-
23
16
## The problem:
24
17
This program solves the problem, for $\Omega \in \mathbb{R^2}$
25
18
@@ -40,8 +33,6 @@ for some positive $\gamma_0$ so that we have coercivity in $L^2$ at the continuo
40
33
41
34
## The weak formulation:
42
35
43
-
44
-
45
36
As trial space we choose $V_h = \{ v_h \in L^2(\Omega): v_h \in P^1(\mathbb{T_h})\} \notin H^1(\Omega)$. If we integrate by parts and sum over all cells
46
37
47
38
@f[
@@ -95,8 +86,6 @@ valid for $u \in H^{k+1}(\Omega)$.
95
86
96
87
See Brezzi-Marini-Süli [3] for more details.
97
88
98
-
99
-
100
89
## A-posteriori error estimator:
101
90
102
91
The estimator is the one proposed by Georgoulis, Edward Hall and Charalambos Makridakis in [3]. This approach is quite different with respect to other works in the field, as the authors are trying to develop an estimator for the original hyperbolic problem, rather than taking the hyperbolic regime as the vanishing diffusivity limit.
@@ -121,8 +110,6 @@ where:
121
110
122
111
- $u_h^+$ is the interior trace from the current cell $T$ of a the finite element function $u_h$.
123
112
124
-
125
-
126
113
## Test case:
127
114
128
115
The following test case has been taken from [3]. Consider:
@@ -138,7 +125,6 @@ The next image is the 3D view of the numerical solution:
138
125
More interestingly, we see that the estimator has been able to capture the layer. Here a bulk-chasing criterion is used, with bottom fraction ´0.5´ and no coarsening. This mesh is obtained after 12 refinement cycles.
If we look at the decrease of the energy norm of the error in the globally refined case and in the adaptively case, with respect to the DoFs, we obtain:
0 commit comments