CompositionalMultiphaseWell.hpp

Go to the documentation of this file.

/*
 * ------------------------------------------------------------------------------------------------------------
 * SPDX-License-Identifier: LGPL-2.1-only
 *
 * Copyright (c) 2016-2024 Lawrence Livermore National Security LLC
 * Copyright (c) 2018-2024 TotalEnergies
 * Copyright (c) 2018-2024 The Board of Trustees of the Leland Stanford Junior University
 * Copyright (c) 2023-2024 Chevron
 * Copyright (c) 2019-     GEOS/GEOSX Contributors
 * All rights reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */
 
#ifndef GEOS_PHYSICSSOLVERS_FLUIDFLOW_WELLS_COMPOSITIONALMULTIPHASEWELL_HPP_
#define GEOS_PHYSICSSOLVERS_FLUIDFLOW_WELLS_COMPOSITIONALMULTIPHASEWELL_HPP_
 
#include "constitutive/fluid/multifluid/Layouts.hpp"
#include "constitutive/relativePermeability/Layouts.hpp"
#include "physicsSolvers/fluidFlow/wells/WellSolverBase.hpp"
#include "physicsSolvers/fluidFlow/CompositionalMultiphaseBase.hpp"
 
#include "physicsSolvers/fluidFlow/wells/WellConstraintsBase.hpp"
 
namespace geos
{
 
namespace constitutive
{
class ConstitutiveManager;
class MultiFluidBase;
}
 
class CompositionalMultiphaseWell : public WellSolverBase
{
public:
 
  CompositionalMultiphaseWell( const string & name,
                               Group * const parent );
 
  CompositionalMultiphaseWell() = delete;
 
  CompositionalMultiphaseWell( CompositionalMultiphaseWell const & ) = delete;
 
  CompositionalMultiphaseWell( CompositionalMultiphaseWell && ) = default;
 
  CompositionalMultiphaseWell & operator=( CompositionalMultiphaseWell const & ) = delete;
 
  CompositionalMultiphaseWell & operator=( CompositionalMultiphaseWell && ) = delete;
 
  virtual ~CompositionalMultiphaseWell() override = default;
 
  static string catalogName() { return "CompositionalMultiphaseWell"; }
  string getCatalogName() const override { return catalogName(); }
 
  virtual void registerDataOnMesh( Group & meshBodies ) override;
 
 
  virtual real64
  calculateWellResidualNorm( real64 const & time_n,
                             real64 const & dt,
                             WellElementSubRegion const & subRegion,
                             DofManager const & dofManager,
                             arrayView1d< real64 const > const & localRhs ) override;
 
  virtual real64
  calculateResidualNorm( real64 const & time_n,
                         real64 const & dt,
                         DomainPartition const & domain,
                         DofManager const & dofManager,
                         arrayView1d< real64 const > const & localRhs ) override;
 
  virtual real64
  scalingForWellSystemSolution( ElementSubRegionBase & subRegion,
                                DofManager const & dofManager,
                                arrayView1d< real64 const > const & localSolution ) override;
 
  virtual real64
  scalingForSystemSolution( DomainPartition & domain,
                            DofManager const & dofManager,
                            arrayView1d< real64 const > const & localSolution ) override;
 
  virtual bool
  checkWellSystemSolution( ElementSubRegionBase & subRegion,
                           DofManager const & dofManager,
                           arrayView1d< real64 const > const & localSolution,
                           real64 const scalingFactor ) override;
 
  virtual bool
  checkSystemSolution( DomainPartition & domain,
                       DofManager const & dofManager,
                       arrayView1d< real64 const > const & localSolution,
                       real64 const scalingFactor ) override;
 
  virtual void
  applyWellSystemSolution( DofManager const & dofManager,
                           arrayView1d< real64 const > const & localSolution,
                           real64 const scalingFactor,
                           real64 const dt,
                           DomainPartition & domain,
                           MeshLevel & mesh,
                           WellElementSubRegion & subRegion ) override;
 
  virtual void
  applySystemSolution( DofManager const & dofManager,
                       arrayView1d< real64 const > const & localSolution,
                       real64 const scalingFactor,
                       real64 const dt,
                       DomainPartition & domain ) override;
 
  virtual void
  resetStateToBeginningOfStep( DomainPartition & domain ) override;
 
  virtual void
  implicitStepSetup( real64 const & time,
                     real64 const & dt,
                     DomainPartition & domain ) override;
 
  virtual void
  implicitStepComplete( real64 const & time,
                        real64 const & dt,
                        DomainPartition & domain ) override;
 
  void updateGlobalComponentFraction( WellElementSubRegion & subRegion ) const;
 
  void updateVolRatesForConstraint( WellElementSubRegion & subRegion );
 
  void updateBHPForConstraint( WellElementSubRegion & subRegion );
 
  void updateFluidModel( WellElementSubRegion & subRegion );
 
  void updateSeparator( WellElementSubRegion & subRegion );
 
  void calculateReferenceElementRates( WellElementSubRegion & subRegion );
 
  real64 updatePhaseVolumeFraction( WellElementSubRegion & subRegion ) const;
 
  void updateTotalMassDensity( WellElementSubRegion & subRegion ) const;
 
  virtual void computeWellPerforationRates( real64 const & time_n,
                                            real64 const & GEOS_UNUSED_PARAM( dt ),
                                            ElementRegionManager const & elemManager,
                                            WellElementSubRegion & subRegion )override;
 
 
  virtual void computePerforationRates( real64 const & time_n,
                                        real64 const & dt, DomainPartition & domain ) override;
 
  virtual real64 updateWellState( WellElementSubRegion & subRegion ) override;
  virtual void updateState( DomainPartition & domain ) override;
 
  virtual real64 updateSubRegionState( WellElementSubRegion & subRegion ) override;
 
  virtual string wellElementDofName() const override { return viewKeyStruct::dofFieldString(); }
 
  virtual string resElementDofName() const override { return CompositionalMultiphaseBase::viewKeyStruct::elemDofFieldString(); }
 
  virtual localIndex numFluidComponents() const override { return m_numComponents; }
 
  virtual localIndex numFluidPhases() const override { return m_numPhases; }
 
  integer useTotalMassEquation() const { return m_useTotalMassEquation; }
 
  virtual void assembleWellFluxTerms( real64 const & time,
                                      real64 const & dt,
                                      WellElementSubRegion const & subRegion,
                                      DofManager const & dofManager,
                                      CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                      arrayView1d< real64 > const & localRhs ) override;
 
  virtual void assembleFluxTerms( real64 const & time_n,
                                  real64 const & dt,
                                  DomainPartition & domain,
                                  DofManager const & dofManager,
                                  CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                  arrayView1d< real64 > const & localRhs )override;
 
  virtual void assembleWellAccumulationTerms( real64 const & time,
                                              real64 const & dt,
                                              WellElementSubRegion & subRegion,
                                              DofManager const & dofManager,
                                              CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                              arrayView1d< real64 > const & localRhs ) override;
 
  virtual void assembleAccumulationTerms( real64 const & time_n,
                                          real64 const & dt,
                                          DomainPartition & domain,
                                          DofManager const & dofManager,
                                          CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                          arrayView1d< real64 > const & localRhs ) override;
 
  virtual void assembleWellConstraintTerms( real64 const & time_n,
                                            real64 const & dt,
                                            WellElementSubRegion const & subRegion,
                                            DofManager const & dofManager,
                                            CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                            arrayView1d< real64 > const & localRhs ) override;
 
  virtual void applyWellBoundaryConditions( real64 const time_n,
                                            real64 const dt,
                                            ElementRegionManager & elemManager,
                                            WellElementSubRegion & subRegion,
                                            DofManager const & dofManager,
                                            arrayView1d< real64 > const & localRhs,
                                            CRSMatrixView< real64, globalIndex const > const & localMatrix ) override;
 
  virtual void applyBoundaryConditions( real64 const GEOS_UNUSED_PARAM( time_n ),
                                        real64 const GEOS_UNUSED_PARAM( dt ),
                                        DomainPartition & GEOS_UNUSED_PARAM( domain ),
                                        DofManager const & GEOS_UNUSED_PARAM( dofManager ),
                                        CRSMatrixView< real64, globalIndex const > const & GEOS_UNUSED_PARAM( localMatrix ),
                                        arrayView1d< real64 > const & GEOS_UNUSED_PARAM( localRhs ) ) override;
 
 
    virtual void assembleWellPressureRelations( real64 const & time_n,
                                              real64 const & dt,
                                              WellElementSubRegion const & subRegion,
                                              DofManager const & dofManager,
                                              CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                              arrayView1d< real64 > const & localRhs ) override;
 
  virtual void assemblePressureRelations( real64 const & time_n,
                                          real64 const & dt,
                                          DomainPartition const & domain,
                                          DofManager const & dofManager,
                                          CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                          arrayView1d< real64 > const & localRhs ) override;
 
 
  void chopNegativeDensities( DomainPartition & domain );
 
  void chopNegativeDensities( WellElementSubRegion & subRegion );
 
  struct viewKeyStruct : WellSolverBase::viewKeyStruct
  {
    static constexpr char const * dofFieldString() { return "compositionalWellVars"; }
 
    // inputs
 
    static constexpr char const * useMassFlagString() { return CompositionalMultiphaseBase::viewKeyStruct::useMassFlagString(); }
 
    static constexpr char const * useTotalMassEquationString() { return CompositionalMultiphaseBase::viewKeyStruct::useTotalMassEquationString(); }
 
    static constexpr char const * maxCompFracChangeString() { return CompositionalMultiphaseBase::viewKeyStruct::maxCompFracChangeString(); }
 
    static constexpr char const * maxRelativePresChangeString() { return "maxRelativePressureChange"; }
 
    static constexpr char const * maxAbsolutePresChangeString() { return "maxAbsolutePressureChange"; }
 
    static constexpr char const * maxRelativeCompDensChangeString() { return "maxRelativeCompDensChange"; }
 
    static constexpr char const * maxRelativeTempChangeString() { return "maxRelativeTemperatureChange"; }
 
    static constexpr char const * allowLocalCompDensChoppingString() { return CompositionalMultiphaseBase::viewKeyStruct::allowLocalCompDensChoppingString(); }
 
    // control data (not registered on the mesh)
 
    static constexpr char const * massDensityString() { return "massDensity";}
 
    static constexpr char const * currentBHPString() { return "currentBHP"; }
    static constexpr char const * dCurrentBHPString() { return "dCurrentBHP"; }
 
    static constexpr char const * dCurrentBHP_dPresString() { return "dCurrentBHP_dPres"; }
    static constexpr char const * dCurrentBHP_dCompDensString() { return "dCurrentBHP_dCompDens"; }
 
    static constexpr char const * currentPhaseVolRateString() { return "currentPhaseVolumetricRate"; }
    static constexpr char const * dCurrentPhaseVolRateString() { return "dCurrentPhaseVolumetricRate"; }
 
 
    static constexpr char const * dCurrentPhaseVolRate_dPresString() { return "dCurrentPhaseVolumetricRate_dPres"; }
 
    static constexpr char const * dCurrentPhaseVolRate_dCompDensString() { return "dCurrentPhaseVolumetricRate_dCompDens"; }
 
    static constexpr char const * dCurrentPhaseVolRate_dRateString() { return "dCurrentPhaseVolumetricRate_dRate"; }
 
    static constexpr char const * currentTotalVolRateString() { return "currentTotalVolumetricRate"; }
    static constexpr char const * dCurrentTotalVolRateString() { return "dCurrentTotalVolumetricRate"; }
 
    static constexpr char const * currentMassRateString() { return "currentMassRate"; }
 
    static constexpr char const * dCurrentTotalVolRate_dPresString() { return "dCurrentTotalVolumetricRate_dPres"; }
 
    static constexpr char const * dCurrentTotalVolRate_dCompDensString() { return "dCurrentTotalVolumetricRate_dCompDens"; }
 
    static constexpr char const * dCurrentTotalVolRate_dRateString() { return "dCurrentTotalVolumetricRate_dRate"; }
 
  } viewKeysCompMultiphaseWell;
 
protected:
 
  virtual void postInputInitialization() override;
 
  virtual void initializePostSubGroups() override;
 
  virtual void initializePostInitialConditionsPreSubGroups() override;
 
  void saveState( WellElementSubRegion & subRegion );
  virtual void postRestartInitialization() override final;
  /*
   * @brief Utility function that checks the consistency of the constitutive models
   * @param[in] domain the domain partition
   * @detail
   * This function will produce an error if one of the well constitutive models
   * is incompatible with the corresponding models in reservoir
   * regions connected to that particular well.
   */
  void validateConstitutiveModels( DomainPartition const & domain ) const;
 
  void validateWellControlsForFluid( WellControls const & wellControls,
                                     constitutive::MultiFluidBase const & fluid ) const;
 
 
  virtual void validateWellConstraints( real64 const & time_n,
                                        real64 const & dt,
                                        WellElementSubRegion const & subRegion,
                                        ElementRegionManager const & elemManager ) override;
 
  void createSeparator();
 
  void printRates( real64 const & time_n,
                   real64 const & dt,
                   DomainPartition & domain ) override;
 
  virtual bool evaluateConstraints( real64 const & time_n,
                                    real64 const & stepDt,
                                    integer const cycleNumber,
                                    integer const coupledIterationNumber,
                                    DomainPartition & domain,
                                    MeshLevel & mesh,
                                    ElementRegionManager & elemManager,
                                    WellElementSubRegion & subRegion,
                                    DofManager const & dofManager ) override;
 
private:
 
  void initializeWells( DomainPartition & domain, real64 const & time_n ) override;
 
  void initializeWell( DomainPartition & domain, MeshLevel & mesh, WellElementSubRegion & subRegion, real64 const & time_n ) override;
 
  virtual void setConstitutiveNames( ElementSubRegionBase & subRegion ) const override;
 
  template< typename ... GROUPTYPES >
  void selectLimitingConstraint( real64 const & time_n, integer const coupledIterationNumber, WellElementSubRegion & subRegion );
 
  void solveConstraint(  WellConstraintBase * constraint,
                        real64 const & time_n,
                        real64 const & dt,
                        integer const cycleNumber,
                        integer const coupledIterationNumber,
                        DomainPartition & domain,
                        MeshLevel & mesh,
                        ElementRegionManager & elemManager,
                        WellElementSubRegion & subRegion,
                        DofManager const & dofManager );
 
 
  integer m_useMass;
 
  integer m_useTotalMassEquation;
 
  real64 m_maxCompFracChange;
 
  real64 m_maxRelativePresChange;
 
  real64 m_maxAbsolutePresChange;
 
  real64 m_maxRelativeCompDensChange;
 
  real64 m_maxRelativeTempChange;
 
  real64 m_minScalingFactor;
 
  integer m_allowCompDensChopping;
 
  localIndex m_targetPhaseIndex;
 
 
 
};
 
} // namespace geos
 
 
#endif //GEOS_PHYSICSSOLVERS_FLUIDFLOW_WELLS_COMPOSITIONALMULTIPHASEWELL_HPP_