[hsingtzu]

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

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

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

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

[jferrari]

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.

[hsingtzu]

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.

[jferrari]

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].