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

Implemented Gauss-Seidel solver for polymer. Not tested!

Browse Source
This commit is contained in:
Xavier Raynaud
2012-05-16 15:06:14 +02:00
parent 7d54d1f31b
commit 60599b9dbf
3 changed files with 341 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
examples/polymer_reorder.cpp
View File

@@ -448,11 +448,7 @@ main(int argc, char** argv)
if (use_segregation_split) {
use_column_solver = param.getDefault("use_column_solver", use_column_solver);
if (use_column_solver) {
// use_gauss_seidel_gravity is not implemented for polymer.
use_gauss_seidel_gravity = param.getDefault("use_gauss_seidel_gravity", use_gauss_seidel_gravity);
if (use_gauss_seidel_gravity) {
THROW("gauss_seidel_gravity is not implemented for polymer");
}
}
}
}
@@ -533,8 +529,7 @@ main(int argc, char** argv)
method, nl_tolerance, nl_maxiter);
if (use_gauss_seidel_gravity) {
THROW("use_gauss_seidel_gravity option not yet implemented for polymer case.");
// reorder_model.initGravity(grav);
reorder_model.initGravity(grav);
}
// Non-reordering solver.
NewtonPolymerTransportModel model(fluid, *grid->c_grid(), porevol, grav, guess_old_solution);
@@ -750,9 +745,9 @@ main(int argc, char** argv)
if (use_segregation_split) {
if (use_column_solver) {
if (use_gauss_seidel_gravity) {
THROW("use_gauss_seidel_gravity option not implemented for polymer.");
// reorder_model.solveGravity(columns, stepsize, reorder_sat);
// Opm::toBothSat(reorder_sat, state.saturation());
reorder_model.solveGravity(columns, &porevol[0], stepsize, reorder_sat,
state.concentration(), state.maxconcentration());
Opm::toBothSat(reorder_sat, state.saturation());
} else {
colsolver.solve(columns, stepsize, state.saturation(), state.concentration(),
state.maxconcentration());