groovyBC issue - k and epsilon
 

Hi, I am trying to implement groovyBC for k and epsilon for a domain inlet located in between Xmin = -0.254 to Xmax = 3.556. The relevant k and epsilon files are attached for reference. The k and epsilon values are always positive for the all pos().x. The groovyBC works for Ux and Uy velocities for that inlet. However its the k and epsilon that's creating troubles. When I use a fixedValue for epsilon (value uniform 0.01) at the boundary with groovyBC function for k, then it gives the following error:

--------------------------------------------------------------

Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model realizableKE
bounding k, min: 0 max: 0.59375 average: 0.1593749999999999
realizableKECoeffs
{
Cmu 0.09;
A0 4;
C2 1.9;
sigmak 1;
sigmaEps 1.2;
}

No field sources present


SIMPLE: convergence criteria
field p tolerance 1e-6
field U tolerance 1e-6
field "(k|epsilon|omega)" tolerance 1e-6


Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.01221319187944816, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.000175166776251753, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.9999999993231693, Final residual = 0.01273879170881757, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.9999999960428109, Final residual = 0.09699403429229032, No Iterations 154
time step continuity errors : sum local = 0.05527557721138961, global = -6.764283616009572e-11, cumulative = -6.764283616009572e-11
DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 0.001747110343310934, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::realizableKE::cor rect() in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#7
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/simpleFoam"
#8 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#9
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)

--------------------------------------------------------------------------------


When I use a fixedValue for k (value uniform 0.1) at the boundary with the groovyBC function for epsilon, then it gives the following error:


------------------------------------------------------------------------------------------

Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model realizableKE
bounding epsilon, min: 0 max: 5.88 average: 5.879999999999997
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::realizableKE::rCm u(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#7 Foam::incompressible::RASModels::realizableKE::rCm u(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#8 Foam::incompressible::RASModels::realizableKE::rea lizableKE(Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#9 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::rea lizableKE>::New(Foam::GeometricField<Foam::Vector< double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#10 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#11
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/simpleFoam"
#12 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#13
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)

----------------------------------------------------------------------------------------------------


When I used fixedValue for k and epsilon at the boundaries, the simulation runs fine.


I have read forum discussions like "http://www.cfd-online.com/Forums/openfoam-solving/71495-groovybc-k-profile.html" but still could not figure out the issue. Even using a groovyBC with constant values does not help for k and epsilon like:

type groovyBC;
value uniform 0.01;
variables "a1=0.01;kcalc=a1;";
valueExpression "kcalc";


Any inputs would be of great help. Thanks a lot.

Sagnik

Looks like a division by 0

The line "bounding epsilon" gives a very clear hint here: somewhere epsilon is 0 and therefor the problem with a "division by 0"

The attached BCs make me think of John von Neumann "With four parameters I can fit an elephant, and with five I can make him wiggle his trunk.". Anyway: have you checked that with your parameters and with the range of x these polynomials don't become smaller than 0?

The simple function you use at last should work. Which version of swak/groovyBC do you use?

Bernhard

We are using the latest version: swakVersion: 0.2.1 (Release date: 2012-10-18)

OpenFOAM version is 2.1.1

The controlDict file calls the following libs:
"libOpenFOAM.so"
"libsimpleSwakFunctionObjects.so"
"libswakFunctionObjects.so"
"libgroovyBC.so"

All k and epsilon values are positive. In fact

type groovyBC;
value uniform 0.01;
variables "a1=0.01;kcalc=a1;";
valueExpression "kcalc";

also doesn't work. I believe it should. Any other suggestions on things that we should check/verify would be of great help !

That is really strange. I don't see why this shouldn't work. In fact for instance the Example/groovyBC/pulsedPitzDaily uses a groovyBC for k at the inlet. Could you set up a small case that demonstrates the problem?

Hi, I have put 2 cases to demonstrate the problem. The cases can be downloaded from:

http://dl.dropbox.com/u/58859393/Case_Demos.zip


For Test_groovyBC_1, I use a

type fixedValue;
value uniform 0.1;

condition for "inletFront" and the case runs.


For Test_groovyBC_2, I use a

type groovyBC;
variables "a1k=0.1;kcalc=a1k;";
valueExpression "kcalc";
value uniform 0.01;

condition for "inletFront" and the case throws out the error I mentioned in the very first iteration. All other parameters are the same for the test cases.

Note that I implemented groovyBC for the velocity as well which seems to work fine.

The geometry, mesh etc. doesn't make too much sense. Its just for testing.

It would be great if you can suggest a way out.

Thanks for all the help.

Sagnik

After fixing a bug (but a different one - related to some work I did recently) both of your cases run OK. I'm not sure what might have been the problem but between the last release and now I did some refactoring and that might have squashed the bug.

The working version is in the public development repository. How to download that is described here: http://openfoamwiki.net/index.php/Co...am#Development

Thanks a lot. We have tested quite a few cases in the past few days and it works perfectly fine.

groovy bc
 

hey every one
I am a new user of the code OpenFoam
I found that groovy bc is an important tool but I still have problems with it
here is my example
I need to program an epsilon profile in my inlet
and the equation of this profil is the following:

epsilon=Ustar^3/(k(y+y0))

how can I do it?
I really need help to have some progress in my thesis..

best regards

Lamia

Not sure what you're trying to do here. For instance: with "an epsilon profile in my inlet" you mean "set epsilon at the inlet as a function of the position"? If that is the case then have a look at the pulsedPitzDaily and other example cases where a quantity is set as a function of the position

Otherwise: please explain the terms in your expression. Not everyone is familiar with your nomenclatur

thank's for your answer
I can make things more clear
I simply need to programm the privious equation with groovy Bc
I tried this:
inlet
{
type groovyBC;
variables "u_f=0.43;y0=0.0007,y=5;epsi=pow(U_f,3)/k(pos( ).y+y0);"
valueExpression "epsi";
}
but it gives me an error message
I don't know why!
and where can I find a tutorial about groovyBC?

Error message would have been helpful. But there are two things problematic here that I see at the first glance:
- the "," before y=5 (which is not needed anyway)
- Case is important: either write U_f or u_f

Tutorial: have you checked out my presentation from the Workshop at PSU in 2011 (it is linked from the swak-Wiki-page)? That is the closest thing to a tutorial that I wrote (don't know if there are any others out there)

Hi Bernhard,
I am interested to use groovyBC within the "immersed boundary method" which was released during the seventh openfoam workshop. I want to apply groovyBC for one of my variables at this boundary. Since the immersed boundary was not defined as a patch, the groovyBC is not able to understand it, do you have any suggestions if I am able to use groovyBC for this case?

Thanks

Haven't looked at the IB-stuff yet so I can't say anything clever here. If it is implemented by fixing cell values and these values are stored in a field then the manipulateField-functionObject might be the way to go (although ugly). Which still leaves the question how the information where the boundary cells are is stored. If it is in a cellZone or cellSet then swak can access it.

But that is strictly without looking at the source

Thanks Bernhard, it seems the implementation is by fixing cell values with Bi-linear interpolation. I will look in detail to your suggestions. I have used the groovyBC for a boundary condition in suspended sediment transport and it is working well for a fixed and movable body-fitted boundary. I am moving to an Immersed boundary method and I am interested to use it if it can work.

Hi,

I compiled swak4foam as post#6 but still having some errors. I faced similar problem before and following post#6 worked on my system. However, recently I compiled swak4foam over a cluster but the simulation is not working. OF version is 2.1.1.

-----------------------------
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
bounding k, min: 0 max: 28 average: 20
bounding omega, min: 0 max: 10 average: 10
kOmegaSSTCoeffs
{
alphaK1 0.85034;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.85616;
gamma1 0.5532;
gamma2 0.4403;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}


Starting time loop

Time = 0.0001

[40] [41] #0 Foam::error::printStack(Foam::Ostream&)#0 Foam::error::printStack(Foam::Ostream&)--------------------------------------------------------------------------
An MPI process has executed an operation involving a call to the
"fork()" system call to create a child process. Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your MPI job may hang, crash, or produce silent
data corruption. The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

Local host: orc214 (PID 10423)
MPI_COMM_WORLD rank: 40

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
addr2line failed
[41] #1 Foam::sigFpe::sigHandler(int) addr2line failed
[40] #1 Foam::sigFpe::sigHandler(int) addr2line failed
[41] #2 addr2line failed
[40] #2 addr2line failed
[41] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) addr2line failed
[40] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) addr2line failed
[41] #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) addr2line failed[40] #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) addr2line failed[41] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) addr2line failed[40] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) addr2line failed
[41] #6 Foam::incompressible::RASModels::kOmegaSST::F2() const addr2line failed
[40] #6 Foam::incompressible::RASModels::kOmegaSST::F2() const addr2line failed
[41] #7 Foam::incompressible::RASModels::kOmegaSST::F23() const addr2line failed
[40] #7 Foam::incompressible::RASModels::kOmegaSST::F23() const addr2line failed[41] #8 Foam::incompressible::RASModels::kOmegaSST::kOmega SST(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) addr2line failed[40] #8 Foam::incompressible::RASModels::kOmegaSST::kOmega SST(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) addr2line failed[41] #9 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kOm egaSST>::New(Foam::GeometricField<Foam::Vector<dou ble>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[40] #9 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kOm egaSST>::New(Foam::GeometricField<Foam::Vector<dou ble>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[41] #10 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[40] #10 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[41] #11 Foam::incompressible::turbulenceModel::addturbulen ceModelConstructorToTable<Foam::incompressible::RA SModel>::NewturbulenceModel(Foam::GeometricField<F oam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[40] #11 Foam::incompressible::turbulenceModel::addturbulen ceModelConstructorToTable<Foam::incompressible::RA SModel>::NewturbulenceModel(Foam::GeometricField<F oam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[41] #12 Foam::incompressible::turbulenceModel::New(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[40] #12 Foam::incompressible::turbulenceModel::New(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) addr2line failed[41] #13 addr2line failed[40] #13
[41][41] #14 __libc_start_main[40]
[40] #14 __libc_start_main addr2line failed[41] #15 addr2line failed[40] #15
[41][orc214:10424] *** Process received signal ***
[orc214:10424] Signal: Floating point exception (8)[orc214:10424] Signal code: (-6)
[orc214:10424] Failing at address: 0x1fe82000028b8[orc214:10424] [ 0] /lib64/libc.so.6(+0x32920) [0x2b0d072c2920][orc214:10424] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x2b0d072c28a5]
[orc214:10424] [ 2] /lib64/libc.so.6(+0x32920) [0x2b0d072c2920][orc214:10424] [ 3] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_ 5UListIdEES6_+0xc1) [0x2b0d0652c861][orc214:10424] [ 4] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam6divideINS_1 2fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldI dT_T0_EERKS6_S9_+0xd8) [0x2b0d04256cf8][orc214:10424] [ 5] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4FoamdvINS_12fvPa tchFieldENS_7volMeshEEENS_3tmpINS_14GeometricField IdT_T0_EEEERKS8_SA_+0x280) [0x2b0d0426f3a0][orc214:10424] [ 6] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZNK4Foam14incompres sible9RASModels9kOmegaSST2F2Ev+0x141) [0x2b0d04286c41][orc214:10424] [ 7] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZNK4Foam14incompres sible9RASModels9kOmegaSST3F23Ev+0x21) [0x2b0d04286ff1][orc214:10424] [ 8] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible9RASModels9kOmegaSSTC1ERKNS_14GeometricFieldIN S_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS3 _IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_14tra nsportModelERKNS_4wordESK_+0xeef) [0x2b0d0428942f][orc214:10424] [ 9] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible8RASModel31adddictionaryConstructorToTableINS0 _9RASModels9kOmegaSSTEE3NewERKNS_14GeometricFieldI NS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS 6_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_14tr ansportModelERKNS_4wordE+0x59) [0x2b0d0429b9d9][orc214:10424] [10] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible8RASModel3NewERKNS_14GeometricFieldINS_6Vector IdEENS_12fvPatchFieldENS_7volMeshEEERKNS2_IdNS_13f vsPatchFieldENS_11surfaceMeshEEERNS_14transportMod elERKNS_4wordE+0x36e) [0x2b0d0420780e][orc214:10424] [11] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible15turbulenceModel36addturbulenceModelConstruct orToTableINS0_8RASModelEE18NewturbulenceModelERKNS _14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldE NS_7volMeshEEERKNS5_IdNS_13fvsPatchFieldENS_11surf aceMeshEEERNS_14transportModelERKNS_4wordE+0x10) [0x2b0d04215cc0][orc214:10424] [12] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so(_ZN4Foam14inco mpressible15turbulenceModel3NewERKNS_14GeometricFi eldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEE RKNS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_ 14transportModelERKNS_4wordE+0x8a0) [0x2b0d03f3c270][orc214:10424] [13] pisoFoam() [0x416b7b]
[orc214:10424] [14] /lib64/libc.so.6(__libc_start_main+0xfd) [0x2b0d072aecdd][orc214:10424] [15] pisoFoam() [0x416429]
[orc214:10424] *** End of error message ***

[40][orc214:10423] *** Process received signal ***
[orc214:10423] Signal: Floating point exception (8)[orc214:10423] Signal code: (-6)[orc214:10423] Failing at address: 0x1fe82000028b7
[orc214:10423] [ 0] /lib64/libc.so.6(+0x32920) [0x2b5452d3e920][orc214:10423] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x2b5452d3e8a5][orc214:10423] [ 2] /lib64/libc.so.6(+0x32920) [0x2b5452d3e920][orc214:10423] [ 3] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_ 5UListIdEES6_+0xc1) [0x2b5451fa8861][orc214:10423] [ 4] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam6divideINS_1 2fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldI dT_T0_EERKS6_S9_+0xd8) [0x2b544fcd2cf8][orc214:10423] [ 5] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4FoamdvINS_12fvPa tchFieldENS_7volMeshEEENS_3tmpINS_14GeometricField IdT_T0_EEEERKS8_SA_+0x280) [0x2b544fceb3a0][orc214:10423] [ 6] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZNK4Foam14incompres sible9RASModels9kOmegaSST2F2Ev+0x141) [0x2b544fd02c41][orc214:10423] [ 7] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZNK4Foam14incompres sible9RASModels9kOmegaSST3F23Ev+0x21) [0x2b544fd02ff1][orc214:10423] [ 8] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible9RASModels9kOmegaSSTC1ERKNS_14GeometricFieldIN S_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS3 _IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_14tra nsportModelERKNS_4wordESK_+0xeef) [0x2b544fd0542f][orc214:10423] [ 9] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible8RASModel31adddictionaryConstructorToTableINS0 _9RASModels9kOmegaSSTEE3NewERKNS_14GeometricFieldI NS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS 6_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_14tr ansportModelERKNS_4wordE+0x59) [0x2b544fd179d9][orc214:10423] [10] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible8RASModel3NewERKNS_14GeometricFieldINS_6Vector IdEENS_12fvPatchFieldENS_7volMeshEEERKNS2_IdNS_13f vsPatchFieldENS_11surfaceMeshEEERNS_14transportMod elERKNS_4wordE+0x36e) [0x2b544fc8380e][orc214:10423] [11] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompress ible15turbulenceModel36addturbulenceModelConstruct orToTableINS0_8RASModelEE18NewturbulenceModelERKNS _14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldE NS_7volMeshEEERKNS5_IdNS_13fvsPatchFieldENS_11surf aceMeshEEERNS_14transportModelERKNS_4wordE+0x10) [0x2b544fc91cc0][orc214:10423] [12] /work/jubayer/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so(_ZN4Foam14inco mpressible15turbulenceModel3NewERKNS_14GeometricFi eldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEE RKNS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_ 14transportModelERKNS_4wordE+0x8a0) [0x2b544f9b8270][orc214:10423] [13] pisoFoam() [0x416b7b]
[orc214:10423] [14] /lib64/libc.so.6(__libc_start_main+0xfd) [0x2b5452d2acdd][orc214:10423] [15] pisoFoam() [0x416429]
[orc214:10423] *** End of error message ***--------------------------------------------------------------------------



Any help will be greatly appreciated.


Jubayer

Search the forum. I'm pretty sure this is a variation of the old favourite "I set k/epsilon/omega to 0 and now the division blows up" blues

I need help about groovy bc
 

hello everybody
simply I need to program an Omega profil for KOmega sst model as inlet boundary condition
I just ignore how to do it
while I have the equation of the profil of Omega

omega=alpha*u/sqrt(C).z

Is it possible to program that equation by using groovy bc?

Hi LamiaOF2.1,

I am attaching one of my omega file for your reference which uses groovyBC for the inlet profile. Hope this helps.


Jubayer

omega profil
 

thank you so much for your reply
I will try with this one
and let you know about the result

ps: I am a new user of OpenFoam and I am preparing my Ph.D in atmospheric pollution so I really need to contact someone who has some information about the code, I can say I had an interesting progress but I want more

if you are interested here is my personal adress
without_truth@hotmail.fr

groovy bc
 

hello
I am trying to program with groovy BC this equation

w=alpha*Uref/sqrt(Cµ)*Yref

I tried this but it gives me an error message
inlet
{
type groovyBC;
variables "Uref=5;Yref=0.5;alpha=0.33;Cµ=0.09;
omg=((Uref*alpha)/(sqrt(Cµ)*Yref));"
valueExpression "omg";
}

so is that make any sens?