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

Is steady state solver grid dependent?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 21, 2013, 07:43
Default Is steady state solver grid dependent?
  #1
Member
 
M Mallikarjuna Reddy
Join Date: Jul 2012
Posts: 91
Rep Power: 5
mmkr825 is on a distinguished road
Hi Foamers,
I am doing simulations on 2D rectangular channel by using steady state solver. I prepared my solver by modifying "simpleFoam" solver. The dimensions (length*breadth) of the channel are 0.143*0.0017 (m). My geometry has 500*10 grids in x- and y- directions respectively. When i run my solver with 500*10 grids, it is giving good results. But when I refine grids (say 600*15) or increase the dimensions of the channel, after running few iterations it is giving the following error.

Error:
Code:
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3   in "/lib/i386-linux-gnu/libm.so.6"
#4  pow in "/lib/i386-linux-gnu/libm.so.6"
#5  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#6  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) at /opt/openfoam211/src/OpenFOAM/lnInclude/GeometricScalarField.C:273
#7  
 at /opt/openfoam211/src/OpenFOAM/lnInclude/GeometricScalarField.C:349
#8  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#9  
 in "/home/malli_reddy/OpenFOAM/malli_reddy-2.1.1/platforms/linuxGccDPOpt/bin/SBM25Foam"
Floating point exception
The geometry pass the checkMesh test. Could someone please help me to understand the reason for the error. Thanks in advance.

Regards
Reddy
mmkr825 is offline   Reply With Quote

Old   March 22, 2013, 09:33
Default
  #2
Senior Member
 
Fumiya Nozaki
Join Date: Jun 2010
Location: Yokohama, Japan
Posts: 170
Rep Power: 7
fumiya is on a distinguished road
Could you upload your case?

Best regards,
Fumiya
fumiya is offline   Reply With Quote

Old   March 22, 2013, 09:44
Default
  #3
Member
 
M Mallikarjuna Reddy
Join Date: Jul 2012
Posts: 91
Rep Power: 5
mmkr825 is on a distinguished road
Hi Fumiya,
Thanks for reply. My case is simple 2D channel.
My blockMeshDict:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.001;

vertices
(
    (0 0 0)
    (173 0 0)
    (173 1.7 0)
    (0 1.7 0)
    (0 0 0.1)
    (173 0 0.1)
    (173 1.7 0.1)
    (0 1.7 0.1)
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (1000 10 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    leftWall
    {
        type patch;
        faces
        (
            (0 4 7 3)
        );
    }
    fixedWalls
    {
        type wall;
        faces
        (
            (7 6 2 3)
            (1 5 4 0)
        );
    }
    rightWall
    {
    type patch;
    faces
    (
        (2 6 5 1)
    );
    }
    frontAndBack
    {
        type empty;
        faces
        (
            (0 3 2 1)
            (4 5 6 7)
        );
    }
);

mergePatchPairs
(
);

// ************************************************************************* //
Regards
Mallikarjuna
mmkr825 is offline   Reply With Quote

Old   March 22, 2013, 11:56
Default
  #4
Senior Member
 
Fumiya Nozaki
Join Date: Jun 2010
Location: Yokohama, Japan
Posts: 170
Rep Power: 7
fumiya is on a distinguished road
Hi Mallikarjuna,

Your mesh seems to be OK.
Could you share the files in 0 and system directories?

Best regards,
Fumiya
fumiya is offline   Reply With Quote

Old   March 22, 2013, 13:46
Default
  #5
Member
 
M Mallikarjuna Reddy
Join Date: Jul 2012
Posts: 91
Rep Power: 5
mmkr825 is on a distinguished road
Hi Fumiya,
Very much thanks for your kindness. My solver is highly non-linear. In my solver i encounter shear rate(magnitude of strain tensor). So my solver needs good initial guess for velocity. So i am attaching my solver and case file with initial guess for velocity.

Solver:
https://www.dropbox.com/s/tm8i4y2u0a...M25Foam.tar.gz

Case File:
https://www.dropbox.com/s/bz7u10bc16urpdc/case.tar.gz

I am very thankful to you if you guide me in correct path. Thanks in advance.

Regards
Mallikarjuna.
mmkr825 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
Solver for transonic flow? Martin Hegedus OpenFOAM Running, Solving & CFD 21 April 27, 2015 11:02
error message cuteapathy CFX 14 March 20, 2012 07:45
steady state version of twophaseeulerfoam miles_davis OpenFOAM 15 November 2, 2011 04:28
coding a steady multiphase solver twophaseeulerfoam miles_davis OpenFOAM Programming & Development 4 October 15, 2011 19:49
grid adaption for unsteady solver matty FLUENT 3 July 18, 2005 13:16


All times are GMT -4. The time now is 06:47.