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.

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

Implementation of thermal pressurisation #28

Merged
merged 10 commits into from
Aug 21, 2019

Conversation

martijnende
Copy link
Collaborator

@martijnende martijnende commented Jul 30, 2019

An improved implementation for high-velocity thermal pressurisation, based on Noda & Lapusta (2010). The diffusion of fluid and heat normal to the fault plane is solved in a discretised spectral domain, which is numerically stable for large time steps (as opposed to discretisation methods in the spatial domain, such as FEM, FD, etc.). This procedure allows for efficient modelling of frictional heating and associated pressurisation of the pore fluid.

This module is compatible with both the rate-and-state friction and the Chen-Niemeijer-Spiers (CNS) frameworks. When the CNS model is used, dilatancy hardening is accounted for through a dilatancy efficiency coefficient (set to 0 to ignore dilatancy hardening).

Fixed a small bug in pyqdyn.py in which multiple columns in the time series/snapshots output were named identically.

@martijnende martijnende mentioned this pull request Jul 30, 2019
@martijnende martijnende merged commit 2019a69 into release/2.1.0 Aug 21, 2019
martijnende added a commit that referenced this pull request Aug 21, 2019
* New feature: LVFZ in FINITE=1 loading conditions (#26)
1D faults with a low-velocity fault zone are now available for FINITE=1 loading conditions (See documentation).

* New feature: Implementation of thermal pressurisation (#28)
An improved implementation for high-velocity thermal pressurisation, based on Noda & Lapusta (2010). The diffusion of fluid and heat normal to the fault plane is solved in a discretised spectral domain, which is numerically stable for large time steps (as opposed to discretisation methods in the spatial domain, such as FEM, FD, etc.). This procedure allows for efficient modelling of frictional heating and associated pressurisation of the pore fluid.

This module is compatible with both the rate-and-state friction and the Chen-Niemeijer-Spiers (CNS) frameworks. When the CNS model is used, dilatancy hardening is accounted for through a dilatancy efficiency coefficient (set to 0 to ignore dilatancy hardening).

* Bug fig: Fixed a small bug in `pyqdyn.py` in which multiple columns in the time series/snapshots output were named identically.
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.

1 participant