[gschaider]

Originally Posted by egp

I am getting an error when trying to use the fpos functionality in swak4Foam>funkySetFields.

funkySetFields -latestTime
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-ext-943ecc426649
Exec : funkySetFields -latestTime
Date : Mar 25 2011
Time : 04:58:35
Host : egpMBA.local
PID : 70817
Case : /Users/egp/OpenFOAM/egp-1.6-ext/run/RayleighTaylor/coarse
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 1

Part: Alpha
Putting "average((fpos().y+surf(0.5)*fproj().y)>=surf(-0.05)*cos(surf(6.2831)*fpos().x) ? surf(1.0) : surf(0.))" into field alpha1 at t = "1" if condition "true" is true

--> FOAM FATAL ERROR:
Parser Error at "1.9-12" :"syntax error, unexpected TOKEN_fposition"
"average((fpos().y+surf(0.5)*fproj().y)>=surf(-0.05)*cos(surf(6.2831)*fpos().x) ? surf(1.0) : surf(0.))"
" ^^^^ "

From function parsingValue
in file lnInclude/CommonValueExpressionDriverI.H at line 365.

FOAM exiting

This is strange since swak4FoamParsers/FieldValueExpressionDriver.C creates an IOobject named "fpos", and the machinery appears to largely be in place for surface field operations.

Note: this funkySetFieldsDict works just fine with the stand-alone FSF.

Hi Eric!

The error popped up here just yesterday too. If you check you'll see that instead of average the function sum won't work too but integrate and reconstruct will still work like the old-school FSF.

The problem is (haven't verified yet but am quite sure) that I introduced another average/sum into the grammar that reduces the whole field to one number and that obviously "shadows" the old average. To avoid any further confusion I will rename the old average (the one you need. Sorry for the incompatibility) to fvAverage and push the change to the mercurial-repository

Bernhard

[gschaider]

Originally Posted by gschaider

Hi Eric!

The error popped up here just yesterday too. If you check you'll see that instead of average the function sum won't work too but integrate and reconstruct will still work like the old-school FSF.

The problem is (haven't verified yet but am quite sure) that I introduced another average/sum into the grammar that reduces the whole field to one number and that obviously "shadows" the old average. To avoid any further confusion I will rename the old average (the one you need. Sorry for the incompatibility) to fvAverage and push the change to the mercurial-repository

While trying to implement the fix I discovered that it is already there (no need to pull anything): the "old" average is now faceAverage. Old sum is now surfSum. So the only thing that remains is documenting it. That will probably take longer

[egp]

Thanks, Bernhard. That fixed my problem!!

Originally Posted by gschaider

While trying to implement the fix I discovered that it is already there (no need to pull anything): the "old" average is now faceAverage. Old sum is now surfSum. So the only thing that remains is documenting it. That will probably take longer

[mabinty]

Dear all,

I m experiencing a problem with the division operator "/". Any time i use it to define a heat source the simulations crashes. The expression looks like:

variables    "Xmin=-0.15;Ymin=0.0;Zmin=-0.15;Xmax=0.15;Ymax=0.3;Zmax=0.15;power=27000;V=(Xmax-Xmin)*(Ymax-Ymin)*(Zmax-Zmin);powerDens=power/V;";

I tried various possibilities but I always get a similar error message (see below) except when using a numerical value in the denominator, e.g. powerDens=power/0.027. I use OF 1.6.x, bison 2.4 and flex 2.5.35.

Thanks in advance for your help!
Aram

Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 9.83737e-28, Final residual = 9.83737e-28, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 1.12552e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 7.26874e-12, Final residual = 7.26874e-12, No Iterations 0
[1] #0[3]   Foam::error::printStack(Foam::Ostream&)#0  Foam::error::printStack(Foam::Ostream&)[0] [2] #0  #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigFpeHandler(int) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #1  Foam::sigFpe::sigFpeHandler(int) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[2] #1  Foam::sigFpe::sigFpeHandler(int) in  in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigFpeHandler(int)"/home/aa/OpenFOAM/OpenFOA in "M-1.6.x/l/home/aa/ib/linux64OGpenFOAMccDPOpt/l/OpenFOAM-1.6ibOpenFO.x/lib/liAM.so"
[1] #2  nux64GccDPOpt/libOpenFOAM.so"
[2] #2  ?? in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[0] #2   in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #2  ?? in "/lib/libc.so.6"
[1] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib/libc.so.6"
[2] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&)???? in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #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 "/lib/libc.so.6"
[3] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
[1] #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 "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[2] #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 "/home/aa/OpenFO in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"A
[1] #6M/OpenFOAM-1.6.x/  lib/linparserField::FieldValueExpressionParser::parse()ux64GccDPOpt/libOpenFOAM.so"
[3] #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 "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
[2] #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 "/lib/libc.so.6"
[0] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
[2] #6  parserField::FieldValueExpressionParser::parse() in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
[3] #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 "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[1] #7  Foam::FieldValueExpressionDriver::parse(std::string const&) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[2] #7  Foam::FieldValueExpressionDriver::parse(std::string const&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
[3] #6  parserField::FieldValueExpressionParser::parse() in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[2] #8  Foam::CommonValueExpressionDriver::evaluateVariable(Foam::word const&, Foam::string const&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
[0] #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 "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[2] #9  Foam::CommonValueExpressionDriver::addVariables(Foam::string const&, bool) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[1] #8  Foam::CommonValueExpressionDriver::evaluateVariable(Foam::word const&, Foam::string const&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
[0] #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 "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[3] #7  Foam::FieldValueExpressionDriver::parse(std::string const&) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[1] #9  Foam::CommonValueExpressionDriver::addVariables(Foam::string const&, bool) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[2] #10  Foam::CommonValueExpressionDriver::addVariables(Foam::List<Foam::string> const&, bool) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
[0] #6  parserField::FieldValueExpressionParser::parse() in "/home/aa/Ope in "/home/aa/OpenFOAM/aa-1.6.x/lib/nFOAM/aa-1.6.x/linux64GccDPOpt/libswak4FoamParlib/linux64GccDPOpt/libswak4FoamParsers.so"
[3] #8  sers.so"
[1] #10  Foam::CommonValueExpressionDriver::addVariables(Foam::List<Foam::string> const&, bool)Foam::CommonValueExpressionDriver::evaluateVariable(Foam::word const&, Foam::string const&) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[2] #11   in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[1] #11  Foam::expressionSource<double>::operator()()Foam::expressionSource<double>::operator()() in "/h in "/home/aa/OpenFOAM/aa-1.6.x/lib/ in "/holme/aa/OpenFOAM/aa-1.6.x/liib/linux64GccDPOpt/libswak4FoamParsers.so"
[0] #7  Foam::FieldValueExpressionDriver::parse(std::string const&)nux64GccDPOpt/libswakSourceFields.so"
[2] #12  ome/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[3] #9  Foam::CommonValueExpressionDriver::addVariables(Foam::string const&, bool) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswakSourceFields.so"
[1] #12   in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[0] #8  Foam::CommonValueExpressionDriver::evaluateVariable(Foam::word const&, Foam::string const&)main in "/home/aa/OpenFOAM in "/ho/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[3] #me/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[0] #9  Foam::CommonValueExpressionDriver::addVariables(Foam::string const&, bool)10  Foam::CommonValueExpressionDriver::addVariables(Foam::List<Foam::string> const&, bool)main in "/home/aa/OpenFOAM/aa-1.6.x/lib/linu in "/homex/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[0] #10  Foam::CommonValueExpressionDriver::addVariables(Foam::List<Foam::string> const&, bool)64GccDPOpt/libswak4FoamParsers.so"
[3] #11  Foam::expressionSource<double>::operator()() in "/home/aa/OpenFOAM/aa-1.6.x/applications/bin/linux64GccDPOpt/chtMultiRegionFoamMF"
[2] #13  __libc_start_main in "/home/aa/OpenFOAM/aa-1.6.x/applications/bin/linux64GccDPOpt/chtMultiRegionFoamMF"
[1] #13  __libc_start_main in "/home/aa/OpenFOAM/aa-1.6.x/lib/lin in "/home/aa/OpenFOAM/aa-1.6.ux64GccDPOx/lib/linux64GccDPOpt/libswak4FoamParsers.so"
[0] #11  Foam::expressionSource<double>::operator()()pt/libswakSourceFields.so"
[3] #12   in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswakSourceFields.so"
[0] #12   in "/lib/libc.so.6"
[2] #14   in "/lib/libc.so.6"
[1] #14  main_startmain_start in "/home/aa/OpenFOAM/aa-1.6.x/applications/bin/linux64GccDPOpt/chtMultiRegionFoamMF"
[0] #13  __libc_start_main in "/home/aa/OpenFOAM/aa-1.6.x/applications/bin/linux64G at /build/buildd/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
[lws16:21853] *** Process received signal ***
[lws16:21853] Signal: Floating point exception (8)
[lws16:21853] Signal code:  (-6)
[lws16:21853] Failing at address: 0x3e80000555d
[lws16:21853] [ 0] /lib/libc.so.6 [0x7f23f762e040]
[lws16:21853] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7f23f762dfb5]
[lws16:21853] [ 2] /lib/libc.so.6 [0x7f23f762e040]
[lws16:21853] [ 3] /home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0xc1) [0x7f23f85daea1]
[lws16:21853] [ 4] /home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9_+0xd8) [0x7f23fa2615a8]
[lws16:21853] [ 5] /home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_SA_+0x214) [0x7f23fa26de74]
[lws16:21853] [ 6] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN11parserField26FieldValueExpressionParser5parseEv+0x6e5d) [0x7f23f8bc4fcd]
[lws16:21853] [ 7] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam26FieldValueExpressionDriver5parseERKSs+0x5d) [0x7f23f8ccef9d]
[lws16:21853] [ 8] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam27CommonValueExpressionDriver16evaluateVariableERKNS_4wordERKNS_6stringE+0x22) [0x7f23f8cbcb72]
[lws16:21853] [ 9] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam27CommonValueExpressionDriver12addVariablesERKNS_6stringEb+0x481) [0x7f23f8cbe7c1]
[lws16:21853] [10] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam27CommonValueExpressionDriver12addVariablesERKNS_4ListINS_6stringEEEb+0x39) [0x7f23f8cbecd9]
[lws16:21853] [11] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswakSourceFields.so(_ZN4Foam16expressionSourceIdEclEv+0x2d) [0x7f23f893c72d]
[lws16:21853] [12] chtMultiRegionFoamMF [0x44c8e9]
[lws16:21853] [13] /lib/libc.so.6(__libc_start_main+0xe6) [0x7f23f76195a6]
[lws16:21853] [14] chtMultiRegionFoamMF [0x423dd9]
[lws16:21853] *** End of error message ***
 in "/lib/libc.so.6"
[0] #14   at /build/buildd/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
[lws16:21852] *** Process received signal ***
[lws16:21852] Signal: Floating point exception (8)
[lws16:21852] Signal code:  (-6)
[lws16:21852] Failing at address: 0x3e80000555c
[lws16:21852] [ 0] /lib/libc.so.6 [0x7fbdc5113040]
[lws16:21852] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7fbdc5112fb5]
[lws16:21852] [ 2] /lib/libc.so.6 [0x7fbdc5113040]
[lws16:21852] [ 3] /home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0xc1) [0x7fbdc60bfea1]
[lws16:21852] [ 4] /home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9_+0xd8) [0x7fbdc7d465a8]
[lws16:21852] [ 5] /home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_SA_+0x214) [0x7fbdc7d52e74]
[lws16:21852] [ 6] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN11parserField26FieldValueExpressionParser5parseEv+0x6e5d) [0x7fbdc66a9fcd]
[lws16:21852] [ 7] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam26FieldValueExpressionDriver5parseERKSs+0x5d) [0x7fbdc67b3f9d]
[lws16:21852] [ 8] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam27CommonValueExpressionDriver16evaluateVariableERKNS_4wordERKNS_6stringE+0x22) [0x7fbdc67a1b72]
[lws16:21852] [ 9] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam27CommonValueExpressionDriver12addVariablesERKNS_6stringEb+0x481) [0x7fbdc67a37c1]
[lws16:21852] [10] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so(_ZN4Foam27CommonValueExpressionDriver12addVariablesERKNS_4ListINS_6stringEEEb+0x39) [0x7fbdc67a3cd9]
[lws16:21852] [11] /home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswakSourceFields.so(_ZN4Foam16expressionSourceIdEclEv+0x2d) [0x7fbdc642172d]
[lws16:21852] [12] chtMultiRegionFoamMF [0x44c8e9]
[lws16:21852] [13] /lib/libc.so.6(__libc_start_main+0xe6) [0x7fbdc50fe5a6]
[lws16:21852] [14] chtMultiRegionFoamMF [0x423dd9]
[lws16:21852] *** End of error message ***
ccDPOpt/chtMultiRegionFoamMF"
[3] #13  __libc_start_main--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 21853 on node lws16 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

[gschaider]

Originally Posted by mabinty

Dear all,

I m experiencing a problem with the division operator "/". Any time i use it to define a heat source the simulations crashes. The expression looks like:

variables    "Xmin=-0.15;Ymin=0.0;Zmin=-0.15;Xmax=0.15;Ymax=0.3;Zmax=0.15;power=27000;V=(Xmax-Xmin)*(Ymax-Ymin)*(Zmax-Zmin);powerDens=power/V;";

I tried various possibilities but I always get a similar error message (see below) except when using a numerical value in the denominator, e.g. powerDens=power/0.027. I use OF 1.6.x, bison 2.4 and flex 2.5.35.

Thanks in advance for your help!
Aram

Hm. That should work. The stack-trace you posted is from a parallel run. Does this also happen in a serial run (that would help me to tell what the problem is). A little bug-report at the Mantis mentioned in the first post of this thread would also help me to not forget to look into it

Bernhard

[mabinty]

Hi Bernhard!

The problem occurs also in single runs (see error message below). I currently face some troubles with my sourceforge account but as soon as i resolve it i report the bug on Mantis.

Thanks for your help!

Aram

Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 9.40997e-28, Final residual = 9.40997e-28, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 1.10761e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 7.29066e-12, Final residual = 7.29066e-12, No Iterations 0
#0  Foam::error::printStack(Foam::Ostream&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  ?? in "/lib/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/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 "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.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 "/home/aa/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#6  parserField::FieldValueExpressionParser::parse() in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
#7  Foam::FieldValueExpressionDriver::parse(std::string const&) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
#8  Foam::CommonValueExpressionDriver::evaluateVariable(Foam::word const&, Foam::string const&) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
#9  Foam::CommonValueExpressionDriver::addVariables(Foam::string const&, bool) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
#10  Foam::CommonValueExpressionDriver::addVariables(Foam::List<Foam::string> const&, bool) in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswak4FoamParsers.so"
#11  Foam::expressionSource<double>::operator()() in "/home/aa/OpenFOAM/aa-1.6.x/lib/linux64GccDPOpt/libswakSourceFields.so"
#12  main in "/home/aa/OpenFOAM/aa-1.6.x/applications/bin/linux64GccDPOpt/chtMultiRegionFoamMF"
#13  __libc_start_main in "/lib/libc.so.6"
#14  _start at /build/buildd/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116

[pablodecastillo]

Hello

I would like to apply this BC in the inlet of my domain and update every iteration because it is moving.

So;
"average((fpos().y + surf(0.5) * fproj().y) < (surf(0.)) ? surf(1.0) : ((fpos().y - surf(0.5) * fproj().y) > surf(0.0) ? surf(0.0) : (surf(0.5) - fpos().y / (fproj().y + surf(0.00000001)))))"

is it possible to translate to groobybc??, i can not find in groovybc fproj and surf functions.

Advanced thanks

[gschaider]

Originally Posted by pablodecastillo

Hello

I would like to apply this BC in the inlet of my domain and update every iteration because it is moving.

So;
"average((fpos().y + surf(0.5) * fproj().y) < (surf(0.)) ? surf(1.0) : ((fpos().y - surf(0.5) * fproj().y) > surf(0.0) ? surf(0.0) : (surf(0.5) - fpos().y / (fproj().y + surf(0.00000001)))))"

is it possible to translate to groobybc??, i can not find in groovybc fproj and surf functions.

Advanced thanks

Reason is that on a patch you only have faces (no cells). Therefore by definition all values are surf (no need for that function). Also pos IS the fpos. The only thing that might be missing is fproj. Maybe for a first try you can do without it. Or you can use pts and toFace/toPoint to achieve something similar (have a look at the groovyBC-docu)

[pablodecastillo]

Thanks,

It should be something like this;

inlet
{
type groovyBC;
variables "yp=pts().y;minY=min(yp);maxY=max(yp);proj_y=m axY-minY;";
valueExpression "(pos().y + 0.5 * proj_y) < 0.0 ? 1.0 : (pos().y - 0.5 * proj_y) > 0.0 ? 0.0 : (0.5 - pos().y / proj_y)";
value uniform (0 0 0);
}

right?

[gschaider]

Originally Posted by pablodecastillo

Thanks,

It should be something like this;

inlet
{
type groovyBC;
variables "yp=pts().y;minY=min(yp);maxY=max(yp);proj_y=m axY-minY;";
valueExpression "(pos().y + 0.5 * proj_y) < 0.0 ? 1.0 : (pos().y - 0.5 * proj_y) > 0.0 ? 0.0 : (0.5 - pos().y / proj_y)";
value uniform (0 0 0);
}

right?

I'm afraid proj is not equivalent to the projection you'd like to have (the one fproj calculated). To rebuild this using existing functionality is not as trivial as I thought (that's why it was built into funky in the first place) - but I havn'T looked at your expression. Use replayTransientBC to see if it "looks right"

Other suggestions:

a. Live with a jagged interface at the inlet
b. Try combinations of pos() and toFace(pts()) to get a blurred edge
c. Use the expressionField-function object to get a field with the fproj() at every time-step and use the in your expression (not sure whehter this works correctly because the field you get will be a volumeField)
d. Write a bug-report asking for fproj for patches and hope that it gets done in the next few weeks

[pablodecastillo]

Hello Bernhard,

You are right, it is not working. When calculate minY and maxY, it is doing for all the inlet. I thought that would take just points for every face.

I will try to think another thing, if i can not find a easy solution, it will be better create a new BC.

Thanks

[pablodecastillo]

Hi again,

pts is defined like pts Vector field with the vertices.

So how can i known how many vertices has the Face, and his coordinates?

Thanks.

[gschaider]

Originally Posted by pablodecastillo

Hi again,

pts is defined like pts Vector field with the vertices.

So how can i known how many vertices has the Face, and his coordinates?

Thanks.

Not with groovyBC. Sorry. That would add way too much complexity to the parser and if that kind of information is needed then programming the BC in C++ is almost the same effort (its hard to process that kind of information without loops anyway)

[albem]

Could you please post your file?

Originally Posted by xisto

I already use both.

funkySetFields and funkySetBoundaryFilelds

And it works for me.

Excellent tool.

Thanks.

Carlos

[mabinty]

dear all!

I m using swal4Foam in a multi region case where I introduced (like shown in the example InterFoamWithSources) a heat source in one of the regions, i.e air, by defining the heatSourceDict as follows:

variables "Xmin=-0.15;Ymin=0.0;Zmin=-0.15;Xmax=0.15;Ymax=0.3;Zmax=0.15;powerDens=1.0e06;"; //= 27e3/0.3^3 W/m3

expression   "(pos().x>=Xmin && pos().x<=Xmax && pos().y>=Ymin && pos().y<=Ymax && pos().z >= Zmin && pos().z <= Zmax) ? powerDens : 0";

dimensions   [1 -1 -3 0 0 0 0];

That work quite well. Now I want to extend the domain of the heat source by establishing a cellSet "heatSource" in the considered region air and use the following heatSourceDict:

variables "powerDens=1.0e06;"; //= 27e3/0.3^3 W/m3

expression   "set(heatSource) ? powerDens : 0";

dimensions   [1 -1 -3 0 0 0 0];

My main problem now is how to create a cellSet "heatSource" in the region air applying the cellRegionAddressing of the corresponding region. I tried to use:

setSet -batch makeCellSets.heatSource -region air

where makeCellSets.heatSource defines the desired cellSet "heatSource", after establishing the region air. The utility produces a cellSet constant/air/polyMesh/sets/heatSource but when I let the solver run I get the following message in the log file:

Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.000182075, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.000208165, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.000173218, No Iterations 1
No set "heatSource" at t=0.1 falling back to 'constant'

and when I visualize the heatSource in paraview the location doesn t correspond to the one defined by the makeCellSets.heatSource file.

Does anybody have an idea how to solve this problem??

Appreciate your comments!

Aram

[lrunber]

Hi All,
I have been trying to set up swak4Foam on ubuntu 11.04. This is what I have done:
Installed packages:
swak4foam-openfoam-1.6-ext
swak4foam-openfoam-1.6-ext-dev

Next I wanted to compile InterFoamWithSources, I changed the these lines in ./Make/options, according to where the packages where installed
from:

    -I../../../Libraries/swak4FoamParsers/lnInclude \
    -I../../../Libraries/swakFiniteArea/lnInclude

to:

    -I$(LIB_SRC)/swak4Foam/swak4FoamParsers/lnInclude \
    -I$(LIB_SRC)/swak4Foam/swakSourceFields/lnInclude

now when I run wmake:

Making dependency list for source file interFoamWithSources.C
SOURCE=interFoamWithSources.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-40 -I/usr/lib/OpenFOAM-1.6-ext/src/transportModels -I/usr/lib/OpenFOAM-1.6-ext/src/transportModels/incompressible/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/transportModels/interfaceProperties/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/turbulenceModels/incompressible/turbulenceModel -I/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swak4FoamParsers/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swakSourceFields/lnInclude -IlnInclude -I. -I/usr/lib/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/interFoamWithSources.o
In file included from /usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swak4FoamParsers/lnInclude/FieldValueExpressionParser_location.hh:44:0,
                 from FieldValueExpressionParser.tab.hh:65,
                 from /usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swak4FoamParsers/lnInclude/FieldValueExpressionDriver.H:18,
                 from /usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swakSourceFields/lnInclude/expressionSource.H:45,
                 from interFoamWithSources.C:48:
position.hh: In function ‘bool parserField::operator==(const parserField::position&, const parserField::position&)’:
position.hh:136:69: warning: suggest parentheses around ‘&&’ within ‘||’
In file included from interFoamWithSources.C:60:0:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:134:5: error: invalid initialization of reference of type ‘const Foam::dictionary&’ from expression of type ‘Foam::volScalarField’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/findRefCell.H:50:6: error: in passing argument 2 of ‘void Foam::setRefCell(const Foam::volScalarField&, const Foam::dictionary&, Foam::label&, Foam::scalar&, bool)’
In file included from interFoamWithSources.C:91:0:
pEqn.H:34:23: error: ‘struct Foam::volScalarField’ has no member named ‘select’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:11:10: warning: unused variable ‘transonic’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:14:9: warning: unused variable ‘nOuterCorr’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nCorr’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:8:10: warning: unused variable ‘momentumPredictor’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:11:10: warning: unused variable ‘transonic’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:14:9: warning: unused variable ‘nOuterCorr’
make: *** [Make/linux64GccDPOpt/interFoamWithSources.o] Error 1

Could someone please help put this fool out of his misery

[gschaider]

Originally Posted by lrunber

Hi All,
I have been trying to set up swak4Foam on ubuntu 11.04. This is what I have done:
Installed packages:
swak4foam-openfoam-1.6-ext
swak4foam-openfoam-1.6-ext-dev

Next I wanted to compile InterFoamWithSources, I changed the these lines in ./Make/options, according to where the packages where installed
from:

    -I../../../Libraries/swak4FoamParsers/lnInclude \
    -I../../../Libraries/swakFiniteArea/lnInclude

to:

    -I$(LIB_SRC)/swak4Foam/swak4FoamParsers/lnInclude \
    -I$(LIB_SRC)/swak4Foam/swakSourceFields/lnInclude

now when I run wmake:

Making dependency list for source file interFoamWithSources.C
SOURCE=interFoamWithSources.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-40 -I/usr/lib/OpenFOAM-1.6-ext/src/transportModels -I/usr/lib/OpenFOAM-1.6-ext/src/transportModels/incompressible/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/transportModels/interfaceProperties/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/turbulenceModels/incompressible/turbulenceModel -I/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swak4FoamParsers/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swakSourceFields/lnInclude -IlnInclude -I. -I/usr/lib/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude -I/usr/lib/OpenFOAM-1.6-ext/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/interFoamWithSources.o
In file included from /usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swak4FoamParsers/lnInclude/FieldValueExpressionParser_location.hh:44:0,
                 from FieldValueExpressionParser.tab.hh:65,
                 from /usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swak4FoamParsers/lnInclude/FieldValueExpressionDriver.H:18,
                 from /usr/lib/OpenFOAM-1.6-ext/src/swak4Foam/swakSourceFields/lnInclude/expressionSource.H:45,
                 from interFoamWithSources.C:48:
position.hh: In function ‘bool parserField::operator==(const parserField::position&, const parserField::position&)’:
position.hh:136:69: warning: suggest parentheses around ‘&&’ within ‘||’
In file included from interFoamWithSources.C:60:0:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:134:5: error: invalid initialization of reference of type ‘const Foam::dictionary&’ from expression of type ‘Foam::volScalarField’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/findRefCell.H:50:6: error: in passing argument 2 of ‘void Foam::setRefCell(const Foam::volScalarField&, const Foam::dictionary&, Foam::label&, Foam::scalar&, bool)’
In file included from interFoamWithSources.C:91:0:
pEqn.H:34:23: error: ‘struct Foam::volScalarField’ has no member named ‘select’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:11:10: warning: unused variable ‘transonic’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:14:9: warning: unused variable ‘nOuterCorr’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nCorr’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:8:10: warning: unused variable ‘momentumPredictor’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:11:10: warning: unused variable ‘transonic’
/usr/lib/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/readPISOControls.H:14:9: warning: unused variable ‘nOuterCorr’
make: *** [Make/linux64GccDPOpt/interFoamWithSources.o] Error 1

Could someone please help put this fool out of his misery

The problem is that the include paths "-I../../ ..." assume that the example is in the same source tree as the Library sources. With the debian package this is not the case I'm afraid. There are two possible workarounds:
- Use some symbolic links to "emulate" an environment where the Headers ARE in ../..
- Modify Make/options so that -I points to the actual place of the headers (I believe the dev package installs them into $FOAM_SRC)

Sorry for the mixup but I'm not using Debian/Ubuntu myself. If anyone fells like maintaining the Debian-packaging he's more than welcome.

[gschaider]

Originally Posted by gschaider

The problem is that the include paths "-I../../ ..." assume that the example is in the same source tree as the Library sources. With the debian package this is not the case I'm afraid. There are two possible workarounds:
- Use some symbolic links to "emulate" an environment where the Headers ARE in ../..
- Modify Make/options so that -I points to the actual place of the headers (I believe the dev package installs them into $FOAM_SRC)

Sorry for the mixup but I'm not using Debian/Ubuntu myself. If anyone fells like maintaining the Debian-packaging he's more than welcome.

Sorry. My bad. You already did this and the problem you have is something different. InterFoam changed a bit in the last time. The example solver was based on 1.7 which is slighly incompatible with 1.6. The best guess would be to take the sources of the original 1.6 solver and put in the 4 or 5 lines that are needed for the sources.

[gschaider]

Originally Posted by mabinty

dear all!

I m using swal4Foam in a multi region case where I introduced (like shown in the example InterFoamWithSources) a heat source in one of the regions, i.e air, by defining the heatSourceDict as follows:

variables "Xmin=-0.15;Ymin=0.0;Zmin=-0.15;Xmax=0.15;Ymax=0.3;Zmax=0.15;powerDens=1.0e06;"; //= 27e3/0.3^3 W/m3

expression   "(pos().x>=Xmin && pos().x<=Xmax && pos().y>=Ymin && pos().y<=Ymax && pos().z >= Zmin && pos().z <= Zmax) ? powerDens : 0";

dimensions   [1 -1 -3 0 0 0 0];

That work quite well. Now I want to extend the domain of the heat source by establishing a cellSet "heatSource" in the considered region air and use the following heatSourceDict:

variables "powerDens=1.0e06;"; //= 27e3/0.3^3 W/m3

expression   "set(heatSource) ? powerDens : 0";

dimensions   [1 -1 -3 0 0 0 0];

My main problem now is how to create a cellSet "heatSource" in the region air applying the cellRegionAddressing of the corresponding region. I tried to use:

setSet -batch makeCellSets.heatSource -region air

where makeCellSets.heatSource defines the desired cellSet "heatSource", after establishing the region air. The utility produces a cellSet constant/air/polyMesh/sets/heatSource but when I let the solver run I get the following message in the log file:

Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.000182075, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.000208165, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.000173218, No Iterations 1
No set "heatSource" at t=0.1 falling back to 'constant'

and when I visualize the heatSource in paraview the location doesn t correspond to the one defined by the makeCellSets.heatSource file.

Does anybody have an idea how to solve this problem??

Appreciate your comments!

Aram

The message you highlighted is not the problem. It just says "because I'm at t=0.1 I had a look in 0.1/air/polyMesh/sets and the cell set wasn't there so I'm going for the one in constant" which is actually what you want anyway.

About the headSource being in the wrong place I have no idea. Have you tried to visualize it without using swak? Either in paraFoam (I think that supports sets) or if you're using the native reader (I think that only supports zones) by converting it with setsToZones. If THAT is in the right place then it is a swak-problem otherwise I'd guess it is a problem with the setup

[lrunber]

Originally Posted by gschaider

Sorry. My bad. You already did this and the problem you have is something different. InterFoam changed a bit in the last time. The example solver was based on 1.7 which is slighly incompatible with 1.6. The best guess would be to take the sources of the original 1.6 solver and put in the 4 or 5 lines that are needed for the sources.

Thanks, I see the problem was the example solver based on 1.7. I have 1.7 installed so I copied the /usr/lib/OpenFOAM-1.6-ext/src/swak4Foam to /opt/openfoam171/src and ran Allclean & Allwmake.
Now InterFoamWithSources compiles fine, however, the utilities e.g. funkySetFields are only in the 1.6-ext bin, and since the src was not included in the ubuntu packages I proceeded to:

DL the utilities from: https://openfoam-extend.svn.sourcefo...oam/Utilities/
and proceed to paste into: /opt/openfoam171/applications/utilities/swak4Foam
After running wmake on the utilities (had to change /Make/options as in previous post), I can succesfully run mixingThing case.

Thanks for the great programs