/*--------------------------------*- C++ -*----------------------------------*| =========				 |												 |
| \	  /  F ield		 | OpenFOAM: The Open Source CFD Toolbox		   |
|  \	/   O peration	 | Version:  5									 |
|   \  /	A nd		   | Web:	  www.OpenFOAM.org					  |
|	\/	 M anipulation  |												 |
\*---------------------------------------------------------------------------*/

// User-defined parameters

Wflume 0.914;

initialVel 0.5;

SWL 0.15182;

pref 0.0;

// Turbulence Calcs

// p_rghIC  #calc "$pref- (0.5*$rho*$initialVel*$initialVel)"; // [Pa] dynamic pressure, p_rgh=p-1/2*rho*U^2
p_rghIC  0.0; // [Pa] dynamic pressure, p_rgh=p-1/2*rho*U^2

L_REF	0.914; // [m]

U_REF   0.25; // [m/s]

INTENSITY   0.25;

C_mu 0.09; // [unitless] 


kIC #eval "(3*($INTENSITY*$U_REF)*($INTENSITY*$U_REF))/2"; // [m^2/s^2] turbulent kinetic energy, TKE or k

omegaIC #eval "(1/$L_REF)*pow($kIC,0.5)/pow($C_mu,0.25)"; // [s^-1] specific turbulence frequency, omega
epsilonIC #eval "(1/$L_REF)*pow($kIC,1.5)*pow($C_mu,0.75)"; // [m^2/s^3] TKE dissipation rate, epsilon




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