/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      kinematicCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#include        "../setUp"  // include file with all the settings

// NUMERICAL SETTINGS
solution
{
	// Activation of particles and coupling method
    active          true;//false; // activation of particles
    coupled         true;//false; //two-way coupling
    transient       yes;
    cellValueSourceCorrection on;

	// Interpolation method
    interpolationSchemes
    {
        rho         	cell;
        U           	cellPoint;
        mu          	cell;
    }

	// Integration method
    integrationSchemes
    {
        U               analytical; //Euler;
    }
	
	// Relaxation
    sourceTerms
    {
        schemes
        {
           U semiImplicit 0.5; //the number is relaxCoeff for the field
        }
    }
}
// FLOW PROPERTIES
constantProperties
{
    rho0            $rhoPar; 
    alphaMax        1;
}
// PARTICLE SUBMODELS
subModels
{
	// Acting forces
    particleForces
    {
        sphereDrag;
        
        gravity;	
    }
	
	// Injection models
    injectionModels
    {
		/*
		model1
        {
            type            patchInjection;
            parcelBasisType mass;
            patchName       atmosphere;    //patch for injecting particles
            U0             (0 0 -1.544);   //terminal fall velocity
            sizeDistribution               //particle size distribution
            {
                type        normal;
                normalDistribution
                {
                    expectation 0.0002;
                    variance 0.00005;
                    minValue 5e-5;
                    maxValue 0.0005;
                }
            }
            flowRateProfile constant 1;   //flow rate
            massTotal       3543.45;      //kg. Mass injected in "duration"
            SOI             0;
            duration        3600;         //s. Time duration where massTotal is injected. 
            parcelsPerSecond 0.05e6;      //number of particles per s.
        } 
        */
    }
	// Particle dispersion model
    dispersionModel none;

	// Rebound-splash of snow particles
    patchInteractionModel   localInteractionStickReboundSplash; //none;
    
     localInteractionStickReboundSplashCoeffs
    {
	
        patches
        (
            atmosphere
            {
                type 	escape;
                e    	1.0;
                mu   	0.0;
            }
            
            // Snow surface settings
            snowBed  
            {
			// Rebound coefficients
			type 		stickReboundSplash;
			Pm			0.9;	 //(-) max probability of rebound.
			gamma		2.0;	 //(-) rebound
			
			// Splash entrainment coefficients
			epsilonr	0.25;	 //(-) energy balance		
			mur			0.5;	 //(-) momentum balance	
			muf			0.4;	 //(-) momentum balance	
			bEne	   	10e-9;	 //(-) bed cohesion
			corrm		0.0;	 //(-) mass-velocity correlation
			corre		0.0;	 //(-) mass-velocity correlation
			pppMax		$pppMax; //maximum number of particles per parcel
			
			// Particle diameter properties
			dm			$dm; 	 //m. mean particle diameter
			ds			$ds; 	 //m. std deviation of diameter
			d_max		$d_max;  //m. maximum particle diameter
			d_min		$d_min;  //m. minimum particle diameter
            }
        );
        // Particle probability distributions
		sizeDistributionTriple
		{
			type        normalLogNormalExponential;
			normalLogNormalExponentialDistribution
			{
				//nothing 
			}
		}         
    }
    heatTransferModel none;

    surfaceFilmModel none;

    collisionModel none;

    stochasticCollisionModel none;

    radiation off;
    
    // Aerodynamic entrainment of snow particles
    bedAerodynamicLiftInjectionModel  logLawShearStress; //none;
    logLawShearStressCoeffs
    {
		// Particle probability distributions
		sizeDistributionTriple
		{
			type        normalLogNormalExponential;
			normalLogNormalExponentialDistribution
			{
				//nothing 
			}
		};
		// Particle and coefficient settings
		aerodynamicLiftPatch    snowBed;    //patch for entrainement
		tauLogLaw       		false; 		//shear stress computation method
		dm						$dm; 		//m. mean particle diameter.
		ds						$ds; 		//m. std deviation of diameter
		d_max					$d_max;  	//m. maximum particle diameter
		d_min					$d_min;  	//m. minimum particle diameter
		z0          			$Z0; 		//m. aerodynamic roughness.
		pppMin					$pppMin;    //number of particles per parcel.
		SOA             		100; 		//s. start of activation.
		Acst            		0.2;        //(-). fluid threshold coefficient.
	}
}

cloudFunctions
{

}

// ************************************************************************* //