Turbulent inlet of mass flow (from csv file)
 

Hi Foamers. I need a turbulent mass flow inlet for a time varying flow rate.
I tried to apply the following BC for a U profile (Solver:twoLiquidMixingFoam).
in
{
type turbulentInlet;
fluctuationScale (0.02 0.01 0.01);
referenceField flowRateInletVelocity;

massFlowRate table
(
( 0.0000000e+001 0.0000000e+000 )
( 0.0100000e+001 5.2890505e-003 )
( 0.050000e+001 1.0507891e-001 )
( 0.1030000e+001 1.5657144 )
( 0.2040000e+001 2.5737434 )
( 0.30000e+001 3.5749383 )
);
rho phase1;
rhoInlet 842.346;
}
No errors are generated starting the simulation, but seems as the inlet field does not feel the perturbation. Hierarchical errors?


In addiction I discover that using the csv data input as shown in the following BC
in
{
type flowRateInletVelocity;
massFlowRate csvFile;
csvFileCoeffs
{
fileName "FlowRate2.dat";
outOfBounds clamp;
hasHeaderLine false;
refColumn 0;
componentColumns (1);
}
rho phase1;
rhoInlet 842.346;

}
The simulation has no problem to start on 1 processor,
Reading g
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 204 source faces and 204 target faces
AMI: Patch source weights min/max/average = 1, 1, 1
AMI: Patch target weights min/max/average = 1, 1, 1
Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type LESModel
Selecting LES turbulence model Vreman
Selecting LES delta type cubeRootVol
VremanCoeffs
{
ce 1.048;
ck 0.094;
}

Calculating field g.h

Courant Number mean: 0.00132949 max: 0.2065

PIMPLE: Operating solver in PISO mode


Starting time loop

Courant Number mean: 0.00132949 max: 0.2065
Interface Courant Number mean: 0 max: 0
deltaT = 0.0111111
Time = 0.0111111

MULES: Solving for alpha1
Phase-1 volume fraction = 0.00147721 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.00295443 Min(alpha1) = 0 Max(alpha1) = 1
DICPCG: Solving for alpha1, Initial residual = 0.0995769, Final residual = 4.06332e-15, No Iterations 1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 3.11085e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.01e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 2.80578e-07, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 3.81469e-07, No Iterations 36
time step continuity errors : sum local = 1.12702e-09, global = -8.37911e-11, cumulative = -8.37911e-11
GAMG: Solving for p_rgh, Initial residual = 0.00724602, Final residual = 7.52881e-07, No Iterations 13
time step continuity errors : sum local = 4.74598e-07, global = 1.64522e-07, cumulative = 1.64439e-07
ExecutionTime = 1.15 s ClockTime = 1 s

but if i try to run it on more than one processor (4 in this case) i obtain the following error:

Reading g
Reading field p_rgh

AMI: Creating addressing and weights between 204 source faces and 204 target faces
AMI: Patch source weights min/max/average = 1, 1, 1
AMI: Patch target weights min/max/average = 1, 1, 1
Reading field U

[3]
[3]
[3] --> FOAM FATAL IO ERROR:
[3] wrong token type - expected int, found on line 46 the word '^A'
[3]
[3] file: /scratch/villoresi/TEST/processor3/0/U.boundaryField.in.csvFileCoeffs.componentColumns at line 46.
[3]
[3] From function operator>>(Istream&, int&)
[3] in file primitives/ints/int/intIO.C at line 68.

Why??

P.S. Yeah, i would like to use swak4foam but i cannot being that i cannot update bison for admin reasons.