mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
37ce9d0db3
- solver now takes parameter for inflow_c - simulator calls polymerInflowAtTime(double time) for its value - currently starts polymer injection after 4 days (0.4 pv) Also added some (inactive) debug printouts. Still not working when injecting polymer.
85 lines
2.3 KiB
C++
85 lines
2.3 KiB
C++
/* Copyright 2011 (c) Jostein R. Natvig <Jostein.R.Natvig at sintef.no> */
|
|
/* Copyright 2012 (c) SINTEF */
|
|
|
|
#include <stdlib.h>
|
|
#include <stdio.h>
|
|
#include <assert.h>
|
|
|
|
#include <opm/polymer/polymermodel.hpp>
|
|
#include <opm/polymer/polymertransport.hpp>
|
|
#include <opm/core/grid.h>
|
|
#include <opm/core/transport/reorder/reordersequence.h>
|
|
#include <opm/core/transport/reorder/tarjan.h>
|
|
#include <opm/core/transport/reorder/nlsolvers.h>
|
|
|
|
|
|
|
|
void polymertransport(
|
|
const double *porevolume,
|
|
const double *porosity,
|
|
const double *source,
|
|
const double dt,
|
|
const double inflow_c,
|
|
struct UnstructuredGrid *grid,
|
|
const Opm::IncompPropertiesInterface* props,
|
|
const PolymerData* polydata,
|
|
const double *darcyflux,
|
|
double *saturation,
|
|
double *concentration,
|
|
double *cmax)
|
|
{
|
|
int i;
|
|
|
|
int ncomponents;
|
|
int *sequence;
|
|
int *components;
|
|
|
|
PolymerSolverData *data = init_solverdata(grid, props, polydata, darcyflux,
|
|
porevolume, porosity, source, dt,
|
|
inflow_c, saturation,
|
|
concentration, cmax);
|
|
|
|
|
|
struct NonlinearSolverCtrl ctrl;
|
|
|
|
|
|
/* Compute sequence of single-cell problems */
|
|
sequence = (int*) malloc( grid->number_of_cells * sizeof *sequence);
|
|
components = (int*) malloc(( grid->number_of_cells+1) * sizeof *components);
|
|
|
|
compute_sequence(grid, darcyflux, sequence, components, &ncomponents);
|
|
assert(ncomponents == grid->number_of_cells);
|
|
|
|
|
|
|
|
ctrl.maxiterations = 20;
|
|
ctrl.nltolerance = 1e-9;
|
|
ctrl.method = NonlinearSolverCtrl::REGULAFALSI;
|
|
ctrl.initialguess = 0.5;
|
|
|
|
/* Assume all strong components are single-cell domains. */
|
|
for(i=0; i<grid->number_of_cells; ++i)
|
|
{
|
|
#ifdef MATLAB_MEX_FILE
|
|
if (interrupt_signal)
|
|
{
|
|
mexPrintf("Reorder loop interrupted by user: %d of %d "
|
|
"cells finished.\n", i, grid->number_of_cells);
|
|
break;
|
|
}
|
|
#endif
|
|
polymer_solvecell(data, &ctrl, sequence[i]);
|
|
/*print("iterations:%d residual:%20.16f\n",
|
|
* ctrl.iterations, ctrl.residual);*/
|
|
}
|
|
|
|
destroy_solverdata(data);
|
|
|
|
free(sequence);
|
|
free(components);
|
|
}
|
|
|
|
/* Local Variables: */
|
|
/* c-basic-offset:4 */
|
|
/* End: */
|