/*--------------------------------*- 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 @!dropRadius!@ ;
De         @!De!@;
Ca         @!Ca!@;
Re         @!Re!@;
Co         @!Co!@;
eta2       @!eta2!@; //matrix
rho2	   @!rho2!@;
m          @!m!@;

cellsPerRadius @!cpr!@;

X      #eval {10.0 * $dropRadius};
Y      #eval {8.0  * $dropRadius};
Z      #eval {4.0*$dropRadius};

nx     #eval {10.0 * $cellsPerRadius};
ny     #eval {8.0 * $cellsPerRadius};
nz     #eval {4.0 * $cellsPerRadius};

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};
Cz      #eval {$Z/2.0};



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