Skip to content

feat(com1DFA): add muI friction model - #1346

Open
fso42 wants to merge 2 commits into
masterfrom
moOfI
Open

fso42 wants to merge 2 commits into
masterfrom
moOfI

Conversation

@fso42

@fso42 fso42 commented Nov 5, 2026 •

Copy link
Copy Markdown
Contributor

Adds first draft of mu(I) rheology

PR Checklist

Please confirm before requesting review:

  • I ran pytest locally without fails
  • I added/updated tests where needed
  • I updated documentation where needed

Confirm before the final merge/rebase into master

  • Commits are sensibly squashed and rebased onto latest master
  • Standardtest run without difference (with recompiled cython code)

@fso42
fso42 requested a review from PaulaSp3 November 5, 2026 08:04
@fso42 fso42 self-assigned this Nov 5, 2026
@fso42 fso42 added the enhancement New feature or request label Nov 5, 2026
@qltysh

qltysh Bot commented Nov 5, 2026 •

Copy link
Copy Markdown
Contributor

Analysis for project AvaFrame

❌ 2 blocking issues (2 total)

Tool Category Rule Count
black Style Incorrect formatting, autoformat by running qlty fmt. 2 ❌

@qltysh one-click actions:

  • Auto-fix formatting (qlty fmt && git push)

Comment thread avaframe/com1DFA/checkCfg.py Resolved
@JuLa96

JuLa96 commented Nov 6, 2026

Copy link
Copy Markdown
Contributor

closes OpenNHM/DebrisFrame#120

Comment thread docs/theoryCom1DFA.rst
the flow density ``rho`` (probably the bulk density was entered) or if :math:`c_\gamma \neq 5/2`.

.. Note:: With :eq:`muI-inertial` and :eq:`muI-shearRate`, :math:`I / I_0` depends on :math:`d`, :math:`I_0` and
:math:`\rho_g` only through the combination :math:`d / (I_0 \sqrt{\rho / \rho_g})` (:math:`\rho` is the flow

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

For the fit below, do you consider flow density rho?
If not, it's easier to understand (IMO) to not mention flow density here, if it is consider I suggest mention the used flow density values.

Comment thread avaframe/com1DFA/com1DFACfg.ini Resolved
Comment thread docs/theoryCom1DFA.rst Resolved
Comment thread docs/theoryCom1DFA.rst Resolved
Comment thread docs/theoryCom1DFA.rst Resolved
Comment thread docs/theoryCom1DFA.rst

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Maybe it makes sense to clarify the difference between depth and thickness, e.g. that other references use depth-averaged model, but the mu(I) model is transferable to thickness-averaged models, that's why we're talking about thickness?

Comment thread docs/theoryCom1DFA.rst
\overline{u}^* = \frac{2}{5}\,\frac{I^*}{d}\,\sqrt{\frac{\rho}{\rho_g}\,g\cos\theta}\;\overline{h}^{3/2}
:label: muI-steady

which is the Pouliquen flow rule :math:`Fr = \beta\,h / h_{stop}(\theta)` with

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Suggested change
which is the Pouliquen flow rule :math:`Fr = \beta\,h / h_{stop}(\theta)` with
which is the Pouliquen flow rule describing the Froud number as :math:`Fr = \beta\,h / h_{stop}(\theta)` with

Comment thread docs/theoryCom1DFA.rst
:math:`\beta = \frac{2}{5}\,I_0\,\sqrt{(\rho / \rho_g)\cos\theta}` :cite:`GrEd2014`. Equation
:eq:`muI-steady` is used as the analytic reference in the unit tests. Consequences:

- Below :math:`\mu_s` the flow stops (the friction force equals the driving force at rest, as for Coulomb friction).

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Suggested change
- Below :math:`\mu_s` the flow stops (the friction force equals the driving force at rest, as for Coulomb friction).
- Below :math:`\arctan\mu_s` the flow stops (the friction force equals the driving force at rest, as for Coulomb friction).

Comment thread avaframe/com1DFA/DFAfunctionsCython.pyx Resolved

This branch has not been deployed

No deployments
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants