Skip to content

Debris-bed friction slip law#1101

Closed
jeremy-brooks-1 wants to merge 2 commits into
sandialabs:masterfrom
jeremy-brooks-1:debris-friction-test
Closed

Debris-bed friction slip law#1101
jeremy-brooks-1 wants to merge 2 commits into
sandialabs:masterfrom
jeremy-brooks-1:debris-friction-test

Conversation

@jeremy-brooks-1

Copy link
Copy Markdown
Contributor

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.

@jeremy-brooks-1

Copy link
Copy Markdown
Contributor Author

@mperego Here is the branch we mentioned on slack attempting to add debris-bed friction. Let me know if you need additional info, thanks!

@jeremy-brooks-1

Copy link
Copy Markdown
Contributor Author

Corresponding MPAS PR is: MALI-Dev/E3SM#135

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

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

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?)

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

sure, can do

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

Apologies for the delay, but it has now been rebased to Albany master.

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));

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

Were you thinking to uncomment this? As it is, the basal debris field and the bulk friction field are unused...

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

yes, this should be uncommented thanks for catching that!

@mperego

mperego commented Dec 20, 2024

Copy link
Copy Markdown
Collaborator

@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.
Note: N0 and u0 are hard-coded, but need to be definable  via YAML file
(as scalars, not fields)
@mperego

mperego commented Feb 18, 2026

Copy link
Copy Markdown
Collaborator

Closing in favor of newer PR #1184.

@mperego mperego closed this Feb 18, 2026
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