CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   OF 2.x, cyclic boundary and MRFSimpleFoam compatibility (http://www.cfd-online.com/Forums/openfoam-solving/95936-2-x-cyclic-boundary-mrfsimplefoam-compatibility.html)

jmf January 8, 2012 14:34

OF 2.x, cyclic boundary and MRFSimpleFoam compatibility
 
Hello dear forum

Did anybody experiment OF 2.1.0 with cyclicAMI boundaries and MRFSimpleFoam ? Is it supposed to work ?

I used to do that with OF 1.5-dev. I'm testing OF 2.1.0, from RPM packages for SUSE.

I get following error :

Code:

Reading field p

Reading field U

AMI: Creating addressing and weights between 10342 source faces and 13525 target faces
AMI: Patch source weights min/max/average = 6.7415422e-14, 1.2787833, 1.0003003
AMI: Patch target weights min/max/average = 0, 1.2746978, 1.0002249
Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib64/libc.so.6"
#3  Foam::FaceCellWave<Foam::wallPoint, int>::updateFace(int, Foam::wallPoint const&, double, Foam::wallPoint&) in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#4  Foam::FaceCellWave<Foam::wallPoint, int>::handleAMICyclicPatches() in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#5  Foam::FaceCellWave<Foam::wallPoint, int>::iterate(int) in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#6  Foam::patchWave::correct() in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#7  Foam::patchWave::patchWave(Foam::polyMesh const&, Foam::HashSet<int, Foam::Hash<int> > const&, bool) in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#8  Foam::wallDist::correct() in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9  Foam::wallDist::wallDist(Foam::fvMesh const&, bool) in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10  Foam::incompressible::RASModels::kOmegaSST::kOmegaSST(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/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#11  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kOmegaSST>::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/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#12  Foam::incompressible::RASModel::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/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#13 
 in "/opt/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/MRFSimpleFoam"
#14  __libc_start_main in "/lib64/libc.so.6"
#15 
 at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
Exception en point flottant

Does anybody have a suggestion ?

Thanks in advance, and best wishes

Jean-Michel


PS : many thanks and congratulations to OF maintainers team for the last evolutions of the software : AMI, edges in SHM, packaging

danishdude January 9, 2012 15:24

I configured a cyclicAMI case last weekend (using rotationally cyclic boundary conditions), and it seemed to work smoothly. Granted, I'm having convergence issues with MRFSimpleFoam, but I'm fairly certain it's a setup issue and has nothing to do with the AMI implementation (It's looking like the boundary patches need to be further offset from the rotating features).

Can you provide more detail on your case? Perhaps post the case files if they aren't too big?

jmf January 9, 2012 16:22

[SOLVED] cyclicAMI + MRFSimpleFoam crash
 
Hi Mickael


Thanks for your attention. Finally I found the fix on the forum. Running following command did the job :

Code:

unset FOAM_SIGFPE
Don't ask me why

This forum is very good. Thanks to everybody


Regards

jean-Michel

buct11019 May 23, 2012 10:10

Quote:

Originally Posted by jmf (Post 338527)
Hi Mickael


Thanks for your attention. Finally I found the fix on the forum. Running following command did the job :

Code:

unset FOAM_SIGFPE
Don't ask me why

This forum is very good. Thanks to everybody


Regards

jean-Michel

It doesen't work...just there was not the warning ,but can't calculate.

buct11019 May 24, 2012 04:01

Quote:

Originally Posted by buct11019 (Post 362688)
It doesen't work...just there was not the warning ,but can't calculate.

I think it's because you over estimate the boundary values of k and epsilon

vahid.najafi July 29, 2012 01:51

plz tell me it is true???or not???
 
hi Dear Foamers.
I want to add k(kinetic turbulence energy) in one solver.

for this reason added this line in code :

const volScalarField &k=U_.db().lookupObject<volScalarField>("k")
and wmake was Successfully .

but How can I Understand that this k is the same with k(kinetic turbulence energy)???Is it true?????

Because I replase another world in this :

const volScalarField &M=U_.db().lookupObject<volScalarField>("M")
but not error occured!!!!!!??????


All times are GMT -4. The time now is 19:30.