/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


phases (air slurry debris);


     air 
     {
         transportModel Newtonian;
         nu nu [ 0 2 -1 0 0 0 0 ] 1.48e-5; 
         rho rho [ 1 -3 0 0 0 0 0 ] 1;
     }

     slurry
     {

	transportModel HerschelBulkleyDebrisFlow;
		nu             nu [ 0 2 -1 0 0 0 0 ] 0.05;//dummy

		HerschelBulkleyDebrisFlowCoeffs
		{ 
		waterContent waterContent [ 0 0 0 0 0 0 0 ] 0.285;
		Cv Cv [ 0 0 0 0 0 0 0 ] 0.48643; // based on density of mixture 1803 kg/m^3 
		P0v P0v [ 0 0 0 0 0 0 0 ] 0.08743; // based on solid clay density 1700 kg/m^3
		tau00 tau00 [ 1 -1 -2 0 0 0 0 ] 41.333; // calibrated value [Pa]
		rhoF rhoF [ 1 -3 0 0 0 0 0 ] 1803;
		nu0 nu0 [ 0 2 -1 0 0 0 0 ] 2;
		// dummies that will be overwritten but are necessary for unit consistency:
		tau0 tau0 [ 0 2 -2 0 0 0 0 ] 0.015018;// tao00 based on  41.333 Pa at 28.5% water content by weight
		k k [ 0 2 -1 0 0 0 0 ] 0.004505;
		n n [ 0 0 0 0 0 0 0 ] 0.34;
		} 
		rho               rho [ 1 -3 0 0 0 0 0 ] 1803;

     }

     debris
     {
	    transportModel PudasainiCoulombViscoPlastic;
		    nu              nu [0 2 -1 0 0 0 0] 0.05;// dummy, not needed
		    rho             rho [1 -3 0 0 0 0 0] 1803;
		    PudasainiCoulombViscoPlasticCoeffs 
		    {
			delta           delta [0 0 0 0 0 0 0] 0.628319;//36 rough gravel based on tilt table and angle of repose tests
		       
			rhoS            rhoS [1 -3 0 0 0 0 0] 1803; 
		 	nuMin           nuMin [0 2 -1 0 0 0 0] 1.e-6;
			nu0             nu0 [0 2 -1 0 0 0 0] 2; 
			my		my [0 0 1 0 0 0 0] 0.2;
		    }

     }



// Surface tension coefficients
sigma12           sigma12 [1 0 -2 0 0 0 0] 0.07;
sigma13           sigma13 [1 0 -2 0 0 0 0] 0.07;

// Diffusivity between miscible phases
D23               D23   [0 2 -1 0 0 0 0]  0;
// solid gravel density
rhoS            rhoS [1 -3 0 0 0 0 0] 2700; 
// grain diameter (not needed)
dS              dS [0 1 0 0 0 0 0] 0.01; 
