SeismicityRate.hpp

/*
 * ------------------------------------------------------------------------------------------------------------
 * 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_INDUCED_SEISMICITY_SEISMICITYRATE_HPP
#define GEOS_PHYSICSSOLVERS_INDUCED_SEISMICITY_SEISMICITYRATE_HPP
 
#include "physicsSolvers/PhysicsSolverBase.hpp"
 
namespace geos
{
 
class SeismicityRate : public PhysicsSolverBase
{
public:
  SeismicityRate() = delete;
 
  SeismicityRate( const string & name,
                  Group * const parent );
 
  virtual ~SeismicityRate() override;
 
  static string catalogName() { return "SeismicityRate"; }
 
  virtual string getCatalogName() const override { return catalogName(); }
 
  virtual void registerDataOnMesh( Group & meshBodies ) override;
 
  struct viewKeyStruct : public PhysicsSolverBase::viewKeyStruct
  {
    static constexpr char const * stressSolverNameString() { return "stressSolverName"; }
    static constexpr char const * initialFaultNormalTractionString() { return "initialFaultNormalTraction"; }
    static constexpr char const * initialFaultShearTractionString() { return "initialFaultShearTraction"; }
    static constexpr char const * faultNormalDirectionString() { return "faultNormalDirection"; }
    static constexpr char const * faultShearDirectionString() { return "faultShearDirection"; }
    static constexpr char const * directEffectString() { return "directEffect"; }
    static constexpr char const * backgroundStressingRateString() { return "backgroundStressingRate"; }
  };
 
  virtual real64 solverStep( real64 const & time_n,
                             real64 const & dt,
                             integer const cycleNumber,
                             DomainPartition & domain ) override final;
 
  void integralSolverStep( real64 const & time_n,
                           real64 const & dt,
                           ElementSubRegionBase & subRegion );
 
  void updateFaultTraction( ElementSubRegionBase & subRegion ) const;
 
  void saveOldState( ElementSubRegionBase & subRegion ) const;
 
  void computeTotalStressOnFault( arrayView1d< real64 const > const & biotCoefficient,
                                  arrayView1d< real64 const > const & pres,
                                  real64 const (&faultNormalProjectionTensor)[6],
                                  real64 const (&faultShearProjectionTensor)[6],
                                  arrayView1d< real64 > const & sig,
                                  arrayView1d< real64 > const & tau ) const;
 
  void initializeFaultTraction( real64 const time_n,
                                integer const cycleNumber,
                                DomainPartition & domain ) const;
 
 
protected:
 
  real64 updateStresses( real64 const & time_n,
                         real64 const & dt,
                         const int cycleNumber,
                         DomainPartition & domain ) const;
 
 
 
  void constructFaultStressProjectionTensors( real64 ( &faultNormalProjectionTensor )[6],
                                              real64 ( &faultShearProjectionTensor )[6] ) const;
 
  virtual void postInputInitialization() override;
 
  PhysicsSolverBase *m_stressSolver;
 
  string m_stressSolverName;
 
  R1Tensor m_faultNormalDirection;
 
  R1Tensor m_faultShearDirection;
 
  real64 m_directEffect;
 
  real64 m_backgroundStressingRate;
};
 
} /* namespace geos */
 
#endif /* GEOS_PHYSICSSOLVERS_INDUCED_SEISMICITY_SEISMICITYRATE_HPP */