CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   icoFoam does not work on different geometry (http://www.cfd-online.com/Forums/openfoam-solving/94944-icofoam-does-not-work-different-geometry.html)

hsingtzu November 30, 2011 14:43

icoFoam does not work on different geometry
 
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

jferrari December 10, 2011 18:52

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 December 14, 2011 12:23

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 December 14, 2011 19:33

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


All times are GMT -4. The time now is 20:18.