/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      parametersDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// L length of the domain
// Y height of the domain
// 1 droplet fluid
// 2 matrix fluid
//
// Non-dimensional numbers
// Re = rho2* shearRate * dropRadius^2 / eta2
// Ca = eta2 * shearRate * dropRadius / sigma
// m = eta1 / eta2
// shearRate = (Ut - Ub) / Y (Ut, Ub) are the boundary velocity at top and bottom respectively
// eta1 = etaS + etaP
// etaS = etaP
// 
// This all lead to
// U = Re * eta2* Y / (2.0*rho2*dropRadius^2)
// sigma = eta2^2 * Re / (rho2*dropRadius*Ca)
// Reported Time-> endTime = 3*(shearRate^(-1)) -> 3*Y/(2*U)


dropRadius 0.001 ;
De         0.4;
Ca         0.24;
Re         0.0003;
Co         0.18;
eta2       6.6; //matrix
rho2       1000.0;
m          1;

cellsPerRadius 50;

X      #eval {20.0 * $dropRadius};
Y      #eval {8.0  * $dropRadius};
Z      #eval {0.1*$dropRadius};

nx     #eval {20.0 * $cellsPerRadius};
ny     #eval {8.0 * $cellsPerRadius};
nz     #eval {1.0};

eta1  #eval {$m*$eta2};
etaS  #eval {$eta2/2.0};
etaP  #eval {$eta2/2.0};

shearRate  #eval {$Re*$eta2/$rho2/$dropRadius/$dropRadius};

U       #eval {$shearRate*$Y/2};

pLambda #eval {$De/$shearRate};

sigma   #eval {$eta2*$shearRate*$dropRadius/$Ca};

//Drop initial position
Cx      #eval {$X/2.0};

Cy      #eval {$Y/2.0};



//X      #eval {20.0 * 0.001};
//Y      #eval {8.0  * 0.001};
//Z      #eval {0.1 * 0.001};

//nx     #eval {20.0 * 50};
//ny     #eval {8.0 *  50};
//nz     #eval {1.0};

//eta1  #eval {1*6.6};
//etaS  #eval {6.6/2.0};
//etaP  #eval {6.6/2.0};

//shearRate  #eval {0.0003*6.6/1000.0/0.001/0.001};

//U       #eval {(0.0003*6.6/1000.0/0.001/0.001)*(8.0*0.001)/2};

//pLambda #eval {0.4/(0.0003*6.6/1000.0/0.001/0.001)};

//sigma   #eval {6.6*(0.0003*6.6/1000.0/0.001/0.001)*0.001/0.24};

//Drop initial position
//Cx      #eval {(20.0 * 0.001)/2.0};

//Cy      #eval {(8.0 * 0.001)/2.0};


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