96 lines
3.3 KiB
C++
96 lines
3.3 KiB
C++
#include <iostream>
|
|
#include <opm/core/eclipse/EclipseGridParser.hpp>
|
|
#include <opm/core/utility/parameters/ParameterGroup.hpp>
|
|
|
|
|
|
#include "opm/core/utility/initState.hpp"
|
|
#include <opm/core/WellsManager.hpp>
|
|
#include <opm/core/GridManager.hpp>
|
|
#include <opm/core/pressure/IncompTpfa.hpp>
|
|
#include <opm/core/fluid/IncompPropertiesFromDeck.hpp>
|
|
#include <opm/core/linalg/LinearSolverUmfpack.hpp>
|
|
#include <opm/core/newwells.h>
|
|
#include <opm/core/grid.h>
|
|
#include <opm/core/utility/miscUtilities.hpp>
|
|
#include <opm/core/TwophaseState.hpp>
|
|
#include <opm/core/pressure/FlowBCManager.hpp>
|
|
|
|
#ifdef EXPERIMENT_ISTL
|
|
#include <opm/core/linalg/LinearSolverIstl.hpp>
|
|
#endif
|
|
int main(int argc, char** argv) {
|
|
|
|
using namespace Opm::parameter;
|
|
using namespace Opm;
|
|
ParameterGroup parameters( argc, argv, false );
|
|
std::string file_name = parameters.getDefault<std::string>("inputdeck", "data.data");
|
|
|
|
// Read input file
|
|
EclipseGridParser parser(file_name);
|
|
std::cout << "Done!" << std::endl;
|
|
// Setup grid
|
|
GridManager grid(parser);
|
|
|
|
// Finally handle the wells
|
|
WellsManager wells(parser, *grid.c_grid(), NULL);
|
|
|
|
std::vector<int> global_cells(grid.c_grid()->global_cell, grid.c_grid()->global_cell + grid.c_grid()->number_of_cells);
|
|
|
|
double gravity[3] = {0.0, 0.0, parameters.getDefault<double>("gravity", 0.0)};
|
|
IncompPropertiesFromDeck incomp_properties(parser, global_cells);
|
|
|
|
#ifdef EXPERIMENT_ISTL
|
|
Opm::LinearSolverIstl linsolver(parameters);
|
|
#else
|
|
Opm::LinearSolverUmfpack linsolver;
|
|
#endif // EXPERIMENT_ISTL
|
|
IncompTpfa pressure_solver(*grid.c_grid(), incomp_properties.permeability(),
|
|
gravity, linsolver, wells.c_wells());
|
|
|
|
|
|
std::vector<int> all_cells;
|
|
for(int i = 0; i < grid.c_grid()->number_of_cells; i++) {
|
|
all_cells.push_back(i);
|
|
}
|
|
|
|
Opm::TwophaseState state;
|
|
|
|
initStateTwophaseFromDeck(*grid.c_grid(), incomp_properties, parser, gravity[2], state);
|
|
|
|
// Compute total mobility and omega
|
|
std::vector<double> totmob;
|
|
std::vector<double> omega;
|
|
computeTotalMobilityOmega(incomp_properties, all_cells, state.saturation(), totmob, omega);
|
|
|
|
std::vector<double> wdp;
|
|
std::vector<double> densities(incomp_properties.density(), incomp_properties.density() + incomp_properties.numPhases());
|
|
computeWDP(*wells.c_wells(), *grid.c_grid(), state.saturation(), densities, wdp);
|
|
|
|
std::vector<double> src;
|
|
Opm::FlowBCManager bcs;
|
|
|
|
std::vector<double> pressure;
|
|
std::vector<double> face_flux;
|
|
|
|
std::vector<double> well_bhp;
|
|
std::vector<double> well_rate_per_cell;
|
|
pressure_solver.solve(totmob, omega, src, wdp, bcs.c_bcs(), pressure, face_flux, well_bhp, well_rate_per_cell);
|
|
std::cout << "Solved" << std::endl;
|
|
|
|
for(size_t i = 0; i < well_rate_per_cell.size(); i++) {
|
|
std::cout << well_rate_per_cell[i] << std::endl;
|
|
}
|
|
std::vector<double> well_rate;
|
|
|
|
computeFlowRatePerWell(*wells.c_wells(), well_rate_per_cell, well_rate);
|
|
if(wells.wellCollection().conditionsMet(well_bhp, well_rate, *grid.c_grid(), state.saturation() )) {
|
|
std::cout << "Conditions met for wells!" << std::endl;
|
|
}
|
|
else
|
|
{
|
|
std::cout << "Conditions not met for wells!"<<std::endl;
|
|
}
|
|
return 0;
|
|
}
|
|
|