Loading...
Searching...
No Matches
Source.hpp
/*
Copyright 2023 Equinor ASA.
Copyright 2023 Norce.
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_SOURCE_PROP_HPP
#define OPM_SOURCE_PROP_HPP
#include <array>
#include <vector>
#include <cstddef>
#include <optional>
namespace Opm {
class Deck;
class DeckRecord;
enum class SourceComponent {
OIL,
GAS,
WATER,
SOLVENT,
POLYMER,
NONE
};
class Source
{
public:
struct SourceCell
{
std::array<int, 3> ijk{};
SourceComponent component{SourceComponent::NONE};
double rate{};
std::optional<double> hrate{};
std::optional<double> temperature{};
SourceCell() = default;
explicit SourceCell(const DeckRecord& record);
static SourceCell serializationTestObject();
bool operator==(const SourceCell& other) const;
bool isSame(const SourceCell& other) const;
bool isSame(const std::pair<std::array<int, 3>, SourceComponent>& other) const;
template<class Serializer>
void serializeOp(Serializer& serializer)
{
serializer(ijk);
serializer(component);
serializer(rate);
serializer(hrate);
serializer(temperature);
}
};
Source() = default;
static Source serializationTestObject();
std::size_t size() const;
std::vector<SourceCell>::const_iterator begin() const;
std::vector<SourceCell>::const_iterator end() const;
bool operator==(const Source& other) const;
double rate(const std::pair<std::array<int, 3>, SourceComponent>& input ) const;
double hrate(const std::pair<std::array<int, 3>, SourceComponent>& input ) const;
double temperature(const std::pair<std::array<int, 3>, SourceComponent>& input) const;
bool hasHrate(const std::pair<std::array<int, 3>, SourceComponent>& input) const;
bool hasTemperature(const std::pair<std::array<int, 3>, SourceComponent>& input) const;
bool hasSource(const std::array<int, 3>& input) const;
void updateSource(const DeckRecord& record);
void addSourceCell(const SourceCell& cell)
{
m_cells.push_back(cell);
}
template<class Serializer>
void serializeOp(Serializer& serializer)
{
serializer(m_cells);
}
private:
std::vector<SourceCell> m_cells;
};
} // namespace Opm
#endif
Opm::DeckRecord
Definition DeckRecord.hpp:32
Opm::Serializer
Class for (de-)serializing.
Definition Serializer.hpp:91
Opm::Source
Definition Source.hpp:44
Opm::Source::addSourceCell
void addSourceCell(const SourceCell &cell)
Add a source term for a grid cell.
Definition Source.hpp:93
Opm
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30
Opm::Source::SourceCell
Definition Source.hpp:47