Debris-bed friction slip law #1101
Conversation
Add debris-bed friction slip relation. Adding two new friction parameters, could set as constants in albany_input.yaml or read fields from MPAS mesh: *bulkFrictionCoefficient *basalDebrisFactor There are corresponding changes in MPAS interface in MALI-Dev/ES3M repository.
|
@mperego Here is the branch we mentioned on slack attempting to add debris-bed friction. Let me know if you need additional info, thanks! |
|
Corresponding MPAS PR is: MALI-Dev/E3SM#135 |
| @@ -714,7 +699,7 @@ void velocity_solver_extrude_3d_grid(int nLayers, int globalTrianglesStride, | |||
|
|
|||
| discretizationList->set("Workset Size", discretizationList->get("Workset Size", -1)); | |||
|
|
|||
| discretizationList->set("Method", discretizationList->get("Method", "STKExtruded")); //set to STKExtruded is not defined | |||
| discretizationList->set("Method", discretizationList->get("Method", "Extruded")); //set to Extruded is not defined | |||
There was a problem hiding this comment.
It seems you based this PR off of a odler version of Albany, can you please rebase to this branch to Albany master? (or merge master into this branch?)
| switch (flowRate_type) { | ||
| case FLOW_RATE_TYPE::CONSTANT: { | ||
| beta(cell,ipt) /= std::pow ( u_norm(cell,ipt) + bedRoughnessValue*scaling*flowRate_val*std::pow(NVal,n), power); //bedRoughness in km | ||
| // beta(cell,ipt) = ((muValue * NVal * std::pow (u_norm(cell,ipt), power-1.0)) / std::pow ( u_norm(cell,ipt) + bedRoughnessValue*scaling*flowRate_val*std::pow(NVal,n), power)) + (bulkFrictionValue * NVal / u_norm(cell,ipt)) + (basalDebrisValue * u_norm(cell,ipt)); |
There was a problem hiding this comment.
Were you thinking to uncomment this? As it is, the basal debris field and the bulk friction field are unused...
There was a problem hiding this comment.
yes, this should be uncommented thanks for catching that!
|
@jeremy-brooks-1 thanks for working on this. Before this can be merged we should add a test that exercise this capability. I can help you with this. |
Add debris-bed friction slip relation.
Adding two new friction parameters, could set as constants in albany_input.yaml or read fields from MPAS mesh:
*bulkFrictionCoefficient
*basalDebrisFactor
There are corresponding changes in MPAS interface in MALI-Dev/ES3M repository.