Skip to main content

more options


The MATLAB m-files below incorporate the theory of Louge and Keast (2001) predicting flows of grains down flat, frictional inclines; in particular, they calculate the effective friction in terms of parameters (restitution, friction, angle of inclination, etc) that can be measured independently.

The first step is to download all programs listed below.  There are two principal programs:



[Fres] = SFD_flat_chute(tanalphai,holdupi,ewi,muwi,b0wi,eei,muui,b00i,muEi,guess)

You can enter

help plot_flat_chute_mueff
help SFD_flat_chute

to read specifications for each file.  Here are two examples you can try:

plot_flat_chute_mueff(.2,10,50,[3],[.2 .4 .6],.8,.141,.35,.972,.177,.25,.593,1,.4) ;
[Fres] = SFD_flat_chute([.278:.001:.364],5,.8,.141,.35,.972,.177,.25,.593,[2 .5]) ;

The first program, plot_flat_chute_mueff, allows you to calculate the effective friction in terms of Fr for a given value of the solid volume fraction at the base.  (It plots the thin lines in Fig. 8 of Louge and Keast).  It does not assume that the flow is steady and fully-developed (SFD).  Therefore, its curves decrease with Fr at constant solid volume fraction.

Instead, the second program, SFD_flat_chute, assumes that the flow is SFD and uses the steady angular momentum balance to find the solid volume fraction in the shear layer.  This program also graph predictions and compares with our own measurements.