feat: feature/byer3/thermal well #3156
There are no files selected for viewing
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,131 @@ | ||
| /* | ||
| * ------------------------------------------------------------------------------------------------------------ | ||
| * SPDX-License-Identifier: LGPL-2.1-only | ||
| * | ||
| * Copyright (c) 2018-2020 Lawrence Livermore National Security LLC | ||
| * Copyright (c) 2018-2020 The Board of Trustees of the Leland Stanford Junior University | ||
| * Copyright (c) 2018-2020 TotalEnergies | ||
| * Copyright (c) 2019- GEOSX Contributors | ||
| * All rights reserved | ||
| * | ||
| * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details. | ||
| * ------------------------------------------------------------------------------------------------------------ | ||
| */ | ||
|
|
||
| /** | ||
| * @file ThermalCompositionalMultiphaseFVM.hpp | ||
| */ | ||
|
|
||
| #ifndef GEOS_LINEARALGEBRA_INTERFACES_HYPREMGRTHERMALCOMPOSITIONALMULTIPHASERESERVOIRFVM_HPP_ | ||
| #define GEOS_LINEARALGEBRA_INTERFACES_HYPREMGRTHERMALCOMPOSITIONALMULTIPHASERESERVOIRFVM_HPP_ | ||
|
|
||
| #include "linearAlgebra/interfaces/hypre/HypreMGR.hpp" | ||
|
|
||
| namespace geos | ||
| { | ||
|
|
||
| namespace hypre | ||
| { | ||
|
|
||
| namespace mgr | ||
| { | ||
|
|
||
| /** | ||
| * @brief ThermalCompositionalMultiphaseReservoirFVM strategy. | ||
| * | ||
| * Labels description stored in point_marker_array | ||
| * 0 = reservoir pressure | ||
| * 1 = reservoir density (component 1 ) | ||
| * ... = reservoir densities (# components , NC) | ||
| * NC + 1 = reservoir temperature | ||
| * | ||
| * numResLabels - 1 = reservoir temperature | ||
| * numResLabels = well pressure | ||
| * numResLabels + 1 = well density | ||
| * ... = ... (well densities) | ||
| * numResLabels + numWellLabels - 3 = last well density | ||
| * numResLabels + numWellLabels - 2 = well rate | ||
| * numResLabels + numWellLabels - 1 = well temperature | ||
| * | ||
| * 3-level MGR reduction strategy | ||
| * - 1st level: eliminate the well block | ||
| * - 2nd level: eliminate the reservoir density associated with the volume constraint | ||
| * - 3rd level: eliminate the remaining the reservoir densities | ||
| * - The coarse grid (pressure and temperature system) is solved with BoomerAMG with numFunctions==2. | ||
| * | ||
| */ | ||
|
|
||
| class ThermalCompositionalMultiphaseReservoirFVM : public MGRStrategyBase< 3 > | ||
| { | ||
| public: | ||
| /** | ||
| * @brief Constructor. | ||
| * @param numComponentsPerField array with number of components for each field | ||
| */ | ||
| explicit ThermalCompositionalMultiphaseReservoirFVM( arrayView1d< int const > const & numComponentsPerField ) | ||
tjb-ltk marked this conversation as resolved.
Show resolved
Hide resolved
|
||
| : MGRStrategyBase( LvArray::integerConversion< HYPRE_Int >( numComponentsPerField[0] + numComponentsPerField[1] ) ) | ||
|
Check warning on line 66 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
| { | ||
| HYPRE_Int const numResLabels = LvArray::integerConversion< HYPRE_Int >( numComponentsPerField[0] ); | ||
|
Check warning on line 68 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| // Level 0: eliminate the well block | ||
| m_labels[0].resize( numResLabels ); | ||
| std::iota( m_labels[0].begin(), m_labels[0].end(), 0 ); | ||
|
Check warning on line 72 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| // Level 1: eliminate last density which corresponds to the volume constraint equation | ||
| m_labels[1].resize( numResLabels - 2 ); | ||
| std::iota( m_labels[1].begin(), m_labels[1].end(), 0 ); | ||
| m_labels[1].push_back( numResLabels-1 ); // keep temperature | ||
|
Check warning on line 77 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
| // Level 2: eliminate the other densities | ||
| m_labels[2].push_back( 0 ); // keep pressure | ||
| m_labels[2].push_back( numResLabels-1 ); // keep temperature | ||
|
Check warning on line 80 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| setupLabels(); | ||
|
Check warning on line 82 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| // level 0 | ||
| m_levelFRelaxType[0] = MGRFRelaxationType::gsElimWInverse; | ||
| m_levelFRelaxIters[0] = 1; | ||
| m_levelInterpType[0] = MGRInterpolationType::blockJacobi; | ||
| m_levelRestrictType[0] = MGRRestrictionType::injection; | ||
| m_levelCoarseGridMethod[0] = MGRCoarseGridMethod::galerkin; | ||
| m_levelGlobalSmootherType[0] = MGRGlobalSmootherType::none; | ||
|
Check warning on line 90 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| m_levelFRelaxType[1] = MGRFRelaxationType::jacobi; | ||
| m_levelFRelaxIters[1] = 1; | ||
| m_levelInterpType[1] = MGRInterpolationType::jacobi; // Diagonal scaling (Jacobi) | ||
| m_levelRestrictType[1] = MGRRestrictionType::injection; | ||
| m_levelCoarseGridMethod[1] = MGRCoarseGridMethod::galerkin; // Standard Galerkin | ||
| m_levelGlobalSmootherType[1] = MGRGlobalSmootherType::blockGaussSeidel; | ||
| m_levelGlobalSmootherIters[1] = 1; | ||
|
Check warning on line 98 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| m_levelFRelaxType[2] = MGRFRelaxationType::jacobi; | ||
| m_levelFRelaxIters[2] = 1; | ||
| m_levelInterpType[2] = MGRInterpolationType::injection; // Injection | ||
| m_levelRestrictType[2] = MGRRestrictionType::injection; | ||
| m_levelCoarseGridMethod[2] = MGRCoarseGridMethod::cprLikeBlockDiag; // Non-Galerkin Quasi-IMPES CPR | ||
| m_levelGlobalSmootherType[2] = MGRGlobalSmootherType::ilu0; | ||
| m_levelGlobalSmootherIters[2] = 1; | ||
| } | ||
|
Check warning on line 107 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| /** | ||
| * @brief Setup the MGR strategy. | ||
| * @param precond preconditioner wrapper | ||
| * @param mgrData auxiliary MGR data | ||
| */ | ||
| void setup( LinearSolverParameters::MGR const &, | ||
|
Check warning on line 114 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
| HyprePrecWrapper & precond, | ||
| HypreMGRData & mgrData ) | ||
| { | ||
| setReduction( precond, mgrData ); | ||
|
Check warning on line 118 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
|
|
||
| // Configure the BoomerAMG solver used as mgr coarse solver for the pressure/temperature reduced system | ||
| setPressureTemperatureAMG( mgrData.coarseSolver ); | ||
| } | ||
|
Check warning on line 122 in src/coreComponents/linearAlgebra/interfaces/hypre/mgrStrategies/ThermalCompositionalMultiphaseReservoirFVM.hpp
|
||
| }; | ||
|
|
||
| } // namespace mgr | ||
|
|
||
| } // namespace hypre | ||
|
|
||
| } // namespace geos | ||
|
|
||
| #endif /*GEOS_LINEARALGEBRA_INTERFACES_HYPREMGRCOMPOSITIONALMULTIPHASERESERVOIRFVM_HPP_*/ | ||
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
do you need to add one with
trueinThermalCompressibleSinglePhaseFluid?There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@tjb-ltk Is this resolved?