CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

OF 2.x, cyclic boundary and MRFSimpleFoam compatibility

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   January 8, 2012, 13:34
Default OF 2.x, cyclic boundary and MRFSimpleFoam compatibility
  #1
jmf
Member
 
J Michel
Join Date: Aug 2009
Location: Orleans - France
Posts: 32
Rep Power: 6
jmf is on a distinguished road
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
jmf is offline   Reply With Quote

Old   January 9, 2012, 14:24
Default
  #2
New Member
 
Michael Ahlmann
Join Date: Feb 2010
Posts: 27
Rep Power: 6
danishdude is on a distinguished road
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?
danishdude is offline   Reply With Quote

Old   January 9, 2012, 15:22
Default [SOLVED] cyclicAMI + MRFSimpleFoam crash
  #3
jmf
Member
 
J Michel
Join Date: Aug 2009
Location: Orleans - France
Posts: 32
Rep Power: 6
jmf is on a distinguished road
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
jmf is offline   Reply With Quote

Old   May 23, 2012, 10:10
Default
  #4
New Member
 
Forrest
Join Date: May 2012
Posts: 5
Rep Power: 3
buct11019 is on a distinguished road
Quote:
Originally Posted by jmf View Post
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 is offline   Reply With Quote

Old   May 24, 2012, 04:01
Default
  #5
New Member
 
Forrest
Join Date: May 2012
Posts: 5
Rep Power: 3
buct11019 is on a distinguished road
Quote:
Originally Posted by buct11019 View Post
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
buct11019 is offline   Reply With Quote

Old   July 29, 2012, 01:51
Default plz tell me it is true???or not???
  #6
Member
 
vahid
Join Date: Feb 2012
Location: Mashhad-Iran
Posts: 80
Rep Power: 3
vahid.najafi is an unknown quantity at this point
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!!!!!!??????
vahid.najafi is offline   Reply With Quote

Reply

Tags
cyclic boundaries, mrfsimplefoam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
MRFSimpleFoam: wrong boundary conditions on rotating walls cves OpenFOAM Bugs 3 May 3, 2010 03:38
case with cyclic boundary openfoam_user OpenFOAM Running, Solving & CFD 0 March 26, 2010 08:46
cyclic BC - createPatch AirS OpenFOAM 7 November 12, 2009 10:04
Cyclic boundary conditions axisymmetric gabriel OpenFOAM Running, Solving & CFD 0 July 27, 2007 11:24
View cyclic boundaries using MRFSimpleFoam gabriel OpenFOAM Paraview & paraFoam 0 July 27, 2007 08:42


All times are GMT -4. The time now is 21:05.