QuasiDynamicEQ.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_QUASIDYNAMICEQ_HPP
#define GEOS_PHYSICSSOLVERS_INDUCED_QUASIDYNAMICEQ_HPP
 
#include "physicsSolvers/PhysicsSolverBase.hpp"
 
namespace geos
{
 
class QuasiDynamicEQ : public PhysicsSolverBase
{
public:
  QuasiDynamicEQ() = delete;
 
  QuasiDynamicEQ( const string & name,
                  Group * const parent );
 
  virtual ~QuasiDynamicEQ() override;
 
  static string catalogName() { return "QuasiDynamicEQ"; }
 
  virtual string getCatalogName() const override { return catalogName(); }
 
  virtual void registerDataOnMesh( Group & meshBodies ) override;
 
  struct viewKeyStruct : public PhysicsSolverBase::viewKeyStruct
  {
    constexpr static char const * stressSolverNameString() { return "stressSolverName"; }
    constexpr static char const * frictionLawNameString() { return "frictionLawName"; }
    constexpr static char const * shearImpedanceString() { return "shearImpedance"; }
    constexpr static char const * targetSlipIncrementString() { return "targetSlipIncrement"; }
  };
 
  virtual real64 solverStep( real64 const & time_n,
                             real64 const & dt,
                             integer const cycleNumber,
                             DomainPartition & domain ) override final;
 
  virtual real64 setNextDt( real64 const & currentDt,
                            DomainPartition & domain ) override final;
 
  real64 updateStresses( real64 const & time_n,
                         real64 const & dt,
                         const int cycleNumber,
                         DomainPartition & domain ) const;
 
  void saveOldStateAndUpdateSlip( ElementSubRegionBase & subRegion, real64 const dt ) const;
 
 
private:
 
 
 
  virtual void postInputInitialization() override;
 
 
 
  PhysicsSolverBase * m_stressSolver;
 
  string m_stressSolverName;
 
  real64 m_shearImpedance;
 
  real64 m_targetSlipIncrement;
 
  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;
  };
};
 
} /* namespace geos */
 
#endif /* GEOS_PHYSICSSOLVERS_INDUCED_QUASIDYNAMICEQ_HPP */