*** Data corresponding to the publication: "Transient behavior and steady-state rheology of dense frictional suspensions in pressure-driven channel flow." ***
Authors: W. Peerbooms (1), A.E.D.M. van der Heijden (2,3), W.-P. Breugem (1)
(1) Process and Energy Department, Multiphase Systems, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
(2) Process and Energy Department, Complex Fluid Processing, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
(3) Department Energetic Materials, Netherlands Organization for Applied Scientific Research (TNO), Ypenburgse Boslaan 2, 2496 ZA The Hague, The Netherlands

Corresponding author: W.-P. Breugem
Contact information: W.P.Breugem@tudelft.nl

*** General Introduction ***
This dataset contains data underlying the figures presented in the publication: 
Peerbooms, W., Heijden, A.v.d. & Breugem, WP. Transient behavior and steady-state rheology of dense frictional suspensions in pressure-driven channel flow. 
Acta Mech (2024). https://doi.org/10.1007/s00707-024-04106-7

The authors gratefully acknowledge Delft University of Technology, the funding of this project under a TU Delft/TNO agreement, and the support of SURF Cooperative for the use of the 
Dutch national e-infrastructure (Cartesius and Snellius supercomputers) under Computing Grant No. 2023.009 “Numerical simulation of complex flows in complex geometries"
from the Dutch Research Council (NWO).

*** Description of the data ***
The data is divided in .txt files for most figures in the article.
The data in the files is organized as follows:

*Fig3a.txt
Contains the normalized concentration profile phi/phi_b as a function of the wall-normal coordinate for various time instants for the case with phi_b=0.4 and muc=0
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) phi/phi_b for 0.2tdev
3) phi/phi_b for 0.4tdev
4) phi/phi_b for 0.6tdev
5) phi/phi_b for 0.8tdev
6) phi/phi_b for tdev

*Fig3b.txt
Contains the normalized concentration profile phi/phi_b as a function of the wall-normal coordinate for various time instants for the case with phi_b=0.4 and muc=0.5
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) phi/phi_b for 0.2tdev
3) phi/phi_b for 0.4tdev
4) phi/phi_b for 0.6tdev
5) phi/phi_b for 0.8tdev
6) phi/phi_b for tdev

*Fig4a.txt
Contains normalized concentration profile phi/phi_b as a function of time at various wall normal positions for the case with phi_b=0.4 and muc=0
The columns correspond to:
1) tvb/d time
2) phi/phi_b for z/d=0.5
3) phi/phi_b for z/d=5
4) phi/phi_b for z/d=10

*Fig4b.txt
Contains the normalized concentration profile phi/phi_b as a function of time at various wall normal positions for the case with phi_b=0.4 and muc=0.5
The columns correspond to:
1) tvb/d time
2) phi/phi_b for z/d=0.5
3) phi/phi_b for z/d=5
4) phi/phi_b for z/d=10

*Fig5a.txt
Contains the intrinsic fluid velocity profile as a function of the wall-normal coordinate for various time instants for the case with phi_b=0.4 and muc=0
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) intrinsic fluid velocity profile for 0.2tdev
3) intrinsic fluid velocity profile for 0.4tdev
4) intrinsic fluid velocity profile for 0.6tdev
5) intrinsic fluid velocity profile for 0.8tdev
6) intrinsic fluid velocity profile for tdev

*Fig5b.txt
Contains the intrinsic fluid velocity profile as a function of the wall-normal coordinate for various time instants for the case with phi_b=0.4 and muc=0.5
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) intrinsic fluid velocity profile for 0.2tdev
3) intrinsic fluid velocity profile for 0.4tdev
4) intrinsic fluid velocity profile for 0.6tdev
5) intrinsic fluid velocity profile for 0.8tdev
6) intrinsic fluid velocity profile for tdev

*Fig6.txt
Contains the evolution parameter as a function of time for all cases
The columns correspond to:
1)  tvb/d time
2)  evolution parameter for the case with phi_b=0.2 and muc=0
3)  evolution parameter for the case with phi_b=0.2 and muc=0.5
4)  evolution parameter for the case with phi_b=0.3 and muc=0
5)  evolution parameter for the case with phi_b=0.3 and muc=0.5
6)  evolution parameter for the case with phi_b=0.4 and muc=0
7)  evolution parameter for the case with phi_b=0.4 and muc=0.5
8)  evolution parameter for the case with phi_b=0.5 and muc=0
9)  evolution parameter for the case with phi_b=0.5 and muc=0.5
10) evolution parameter for the case with phi_b=0.6 and muc=0
11) evolution parameter for the case with phi_b=0.6 and muc=0.5

*Fig7a.txt
Contains the timescale of particle migration as a function of the bulk solid volume fraction
The columns correspond to:
1) phi_b bulk solid volume fraction
2) the timescale of particle migration for the cases with mu_c = 0
3) the timescale of particle migration for the cases with mu_c = 0.5

*Fig7b.txt
Contains the development length as a function of normalized bulk solid volume fraction
The columns correspond to:
1) phi_b/Phi_R normalized bulk solid volume fraction
2) L_dev d^2/4H^3 development length
3) mu_c Coulomb coefficient of siding friction

*Fig8a.txt
Contains the superficial wall-normal particle stress profile as a function of the wall-normal coordinate for various time instants for the case with phi_b=0.4 and muc=0
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) wall-normal particle stress profile for 0.2tdev
3) wall-normal particle stress profile for 0.4tdev
4) wall-normal particle stress profile for 0.6tdev
5) wall-normal particle stress profile for 0.8tdev
6) wall-normal particle stress profile for tdev

*Fig8b.txt
Contains the superficial wall-normal particle stress profile as a function of the wall-normal coordinate for various time instants for the case with phi_b=0.4 and muc=0.5
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) wall-normal particle stress profile for 0.2tdev
3) wall-normal particle stress profile for 0.4tdev
4) wall-normal particle stress profile for 0.6tdev
5) wall-normal particle stress profile for 0.8tdev
6) wall-normal particle stress profile for tdev

*Fig9a.txt
Contains the steady-state solid volume fraction as a function of the wall-normal coordinate for the cases with muc=0
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) steady-state solid volume fraction for the case with phi_b = 0.4
3) steady-state solid volume fraction for the case with phi_b = 0.5
4) steady-state solid volume fraction for the case with phi_b = 0.6

*Fig9b.txt
Contains the steady-state solid volume fraction as a function of the wall-normal coordinate for the cases with muc=0.5
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) steady-state solid volume fraction for the case with phi_b = 0.3
3) steady-state solid volume fraction for the case with phi_b = 0.4
4) steady-state solid volume fraction for the case with phi_b = 0.5
5) steady-state solid volume fraction for the case with phi_b = 0.6

*Fig10a.txt
Contains the steady-state intrinsic streamwise fluid velocity as a function of the wall-normal coordinate for the cases with muc=0
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) steady-state intrinsic streamwise fluid velocity for the case with phi_b = 0.4
3) steady-state intrinsic streamwise fluid velocity for the case with phi_b = 0.5
4) steady-state intrinsic streamwise fluid velocity for the case with phi_b = 0.6

*Fig10b.txt
Contains the steady-state intrinsic streamwise fluid velocity as a function of the wall-normal coordinate for the cases with muc=0.5
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) steady-state intrinsic streamwise fluid velocity for the case with phi_b = 0.3
3) steady-state intrinsic streamwise fluid velocity for the case with phi_b = 0.4
4) steady-state intrinsic streamwise fluid velocity for the case with phi_b = 0.5
5) steady-state intrinsic streamwise fluid velocity for the case with phi_b = 0.6

*Fig13a.txt
Contains the steady-state mean superficial shear stress as a function of the wall-normal coordinate for the case with phi_b=0.4 and muc=0
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) dp_e/dy
3) sigma_{visc,12}
4) sigma_{p,12}
5) sigma_{p,12}^{hyd}
6) sigma_{p,12}^{cont}
7) sigma_{12}

*Fig13b.txt
Contains the steady-state mean superficial shear stress as a function of the wall-normal coordinate for the case with phi_b=0.4 and muc=0.5
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) dp_e/dy
3) sigma_{visc,12}
4) sigma_{p,12}
5) sigma_{p,12}^{hyd}
6) sigma_{p,12}^{cont}
7) sigma_{12}

*Fig14.txt
Contains the local relative shear viscosity as a function of the local volume fraction
The columns correspond to:
1)  local volume fraction for the case with 		phi_b=0.3 and mu_c=0.5
2)  local relative shear viscosity for the case with 	phi_b=0.3 and mu_c=0.5
3)  local volume fraction for the case with 		phi_b=0.4 and mu_c=0
4)  local relative shear viscosity for the case with 	phi_b=0.4 and mu_c=0
5)  local volume fraction for the case with 		phi_b=0.4 and mu_c=0.5
6)  local relative shear viscosity for the case with 	phi_b=0.4 and mu_c=0.5
7)  local volume fraction for the case with 		phi_b=0.5 and mu_c=0
8)  local relative shear viscosity for the case with 	phi_b=0.5 and mu_c=0
9)  local volume fraction for the case with 		phi_b=0.5 and mu_c=0.5
10) local relative shear viscosity for the case with 	phi_b=0.5 and mu_c=0.5

*Fig15.txt
Contains the hydrodynamic and contact contributions to the local relative shear viscosity as a function of the local volume fraction
The columns correspond to:
1)   local volume fraction for the case with 						phi_b=0.3 and mu_c=0.5
2)   hydrodynamic contribution to the local relative shear viscosity for the case with 	phi_b=0.3 and mu_c=0.5
3)   contact contribution to the local relative shear viscosity for the case with 	phi_b=0.3 and mu_c=0.5
4)   local volume fraction for the case with 						phi_b=0.4 and mu_c=0
5)   hydrodynamic contribution to the local relative shear viscosity for the case with 	phi_b=0.4 and mu_c=0
6)   contact contribution to the local relative shear viscosity for the case with 	phi_b=0.4 and mu_c=0
7)   local volume fraction for the case with 						phi_b=0.4 and mu_c=0.5
8)   hydrodynamic contribution to the local relative shear viscosity for the case with 	phi_b=0.4 and mu_c=0.5
9)   contact contribution to the local relative shear viscosity for the case with 	phi_b=0.4 and mu_c=0.5
10)  local volume fraction for the case with 						phi_b=0.5 and mu_c=0
11)  hydrodynamic contribution to the local relative shear viscosity for the case with 	phi_b=0.5 and mu_c=0
12)  contact contribution to the local relative shear viscosity for the case with 	phi_b=0.5 and mu_c=0
13)  local volume fraction for the case with 						phi_b=0.5 and mu_c=0.5
14)  hydrodynamic contribution to the local relative shear viscosity for the case with 	phi_b=0.5 and mu_c=0.5
15)  contact contribution to the local relative shear viscosity for the case with 	phi_b=0.5 and mu_c=0.5


*Fig16a.txt
Contains the steady-state mean superficial wall-normal stress as a function of the wall-normal coordinate for the case with phi_b=0.4 and muc=0
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) -p
3) sigma_{p,22}
4) sigma_{p,22}^{hyd}
5) sigma_{p,22}^{cont}
6) sigma_{22}
7) sigma_{wall,22}

*Fig16b.txt
Contains the steady-state mean superficial wall-normal stress as a function of the wall-normal coordinate for the case with phi_b=0.4 and muc=0.5
The columns correspond to:
1) z/d wall-normal coordinate normalized by the particle diameter
2) -p
3) sigma_{p,22}
4) sigma_{p,22}^{hyd}
5) sigma_{p,22}^{cont}
6) sigma_{22}
7) sigma_{wall,22}

*Fig17.txt
Contains the normal particle stresses as a function of the normalized local volume fraction
The columns correspond to:
1)   normalized local volume fraction for the case with phi_b=0.3 and mu_c=0.5
2)   sigma_p,11/sigma_12 for the case with 		phi_b=0.3 and mu_c=0.5
3)   sigma_p,22/sigma_12 for the case with  		phi_b=0.3 and mu_c=0.5
4)   sigma_p,33/sigma_12 for the case with  		phi_b=0.3 and mu_c=0.5
5)   normalized local volume fraction for the case with phi_b=0.4 and mu_c=0
6)   sigma_p,11/sigma_12 for the case with 		phi_b=0.4 and mu_c=0
7)   sigma_p,22/sigma_12 for the case with  		phi_b=0.4 and mu_c=0
8)   sigma_p,33/sigma_12 for the case with  		phi_b=0.4 and mu_c=0
9)   normalized local volume fraction for the case with phi_b=0.4 and mu_c=0.5
10)  sigma_p,11/sigma_12 for the case with 		phi_b=0.4 and mu_c=0.5
11)  sigma_p,22/sigma_12 for the case with  		phi_b=0.4 and mu_c=0.5
12)  sigma_p,33/sigma_12 for the case with  		phi_b=0.4 and mu_c=0.5
13)  normalized local volume fraction for the case with phi_b=0.5 and mu_c=0
14)  sigma_p,11/sigma_12 for the case with 		phi_b=0.5 and mu_c=0
15)  sigma_p,22/sigma_12 for the case with  		phi_b=0.5 and mu_c=0
16)  sigma_p,33/sigma_12 for the case with  		phi_b=0.5 and mu_c=0
17)  normalized local volume fraction for the case with	phi_b=0.5 and mu_c=0.5
18)  sigma_p,11/sigma_12 for the case with 		phi_b=0.5 and mu_c=0.5
19)  sigma_p,22/sigma_12 for the case with  		phi_b=0.5 and mu_c=0.5
20)  sigma_p,33/sigma_12 for the case with  		phi_b=0.5 and mu_c=0.5

*Fig18.txt
Contains the normal stress differences as a function of the normalized local volume fraction
The columns correspond to:
1)  normalized local volume fraction for the case with 	phi_b=0.3 and mu_c=0.5
2)  N1/sigma_12 for the case with 			phi_b=0.3 and mu_c=0.5
3)  N2/sigma_12 for the case with  			phi_b=0.3 and mu_c=0.5
4)  normalized local volume fraction for the case with 	phi_b=0.4 and mu_c=0
5)  N1/sigma_12 for the case with 			phi_b=0.4 and mu_c=0
6)  N2/sigma_12 for the case with  			phi_b=0.4 and mu_c=0
7)  normalized local volume fraction for the case with 	phi_b=0.4 and mu_c=0.5
8)  N1/sigma_12 for the case with 			phi_b=0.4 and mu_c=0.5
9)  N2/sigma_12 for the case with  			phi_b=0.4 and mu_c=0.5
10) normalized local volume fraction for the case with 	phi_b=0.5 and mu_c=0
11) N1/sigma_12 for the case with 			phi_b=0.5 and mu_c=0
12) N2/sigma_12 for the case with  			phi_b=0.5 and mu_c=0
13) normalized local volume fraction for the case with 	phi_b=0.5 and mu_c=0.5
14) N1/sigma_12 for the case with 			phi_b=0.5 and mu_c=0.5
15) N2/sigma_12 for the case with  			phi_b=0.5 and mu_c=0.5

*Fig19.txt
Contains the local frictional rheology as a function of the viscous number
The columns correspond to:
1)   normalized local volume fraction for the case with phi_b=0.3 and mu_c=0.5
2)   Iv for the case with 				phi_b=0.3 and mu_c=0.5
3)   mu for the case with  				phi_b=0.3 and mu_c=0.5
4)   mu^h for the case with  				phi_b=0.3 and mu_c=0.5
5)   mu^c for the case with  				phi_b=0.3 and mu_c=0.5
6)   normalized local volume fraction for the case with phi_b=0.4 and mu_c=0
7)   Iv for the case with 				phi_b=0.4 and mu_c=0
8)   mu for the case with  				phi_b=0.4 and mu_c=0
9)   mu^h for the case with  				phi_b=0.4 and mu_c=0
10)  mu^c for the case with  				phi_b=0.4 and mu_c=0
11)  normalized local volume fraction for the case with phi_b=0.4 and mu_c=0.5
12)  Iv for the case with 				phi_b=0.4 and mu_c=0.5
13)  mu for the case with  				phi_b=0.4 and mu_c=0.5
14)  mu^h for the case with  				phi_b=0.4 and mu_c=0.5
15)  mu^c for the case with  				phi_b=0.4 and mu_c=0.5
16)  normalized local volume fraction for the case with phi_b=0.5 and mu_c=0
17)  Iv for the case with 				phi_b=0.5 and mu_c=0
18)  mu for the case with  				phi_b=0.5 and mu_c=0
19)  mu^h for the case with  				phi_b=0.5 and mu_c=0
20)  mu^c for the case with  				phi_b=0.5 and mu_c=0
21)  normalized local volume fraction for the case with phi_b=0.5 and mu_c=0.5
22)  Iv for the case with 				phi_b=0.5 and mu_c=0.5
23)  mu for the case with  				phi_b=0.5 and mu_c=0.5
24)  mu^h for the case with  				phi_b=0.5 and mu_c=0.5
25)  mu^c for the case with  				phi_b=0.5 and mu_c=0.5

*Fig20.txt
Contains the fitting parameters for the modified frictional rheology model
The columns correspond to:
1) phi_b 
2) mu_c 
3) I_0^p
4) mu_2^p
5) phi_p

*Fig21.txt
Contains the rms shear viscous number as a function of the mean shear viscous number
The columns correspond to:
1)  Theta_v for the case with 	phi_b=0.3 and mu_c=0.5
2)  Iv for the case with 	phi_b=0.3 and mu_c=0.5
3)  Theta_v for the case with 	phi_b=0.4 and mu_c=0
4)  Iv for the case with 	phi_b=0.4 and mu_c=0
5)  Theta_v for the case with 	phi_b=0.4 and mu_c=0.5
6)  Iv for the case with 	phi_b=0.4 and mu_c=0.5
7)  Theta_v for the case with 	phi_b=0.5 and mu_c=0
8)  Iv for the case with 	phi_b=0.5 and mu_c=0
9)  Theta_v for the case with 	phi_b=0.5 and mu_c=0.5
10) Iv for the case with 	phi_b=0.5 and mu_c=0.5

*Fig22.txt
Contains I_v^*, I_0^p, and I_v^sub
The columns correspond to:
1) phi_b 
2) mu_c 
3) I_v^*
4) I_0^p
5) I_v^sub

*Fig23.txt
Contains the Local mean volume fraction normalized with the maximum flowable packing fraction as a function of the rms shear viscous number
The columns correspond to:
1)  Theta_v for the case with 	phi_b=0.3 and mu_c=0.5
2)  phi/phi_m for the case with phi_b=0.3 and mu_c=0.5
3)  Theta_v for the case with 	phi_b=0.4 and mu_c=0
4)  phi/phi_m for the case with phi_b=0.4 and mu_c=0
5)  Theta_v for the case with 	phi_b=0.4 and mu_c=0.5
6)  phi/phi_m for the case with phi_b=0.4 and mu_c=0.5
7)  Theta_v for the case with 	phi_b=0.5 and mu_c=0
8)  phi/phi_m for the case with phi_b=0.5 and mu_c=0
9)  Theta_v for the case with 	phi_b=0.5 and mu_c=0.5
10) phi/phi_m for the case with phi_b=0.5 and mu_c=0.5

*Fig24.txt
Contains the Local mean volume fraction normalized with the maximum flowable packing fraction as a function of the mean shear viscous number
The columns correspond to:
1)  Iv for the case with 	phi_b=0.3 and mu_c=0.5
2)  phi/phi_m for the case with phi_b=0.3 and mu_c=0.5
3)  Iv for the case with 	phi_b=0.4 and mu_c=0
4)  phi/phi_m for the case with phi_b=0.4 and mu_c=0
5)  Iv for the case with 	phi_b=0.4 and mu_c=0.5
6)  phi/phi_m for the case with phi_b=0.4 and mu_c=0.5
7)  Iv for the case with 	phi_b=0.5 and mu_c=0
8)  phi/phi_m for the case with phi_b=0.5 and mu_c=0
9)  Iv for the case with 	phi_b=0.5 and mu_c=0.5
10) phi/phi_m for the case with phi_b=0.5 and mu_c=0.5