feat: add MGR for single phase poromechanics with conforming fractures and wells #3731
Conversation
|
Currently, I just borrow the MGR from |
jhuang2601
requested review from
CusiniM,
Guotong-Ren,
frankfeifan,
paveltomin,
rrsettgast and
victorapm
as code owners
|
@castelletto1 @victorapm can you please help? |
|
@jhuang2601 Do we have an input XML to test the new MGR strategy? |
Yes, @victorapm, just push an example to this PR |
|
Code owners' review needed (@castelletto1 or @matteofrigo5 ?) |
|
Hello @victorapm @castelletto1 @matteofrigo5, any update on this MGR recipe, which is required for InSalah case? |
|
@victorapm @castelletto1 @matteofrigo5 : just checking on the status on this. Can you take a look? |
|
We should keep in mind that two strategies are needed, one for the isothermal case and one for the thermal case. |
| * dofLabel: 0 = displacement, x-component | ||
| * dofLabel: 1 = displacement, y-component | ||
| * dofLabel: 2 = displacement, z-component | ||
| * dofLabel: 3 = face-centered lagrange multiplier (tn) | ||
| * dofLabel: 4 = face-centered lagrange multiplier (tt1) | ||
| * dofLabel: 5 = face-centered lagrange multiplier (tt2) | ||
| * dofLabel: 6 = pressure (cell elem + fracture elems) |
There was a problem hiding this comment.
Isothermal case:
| * dofLabel: 0 = displacement, x-component | |
| * dofLabel: 1 = displacement, y-component | |
| * dofLabel: 2 = displacement, z-component | |
| * dofLabel: 3 = face-centered lagrange multiplier (tn) | |
| * dofLabel: 4 = face-centered lagrange multiplier (tt1) | |
| * dofLabel: 5 = face-centered lagrange multiplier (tt2) | |
| * dofLabel: 6 = pressure (cell elem + fracture elems) | |
| * dofLabel: 0 = displacement, x-component | |
| * dofLabel: 1 = displacement, y-component | |
| * dofLabel: 2 = displacement, z-component | |
| * dofLabel: 3 = face-centered lagrange multiplier (tn) | |
| * dofLabel: 4 = face-centered lagrange multiplier (tt1) | |
| * dofLabel: 5 = face-centered lagrange multiplier (tt2) | |
| * dofLabel: 6 = pressure (cell elem + fracture elems) | |
| * dofLabel: 7 = well pressure | |
| * dofLabel: 8 = well rate |
Thermal case (should be another included in another strategy)
| * dofLabel: 0 = displacement, x-component | |
| * dofLabel: 1 = displacement, y-component | |
| * dofLabel: 2 = displacement, z-component | |
| * dofLabel: 3 = face-centered lagrange multiplier (tn) | |
| * dofLabel: 4 = face-centered lagrange multiplier (tt1) | |
| * dofLabel: 5 = face-centered lagrange multiplier (tt2) | |
| * dofLabel: 6 = pressure (cell elem + fracture elems) | |
| * dofLabel: 0 = displacement, x-component | |
| * dofLabel: 1 = displacement, y-component | |
| * dofLabel: 2 = displacement, z-component | |
| * dofLabel: 3 = face-centered lagrange multiplier (tn) | |
| * dofLabel: 4 = face-centered lagrange multiplier (tt1) | |
| * dofLabel: 5 = face-centered lagrange multiplier (tt2) | |
| * dofLabel: 6 = pressure (cell elem + fracture elems) | |
| * dofLabel: 7 = temperature (cell elem + fracture elems) | |
| * dofLabel: 8 = well pressure | |
| * dofLabel: 9 = well termperature | |
| * dofLabel: 10 = well rate |
| * @brief Constructor. | ||
| */ | ||
| explicit SinglePhasePoromechanicsConformingFracturesReservoirFVM( arrayView1d< int const > const & ) | ||
| : MGRStrategyBase( 7 ) |
There was a problem hiding this comment.
| : MGRStrategyBase( 7 ) | |
| : MGRStrategyBase( 9 ) |
| // we keep u and p | ||
| m_labels[0].push_back( 0 ); | ||
| m_labels[0].push_back( 1 ); | ||
| m_labels[0].push_back( 2 ); | ||
| m_labels[0].push_back( 6 ); | ||
| // we keep p | ||
| m_labels[1].push_back( 6 ); | ||
|
|
||
| setupLabels(); | ||
|
|
||
| // Level 0 | ||
| m_levelFRelaxType[0] = MGRFRelaxationType::none; | ||
| m_levelFRelaxIters[0] = 0; | ||
|
|
||
| m_levelGlobalSmootherType[0] = MGRGlobalSmootherType::ilu0; | ||
| m_levelGlobalSmootherIters[0] = 1; | ||
|
|
||
| m_levelInterpType[0] = MGRInterpolationType::blockJacobi; | ||
| m_levelRestrictType[0] = MGRRestrictionType::injection; | ||
| m_levelCoarseGridMethod[0] = MGRCoarseGridMethod::galerkin; | ||
|
|
||
| // Level 1 | ||
| m_levelFRelaxType[1] = MGRFRelaxationType::amgVCycle; | ||
| m_levelFRelaxIters[1] = 1; | ||
| m_levelGlobalSmootherType[1] = MGRGlobalSmootherType::none; | ||
| m_levelInterpType[1] = MGRInterpolationType::jacobi; | ||
| m_levelRestrictType[1] = MGRRestrictionType::injection; | ||
| m_levelCoarseGridMethod[1] = MGRCoarseGridMethod::nonGalerkin; |
There was a problem hiding this comment.
Focusing on the isothermal case, I would first try a 3 level strategy:
- Eliminate Lagrange multipliers as done in the LagrangianContactMechanics strategy
- Follow the same reduction strategy adopted in the SinglePhasePoromechanicsReservoirFVM strategy
This PR adds MGR strategy for
SinglePhasePoromechanicsConformingFracturesReservoir