CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

rhoRotorSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 23, 2009, 09:25
Default rhoRotorSimpleFoam
  #1
New Member
 
Marcin Hinz
Join Date: May 2009
Location: Monheim am Rhein, Germany
Posts: 24
Rep Power: 8
Marcin is on a distinguished road
Hi everyone!

I have tried to make a new solver based on the rhoSimpleFoam with the name below. The idea is to calculate with this solver rotor blades of a turbine, so an additional velocity therm is needed.
I have copied all files and made some modifications. I've added a vector field Uabs in the rhoRotorSimleFoam.C file, added the centrifugal and Coriolis forces in the UEqn.H ( rho*SRF->Su() ) and created a single rotating frame model in the createFields.H
After I have compiled it I didn't get any bugs, but when I wonted to start a simulation of a rotor blade I have got this bug:

Attempt to cast type wall to type lduInterface#0 Foam::error:rintStack(Foam::Ostream&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::Ostream& Foam:perator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rhoRotorSimpleFoam"
#3 Foam::lduInterface const& Foam::refCast<Foam::lduInterface const, Foam::fvPatch const>(Foam::fvPatch const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::cyclicFvPatchField<double>::cyclicFvPatchFie ld(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::cyclicFvPatchField<double> >::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rhoRotorSimpleFoam"
#7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::GeometricB oundaryField(Foam::fvBoundaryMesh const&, Foam:imensionedField<double,Foam::volMesh> const&, Foam::dictionary const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rhoRotorSimpleFoam"
#8 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readField(Foam::Istream&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rhoRotorSimpleFoam"
#9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rhoRotorSimpleFoam"
#10 Foam::basicThermo::basicThermo(Foam::fvMesh const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#11 Foam::hThermo<Foam:ureMixture<Foam::constTranspo rt<Foam::specieThermo<Foam::hConstThermo<Foam:er fectGas> > > > >::hThermo(Foam::fvMesh const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#12 Foam::basicThermo::addfvMeshConstructorToTable<Foa m::hThermo<Foam:ureMixture<Foam::constTransport< Foam::specieThermo<Foam::hConstThermo<Foam:erfec tGas> > > > > >::New(Foam::fvMesh const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#13 Foam::basicThermo::New(Foam::fvMesh const&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#14 main in "/home/marcin/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rhoRotorSimpleFoam"
#15 __libc_start_main in "/lib64/libc.so.6"
#16 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116

From function refCast<To>(From&)
in file /home/marcin/OpenFOAM/OpenFOAM-1.5/src/OpenFOAM/lnInclude/typeInfo.H at line 106.


Does anyone know what it means? Thanks for any help!

Best regards, Marcin.
Marcin is offline   Reply With Quote

Old   June 23, 2009, 13:33
Default
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 248
Rep Power: 9
jaswi is on a distinguished road
Hi Marcin

There seems to be some clash in your boundary conditions.
I am no experts but in your error message it says

Quote:
Attempt to cast type wall to type lduInterface#0 Foam::error:rintStack(Foam::Ostream&) in "/home/marcin/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
Now lduInterface is an abstract base class for implicitly-coupled interfaces e.g. processor and cyclic patches and above message implies that you wall patch is being recasted into lduInterface type. One more question - are you using AMG solver in you fvSolution file. Try changing it to some other solver and see what happens.

There might be problem also if you are using Cylic Boundary condition on your field variable for a patch which is defined as wall in your polyMesh/boundary file.

Please cross check you case setup and if nothing works then try running your case without the modified solver, just to find out whether error is due to your change or due to case setup.

Best Regards
Jaswi
jaswi is offline   Reply With Quote

Old   June 24, 2009, 03:26
Default
  #3
New Member
 
Marcin Hinz
Join Date: May 2009
Location: Monheim am Rhein, Germany
Posts: 24
Rep Power: 8
Marcin is on a distinguished road
Hi Jaswi,

thanks a lot for your reply. Just to anwser your question I'm using only PCG and PBiCG solver, so it isn't caused with the AMG.
In the fact I have verified my boudary conditions and the boundary file and I found the reason thanks to you. The problem was that by mesh import from icem the cyclic boudary condition was lost. I have fixed it in the polyMesh/boudary file and now it seems to work.

Thanks a lot once more!

Best regards, Marcin.
Marcin is offline   Reply With Quote

Reply

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



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