GEOS-DEV / GEOS-DEV/GEOS

A radial frac case cannot be running with ElasticTransverseIsotropic

Open
#1,825 8 comments 0 reactions 2 assignees View on GitHub

@jhuang2601 is already working on this.

Since Mar 24, 2022.

type: bug
Dominant language
C++
Stars
287
Forks
109
Avg merge
4d 41m
Merged PRs (30d)
5

Description

Describe the bug
Cannot run pennyShapedViscosityDominated case with ElasticTransverseIsotropic and unexpected segmentation fault is reported.

To Reproduce

The xml files are attached: case.zip

Run pennyShapedViscosityDominated_TransverseIsotropic_smoke.xml with the develop repo to reproduce the reported issue

Screenshots
Below is the running log:

Max threads: 2
MKL max threads: 1
GEOSX version 0.2.0 (HEAD, sha1: 4312248)
Adding Solver of type Hydrofracture, named hydrofracture
Adding Solver of type SolidMechanicsLagrangianSSLE, named lagsolve
Adding Solver of type SinglePhaseFVM, named SinglePhaseFlow
Adding Solver of type SurfaceGenerator, named SurfaceGen
Adding Mesh: InternalMesh, mesh1
Adding Geometric Object: Box, fracture
Adding Geometric Object: Box, source
Adding Geometric Object: Box, core
Adding Event: SoloEvent, preFracture
Adding Event: PeriodicEvent, solverApplications1
Adding Event: PeriodicEvent, outputs
Adding Event: PeriodicEvent, restarts
   TableFunction: apertureTable
Adding Output: Silo, siloOutput
Adding Output: Restart, restartOutput
Adding Object CellElementRegion named Domain from ObjectManager::Catalog.
Adding Object SurfaceElementRegion named Fracture from ObjectManager::Catalog.
Total number of nodes:1372
Total number of elems:1014
Rank 0: Total number of nodes:1372
  Domain/cb1/water is allocated with 8 quadrature points.
  Domain/cb1/rock is allocated with 8 quadrature points.
  Fracture/faceElementSubRegion/water is allocated with 8 quadrature points.
  Fracture/faceElementSubRegion/rock is allocated with 8 quadrature points.
  Fracture/faceElementSubRegion/fractureFilling is allocated with 8 quadrature points.
  Fracture/faceElementSubRegion/fractureContact is allocated with 8 quadrature points.
Time: 0s, dt:1s, Cycle: 0
    Attempt:  0, NewtonIter:  0
    ( Rsolid, Rfluid ) = ( 0.00e+00, 4.42e+00 )
    ( R ) = ( 4.42e+00 ) ; 
    Attempt:  0, NewtonIter:  1
    ( Rsolid, Rfluid ) = ( 7.66e-16, 4.36e+02 )
    ( R ) = ( 4.36e+02 ) ; 
    Last LinSolve(iter,res) = (   1, 7.44e-08 ) ; 
        Line search @ 0.500:          ( Rsolid, Rfluid ) = ( 7.66e-16, 2.15e+02 )

        Line search @ 0.250:          ( Rsolid, Rfluid ) = ( 7.66e-16, 1.05e+02 )

        Line search @ 0.125:          ( Rsolid, Rfluid ) = ( 7.66e-16, 4.95e+01 )

        Line search @ 0.062:          ( Rsolid, Rfluid ) = ( 7.66e-16, 2.19e+01 )

        Line search failed to produce reduced residual. Accepting iteration.
    Attempt:  0, NewtonIter:  2
    ( Rsolid, Rfluid ) = ( 3.51e-15, 2.00e-04 )
    ( R ) = ( 2.00e-04 ) ; 
    Last LinSolve(iter,res) = (   1, 6.68e-10 ) ; 
    Attempt:  0, NewtonIter:  3
    ( Rsolid, Rfluid ) = ( 4.51e-16, 1.13e-07 )
    ( R ) = ( 1.13e-07 ) ; 
    Last LinSolve(iter,res) = (   1, 6.71e-10 ) ; 
Received signal 11: Segmentation fault

** StackTrace of 12 frames **
Frame 0: /lib64/libc.so.6 
Frame 1: geosx::SurfaceGenerator::calculateNodeAndFaceSif(geosx::DomainPartition const&, geosx::NodeManager&, geosx::EdgeManager const&, geosx::FaceManager&, geosx::ElementRegionManager const&) 
Frame 2: geosx::SurfaceGenerator::identifyRupturedFaces(geosx::DomainPartition const&, geosx::NodeManager&, geosx::EdgeManager&, geosx::FaceManager&, geosx::ElementRegionManager const&, bool) 
Frame 3: geosx::SurfaceGenerator::separationDriver(geosx::DomainPartition&, geosx::MeshLevel&, std::vector<geosx::NeighborCommunicator, std::allocator<geosx::NeighborCommunicator> >&, int, int, bool, double) 
Frame 4: geosx::SurfaceGenerator::solverStep(double const&, double const&, int, geosx::DomainPartition&) 
Frame 5: geosx::HydrofractureSolver::solverStep(double const&, double const&, int, geosx::DomainPartition&) 
Frame 6: geosx::SolverBase::execute(double, double, int, int, double, geosx::DomainPartition&) 
Frame 7: geosx::EventBase::execute(double, double, int, int, double, geosx::DomainPartition&) 
Frame 8: geosx::EventManager::run(geosx::DomainPartition&) 
Frame 9: geosx::GeosxState::run() 
Frame 10: main 
Frame 11: __libc_start_main 
Frame 12: /data/gpfs/Users/j0436735/travis-deployments/CPU/GEOSX-4312248/bin/geosx 
=====

Expected behavior
In addition to ElasticIsotropic, this hydrofrac case should be running with different solid models

Platform (please complete the following information):

  • Machine: Pecan
  • GEOSX Version: develop 4312248

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.