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

[cfMesh] How can i get the proper result by cfMesh and icoFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 17, 2015, 16:45
Default How can i get the proper result by cfMesh and icoFoam
  #1
New Member
 
Choi Hak Kyu
Join Date: Apr 2015
Posts: 12
Rep Power: 7
prodo is on a distinguished road
Hello,

I made a mesh by cfMesh, like attached picture.

and then solve this case by icoFoam

initial conditon is v=0.1m/s

but, you can see there is non-physical result on the mesh grading region.

so, how can i get proper result in this case?

this is my fvSolution, fvSchemes file.

thank you

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
}

divSchemes
{
    default                 Gauss upwind;
    div(phi,U)                  Gauss linearUpwindV Gauss linear;
    div((nuEff*dev(grad(U).T())))     Gauss linear;

}

laplacianSchemes
{
    default         Gauss linear orthogonal;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default     corrected;
 }

fluxRequired
{
    default         no;
    p               ;
}


// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver           GAMG;
        tolerance        1e-7;
        relTol           0;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10;
        mergeLevels      1;

  /*    solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0;             */
    }

    U
    {


solver    smoothSolver;
        smoother    GaussSeidel;
        tolerance    1e-8;
        relTol        0;
        nSweeps        1;
/*        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0;            */
    }
}

PISO
{
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       0;
}
relaxationFactors
{
    p               0.3;
    U               0.7;
    k               0.7;
    epsilon         0.7;
    omega            0.7;
}

// ************************************************************************* //




Attached Images
File Type: jpg 1.jpg (61.5 KB, 17 views)
File Type: jpg 2.jpg (94.4 KB, 23 views)
File Type: jpg 3.jpg (11.3 KB, 22 views)
prodo is offline   Reply With Quote

Old   May 18, 2015, 08:24
Default
  #2
Senior Member
 
Franjo Juretic
Join Date: Aug 2011
Location: Velika Gorica, Croatia
Posts: 124
Rep Power: 13
franjo_j is on a distinguished road
Send a message via Skype™ to franjo_j
Hello,

What happens when you replace linearUpwindV or the momentum equation with vanLeer, Gamma or limitedLinear?
In addition, what happens when you increase the number of sweeps in solver?

Regards,

Franjo
__________________
Principal Developer of cfMesh and CF-MESH+
www.cfmesh.com
Social media: LinkedIn, Twitter, YouTube, Facebook, Pinterest, Instagram
franjo_j is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 15:36.