mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-01-18 21:02:57 -06:00
Added IncompTpfaSinglePhase class.
This commit is contained in:
Atgeirr Flø Rasmussen
parent
fa16d8f616
commit
998d5f9d5a
155
opm/core/pressure/IncompTpfaSinglePhase.cpp
Normal file
155
opm/core/pressure/IncompTpfaSinglePhase.cpp
Normal file
@ -0,0 +1,155 @@
|
||||
/*
|
||||
Copyright 2015 SINTEF ICT, Applied Mathematics.
|
||||
|
||||
This file is part of the Open Porous Media project (OPM).
|
||||
|
||||
OPM is free software: you can redistribute it and/or modify
|
||||
it under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OPM is distributed in the hope that it will be useful,
|
||||
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
||||
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
||||
GNU General Public License for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
||||
*/
|
||||
|
||||
|
||||
#include "config.h"
|
||||
#include <opm/core/pressure/IncompTpfaSinglePhase.hpp>
|
||||
|
||||
#include <opm/core/props/IncompPropertiesSinglePhase.hpp>
|
||||
#include <opm/core/pressure/tpfa/ifs_tpfa.h>
|
||||
#include <opm/core/pressure/tpfa/trans_tpfa.h>
|
||||
// #include <opm/core/pressure/mimetic/mimetic.h>
|
||||
// #include <opm/core/pressure/flow_bc.h>
|
||||
#include <opm/core/linalg/LinearSolverInterface.hpp>
|
||||
#include <opm/core/linalg/sparse_sys.h>
|
||||
// #include <opm/core/simulator/TwophaseState.hpp>
|
||||
// #include <opm/core/simulator/WellState.hpp>
|
||||
#include <opm/core/utility/ErrorMacros.hpp>
|
||||
// #include <opm/core/utility/miscUtilities.hpp>
|
||||
#include <opm/core/wells.h>
|
||||
// #include <iostream>
|
||||
// #include <iomanip>
|
||||
// #include <cmath>
|
||||
// #include <algorithm>
|
||||
|
||||
namespace Opm
|
||||
{
|
||||
|
||||
|
||||
|
||||
/// Construct solver for incompressible case.
|
||||
/// \param[in] grid A 2d or 3d grid.
|
||||
/// \param[in] props Rock and fluid properties.
|
||||
/// \param[in] linsolver Linear solver to use.
|
||||
/// \param[in] wells The wells used as driving forces.
|
||||
IncompTpfaSinglePhase::IncompTpfaSinglePhase(const UnstructuredGrid& grid,
|
||||
const IncompPropertiesSinglePhase& props,
|
||||
const LinearSolverInterface& linsolver,
|
||||
const Wells& wells)
|
||||
: grid_(grid),
|
||||
props_(props),
|
||||
linsolver_(linsolver),
|
||||
wells_(wells),
|
||||
htrans_(grid.cell_facepos[ grid.number_of_cells ]),
|
||||
trans_ (grid.number_of_faces),
|
||||
zeros_(grid.number_of_faces)
|
||||
{
|
||||
computeStaticData();
|
||||
}
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
/// Destructor.
|
||||
IncompTpfaSinglePhase::~IncompTpfaSinglePhase()
|
||||
{
|
||||
ifs_tpfa_destroy(h_);
|
||||
}
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
/// Solve the pressure equation.
|
||||
void IncompTpfaSinglePhase::solve(std::vector<double>& press,
|
||||
std::vector<double>& flux,
|
||||
std::vector<double>& bhp,
|
||||
std::vector<double>& wellrates)
|
||||
{
|
||||
// Set up properties.
|
||||
computePerSolveDynamicData();
|
||||
|
||||
// Assemble.
|
||||
UnstructuredGrid* gg = const_cast<UnstructuredGrid*>(&grid_);
|
||||
int ok = ifs_tpfa_assemble(gg, &forces_, trans_.data(), zeros_.data(), h_);
|
||||
if (!ok) {
|
||||
OPM_THROW(std::runtime_error, "Failed assembling pressure system.");
|
||||
}
|
||||
|
||||
// Solve.
|
||||
linsolver_.solve(h_->A, h_->b, h_->x);
|
||||
|
||||
// Obtain solution.
|
||||
press.resize(grid_.number_of_cells);
|
||||
flux.resize(grid_.number_of_faces);
|
||||
wellrates.resize(wells_.well_connpos[ wells_.number_of_wells ]);
|
||||
bhp.resize(wells_.number_of_wells);
|
||||
ifs_tpfa_solution soln = { NULL, NULL, NULL, NULL };
|
||||
soln.cell_press = press.data();
|
||||
soln.face_flux = flux.data();
|
||||
soln.well_press = bhp.data();
|
||||
soln.well_flux = wellrates.data();
|
||||
ifs_tpfa_press_flux(gg, &forces_, &trans_[0], h_, &soln);
|
||||
}
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
/// Compute data that never changes (after construction).
|
||||
void IncompTpfaSinglePhase::computeStaticData()
|
||||
{
|
||||
UnstructuredGrid* gg = const_cast<UnstructuredGrid*>(&grid_);
|
||||
tpfa_htrans_compute(gg, props_.permeability(), &htrans_[0]);
|
||||
h_ = ifs_tpfa_construct(gg, const_cast<struct Wells*>(&wells_));
|
||||
}
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
/// Compute per-solve dynamic properties.
|
||||
void IncompTpfaSinglePhase::computePerSolveDynamicData()
|
||||
{
|
||||
// Computed here:
|
||||
//
|
||||
// std::vector<double> totmob_;
|
||||
// std::vector<double> trans_;
|
||||
// ifs_tpfa_forces forces_;
|
||||
|
||||
// totmob_
|
||||
totmob_.clear();
|
||||
totmob_.resize(grid_.number_of_cells, 1.0/(*props_.viscosity()));
|
||||
// trans_
|
||||
tpfa_eff_trans_compute(const_cast<UnstructuredGrid*>(&grid_), totmob_.data(), htrans_.data(), trans_.data());
|
||||
// forces_
|
||||
forces_.src = NULL;
|
||||
forces_.bc = NULL;
|
||||
forces_.W = &wells_;
|
||||
forces_.totmob = totmob_.data();
|
||||
forces_.wdp = zeros_.data();
|
||||
}
|
||||
|
||||
|
||||
} // namespace Opm
|
||||
88
opm/core/pressure/IncompTpfaSinglePhase.hpp
Normal file
88
opm/core/pressure/IncompTpfaSinglePhase.hpp
Normal file
@ -0,0 +1,88 @@
|
||||
/*
|
||||
Copyright 2015 SINTEF ICT, Applied Mathematics.
|
||||
|
||||
This file is part of the Open Porous Media project (OPM).
|
||||
|
||||
OPM is free software: you can redistribute it and/or modify
|
||||
it under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OPM is distributed in the hope that it will be useful,
|
||||
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
||||
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
||||
GNU General Public License for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
||||
*/
|
||||
|
||||
#ifndef OPM_INCOMPTPFASINGLEPHASE_HEADER_INCLUDED
|
||||
#define OPM_INCOMPTPFASINGLEPHASE_HEADER_INCLUDED
|
||||
|
||||
|
||||
#include <opm/core/pressure/tpfa/ifs_tpfa.h>
|
||||
#include <vector>
|
||||
|
||||
struct UnstructuredGrid;
|
||||
struct Wells;
|
||||
|
||||
namespace Opm
|
||||
{
|
||||
|
||||
class IncompPropertiesSinglePhase;
|
||||
class LinearSolverInterface;
|
||||
|
||||
/// Encapsulating a tpfa pressure solver for the incompressible-fluid case.
|
||||
/// Supports gravity, wells controlled by bhp or reservoir rates,
|
||||
/// boundary conditions and simple sources as driving forces.
|
||||
/// Rock compressibility can be included, and necessary nonlinear
|
||||
/// iterations are handled.
|
||||
/// Below we use the shortcuts D for the number of dimensions, N
|
||||
/// for the number of cells and F for the number of faces.
|
||||
class IncompTpfaSinglePhase
|
||||
{
|
||||
public:
|
||||
/// Construct solver for incompressible case.
|
||||
/// \param[in] grid A 2d or 3d grid.
|
||||
/// \param[in] props Rock and fluid properties.
|
||||
/// \param[in] linsolver Linear solver to use.
|
||||
/// \param[in] wells The wells used as driving forces.
|
||||
IncompTpfaSinglePhase(const UnstructuredGrid& grid,
|
||||
const IncompPropertiesSinglePhase& props,
|
||||
const LinearSolverInterface& linsolver,
|
||||
const Wells& wells);
|
||||
|
||||
/// Destructor.
|
||||
~IncompTpfaSinglePhase();
|
||||
|
||||
/// Solve the pressure equation.
|
||||
void solve(std::vector<double>& press,
|
||||
std::vector<double>& flux,
|
||||
std::vector<double>& bhp,
|
||||
std::vector<double>& wellrates);
|
||||
|
||||
private:
|
||||
// Helper functions.
|
||||
void computeStaticData();
|
||||
void computePerSolveDynamicData();
|
||||
|
||||
protected:
|
||||
// ------ Data that will remain unmodified after construction. ------
|
||||
const UnstructuredGrid& grid_;
|
||||
const IncompPropertiesSinglePhase& props_;
|
||||
const LinearSolverInterface& linsolver_;
|
||||
const Wells& wells_;
|
||||
std::vector<double> htrans_;
|
||||
std::vector<double> trans_ ;
|
||||
std::vector<double> zeros_;
|
||||
std::vector<double> totmob_;
|
||||
struct ifs_tpfa_forces forces_;
|
||||
|
||||
// ------ Internal data for the ifs_tpfa solver. ------
|
||||
struct ifs_tpfa_data* h_;
|
||||
};
|
||||
|
||||
} // namespace Opm
|
||||
|
||||
#endif // OPM_INCOMPTPFASINGLEPHASE_HEADER_INCLUDED
|
||||
Loading…
Reference in New Issue
Block a user