Add debris-bed friction slip law #135
Conversation
jeremy-brooks-1
added
the
Albany
|
Requires Albany changes in: sandialabs/Albany#1101 |
jeremy-brooks-1
force-pushed
the
jeremy-brooks-1/mali/debris-friction-full
branch
from
db178a3 to
c924550
Compare
This commit adds a new edRoughnessRC variable to MALI on the MPAS side. Corresponding changes are needed on the Albany side.
Adding two fields (bulkFriction, basalDebris) to use with debris-bed friction slip law.
jeremy-brooks-1
force-pushed
the
jeremy-brooks-1/mali/debris-friction-full
branch
from
c924550 to
079213d
Compare
Updated Registry and MPAS velocity interface
| const std::vector<double>& smbData, | ||
| const std::vector<double>& stiffnessFactorData, | ||
| const std::vector<double>& effecPressData, | ||
| std::vector<double>& effecPressData, |
There was a problem hiding this comment.
@mperego , was this change intentional?
There was a problem hiding this comment.
Yes. Now Albany can return the effective pressure used (e.g., the open-connection one).
| <var name="bulkFriction" type="real" dimensions="nCells Time" default_value="0.0" | ||
| description="input value of bulk friction coefficient in debris-bed friction sliding law." | ||
| /> | ||
| <var name="basalDebrisFactor" type="real" dimensions="nCells Time" units="Pa yr m^-1" default_value="0.0" |
There was a problem hiding this comment.
@mperego , are these the correct units we should use here? Do we need to account for any scaling that Albany does here? (noting that I'm not sure if Jeremy had implemented scaling on these fields on the Albany side)
There was a problem hiding this comment.
Based on how they are used in Albany, I think bedRoughnessRC is fine, but basalDebrisFactor should be in kPa yr m^{-1}, like muFriction. We can convert the unis in the interface_velocity_solver.cpp file.
TBD