A radial frac case cannot be running with ElasticTransverseIsotropic
Open
@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
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Assessment
This issue has not been assessed yet.