TractionBoundaryCondition.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_FIELDSPECIFICATION_TRACTIONBOUNDARYCONDITION_HPP
#define GEOS_FIELDSPECIFICATION_TRACTIONBOUNDARYCONDITION_HPP
 
#include "FieldSpecificationBase.hpp"
#include "mesh/FaceManager.hpp"
#include "mesh/NodeManager.hpp"
 
namespace geos
{
 
class TableFunction;
 
class TractionBoundaryCondition : public FieldSpecificationBase
{
public:
  TractionBoundaryCondition( string const & name, Group * parent );
 
  TractionBoundaryCondition() = delete;
 
  virtual ~TractionBoundaryCondition() = default;
 
  TractionBoundaryCondition( TractionBoundaryCondition const & ) = delete;
 
  TractionBoundaryCondition( TractionBoundaryCondition && ) = default;
 
  TractionBoundaryCondition & operator=( TractionBoundaryCondition const & ) = delete;
 
  TractionBoundaryCondition & operator=( TractionBoundaryCondition && ) = delete;
 
 
  static string catalogName() { return "Traction"; }
 
 
  void launch( real64 const time,
               arrayView1d< globalIndex const > const blockLocalDofNumber,
               globalIndex const dofRankOffset,
               FaceManager const & faceManager,
               arrayView2d< real64 const, nodes::REFERENCE_POSITION_USD > const nodePositions,
               SortedArrayView< localIndex const > const & targetSet,
               arrayView1d< real64 > const & localRhs ) const;
 
  void reinitScaleSet( FaceManager const & faceManager,
                       SortedArrayView< localIndex const > const & targetSet,
                       arrayView1d< real64 const > const nodalScaleSet );
 
  struct viewKeyStruct : public FieldSpecificationBase::viewKeyStruct
  {
    constexpr static char const * tractionTypeString() { return "tractionType"; }
 
    constexpr static char const * inputStressString() { return "inputStress"; }
 
//    /// @return The key for the function describing the components of stress.
//    constexpr static char const * stressFunctionString() { return "stressFunctions"; }
 
    constexpr static char const * scaleSetString() { return "scaleSet"; }
 
    constexpr static char const * nodalScaleFlagString() { return "nodalScaleFlag"; }
 
  };
 
  enum class TractionType : integer
  {
    vector, 
    normal, 
    stress  
  };
 
protected:
 
  virtual void postInputInitialization() override final;
 
  virtual void initializePreSubGroups() override final;
 
  TractionType m_tractionType;
 
  R2SymTensor m_inputStress;
 
  array1d< real64 > m_scaleSet;
 
  integer m_nodalScaleFlag;
 
//  /// names of the functions used to specify stress for the generation of tractions.
//  string_array m_stressFunctionNames;
//
//  bool m_useStressFunctions;
//
//  TableFunction const * m_stressFunctions[6]; // if this line is re-enabled, ensure that those pointers are initialized
 
};
 
ENUM_STRINGS( TractionBoundaryCondition::TractionType,
              "vector",
              "normal",
              "stress" );
 
} /* namespace geos */
 
#endif /* GEOS_FIELDSPECIFICATION_TRACTIONBOUNDARYCONDITION_HPP */