SpringSlider.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_INDUCEDSEISMICITY_SPRINGSLIDER_HPP
#define GEOS_PHYSICSSOLVERS_INDUCEDSEISMICITY_SPRINGSLIDER_HPP
 
#include "physicsSolvers/inducedSeismicity/ImplicitQDRateAndState.hpp"
#include "physicsSolvers/solidMechanics/contact/ContactFields.hpp"
#include "physicsSolvers/inducedSeismicity/rateAndStateFields.hpp"
 
namespace geos
{
 
template< typename RSSOLVER_TYPE = ImplicitQDRateAndState >
class SpringSlider : public RSSOLVER_TYPE
{
public:
 
  SpringSlider() = delete;
 
  SpringSlider( const string & name,
                dataRepository::Group * const parent );
 
  virtual ~SpringSlider() override;
 
  static string catalogName() { return RSSOLVER_TYPE::derivedSolverPrefix() + "SpringSlider"; }
 
  virtual string getCatalogName() const override { return catalogName(); }
 
  virtual void registerDataOnMesh( dataRepository::Group & meshBodies ) override;
 
  struct viewKeyStruct : public RSSOLVER_TYPE::viewKeyStruct
  {};
 
  virtual real64 updateStresses( real64 const & time_n,
                                 real64 const & dt,
                                 const int cycleNumber,
                                 DomainPartition & domain ) const override final;
 
  template< typename FRICTION_TYPE >
  void updateShearTraction( SurfaceElementSubRegion & subRegion,
                            FRICTION_TYPE & frictionLaw,
                            real64 const & dt ) const;
 
private:
 
  class SpringSliderParameters
  {
public:
 
    GEOS_HOST_DEVICE
    SpringSliderParameters( real64 const normalTraction, real64 const a, real64 const b, real64 const Dc ):
      tauRate( 1e-4 ),
      springStiffness( 0.0 )
    {
      real64 const criticalStiffness = normalTraction * (b - a) / Dc;
      springStiffness = 0.9 * criticalStiffness;
    }
 
    SpringSliderParameters( SpringSliderParameters const & ) = default;
 
    SpringSliderParameters( SpringSliderParameters && ) = default;
 
    SpringSliderParameters() = delete;
 
    SpringSliderParameters & operator=( SpringSliderParameters const & ) = delete;
 
    SpringSliderParameters & operator=( SpringSliderParameters && ) =  delete;
 
    real64 tauRate;
 
    real64 springStiffness;
  };
};
 
template< typename RSSOLVER_TYPE >
template< typename FRICTION_TYPE >
void SpringSlider< RSSOLVER_TYPE >::updateShearTraction( SurfaceElementSubRegion & subRegion,
                                                         FRICTION_TYPE & frictionLaw,
                                                         real64 const & dt ) const
{
  arrayView2d< real64 const > const deltaSlip = subRegion.getField< fields::contact::deltaSlip >();
  arrayView2d< real64 > const shearTraction   = subRegion.getField< fields::rateAndState::shearTraction >();
  arrayView2d< real64 > const shearTraction_n = subRegion.getField< fields::rateAndState::shearTraction_n >();
 
  arrayView1d< real64 > const normalTraction  = subRegion.getField< fields::rateAndState::normalTraction >();
  forAll< parallelDevicePolicy<> >( subRegion.size(), [=] GEOS_DEVICE ( localIndex const k )
  {
    SpringSliderParameters springSliderParameters = SpringSliderParameters( normalTraction[k],
                                                                            frictionLaw.getACoefficient( k ),
                                                                            frictionLaw.getBCoefficient( k ),
                                                                            frictionLaw.getDcCoefficient( k ) );
 
 
 
    shearTraction[k][0] = shearTraction_n[k][0] + springSliderParameters.tauRate * dt
                          - springSliderParameters.springStiffness * deltaSlip[k][0];
    shearTraction[k][1] = shearTraction_n[k][1] + springSliderParameters.tauRate * dt
                          - springSliderParameters.springStiffness * deltaSlip[k][1];
  } );
}
 
} /* namespace geos */
 
#endif /* GEOS_PHYSICSSOLVERS_INDUCEDSEISMICITY_SPRINGSLIDER_HPP */