Pressure calculation issue - granular collapse on slopes

45 views
Skip to first unread message

Alexis Bougouin

unread,
Dec 9, 2025, 4:07:35 AM (10 days ago) Dec 9
to basilisk-fr
Dear community,

I'm currently performing Navier-Stokes simulations of granular collapses on different slopes, based on Pierre-Yves Lagrée's example entitled 'Continuum collapse of a column of grains real 2D, no slip at the base'. The general idea is to vary the slope angle up to 36° in order to compare the results with large-scale experiments already carried out, using Coulomb and mu(I)-rheology.

I have several issues, all of which are likely related to pressure calculation, including the slope effect.

- Regardless of rheology, the calculated pressure appears unstable as the slope increases without affecting the velocity (see attached videos). I included the slope only in gravity and pressure initiation (see attached script).
- I have no deposits along the plane above theta ~ 25°, while the angle of repose is set to 32°.
- Inherent to the numerical method, bubbles are trapped at the front, which significantly affect the flow dynamics at steep slopes. There is no solution to this problem, and therefore this numerical method is ineffective at steep slopes?

Thank you for your help,
velocity.mpg
UnsteadyGranularDamBreak.c
pressure.mpg

"Pierre-Yves Lagrée (PYL)"

unread,
Dec 9, 2025, 4:15:41 AM (10 days ago) Dec 9
to Alexis Bougouin, basilisk-fr
Dear Alexis 

Maybe the instability of the m(I) is acting there, we should add the Barker Gray regularisation
For the "bubbles" it is inherent to the Boundary conditions u=0.  Maybe we should add a special BC at the tip (where the mu(I) rheology fails by the way)


Cheers 

PS: on peut continuer en privé en français 

--
You received this message because you are subscribed to the Google Groups "basilisk-fr" group.
To unsubscribe from this group and stop receiving emails from it, send an email to basilisk-fr...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/basilisk-fr/8f4cb685-bfa3-453d-827c-ecdb3f9c51ean%40googlegroups.com.
<velocity.mpg><UnsteadyGranularDamBreak.c><pressure.mpg>

----> ----> ----> ----> ----> ----> ----> ----> 
--->  ---> Pierre-Yves Lagrée (PYL) ----> ---->
//////////////////////////////////////////////

Reply all
Reply to author
Forward
0 new messages