vitorspadeto |
June 18, 2014 16:04 |
groovyBC error: velocity profile (2D) >> what's wrong?
groovyBC error
Hi everybody!
I'm trying to implement a parabolic velocity profile to my inlet (2D), cordinates: (x,z) . Inlet is located between the coordinates z= 0 and z=0.001. the parabolic equation for this profile is
Velocity_x_direction = 99*z - 98072*z².
This equation in this range of points (between z= 0 and z=0.001) give the exact velocity profile that I need.
I already add "libs" lines in the end of ControlDict file (I'm using swak4Foam to use GroovyBC):
Code:
/*--------------------------------*- 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 "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application icoFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 10;
deltaT 0.1;
writeControl adjustableRunTime;
writeInterval 1;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
adjustTimeStep yes;
maxCo 1;
maxAlphaCo 1;
maxDeltaT 1;
libs (
"libOpenFOAM.so"
"libsimpleSwakFunctionObjects.so"
"libswakFunctionObjects.so"
"libgroovyBC.so"
);
// ************************************************************************* //
And this is my U file, with the profile velocity:
Code:
/*--------------------------------*- 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 volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type groovyBC;
variables "zp=pts().z; velocity=(99.036*zp-98072.289*pow(zp,2));";
valueExpression "velocity";
value uniform (0.02 0 0);
}
outlet
{
type zeroGradient; // At the outlet the velocity is unknown so we take the gradient of velocity as zero here
}
wall
{
type fixedValue;
value uniform (0 0 0); // we use the no slip condition at the walls so u = v = 0
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //
This is the error that appeared on the terminal:
HTML Code:
sp@sp-945GCM-S2C:~/PGSPAD/exemplos_para_o_PG/blood/Newtonian$ icoFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.0-5be49240882f
Exec : icoFoam
Date : Jun 18 2014
Time : 17:01:52
Host : "sp-945GCM-S2C"
PID : 4579
Case : /home/sp/PGSPAD/exemplos_para_o_PG/blood/Newtonian
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading transportProperties
Reading field p
Reading field U
Reading/calculating face flux field phi
Starting time loop
Time = 0.1
Courant Number mean: 0.00555556 max: 10
word::stripInvalid() called for word velocity
For debug level (= 2) > 1 this is considered fatal
Abortado (imagem do núcleo gravada)
What Should I do to correct this error?
THANKS FRIENDS!
|