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

cyclicAMI - kOmegaSST - divergence issue

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 2, 2018, 16:43
Default cyclicAMI - kOmegaSST - divergence issue
  #1
Senior Member
 
cyln
Join Date: Jul 2016
Posts: 102
Rep Power: 10
cyln is on a distinguished road
Hello, I am running a transient RANS simulation using pimpleFoam on OF4.1. My grid is a combination of non-comformal unstructured and structured grids. I defined an interface between the two grids using cyclicAMI (noOrdering).

I have completed my laminar simulation as it had been running without any issues. However, whenever I switch to kOmegaSST turbulence model, my simulation diverges whenever turbulence transport equations are solved for the first time. Before, I was using kOmegaSST for the same case with fully unstructured grid (without cyclicAMI) and it had been running for several thousand time steps. So, I know the divergence is related to cyclicAMI (noOrdering). I tried kEpsilon and kOmega, and I still got divergence.

As for discretization schemes, I am using linearUpwind for U, k and omega.

My boundary conditions are given below. Can anybody suggest me why my turbulent simulation (with cyclicAMI) is diverging?

polyMesh/boundary
Code:
7
(
    FARFIELD
    {   
        type            patch;
        nFaces          65021;
        startFace       67971616;
    }
    OUTLET
    {   
        type            patch;
        nFaces          6249;
        startFace       68036637;
    }
    BODY
    {   
        type            wall;
        inGroups        1(wall);
        nFaces          158614;
        startFace       68042886;
    }
    INTEROUTT2
    {   
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          40559;
        startFace       68201500;
        matchTolerance  0.001;
        transform       noOrdering;
        neighbourPatch  INTEROUT2;
    }
    INTEROUT2
    {   
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          53599;
        startFace       68242059;
        matchTolerance  0.001;
        transform       noOrdering;
        neighbourPatch  INTEROUTT2;
    }
    INTERFACET2
    {   
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          249569;
        startFace       68295658;
        matchTolerance  0.001;
        transform       noOrdering;
        neighbourPatch  INTERFACE2;
    }
    INTERFACE2
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          308799;
        startFace       68545227;
        matchTolerance  0.001;
        transform       noOrdering;
        neighbourPatch  INTERFACET2;
    }
0/U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (-40 0 0);

boundaryField
{
    FARFIELD
    {
        type fixedValue;
        value uniform (-40 0 0);
    }

    BODY
    {
         type            noSlip;
    }

    OUTLET
    {
        type            zeroGradient;
    }

        INTEROUTT2
    {
        type cyclicAMI;
    }

    INTEROUT2
    {
        type cyclicAMI;
    }

    INTERFACET2
    {
        type cyclicAMI;
    }

    INTERFACE2
    {
        type cyclicAMI;
    }
}
0/p
Code:
internalField   uniform 0.0;

boundaryField
{
    FARFIELD
    {
             type           zeroGradient;
    }

    BODY
    {
         type            zeroGradient;
    }

    OUTLET
    {
         type            fixedValue;
         value           uniform 0.0;
    }

         INTEROUTT2
    {
         type            cyclicAMI;
    }
    INTEROUT2
    {
         type            cyclicAMI;
    }


    INTERFACET2
    {
         type            cyclicAMI;
    }
    INTERFACE2
    {
         type            cyclicAMI;
    }
}
0/k
Code:
internalField   uniform 1e-6; 

boundaryField
{

    FARFIELD
    {
        type fixedValue;
        value uniform 1e-6;
    }

    BODY
    {
        type           fixedValue;
        value          uniform 0.0;
    }

    OUTLET
    {
        type            zeroGradient;
    }

                INTEROUTT2
    {
        type cyclicAMI;
    }

    INTEROUT2
    {
        type cyclicAMI;
    }

    INTERFACET2
    {
        type cyclicAMI;
    }

    INTERFACE2
    {
        type cyclicAMI;
    }
0/omega
Code:
internalField   uniform 0.55 ; 


boundaryField
{

    FARFIELD
    {
        type fixedValue;
        value uniform 0.55;     //  0.0017; // 22.36;
    }

    BODY
    {
        type           fixedValue;
        value          uniform 241632653;
    }

    OUTLET
    {
        type            zeroGradient;

        }

        INTEROUTT2
    {
        type cyclicAMI;
    }

    INTEROUT2
    {
        type cyclicAMI;
    }

    INTERFACET2
    {
        type cyclicAMI;
    }

    INTERFACE2
    {
        type cyclicAMI;
    }
0/nut
Code:
internalField   uniform 0;

boundaryField
{

    FARFIELD
    {
      type            calculated;
      value           uniform 0;
    }

    OUTLET
    {
         type            calculated;
        value           uniform 0;
    }

    BODY
    {
        type            fixedValue; 
        value           uniform 0;
    }

                INTEROUTT2
    {
        type cyclicAMI;
   }

    INTEROUT2
    {
    type cyclicAMI;


}

    INTERFACET2
    {
       type cyclicAMI;


 }

    INTERFACE2
    {
        type cyclicAMI;

}
cyln is offline   Reply With Quote

Old   October 5, 2018, 11:06
Default
  #2
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 647
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi,

I am guessing (as you did not show the exact output) that you have a divide by zero at the cyclicAMI interfaces. The reason is that you have no value entry and this may mean that a default of 0 is used. Just give a value equal to your internalField and rerun?

Regards,
Tom
tomf 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


Similar Threads
Thread Thread Starter Forum Replies Last Post
time step continuity problem in VAWT simulation lpz_michele OpenFOAM Running, Solving & CFD 5 February 22, 2018 20:50
Divergence in non-Newtonian fluid UDF moabdi FLUENT 0 June 23, 2016 12:30
Divergence issue with SST Turb Model jehanz SU2 3 May 19, 2015 19:01
simpleFoam kOmegaSST LowRe pressure divergence Pat84 OpenFOAM Running, Solving & CFD 2 August 12, 2013 18:42
Floating point error and divergence detected aannjj FLUENT 0 July 2, 2013 04:44


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