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

icoFoam does not work on different geometry

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 30, 2011, 14:43
Default icoFoam does not work on different geometry
  #1
New Member
 
Hsingtzu Wu
Join Date: Feb 2011
Posts: 21
Rep Power: 5
hsingtzu is on a distinguished road
Hello

I was trying to apply icoFoam on different geometry, but it seems not working.
I borrowed the "cylinder" polyMesh from tutorials/basic/potentialFoam and did some changes. I changed the type of cylinder to patch and changed the 0 folder accordingly.

my p file looks like

Code:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    down            
    {
        type            symmetryPlane;
    }

    right           
    {
        type            fixedValue;
        value           uniform 0;
    }

    up              
    {
        type            symmetryPlane;
    }

    left            
    {
        type            zeroGradient;
    }

    cylinder        
    {
        type            zeroGradient;
    }

    defaultFaces    
    {
        type            empty;
    }
}
and my U file looks like

Code:
internalField   uniform (1 0 0); 

boundaryField
{
    down
    {
        type            symmetryPlane;
    }
    right
    {
        type            zeroGradient;
    }
    up
    {
        type            symmetryPlane;
    }
    left
    {
        type            fixedValue;
        value           uniform (1 0 0);
    }
    cylinder
    {
        type            zeroGradient;
    }
    defaultFaces
    {
        type            empty;
    }
}
I borrowed the transportProperties file and the system folder from tutorial/incompressible/icoFoam/cavity
Then I ran it and got the following error message

Code:
Time = 0.005

Courant Number mean: 7.1527 max: 5.99709e+298
#0 Foam::error::printStack(Foam::Ostream&) in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/lib/libOpenFOAM.so”
#1 Foam::sigFpe::sigHandler(int) in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/lib/libOpenFOAM.so”
#2 Uninterpreted:
#3 void Foam::fvc::surfaceIntegrate<Foam::Vector<double> >(Foam::Field<Foam::Vector<double> >&, Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/lib/libfiniteVolume.so”
#4 Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/lib/libfiniteVolume.so”
#5 Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/lib/libfiniteVolume.so”
#6 Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<Foam::Vector<double> >(Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/lib/libfiniteVolume.so”
#7 Foam::fv::gaussLaplacianScheme<Foam::Vector<double>, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/lib/libfiniteVolume.so”
#8
in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/bin/icoFoam”
#9
in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/bin/icoFoam”
#10
in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/bin/icoFoam”
#11
in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/bin/icoFoam”
#12 __libc_start_main in “/lib/tls/i686/cmov/libc.so.6″
#13
in “/home/uiuc/OpenFOAM/uiuc-2.0.1/platforms/linuxGccDPOpt/bin/icoFoam”
Floating point exception
I was wondering if anyone could kindly help me with this issue.
Thanks
hsingtzu is offline   Reply With Quote

Old   December 10, 2011, 18:52
Default
  #2
Member
 
Joe
Join Date: Dec 2011
Location: Groton, CT
Posts: 67
Rep Power: 4
jferrari is on a distinguished road
A few things I noticed about your case. First, the units for pressure seem incorrect, looks like meters^2/seconds^2. It should be [1 -1 -2 0 0 0 0] for pressure. Second, viscous flow over a cylinder isn't symmetric. Third, and I'm not sure what's causing this, but your courant number is exploding.

Sorry I can't help more - I have about two weeks of experience with OpenFOAM and can't get my own cylinder case to run. Hopefully someone with more experience than me will come along and help the both of us.
jferrari is offline   Reply With Quote

Old   December 14, 2011, 12:23
Default
  #3
New Member
 
Hsingtzu Wu
Join Date: Feb 2011
Posts: 21
Rep Power: 5
hsingtzu is on a distinguished road
Thanks for the comment. I really appreciate your time and help.
I am working on my final project which is due this Sun. This means that I will work on this problem next week and reply to your comment.
hsingtzu is offline   Reply With Quote

Old   December 14, 2011, 19:33
Default
  #4
Member
 
Joe
Join Date: Dec 2011
Location: Groton, CT
Posts: 67
Rep Power: 4
jferrari is on a distinguished road
After learning a bit more I realize that I was incorrect about the dimensions. Apparently for pressure, OpenFOAM uses kinematic pressure, or pressure divided by density. That would have dimensions [0 2 -2 0 0 0 0].
jferrari 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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Inclined Geometry TommyLui FLUENT 0 November 7, 2011 05:32
[ICEM] O-grid around wing geometry Anorky ANSYS Meshing & Geometry 12 August 17, 2011 22:28
Simulation of Flow through Complex 3D Geometry EmersonKB CFX 5 July 2, 2009 09:17
[GAMBIT] import terrain geometry cfdworker ANSYS Meshing & Geometry 2 May 24, 2009 22:40
Rotor 67 full compressor geometry jf Main CFD Forum 0 January 24, 2005 07:59


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