OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Increased quadrature order in a term.

Browse Source 
This is for the benefit of bi/tri-linear basis functions, as stated
in the comments.
This commit is contained in:
Atgeirr Flø Rasmussen 2013-01-16 13:45:15 +01:00
parent 33cc611ced
commit 564db5df87
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
opm/core/transport/reorder/TransportModelTracerTofDiscGal.cpp
View File

@ -510,8 +510,15 @@ namespace Opm
// for jac_, i.e. rows cycling fastest. // for jac_, i.e. rows cycling fastest.
{ {
// Even with ECVI velocity interpolation, degree of precision 1 // Even with ECVI velocity interpolation, degree of precision 1
// is sufficient for optimal convergence order for DG1. // is sufficient for optimal convergence order for DG1 when we
const int deg_needed = 2*degree_ - 1; // use linear (total degree 1) basis functions.
// With bi(tri)-linear basis functions, it still seems sufficient
// for convergence order 2, but the solution looks much better and
// has significantly lower error with degree of precision 2.
// For now, we err on the side of caution, and use 2*degree, even
// though this is wasteful for the pure linear basis functions.
// const int deg_needed = 2*degree_ - 1;
const int deg_needed = 2*degree_;
CellQuadrature quad(grid_, cell, deg_needed); CellQuadrature quad(grid_, cell, deg_needed);
for (int quad_pt = 0; quad_pt < quad.numQuadPts(); ++quad_pt) { for (int quad_pt = 0; quad_pt < quad.numQuadPts(); ++quad_pt) {
// b_i (v \cdot \grad b_j) // b_i (v \cdot \grad b_j)