search

geodynamics / aspect

A parallel, extensible finite element code to simulate convection in both 2D and 3D models.
https://aspect.geodynamics.org/
Other
223 stars 235 forks source link

Add cutoff for dry/wet creep in water fugacity calculation #5808

Closed danieldouglas92 closed 4 months ago

danieldouglas92 commented 4 months ago

Note that this PR is an extension of PR #5805, it takes the test implemented in that PR and modifies it slightly to serve as a test for the functionality implemented in this PR.

If there is 0 bound water in the model domain and the user tries using the HK04 olivine hydration viscosity prefactor scheme there is a divide by 0 error. In Hirth & Kohlstaedt 2004, they discuss that when the hydration state of olivine is below 50 H/Si ppm (6.15e-4 wt% water), olivine deformation is consistent with dry conditions, and above 6.15e-4 wt% water olivine deformation is consistent with wet conditions. This PR adds a user parameter that defines this cutoff between dry and wet deformation, with the default value representing the value discussed in Hirth & Kohlstaedt 2004.

bobmyhill commented 4 months ago

/rebuild

naliboff commented 4 months ago

@bobmyhill - this looks good from my end if you would like to take a final look before merging.