perf: Improve MGR recipes for thermal#3939
Conversation
| cprLike = 13, //!< CPR-like restriction | ||
| blockColLumped = 14 //!< Block column-lumped approximation | ||
| blockColLumped = 14, //!< Block column-lumped approximation | ||
| partialColLumped = 15 //!< Partial column-lumped approximation |
There was a problem hiding this comment.
What does "Partial column-lumped approximation" mean?
There was a problem hiding this comment.
It means W_{TS} = 0 instead of doing colsum^{-1}(A_{TS}) * colsum^{-1}(A_{SS}) as done by blockColLumped. W_{TS} being the submatrix of the restriction operator coupling temperature and densities.
The pressure part remains the same: W_{PS} = colsum^{-1}(A_{PS}) * colsum^{-1}(A_{SS})
So, I called it a "partial" column-lumped. Let me know if you have other names in mind!
There was a problem hiding this comment.
So for ThermalCompositionalMultiphaseFVM we have:
| A_SS A_SP A_ST |
| A_PS A_PP A_PT |
| A_TS A_TP A_TT |
With "Partial column-lumped approximation", the correction to block:
| A_PP A_PT |
| A_TP A_TT |
is
- | colsum(A_{PS}) | * colsum^{-1}(A_{SS}) * | A_SP A_ST |
| 0 |
Correct?
There was a problem hiding this comment.
That's correct, so essentially the energy equation isn't perturbed
| m_levelRestrictType[1] = MGRRestrictionType::partialColLumped; // True-IMPES for thermal | ||
| m_levelCoarseGridMethod[1] = MGRCoarseGridMethod::galerkin; |
There was a problem hiding this comment.
I like the use of the Galerkin projection instead of the non-Galerkin. Again, I am not sure I fully understand what partial column lumping is.
There was a problem hiding this comment.
Same comment applies also to the the other two modified strategies
paveltomin
left a comment
There was a problem hiding this comment.
would be interesting to discuss this during our solvers meeting on Thursdays
|
Note: This PR is dependent on #3923 . |
Codecov Report❌ Patch coverage is Additional details and impacted files@@ Coverage Diff @@
## develop #3939 +/- ##
========================================
Coverage 58.87% 58.87%
========================================
Files 1353 1353
Lines 112954 112953 -1
========================================
Hits 66504 66504
+ Misses 46450 46449 -1 ☔ View full report in Codecov by Sentry. 🚀 New features to boost your workflow:
|
…EV/GEOS into feature/paludettomag1/thermal-mgr
Depends on GEOS-DEV/thirdPartyLibs#332