TransportSolverCompressibleTwophaseReorder.hpp

Go to the documentation of this file.

/*
  Copyright 2012 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_TRANSPORTSOLVERCOMPRESSIBLETWOPHASEREORDER_HEADER_INCLUDED
#define OPM_TRANSPORTSOLVERCOMPRESSIBLETWOPHASEREORDER_HEADER_INCLUDED
 
#include <opm/core/transport/reorder/ReorderSolverInterface.hpp>
#include <vector>
 
struct UnstructuredGrid;
 
namespace Opm
{
 
    class BlackoilPropertiesInterface;
 
    class TransportSolverCompressibleTwophaseReorder : public ReorderSolverInterface
    {
    public:
        TransportSolverCompressibleTwophaseReorder(const UnstructuredGrid& grid,
                                           const Opm::BlackoilPropertiesInterface& props,
                                           const double tol,
                                           const int maxit);
 
        void solve(const double* darcyflux,
                   const double* pressure,
                   const double* temperature,
                   const double* porevolume0,
                   const double* porevolume,
                   const double* source,
                   const double dt,
                   std::vector<double>& saturation,
                   std::vector<double>& surfacevol);
 
        void initGravity(const double* grav);
 
        void solveGravity(const std::vector<std::vector<int> >& columns,
                          const double dt,
                          std::vector<double>& saturation,
                          std::vector<double>& surfacevol);
 
    private:
        virtual void solveSingleCell(const int cell);
        virtual void solveMultiCell(const int num_cells, const int* cells);
        void solveSingleCellGravity(const std::vector<int>& cells,
                                    const int pos,
                                    const double* gravflux);
        int solveGravityColumn(const std::vector<int>& cells);
        void initGravityDynamic();
 
    private:
        const UnstructuredGrid& grid_;
        const BlackoilPropertiesInterface& props_;
        std::vector<int> allcells_;
        std::vector<double> visc_;
        std::vector<double> A_;
        std::vector<double> smin_;
        std::vector<double> smax_;
        double tol_;
        int maxit_;
 
        const double* darcyflux_;   // one flux per grid face
        const double* surfacevol0_; // one per phase per cell
        const double* porevolume0_; // one volume per cell
        const double* porevolume_;  // one volume per cell
        const double* source_;      // one source per cell
        double dt_;
        std::vector<double> saturation_;        // P (= num. phases) per cell
        std::vector<double> fractionalflow_;  // = m[0]/(m[0] + m[1]) per cell
        // For gravity segregation.
        const double* gravity_;
        std::vector<double> trans_;
        std::vector<double> density_;
        std::vector<double> gravflux_;
        std::vector<double> mob_;
        std::vector<double> s0_;
 
        // Storing the upwind and downwind graphs for experiments.
        std::vector<int> ia_upw_;
        std::vector<int> ja_upw_;
        std::vector<int> ia_downw_;
        std::vector<int> ja_downw_;
 
        struct Residual;
        double fracFlow(double s, int cell) const;
 
        struct GravityResidual;
        void mobility(double s, int cell, double* mob) const;
    };
 
} // namespace Opm
 
#endif // OPM_TRANSPORTMODELCOMPRESSIBLETWOPHASE_HEADER_INCLUDED