OpenCAEPoro  v0.5.0
A simulator for multicomponent porous media flow
OCPFluidMethod.hpp
Go to the documentation of this file.
1 
12 #ifndef __OCPFLUIDMETHOD_HEADER__
13 #define __OCPFLUIDMETHOD_HEADER__
14 
15 // OpenCAEPoro header files
16 #include "LinearSystem.hpp"
17 #include "OCPControl.hpp"
18 #include "Reservoir.hpp"
19 #include "UtilOutput.hpp"
20 #include "UtilTiming.hpp"
21 
22 class IsothermalMethod
23 {
24 public:
25  void InitRock(Bulk& bk) const;
26  void CalRock(Bulk& bk) const;
27 };
28 
30 class IsoT_IMPEC : virtual public IsothermalMethod
31 {
32 public:
34  void Setup(Reservoir& rs, LinearSystem& ls, const OCPControl& ctrl);
36  void InitReservoir(Reservoir& rs) const;
38  void Prepare(Reservoir& rs, OCPControl& ctrl);
40  void AssembleMat(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
42  void SolveLinearSystem(LinearSystem& ls, Reservoir& rs, OCPControl& ctrl);
44  OCP_BOOL UpdateProperty(Reservoir& rs, OCPControl& ctrl);
46  OCP_BOOL FinishNR(const Reservoir& rs);
47  void FinishStep(Reservoir& rs, OCPControl& ctrl);
48 
49 protected:
51  void InitFlash(Bulk& bk) const;
53  void CalKrPc(Bulk& bk) const;
55  void PassFlashValue(Bulk& bk, const OCP_USI& n) const;
56 
57 private:
59  void AllocateReservoir(Reservoir& rs);
61  void
62  AllocateLinearSystem(LinearSystem& ls, const Reservoir& rs, const OCPControl& ctrl);
64  void CalFlash(Bulk& bk);
66  void CalFlux(Reservoir& rs) const;
68  void CalBulkFlux(Reservoir& rs) const;
70  void MassConserve(Reservoir& rs, const OCP_DBL& dt) const;
72  void
73  AssembleMatBulks(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
75  void
76  AssembleMatWells(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
78  void AssembleMatInjWells(LinearSystem& ls,
79  const Bulk& bk,
80  const Well& wl,
81  const OCP_DBL& dt) const;
83  void AssembleMatProdWells(LinearSystem& ls,
84  const Bulk& bk,
85  const Well& wl,
86  const OCP_DBL& dt) const;
88  void GetSolution(Reservoir& rs, const vector<OCP_DBL>& u);
90  void ResetToLastTimeStep01(Reservoir& rs, OCPControl& ctrl);
91  void ResetToLastTimeStep02(Reservoir& rs, OCPControl& ctrl);
92  void ResetToLastTimeStep03(Reservoir& rs, OCPControl& ctrl);
94  void UpdateLastTimeStep(Reservoir& rs) const;
95 };
96 
98 class IsoT_FIM : virtual public IsothermalMethod
99 {
100 public:
102  void Setup(Reservoir& rs, LinearSystem& ls, const OCPControl& ctrl);
104  void InitReservoir(Reservoir& rs) const;
106  void Prepare(Reservoir& rs, const OCP_DBL& dt);
108  void AssembleMat(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
110  void SolveLinearSystem(LinearSystem& ls, Reservoir& rs, OCPControl& ctrl) const;
112  OCP_BOOL UpdateProperty(Reservoir& rs, OCPControl& ctrl);
114  OCP_BOOL FinishNR(Reservoir& rs, OCPControl& ctrl);
116  void FinishStep(Reservoir& rs, OCPControl& ctrl);
117 
118 protected:
120  void AllocateReservoir(Reservoir& rs);
122  void
123  AllocateLinearSystem(LinearSystem& ls, const Reservoir& rs, const OCPControl& ctrl);
125  void PassFlashValue(Bulk& bk, const OCP_USI& n) const;
127  void CalKrPc(Bulk& bk) const;
129  void CalRes(Reservoir& rs, const OCP_DBL& dt, const OCP_BOOL& resetRes0) const;
131  void
132  AssembleMatWells(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
133  void
134  AssembleMatWellsNew(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
136  void ResetToLastTimeStep(Reservoir& rs, OCPControl& ctrl);
138  void UpdateLastTimeStep(Reservoir& rs) const;
139 
140 private:
142  void InitFlash(Bulk& bk) const;
144  void CalFlash(Bulk& bk);
146  void
147  AssembleMatBulks(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
149  void AssembleMatInjWells(LinearSystem& ls,
150  const Bulk& bk,
151  const Well& wl,
152  const OCP_DBL& dt) const;
154  void AssembleMatProdWells(LinearSystem& ls,
155  const Bulk& bk,
156  const Well& wl,
157  const OCP_DBL& dt) const;
159  void
160  AssembleMatBulksNew(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
161  void AssembleMatBulksNewS(LinearSystem& ls,
162  const Reservoir& rs,
163  const OCP_DBL& dt) const;
165  void AssembleMatInjWellsNew(LinearSystem& ls,
166  const Bulk& bk,
167  const Well& wl,
168  const OCP_DBL& dt) const;
170  void AssembleMatProdWellsNew(LinearSystem& ls,
171  const Bulk& bk,
172  const Well& wl,
173  const OCP_DBL& dt) const;
175  void
176  GetSolution(Reservoir& rs, const vector<OCP_DBL>& u, const OCPControl& ctrl) const;
177 };
178 
179 class IsoT_FIMn : protected IsoT_FIM
180 {
181 public:
183  void Setup(Reservoir& rs, LinearSystem& ls, const OCPControl& ctrl);
185  void InitReservoir(Reservoir& rs) const;
187  void Prepare(Reservoir& rs, const OCP_DBL& dt);
189  void AssembleMat(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
191  void SolveLinearSystem(LinearSystem& ls, Reservoir& rs, OCPControl& ctrl) const;
193  OCP_BOOL UpdateProperty(Reservoir& rs, OCPControl& ctrl);
195  OCP_BOOL FinishNR(Reservoir& rs, OCPControl& ctrl);
197  void FinishStep(Reservoir& rs, OCPControl& ctrl) const;
198 
199 protected:
201  void AllocateReservoir(Reservoir& rs);
203  void InitFlash(Bulk& bk) const;
205  void CalFlash(Bulk& bk);
207  void PassFlashValue(Bulk& bk, const OCP_USI& n) const;
209  void
210  AssembleMatBulksNew(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
212  void
213  AssembleMatWellsNew(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
215  void AssembleMatInjWellsNew(LinearSystem& ls,
216  const Bulk& bk,
217  const Well& wl,
218  const OCP_DBL& dt) const;
220  void AssembleMatProdWellsNew(LinearSystem& ls,
221  const Bulk& bk,
222  const Well& wl,
223  const OCP_DBL& dt) const;
225  void
226  GetSolution(Reservoir& rs, const vector<OCP_DBL>& u, const OCPControl& ctrl) const;
228  void ResetToLastTimeStep(Reservoir& rs, OCPControl& ctrl);
230  void UpdateLastTimeStep(Reservoir& rs) const;
231 };
232 
233 class IsoT_AIMc : protected IsoT_IMPEC, protected IsoT_FIM
234 {
235 public:
237  void Setup(Reservoir& rs, LinearSystem& ls, const OCPControl& ctrl);
239  void InitReservoir(Reservoir& rs) const;
241  void Prepare(Reservoir& rs, const OCP_DBL& dt);
243  void AssembleMat(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
245  void SolveLinearSystem(LinearSystem& ls, Reservoir& rs, OCPControl& ctrl);
247  OCP_BOOL UpdateProperty(Reservoir& rs, OCPControl& ctrl);
249  OCP_BOOL FinishNR(Reservoir& rs, OCPControl& ctrl);
251  void FinishStep(Reservoir& rs, OCPControl& ctrl) const;
252 
253 protected:
255  void AllocateReservoir(Reservoir& rs);
257  void SetFIMBulk(Reservoir& rs);
259  void CalFlashEp(Bulk& bk);
261  void CalFlashEa(Bulk& bk);
263  void CalFlashI(Bulk& bk);
265  void PassFlashValueEp(Bulk& bk, const OCP_USI& n);
267  void CalKrPcE(Bulk& bk);
269  void CalKrPcI(Bulk& bk);
271  void
272  AssembleMatBulks(LinearSystem& ls, const Reservoir& rs, const OCP_DBL& dt) const;
274  void
275  GetSolution(Reservoir& rs, const vector<OCP_DBL>& u, const OCPControl& ctrl) const;
277  void ResetToLastTimeStep(Reservoir& rs, OCPControl& ctrl);
279  void UpdateLastTimeStep(Reservoir& rs) const;
280 };
281 
282 #endif /* end if __OCPFLUIDMETHOD_HEADER__ */
283 
284 /*----------------------------------------------------------------------------*/
285 /* Brief Change History of This File */
286 /*----------------------------------------------------------------------------*/
287 /* Author Date Actions */
288 /*----------------------------------------------------------------------------*/
289 /* Shizhe Li Oct/01/2021 Create file */
290 /* Chensong Zhang Oct/15/2021 Format file */
291 /*----------------------------------------------------------------------------*/
LinearSystem.hpp
Linear solver class declaration.
OCP_DBL
double OCP_DBL
Double precision.
Definition: OCPConst.hpp:27
OCP_USI
unsigned int OCP_USI
Long unsigned integer.
Definition: OCPConst.hpp:25
OCP_BOOL
unsigned int OCP_BOOL
OCP_BOOL in OCP.
Definition: OCPConst.hpp:29
OCPControl.hpp
OCPControl class declaration.
Reservoir.hpp
Reservoir class declaration.
UtilOutput.hpp
Supply basic tools used to output files.
UtilTiming.hpp
Elapsed wall-time and CPU-cycles declaration.
Bulk
Physical information of each active reservoir bulk.
Definition: Bulk.hpp:104
IsoT_AIMc
Definition: OCPFluidMethod.hpp:234
IsoT_AIMc::GetSolution
void GetSolution(Reservoir &rs, const vector< OCP_DBL > &u, const OCPControl &ctrl) const
Update P, Ni, BHP after linear system is solved.
Definition: OCPFluidMethod.cpp:4843
IsoT_AIMc::AssembleMatBulks
void AssembleMatBulks(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble linear system for bulks.
Definition: OCPFluidMethod.cpp:4380
IsoT_AIMc::CalFlashEa
void CalFlashEa(Bulk &bk)
Perform flash calculation with Ni for Explicit bulk – Update all properties.
Definition: OCPFluidMethod.cpp:4276
IsoT_AIMc::SolveLinearSystem
void SolveLinearSystem(LinearSystem &ls, Reservoir &rs, OCPControl &ctrl)
Solve the linear system.
Definition: OCPFluidMethod.cpp:4069
IsoT_AIMc::CalFlashEp
void CalFlashEp(Bulk &bk)
Perform flash calculation with Ni for Explicit bulk – Update partial properties.
Definition: OCPFluidMethod.cpp:4257
IsoT_AIMc::CalKrPcI
void CalKrPcI(Bulk &bk)
Calculate relative permeability and capillary pressure for Implicit bulk.
Definition: OCPFluidMethod.cpp:4363
IsoT_AIMc::Setup
void Setup(Reservoir &rs, LinearSystem &ls, const OCPControl &ctrl)
Setup AIMc.
Definition: OCPFluidMethod.cpp:4020
IsoT_AIMc::AssembleMat
void AssembleMat(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble Matrix.
Definition: OCPFluidMethod.cpp:4060
IsoT_AIMc::FinishNR
OCP_BOOL FinishNR(Reservoir &rs, OCPControl &ctrl)
Finish a Newton-Raphson iteration.
Definition: OCPFluidMethod.cpp:4123
IsoT_AIMc::CalFlashI
void CalFlashI(Bulk &bk)
Perform flash calculation with Ni for Implicit bulk.
Definition: OCPFluidMethod.cpp:4296
IsoT_AIMc::SetFIMBulk
void SetFIMBulk(Reservoir &rs)
Determine which bulk are treated Implicit.
Definition: OCPFluidMethod.cpp:4191
IsoT_AIMc::UpdateProperty
OCP_BOOL UpdateProperty(Reservoir &rs, OCPControl &ctrl)
Update properties of fluids.
Definition: OCPFluidMethod.cpp:4098
IsoT_AIMc::UpdateLastTimeStep
void UpdateLastTimeStep(Reservoir &rs) const
Update values of last step for AIMc.
Definition: OCPFluidMethod.cpp:4983
IsoT_AIMc::PassFlashValueEp
void PassFlashValueEp(Bulk &bk, const OCP_USI &n)
Pass flash value needed for Explicit bulk – Update partial properties.
Definition: OCPFluidMethod.cpp:4317
IsoT_AIMc::CalKrPcE
void CalKrPcE(Bulk &bk)
Calculate relative permeability and capillary pressure for Explicit bulk.
Definition: OCPFluidMethod.cpp:4348
IsoT_AIMc::ResetToLastTimeStep
void ResetToLastTimeStep(Reservoir &rs, OCPControl &ctrl)
Reset variables to last time step.
Definition: OCPFluidMethod.cpp:4976
IsoT_AIMc::Prepare
void Prepare(Reservoir &rs, const OCP_DBL &dt)
Prepare for Assembling matrix.
Definition: OCPFluidMethod.cpp:4043
IsoT_AIMc::InitReservoir
void InitReservoir(Reservoir &rs) const
Init.
Definition: OCPFluidMethod.cpp:4028
IsoT_AIMc::FinishStep
void FinishStep(Reservoir &rs, OCPControl &ctrl) const
Finish a time step.
Definition: OCPFluidMethod.cpp:4165
IsoT_AIMc::AllocateReservoir
void AllocateReservoir(Reservoir &rs)
Allocate memory for reservoir.
Definition: OCPFluidMethod.cpp:4174
IsoT_FIM
IsoT_FIM is FIM (Fully Implicit Method).
Definition: OCPFluidMethod.hpp:99
IsoT_FIM::Setup
void Setup(Reservoir &rs, LinearSystem &ls, const OCPControl &ctrl)
Setup FIM.
Definition: OCPFluidMethod.cpp:828
IsoT_FIM::AssembleMat
void AssembleMat(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble Matrix.
Definition: OCPFluidMethod.cpp:860
IsoT_FIM::UpdateLastTimeStep
void UpdateLastTimeStep(Reservoir &rs) const
Update values of last step for FIM.
Definition: OCPFluidMethod.cpp:2915
IsoT_FIM::InitReservoir
void InitReservoir(Reservoir &rs) const
Init.
Definition: OCPFluidMethod.cpp:836
IsoT_FIM::FinishNR
OCP_BOOL FinishNR(Reservoir &rs, OCPControl &ctrl)
Finish a Newton-Raphson iteration.
Definition: OCPFluidMethod.cpp:939
IsoT_FIM::AllocateLinearSystem
void AllocateLinearSystem(LinearSystem &ls, const Reservoir &rs, const OCPControl &ctrl)
Allocate memory for linear system.
Definition: OCPFluidMethod.cpp:1086
IsoT_FIM::Prepare
void Prepare(Reservoir &rs, const OCP_DBL &dt)
Prepare for Assembling matrix.
Definition: OCPFluidMethod.cpp:852
IsoT_FIM::FinishStep
void FinishStep(Reservoir &rs, OCPControl &ctrl)
Finish a time step.
Definition: OCPFluidMethod.cpp:972
IsoT_FIM::PassFlashValue
void PassFlashValue(Bulk &bk, const OCP_USI &n) const
Pass value needed for FIM from flash to bulk.
Definition: OCPFluidMethod.cpp:1122
IsoT_FIM::SolveLinearSystem
void SolveLinearSystem(LinearSystem &ls, Reservoir &rs, OCPControl &ctrl) const
Solve the linear system.
Definition: OCPFluidMethod.cpp:877
IsoT_FIM::UpdateProperty
OCP_BOOL UpdateProperty(Reservoir &rs, OCPControl &ctrl)
Update properties of fluids.
Definition: OCPFluidMethod.cpp:914
IsoT_FIM::CalRes
void CalRes(Reservoir &rs, const OCP_DBL &dt, const OCP_BOOL &resetRes0) const
Calculate residual.
Definition: OCPFluidMethod.cpp:1201
IsoT_FIM::CalKrPc
void CalKrPc(Bulk &bk) const
Calculate relative permeability and capillary pressure needed for FIM.
Definition: OCPFluidMethod.cpp:1189
IsoT_FIM::ResetToLastTimeStep
void ResetToLastTimeStep(Reservoir &rs, OCPControl &ctrl)
Reset variables to last time step.
Definition: OCPFluidMethod.cpp:2859
IsoT_FIM::AllocateReservoir
void AllocateReservoir(Reservoir &rs)
Allocate memory for reservoir.
Definition: OCPFluidMethod.cpp:980
IsoT_FIM::AssembleMatWells
void AssembleMatWells(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble linear system for wells.
Definition: OCPFluidMethod.cpp:2135
IsoT_FIMn
Definition: OCPFluidMethod.hpp:180
IsoT_FIMn::AssembleMatBulksNew
void AssembleMatBulksNew(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble linear system for bulks.
Definition: OCPFluidMethod.cpp:3198
IsoT_FIMn::FinishStep
void FinishStep(Reservoir &rs, OCPControl &ctrl) const
Finish a time step.
Definition: OCPFluidMethod.cpp:3103
IsoT_FIMn::Prepare
void Prepare(Reservoir &rs, const OCP_DBL &dt)
Prepare for Assembling matrix.
Definition: OCPFluidMethod.cpp:2991
IsoT_FIMn::AssembleMatProdWellsNew
void AssembleMatProdWellsNew(LinearSystem &ls, const Bulk &bk, const Well &wl, const OCP_DBL &dt) const
Assemble linear system for production wells.
Definition: OCPFluidMethod.cpp:3679
IsoT_FIMn::GetSolution
void GetSolution(Reservoir &rs, const vector< OCP_DBL > &u, const OCPControl &ctrl) const
Update P, Ni, BHP after linear system is solved.
Definition: OCPFluidMethod.cpp:3847
IsoT_FIMn::AllocateReservoir
void AllocateReservoir(Reservoir &rs)
Allocate memory for reservoir.
Definition: OCPFluidMethod.cpp:3111
IsoT_FIMn::ResetToLastTimeStep
void ResetToLastTimeStep(Reservoir &rs, OCPControl &ctrl)
Reset variables to last time step.
Definition: OCPFluidMethod.cpp:3999
IsoT_FIMn::UpdateProperty
OCP_BOOL UpdateProperty(Reservoir &rs, OCPControl &ctrl)
Update properties of fluids.
Definition: OCPFluidMethod.cpp:3042
IsoT_FIMn::AssembleMatInjWellsNew
void AssembleMatInjWellsNew(LinearSystem &ls, const Bulk &bk, const Well &wl, const OCP_DBL &dt) const
Assemble linear system for injection wells.
Definition: OCPFluidMethod.cpp:3524
IsoT_FIMn::AssembleMatWellsNew
void AssembleMatWellsNew(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble linear system for wells.
Definition: OCPFluidMethod.cpp:3503
IsoT_FIMn::InitFlash
void InitFlash(Bulk &bk) const
Perform Flash with Sj and calculate values needed for FIMn.
Definition: OCPFluidMethod.cpp:3129
IsoT_FIMn::PassFlashValue
void PassFlashValue(Bulk &bk, const OCP_USI &n) const
Pass value needed for FIMn from flash to bulk.
Definition: OCPFluidMethod.cpp:3179
IsoT_FIMn::AssembleMat
void AssembleMat(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble Matrix.
Definition: OCPFluidMethod.cpp:2998
IsoT_FIMn::Setup
void Setup(Reservoir &rs, LinearSystem &ls, const OCPControl &ctrl)
Setup FIM.
Definition: OCPFluidMethod.cpp:2967
IsoT_FIMn::InitReservoir
void InitReservoir(Reservoir &rs) const
Init.
Definition: OCPFluidMethod.cpp:2975
IsoT_FIMn::CalFlash
void CalFlash(Bulk &bk)
Perform Flash with Ni and calculate values needed for FIMn.
Definition: OCPFluidMethod.cpp:3152
IsoT_FIMn::SolveLinearSystem
void SolveLinearSystem(LinearSystem &ls, Reservoir &rs, OCPControl &ctrl) const
Solve the linear system.
Definition: OCPFluidMethod.cpp:3008
IsoT_FIMn::FinishNR
OCP_BOOL FinishNR(Reservoir &rs, OCPControl &ctrl)
Finish a Newton-Raphson iteration.
Definition: OCPFluidMethod.cpp:3068
IsoT_FIMn::UpdateLastTimeStep
void UpdateLastTimeStep(Reservoir &rs) const
Update values of last step for FIMn.
Definition: OCPFluidMethod.cpp:4008
IsoT_IMPEC
IsoT_IMPEC is IMPEC (implicit pressure explict saturation) method.
Definition: OCPFluidMethod.hpp:31
IsoT_IMPEC::Prepare
void Prepare(Reservoir &rs, OCPControl &ctrl)
Prepare for Assembling matrix.
Definition: OCPFluidMethod.cpp:64
IsoT_IMPEC::SolveLinearSystem
void SolveLinearSystem(LinearSystem &ls, Reservoir &rs, OCPControl &ctrl)
Solve the linear system.
Definition: OCPFluidMethod.cpp:79
IsoT_IMPEC::AssembleMat
void AssembleMat(LinearSystem &ls, const Reservoir &rs, const OCP_DBL &dt) const
Assemble Matrix.
Definition: OCPFluidMethod.cpp:71
IsoT_IMPEC::InitFlash
void InitFlash(Bulk &bk) const
Perform Flash with Sj and calculate values needed for FIM.
Definition: OCPFluidMethod.cpp:251
IsoT_IMPEC::Setup
void Setup(Reservoir &rs, LinearSystem &ls, const OCPControl &ctrl)
Setup IMPEC.
Definition: OCPFluidMethod.cpp:38
IsoT_IMPEC::UpdateProperty
OCP_BOOL UpdateProperty(Reservoir &rs, OCPControl &ctrl)
Update properties of fluids.
Definition: OCPFluidMethod.cpp:111
IsoT_IMPEC::PassFlashValue
void PassFlashValue(Bulk &bk, const OCP_USI &n) const
Pass value needed for FIM from flash to bulk.
Definition: OCPFluidMethod.cpp:275
IsoT_IMPEC::CalKrPc
void CalKrPc(Bulk &bk) const
Calculate relative permeability and capillary pressure needed for FIM.
Definition: OCPFluidMethod.cpp:311
IsoT_IMPEC::FinishNR
OCP_BOOL FinishNR(const Reservoir &rs)
Determine if NR iteration finishes.
Definition: OCPFluidMethod.cpp:154
IsoT_IMPEC::InitReservoir
void InitReservoir(Reservoir &rs) const
Init.
Definition: OCPFluidMethod.cpp:47
LinearSystem
Linear solvers for discrete systems.
Definition: LinearSystem.hpp:31
OCPControl
All control parameters except for well controllers.
Definition: OCPControl.hpp:94
Well
Definition: Well.hpp:45
Generated on Tue Feb 7 2023 01:49:51 for OpenCAEPoro by doxygen 1.9.1