Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

feat: Make the rock toughness scalable with the fracture size and test anisotropic toughness #3310

Merged
merged 15 commits into from
Jan 21, 2025
Merged

feat: Make the rock toughness scalable with the fracture size and test anisotropic toughness #3310

Show file tree
Hide file tree
Changes from 4 commits
Commits
Show all changes
15 commits
Select commit Hold shift + click to select a range
28849b0
make the toughness scalable with the fracture size
frankfeifan Aug 27, 2024
d5049c3
allow to plot K_IC
frankfeifan Aug 28, 2024
952b9d1
fixed xmls
frankfeifan Aug 28, 2024
4e174f1
Merge remote-tracking branch 'origin/develop' into feat/frankfei/scal…
frankfeifan Aug 28, 2024
2281f88
Merge remote-tracking branch 'origin/develop' into feat/frankfei/scal…
frankfeifan Sep 16, 2024
75a2330
changed baseRockToughness to initialRockToughness and removed unneces…
frankfeifan Sep 17, 2024
af20655
added an example for anisotropic toughness
frankfeifan Sep 17, 2024
79f5c28
Merge remote-tracking branch 'origin/develop' into feat/frankfei/scal…
frankfeifan Jan 11, 2025
6f88d62
added an example with anisotropic size dependent toughness
frankfeifan Jan 11, 2025
4b4d2c8
Merge remote-tracking branch 'origin/develop' into feat/frankfei/scal…
frankfeifan Jan 13, 2025
6d988e3
removed unused faceManager
frankfeifan Jan 13, 2025
35af17b
uncrustify
frankfeifan Jan 13, 2025
9a0bd40
rebaseline.
CusiniM Jan 21, 2025
f7b5045
Merge branch 'develop' into feat/frankfei/scalingToughness
CusiniM Jan 21, 2025
1eb3f94
rebaseline again.
CusiniM Jan 21, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion examples/pygeosxExamples/modifyBoundaryCondition/pkn_example.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -62,7 +62,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Domain }"
rockToughness="1e6"
baseRockToughness="1e6"
nodeBasedSIF="1"
mpiCommOrder="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/almContactMechanics/ALM_stickFault_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@
logLevel="0"
fractureRegion="Fracture"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/Sneddon_hydroFrac_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,7 +37,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="10.0e6"
baseRockToughness="10.0e6"
mpiCommOrder="1"/>
<!-- SPHINX_SNEDDON_SOLVER_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/heterogeneousInSitu_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -53,7 +53,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Domain }"
rockToughness="$K_upscaled$"
baseRockToughness="$K_upscaled$"
nodeBasedSIF="1"
mpiCommOrder="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/heterogeneousInSitu_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -71,7 +71,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Domain }"
rockToughness="2.236068e6"
baseRockToughness="2.236068e6"
nodeBasedSIF="1"
mpiCommOrder="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/hydrofractureSinglePhase2d.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -34,7 +34,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Domain }"
rockToughness="1.0e6"/>
baseRockToughness="1.0e6"/>
</Solvers>

<Mesh>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/kgdEdgeBased_C3D6_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +49,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Domain }"
rockToughness="0.707e7"
baseRockToughness="0.707e7"
mpiCommOrder="1"/>
</Solvers>
</Problem>
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/kgdNodeBased_C3D6_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Domain }"
rockToughness="0.707e7"
baseRockToughness="0.707e7"
nodeBasedSIF="1"
mpiCommOrder="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/kgdToughnessDominated_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="1e6"
baseRockToughness="1e6"
mpiCommOrder="1"/>
<!-- Sphinx_Solvers_SurfaceGenerator_End -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/kgdToughnessDominated_poroelastic_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="1e6"
baseRockToughness="1e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/kgdValidation_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@
logLevel="1"
targetRegions="{ Domain }"
nodeBasedSIF="0"
rockToughness="1.2e6"
baseRockToughness="1.2e6"
mpiCommOrder="1"/>
<!-- SPHINX_SURFACEGENERATOR_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/kgdViscosityDominated_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="1e4"
baseRockToughness="1e4"
mpiCommOrder="1"/>
<!-- Sphinx_Solvers_SurfaceGenerator_End -->

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/kgdViscosityDominated_poroelastic_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="1e4"
baseRockToughness="1e4"
mpiCommOrder="1"
isPoroelastic="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedToughnessDominated_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -67,7 +67,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="3.0e6"
baseRockToughness="3.0e6"
mpiCommOrder="1"/>
<!-- SPHINX_SURFACEGENERATOR_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedToughnessDominated_poroelastic_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -68,7 +68,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="3.0e6"
baseRockToughness="3.0e6"
mpiCommOrder="1"
isPoroelastic="1"/>
<!-- SPHINX_SURFACEGENERATOR_END -->
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedToughnessDominated_poroelastic_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="3.0e6"
baseRockToughness="3.0e6"
mpiCommOrder="1"
isPoroelastic="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedToughnessDominated_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,7 +43,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="3.0e6"
baseRockToughness="3.0e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedViscosityDominated_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -68,7 +68,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="0.3e6"
baseRockToughness="0.3e6"
mpiCommOrder="1"/>
<!-- SPHINX_SURFACEGENERATOR_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedViscosityDominated_poroelastic_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -69,7 +69,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="0.3e6"
baseRockToughness="0.3e6"
mpiCommOrder="1"/>
<!-- SPHINX_SURFACEGENERATOR_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedViscosityDominated_poroelastic_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="0.3e6"
baseRockToughness="0.3e6"
mpiCommOrder="1"
isPoroelastic="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pennyShapedViscosityDominated_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,7 +43,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="0.3e6"
baseRockToughness="0.3e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pknViscosityDominated_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -66,7 +66,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="0.1e6"
baseRockToughness="0.1e6"
mpiCommOrder="1"/>
<!-- SPHINX_SURFACEGENERATOR_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pknViscosityDominated_poroelastic_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -66,7 +66,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="0.1e6"
baseRockToughness="0.1e6"
mpiCommOrder="1"/>
<!-- SPHINX_SURFACEGENERATOR_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pknViscosityDominated_poroelastic_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="0.1e6"
baseRockToughness="0.1e6"
mpiCommOrder="1"
isPoroelastic="1"/>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/pknViscosityDominated_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,7 +43,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
nodeBasedSIF="1"
rockToughness="1e4"
baseRockToughness="1e4"
mpiCommOrder="1"/>
</Solvers>

Expand Down
6 changes: 3 additions & 3 deletions inputFiles/hydraulicFracturing/scripts/hydrofractureFigure.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -36,7 +36,7 @@ def getParametersFromXML(xmlFilePath):
poissonRatio = float(elasticParam.get('defaultPoissonRatio'))
injectionRate = -2.0 * float(tree.find('FieldSpecifications/SourceFlux').get('scale')) / fluidDensity
tree = ElementTree.parse(prefix + xmlFilePath + "_base.xml")
toughness = float(tree.find('Solvers/SurfaceGenerator').get('rockToughness'))
toughness = float(tree.find('Solvers/SurfaceGenerator').get('baseRockToughness'))

elif xmlFilePath == 'pennyShapedToughnessDominated' or xmlFilePath == 'pennyShapedViscosityDominated':
K = float(elasticParam.get('defaultBulkModulus'))
Expand All @@ -45,7 +45,7 @@ def getParametersFromXML(xmlFilePath):
poissonRatio = youngModulus / (2.0 * G) - 1.0
injectionRate = -4.0 * float(tree.find('FieldSpecifications/SourceFlux').get('scale')) / fluidDensity
tree = ElementTree.parse(prefix + xmlFilePath + "_benchmark.xml")
toughness = float(tree.find('Solvers/SurfaceGenerator').get('rockToughness'))
toughness = float(tree.find('Solvers/SurfaceGenerator').get('baseRockToughness'))

elif xmlFilePath == 'pknViscosityDominated':
K = float(elasticParam.get('defaultBulkModulus'))
Expand All @@ -54,7 +54,7 @@ def getParametersFromXML(xmlFilePath):
poissonRatio = youngModulus / (2.0 * G) - 1.0
injectionRate = -4.0 * float(tree.find('FieldSpecifications/SourceFlux').get('scale')) / fluidDensity
tree = ElementTree.parse(prefix + xmlFilePath + "_benchmark.xml")
toughness = float(tree.find('Solvers/SurfaceGenerator').get('rockToughness'))
toughness = float(tree.find('Solvers/SurfaceGenerator').get('baseRockToughness'))
found_core = False
for elem in param:
if elem.get("name") == "core":
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/walshQuarterNoChombo_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,7 +79,7 @@
fractureRegion="Fracture"
targetRegions="{Domain, void}"
nodeBasedSIF="1"
rockToughness="3.0e6">
baseRockToughness="3.0e6">
</SurfaceGenerator>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/hydraulicFracturing/walshQuarterNoChombo_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -64,7 +64,7 @@
fractureRegion="Fracture"
targetRegions="{Domain, void}"
nodeBasedSIF="1"
rockToughness="3.0e6"
baseRockToughness="3.0e6"
mpiCommOrder="1">
</SurfaceGenerator>
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/ContactMechanics_PassingCrack_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/ContactMechanics_SimpleCubes_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/ContactMechanics_SingleFracCompression_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,7 +30,7 @@
logLevel="0"
fractureRegion="Fracture"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>
<!-- SPHINX_SOLVER_END -->
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/ContactMechanics_TFrac_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -34,7 +34,7 @@
logLevel="0"
fractureRegion="Fracture"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
<!-- SPHINX_SOLVER_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/ContactMechanics_TFrac_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -34,7 +34,7 @@
logLevel="0"
fractureRegion="Fracture"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
<!-- SPHINX_SOLVER_END -->
</Solvers>
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/ContactMechanics_UnstructuredCrack_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/ContactMechanics_slippingFault_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,7 +30,7 @@
logLevel="0"
fractureRegion="Fracture"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/Sneddon_contactMechanics_benchmark.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,7 @@
logLevel="0"
fractureRegion="Fracture"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>
<!-- SPHINX_SNEDDON_SOLVER_END -->
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/lagrangianContactMechanics/Sneddon_contactMechanics_smoke.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,7 @@
logLevel="0"
fractureRegion="Fracture"
targetRegions="{ Region }"
rockToughness="1.0e6"
baseRockToughness="1.0e6"
mpiCommOrder="1"/>
</Solvers>
<!-- SPHINX_SNEDDON_SOLVER_END -->
Expand Down
2 changes: 1 addition & 1 deletion inputFiles/phaseField/PhaseFieldFracture_Nucleation_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -52,7 +52,7 @@
<SurfaceGenerator
name="SurfaceGen"
targetRegions="{ Region1 }"
rockToughness="0"
baseRockToughness="0"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion inputFiles/poromechanicsFractures/ExponentialDecayPermeability_conformingFracture_base.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
fractureRegion="Fracture"
rockToughness="1e6"
baseRockToughness="1e6"
mpiCommOrder="1"/>
</Solvers>

Expand Down
2 changes: 1 addition & 1 deletion ...poromechanicsFractures/PoroElastic_conformingFracture_2d_faultSlip_sequential_solvers.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -53,7 +53,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
fractureRegion="Fracture"
rockToughness="1e6"
baseRockToughness="1e6"
mpiCommOrder="1"/>
</Solvers>
</Problem>
2 changes: 1 addition & 1 deletion inputFiles/poromechanicsFractures/PoroElastic_conformingFracture_2d_faultSlip_solvers.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
fractureRegion="Fracture"
rockToughness="1e6"
baseRockToughness="1e6"
mpiCommOrder="1"/>
</Solvers>
</Problem>
2 changes: 1 addition & 1 deletion ...romechanicsFractures/PoroElastic_conformingFracture_2d_openingFrac_sequential_solvers.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -52,7 +52,7 @@
name="SurfaceGen"
targetRegions="{ Domain }"
fractureRegion="Fracture"
rockToughness="1e6"
baseRockToughness="1e6"
mpiCommOrder="1"/>
</Solvers>
</Problem>
Loading
Loading