SinglePhasePoromechanicsFractures.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_MULTIPHYSICS_POROMECHANICSKERNELS_SINGLEPHASEPOROMECHANICSFRACTURES_HPP_
#define GEOS_PHYSICSSOLVERS_MULTIPHYSICS_POROMECHANICSKERNELS_SINGLEPHASEPOROMECHANICSFRACTURES_HPP_
#include "physicsSolvers/solidMechanics/contact/FractureState.hpp"
 
namespace geos
{
 
namespace poromechanicsFracturesKernels
{
 
struct StateUpdateKernel
{
 
  template< typename POLICY, typename POROUS_WRAPPER, typename HYDRAULIC_APERTURE_WRAPPER >
  static void
  launch( localIndex const size,
          POROUS_WRAPPER const & porousMaterialWrapper,
          HYDRAULIC_APERTURE_WRAPPER const & contactWrapper,
          arrayView2d< real64 const > const & dispJump,
          arrayView1d< real64 const > const & pressure,
          arrayView1d< real64 const > const & area,
          arrayView1d< real64 const > const & volume,
          arrayView1d< real64 > const & deltaVolume,
          arrayView1d< real64 > const & aperture,
          arrayView1d< real64 const > const & oldHydraulicAperture,
          arrayView1d< real64 > const & hydraulicAperture,
          arrayView2d< real64 const > const & fractureTraction,
          arrayView1d< integer > const & fractureState )
  {
 
    forAll< POLICY >( size, [=] GEOS_HOST_DEVICE ( localIndex const k )
    {
      // update aperture to be equal to the normal displacement jump
      aperture[k] = dispJump[k][0]; // the first component of the jump is the normal one.
 
      real64 dHydraulicAperture_dNormalJump = 0.0;
      real64 dHydraulicAperture_dNormalTraction = 0.0;
      hydraulicAperture[k] = contactWrapper.computeHydraulicAperture( aperture[k],
                                                                      fractureTraction[k][0],
                                                                      fractureState[k],
                                                                      dHydraulicAperture_dNormalJump,
                                                                      dHydraulicAperture_dNormalTraction );
 
      deltaVolume[k] = hydraulicAperture[k] * area[k] - volume[k];
 
      real64 const jump[3] = LVARRAY_TENSOROPS_INIT_LOCAL_3 ( dispJump[k] );
      real64 const traction[3] = LVARRAY_TENSOROPS_INIT_LOCAL_3 ( fractureTraction[k] );
 
      porousMaterialWrapper.updateStateFromPressureApertureJumpAndTraction( k, 0, pressure[k],
                                                                            oldHydraulicAperture[k],
                                                                            hydraulicAperture[k],
                                                                            dHydraulicAperture_dNormalJump,
                                                                            jump,
                                                                            traction );
 
    } );
  }
 
};
 
} // namespace poromechanicsFracturesKernels
 
} /* namespace geos */
 
#endif // GEOS_PHYSICSSOLVERS_MULTIPHYSICS_POROMECHANICSKERNELS_SINGLEPHASEPOROMECHANICSFRACTURES_HPP_