OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Fix: forgotten multiply by cell volume in a quadrature.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-09-27 09:49:36 +02:00
parent 04a76988d9
commit d83ab5856d
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -325,6 +325,7 @@ namespace Opm
// Compute cell jacobian contribution. We use Fortran ordering // Compute cell jacobian contribution. We use Fortran ordering
// for jac_, i.e. rows cycling fastest. // for jac_, i.e. rows cycling fastest.
{ {
const double cvol = grid_.cell_volumes[cell];
CellQuadrature quad(grid_, cell, degree_); CellQuadrature quad(grid_, cell, degree_);
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)
@ -336,7 +337,7 @@ namespace Opm
for (int j = 0; j < num_basis; ++j) { for (int j = 0; j < num_basis; ++j) {
for (int i = 0; i < num_basis; ++i) { for (int i = 0; i < num_basis; ++i) {
for (int dd = 0; dd < dim; ++dd) { for (int dd = 0; dd < dim; ++dd) {
jac_[j*num_basis + i] += w * basis_[i] * grad_basis_[dim*j + dd] * velocity_[dd]; jac_[j*num_basis + i] += w * basis_[i] * grad_basis_[dim*j + dd] * velocity_[dd] * cvol;
} }
} }
} }