[DuMuX] Calculation of Original Oil in Place(OOIP)

[Ranjeet kumar]

Hi all,
I want to calculate original oil in place for 2p problem where
OOIP = 10^4 * element.volume()*phi*(1-Sw)/Bo.

I have implemented as follow in the constructor of spatialparams.hh file:

        Sw_ = 0.4;
        for (const auto& element : elements(gridView)) {
            auto geo = element.geometry();
            //get cell index
             int eIdx = gridView.indexSet().index(element);
            OOIP_ += geo.volume()*porosity_[eIdx] *(1-Sw);  //assumed Bo =
1;
        }
        OOIP_ = OOIP_*10000;
        std::cout << "OOIP =" << OOIP_ << " m^3." << std::endl;
Is this correct?

_______________________________________________
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux