Skip to content

feat: Enable aperture modifications without a full mechanics solver#3906

Draft
LilianeRA wants to merge 32 commits into
developfrom
feature/prescribedStressPath_constitutive
Draft

feat: Enable aperture modifications without a full mechanics solver#3906
LilianeRA wants to merge 32 commits into
developfrom
feature/prescribedStressPath_constitutive

Conversation

@LilianeRA

@LilianeRA LilianeRA commented Nov 7, 2025

Copy link
Copy Markdown

This PR is the first step for implementing a simplified hydromechanical coupling based on the Barton–Bandis aperture model and a prescribed stress path.

The main change introduced here is the new permeability constitutive law (BartonBandisPermeability) that internally computes fracture aperture, according to the current pressure and normal stress, before updating its permeability.


EDIT: as advised, we are reformulating the code design to better encapsulate what is new: the prescribed stress path implementation. Please, don't review the current code until further notice.

…ult. A problem sample is defined in 'stressPathDrivenGeomechanics/'. The pressure and aperture data will be saved in the files 'fracPressure_history.hdf5' and 'fracAperture_history.hdf5'.
…ic trajectory: BartonBandisStressPathDriven.
…ePassThru function for our new constitutive law; Deleting commented function in FlowSolverBase; Using previous implementation of the Barton Bandis law, with Kni
… the RST documentation; Removing flag for updating stencils; Minor documentation update.
…he new fracture aperture is calculated internally, and it is not propagated to the aperture field. For now, the test requires that the pressure difference to be very small.
…. Replacing it by a constitutive relation for permeability. Updated aperture computed internally. Results still validated by analytical solution. Sample file updated.
@jafranc jafranc self-requested a review April 15, 2026 06:55
Comment thread src/CMakeLists.txt Outdated
Comment thread .gitignore Outdated
Comment thread src/coreComponents/physicsSolvers/fluidFlow/SinglePhaseBase.cpp Outdated
Comment thread src/coreComponents/physicsSolvers/fluidFlow/FlowSolverBase.hpp Outdated
Comment thread src/coreComponents/physicsSolvers/fluidFlow/FlowSolverBase.cpp
Comment thread src/coreComponents/constitutive/permeability/BartonBandisPermeability.hpp Outdated
Comment thread src/coreComponents/constitutive/permeability/BartonBandisPermeability.hpp Outdated
Comment thread src/coreComponents/constitutive/permeability/BartonBandisPermeability.hpp Outdated
Comment thread src/coreComponents/constitutive/permeability/BartonBandisPermeability.hpp Outdated

@dkachuma dkachuma left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Can you please remove any debug/dead code and sync to develop. Things have moved quite a bit making this a little harder to review.

Comment thread src/coreComponents/constitutive/permeability/BartonBandisPermeability.hpp Outdated
Comment thread src/coreComponents/physicsSolvers/fluidFlow/FlowSolverBase.hpp Outdated
Comment thread src/coreComponents/physicsSolvers/fluidFlow/FlowSolverBase.cpp Outdated
Comment thread src/coreComponents/constitutive/solid/CompressibleSolid.cpp
LilianeRA and others added 15 commits May 22, 2026 16:05
Removing unused optional flags for input file. Remnants of when the proposed method was modeled as a contact constituive class.

Co-authored-by: Dickson Kachuma <81433670+dkachuma@users.noreply.github.com>
…ability.hpp


Better variable name description

Co-authored-by: Dickson Kachuma <81433670+dkachuma@users.noreply.github.com>
…meability.hpp


Suggestion by Kachuma: "I would still use the quadrature index although it's always 0. Just for posterity."

Co-authored-by: Dickson Kachuma <81433670+dkachuma@users.noreply.github.com>
…ult. A problem sample is defined in 'stressPathDrivenGeomechanics/'. The pressure and aperture data will be saved in the files 'fracPressure_history.hdf5' and 'fracAperture_history.hdf5'.
…ePassThru function for our new constitutive law; Deleting commented function in FlowSolverBase; Using previous implementation of the Barton Bandis law, with Kni
…he new fracture aperture is calculated internally, and it is not propagated to the aperture field. For now, the test requires that the pressure difference to be very small.
…. Replacing it by a constitutive relation for permeability. Updated aperture computed internally. Results still validated by analytical solution. Sample file updated.
…lts change slightly in the first few iterations.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants