feat: phase-field nucleation modeling of hydraulic fracturing #2125

frankfeifan · 2022-10-15T20:29:28Z

In this PR, we have added a phase-field poromechanics solver that can simulate the hydraulic fracturing nucleation and propagation. Detailed formulations can be seen in Fei et al. (2023). We have added and used following features in this PR:

On a new solver

A new phase-field poromechanics solver that enables sequential solve of the fully coupled poromechanics and phase-field damage evolution
Added a "PorousDamageSolid" where we can compute porosity and permeability specifically for porous solid with damage
Added a DamagePermeability to simply calculate the phase-field fracture permeability
Added a number of new examples for the phase-field poromechanics solver

On existing solvers and/or functions

Use the pressure gradient calculation in the SinglePhaseHybridFVM
Applied a scaling factor to the nodal traction magnitude so that we can degrade the traction on the surface when there is damage/fracture generated
Added an ImplicitSmallStrainQuasiStaticPressurizedDamage kernel for solid mechanics solver so that we can solve for pressurized fracture nucleation and propagation in a dry solid without using poromechanics (this is for calibration of some phase-field nucleation parameters) -> removed here, and can be added in another PR

…ction to output it

…tion function

… value non-negative

…to the Damage base code

… formulation is called

… the damage base code

…ation To keep the branch up-to-date

…ation

CusiniM

It's a pretty big PR and we may need to iterate a few times but I don't really see anything that should prevent us from merging this quickly.

src/coreComponents/schema/docs/PhaseFieldPoromechanics_other.rst

CusiniM · 2025-01-03T15:25:25Z

src/coreComponents/fieldSpecification/TractionBoundaryCondition.cpp

+      faceScale  += nodalScaleSet[ faceToNodeMap( kf, a ) ];
+    }
+
+    m_scaleSet[i] = std::min( 1.0, (1.0 - faceScale/numNodes)*(1.0 - faceScale/numNodes) );


Suggested change

m_scaleSet[i] = std::min( 1.0, (1.0 - faceScale/numNodes)*(1.0 - faceScale/numNodes) );

m_scaleSet[i] = LvArray::math::min( 1.0, (1.0 - faceScale/numNodes)*(1.0 - faceScale/numNodes) );

CusiniM · 2025-01-03T15:26:26Z

src/coreComponents/fieldSpecification/TractionBoundaryCondition.hpp

@@ -81,6 +81,10 @@ class TractionBoundaryCondition : public FieldSpecificationBase
               SortedArrayView< localIndex const > const & targetSet,
               arrayView1d< real64 > const & localRhs ) const;

+  void reinitScaleSet( FaceManager const & faceManager,


write doxygen.

write doxygen.

The doxygen is added

CusiniM · 2025-01-03T15:28:36Z

src/coreComponents/fieldSpecification/TractionBoundaryCondition.cpp

@@ -220,6 +234,33 @@ void TractionBoundaryCondition::launch( real64 const time,
  }
 }

+void TractionBoundaryCondition::reinitScaleSet( FaceManager const & faceManager,


I am not sure I understand what's happening in here and why you need it.

I am not sure I understand what's happening in here and why you need it.

I added this function to account for the reduction of traction on the boundary where phase-field damage is generated. Maybe this function can be used for other purposes when the traction magnitude depends on any nodal variables.

…efield_nucleation_poromech

...mponents/physicsSolvers/multiphysics/poromechanicsKernels/SinglePhasePoromechanicsDamage.hpp

…chanics kernel and updated copyright content

…efield_nucleation_poromech

…this branch

codecov · 2025-01-14T20:26:56Z

Codecov Report

Attention: Patch coverage is 18.28666% with 496 lines in your changes missing coverage. Please review.

Project coverage is 56.58%. Comparing base (fd36d2a) to head (9b3a5f0).
Report is 1 commits behind head on develop.

Files with missing lines	Patch %	Lines
...icsKernels/SinglePhasePoromechanicsDamage_impl.hpp	0.00%	93 Missing ⚠️
...omponents/constitutive/solid/PorousDamageSolid.hpp	1.16%	85 Missing ⚠️
...ers/multiphysics/PhaseFieldPoromechanicsSolver.cpp	0.00%	73 Missing ⚠️
src/coreComponents/constitutive/solid/Damage.hpp	24.65%	55 Missing ⚠️
...implePDE/PhaseFieldPressurizedDamageFEMKernels.hpp	0.00%	37 Missing ⚠️
...ers/multiphysics/PhaseFieldPoromechanicsSolver.hpp	6.45%	29 Missing ⚠️
...s/fieldSpecification/TractionBoundaryCondition.cpp	0.00%	26 Missing ⚠️
...s/constitutive/permeability/DamagePermeability.cpp	0.00%	18 Missing ⚠️
...s/physicsSolvers/simplePDE/PhaseFieldDamageFEM.cpp	0.00%	18 Missing ⚠️
...s/constitutive/permeability/DamagePermeability.hpp	5.88%	16 Missing ⚠️
... and 13 more

Additional details and impacted files

@@             Coverage Diff             @@
##           develop    #2125      +/-   ##
===========================================
- Coverage    56.79%   56.58%   -0.22%     
===========================================
  Files         1154     1163       +9     
  Lines       100074   100606     +532     
===========================================
+ Hits         56836    56923      +87     
- Misses       43238    43683     +445

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

frankfeifan and others added 30 commits April 5, 2022 11:43

add constitutive files for the phase-field nucleation formulation

456663b

add constitutive files for the phase-field nucleation formulation

5bed9c3

modify the pf formulation by adding an external driving force

0ee51b7

add a member variable for the external driving force and a member fun…

ec02d8c

…ction to output it

add the new constitutive for pf nucleation into the passthru

f2a52a1

Add DamageExtDrivingForce files into CMakeLists

f5a5487

solve the segmentation error

6f19f91

Add an input file for the phase-field nucleation formulation

d534626

correct the calculation of the external driving force

94aeed0

adding TableFunction files for the surfing boundary condition

8cde800

adding an option variable for the external driving force term

400ffcb

debug

5c69cf2

adding the initial fracture and surfing boundary condition

3cb9db3

set upper and lower bounds of the damage when calculating the degrada…

cea5b2b

…tion function

correct the typo in the formulation

16ffa58

clean the input file

bc085c8

set a threshold for the strain energy based on c_e to keep the damage…

b89db28

… value non-negative

clean the phase-field code and add the external driving force term in…

ebdfae2

…to the Damage base code

input files for the phase-field nucleation formulation

c0995ff

Add a switch of the external driving force to check if the nucleation…

dda9988

… formulation is called

clean the phase-field code by merging the nucleation formulation into…

c26c70f

… the damage base code

create the input file and mesh for the three point test

333b5a1

Manually delete the unneeded input file

0480347

correct the sign of higher order terms for beta1

0e1954d

add new mesh files for the 3 points loading test

913cbec

change the upper limit of d to 1 in the degradation function

890a7ed

Merge branch 'develop' into feature/frank/phasefeld_nucleation_formul…

495846f

…ation To keep the branch up-to-date

Merge branch 'develop' into feature/frank/phasefeld_nucleation_formul…

3f339d9

…ation

Merge branch 'develop' into feature/frank/phasefeld_nucleation_formul…

8d66f15

…ation

add crack irreversibility function

af901ec

frankfeifan changed the title ~~Added a new solver for the phase-field nucleation modeling of hydraulic fracturing~~ feat: phase-field nucleation modeling of hydraulic fracturing Dec 4, 2024

frankfeifan marked this pull request as ready for review December 4, 2024 18:15

frankfeifan requested review from castelletto1, jhuang2601, paveltomin, ryar9534, corbett5, cssherman and dkachuma as code owners December 4, 2024 18:15

CusiniM reviewed Jan 3, 2025

View reviewed changes

frankfeifan added 3 commits January 3, 2025 14:18

removed scheme docs files

d3142a5

addressed Matteo's comments

adfff8a

Merge remote-tracking branch 'origin/develop' into feature/frank/phas…

42d3e88

…efield_nucleation_poromech

frankfeifan self-assigned this Jan 4, 2025

frankfeifan commented Jan 4, 2025

View reviewed changes

...mponents/physicsSolvers/multiphysics/poromechanicsKernels/SinglePhasePoromechanicsDamage.hpp Outdated Show resolved Hide resolved

frankfeifan added 8 commits January 10, 2025 11:41

made SinglePhasePoromechanicsDamage a derivative of SinglePhasePorome…

2071f1b

…chanics kernel and updated copyright content

made the damage dependent constant an input parameter

a84b3f3

cleaned up xmls

7f1ba07

cleaned up phase-field damage kernels

4383a0e

uncrustify

722c8ac

Merge remote-tracking branch 'origin/develop' into feature/frank/phas…

4221b78

…efield_nucleation_poromech

removed pressurized damage kernel for nonporous impermeable solid in …

0016489

…this branch

modified xmls and added a smoke test

9b3a5f0

frankfeifan requested a review from wrtobin as a code owner January 13, 2025 19:31

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

feat: phase-field nucleation modeling of hydraulic fracturing #2125

feat: phase-field nucleation modeling of hydraulic fracturing #2125

frankfeifan commented Oct 15, 2022 •

edited

Loading

CusiniM left a comment

CusiniM Jan 3, 2025

CusiniM Jan 3, 2025

frankfeifan Jan 4, 2025

CusiniM Jan 3, 2025

frankfeifan Jan 4, 2025

codecov bot commented Jan 14, 2025

	m_scaleSet[i] = std::min( 1.0, (1.0 - faceScale/numNodes)*(1.0 - faceScale/numNodes) );
	m_scaleSet[i] = LvArray::math::min( 1.0, (1.0 - faceScale/numNodes)*(1.0 - faceScale/numNodes) );

feat: phase-field nucleation modeling of hydraulic fracturing #2125

Are you sure you want to change the base?

feat: phase-field nucleation modeling of hydraulic fracturing #2125

Conversation

frankfeifan commented Oct 15, 2022 • edited Loading

CusiniM left a comment

Choose a reason for hiding this comment

CusiniM Jan 3, 2025

Choose a reason for hiding this comment

CusiniM Jan 3, 2025

Choose a reason for hiding this comment

frankfeifan Jan 4, 2025

Choose a reason for hiding this comment

CusiniM Jan 3, 2025

Choose a reason for hiding this comment

frankfeifan Jan 4, 2025

Choose a reason for hiding this comment

codecov bot commented Jan 14, 2025

Codecov Report

frankfeifan commented Oct 15, 2022 •

edited

Loading