[jmf]

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 :

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]

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]

Hi Mickael


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

unset FOAM_SIGFPE

Don't ask me why

This forum is very good. Thanks to everybody


Regards

jean-Michel

[buct11019]

Originally Posted by jmf

Hi Mickael


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

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]

Originally Posted by buct11019

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]

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!!!!!!??????