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

TurbDyMFoam-Ggi

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 30, 2011, 05:31
Default TurbDyMFoam-Ggi
  #1
Member
 
Jason Eason
Join Date: Jan 2010
Location: Portage, Michigan
Posts: 44
Rep Power: 7
JulytoNovember is on a distinguished road
My 3D helicopter rotor has a max skewness of .24(ggi included). I've changed the cellZones name to movingCells. I have all my sets and zones, everything is in order. When I run the simulation this is the error I get
Code:
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh mixerGgiFvMesh
void mixerGgiFvMesh::addZonesAndModifiers() : Zones and modifiers already present.  Skipping.
Mixer mesh:
    origin: (0 0 0)
    axis  : (0 1 0)
    rpm   : 600
Reading field p

Reading field U

Reading/calculating face flux field phi

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.00027763051 average: 0.00027760734
  Largest master weighting factor correction: 0 average: 0

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    alphaEps        0.76923;
}

Reading field rAU if present


Starting time loop

Courant Number mean: 0 max: 0 velocity magnitude: 0
deltaT = 1.2e-07
--> FOAM Warning : 
    From function dlLibraryTable::open(const dictionary& dict, const word& libsEntry, const TablePtr tablePtr)
    in file lnInclude/dlLibraryTableTemplates.C at line 68
    library "libsampling.so" did not introduce any new entries

Creating ggi check
Time = 1.2e-07

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.00027758037 average: 0.00027756526
  Largest master weighting factor correction: 7.6605389e-14 average: 5.3884092e-15

BiCGStab:  Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
BiCGStab:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
BiCGStab:  Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
BiCGStab:  Solving for p, Initial residual = 1, Final residual = 4.5742132e+29, No Iterations 1000
BiCGStab:  Solving for p, Initial residual = 1.3265969e-10, Final residual = 1.3265969e-10, No Iterations 0
time step continuity errors : sum local = 2.8490162e+16, global = -7.4676087e+10, cumulative = -7.4676087e+10
BiCGStab:  Solving for p, Initial residual = 1.3302648e-10, Final residual = 1.3302648e-10, No Iterations 0
BiCGStab:  Solving for p, Initial residual = 1.3302648e-10, Final residual = 1.3302648e-10, No Iterations 0
time step continuity errors : sum local = 2.8568241e+16, global = -7.4676087e+10, cumulative = -1.4935217e+11
BiCGStab:  Solving for epsilon, Initial residual = 1, Final residual = 0.10807524, No Iterations 1000
bounding epsilon, min: -7.636971e+21 max: 5.0378494e+23 average: 6.2606837e+21
BiCGStab:  Solving for k, Initial residual = 1, Final residual = 21.857977, No Iterations 1000
ExecutionTime = 1431.23 s  ClockTime = 1433 s

Courant Number mean: 7.5378835e+16 max: 2.3315945e+18 velocity magnitude: 1.4461871e+22
deltaT = 1.0293385e-26
GGI pair (insideSlider, outsideSlider) : 8.0916216e+16 8.0938583e+16 Diff = -1.4025123e+13 or 0.017332896 %
Time = 1.2e-07

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.00027758037 average: 0.00027756526
  Largest master weighting factor correction: 2.7844393e-13 average: 5.0794979e-15

BiCGStab:  Solving for Ux, Initial residual = 0.94560656, Final residual = 0.00018205648, No Iterations 1000
BiCGStab:  Solving for Uy, Initial residual = 0.99931858, Final residual = 0.00010471905, No Iterations 1000
BiCGStab:  Solving for Uz, Initial residual = 0.94082729, Final residual = 9.5636703e-05, No Iterations 1000
BiCGStab:  Solving for p, Initial residual = 1, Final residual = 49.154867, No Iterations 1000
BiCGStab:  Solving for p, Initial residual = 0.00024438771, Final residual = 0.00069014118, No Iterations 1000
time step continuity errors : sum local = 0.61264003, global = -0.0062865807, cumulative = -1.4935217e+11
BiCGStab:  Solving for p, Initial residual = 0.00087427042, Final residual = 7.2468904e-05, No Iterations 1000
BiCGStab:  Solving for p, Initial residual = 0.0010638005, Final residual = 0.00010481927, No Iterations 1000
time step continuity errors : sum local = 0.10240109, global = -0.0012875196, cumulative = -1.4935217e+11
BiCGStab:  Solving for epsilon, Initial residual = 1, Final residual = 4.0271798e+60, No Iterations 1000
bounding epsilon, min: -1.2443686e+113 max: 7.7686314e+113 average: 5.829649e+109
BiCGStab:  Solving for k, Initial residual = 0.0018094685, Final residual = 1.9350511e-32, No Iterations 1
bounding k, min: -4194304 max: 3.3802498e+54 average: 2.7688591e+50
ExecutionTime = 2991.93 s  ClockTime = 2997 s

Courant Number mean: 38.087887 max: 55509.304 velocity magnitude: 1.3151239e+28
deltaT = 3.7087063e-32
GGI pair (insideSlider, outsideSlider) : 7.810657e+20 8.1920659e+20 Diff = -2.8190449e+18 or 0.36092289 %
Time = 1.2e-07

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.00027758037 average: 0.00027756526
  Largest master weighting factor correction: 2.7844393e-13 average: 5.2569032e-15

BiCGStab:  Solving for Ux, Initial residual = 0.99999844, Final residual = 1.7850849e-05, No Iterations 1000
Floating point exception
Any thoughts of what the problem is?
__________________
Debian Squeeze - OpenFOAM-2.1.x, Paraview-3.12.0
JulytoNovember is offline   Reply With Quote

Old   January 31, 2011, 00:45
Default
  #2
Senior Member
 
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 369
Rep Power: 12
cnsidero is on a distinguished road
You might want to upload the main files of the case so we can take a peek but upon quick inspection there appears to be something wrong with your GGI interface - there's a large difference between the magnitudes of the weighting factors on the master and shadow patches.

What's the mesh resolution of the interface meshes on innerSlider and outerSlider patches? Are they similar?
cnsidero is offline   Reply With Quote

Old   January 31, 2011, 02:36
Default
  #3
Member
 
Jason Eason
Join Date: Jan 2010
Location: Portage, Michigan
Posts: 44
Rep Power: 7
JulytoNovember is on a distinguished road
Code:
boundary file
4
(
    rotor
    {
        type            patch;
        nFaces          1040;
        startFace       201040;
    }
    insideSlider
    {
        type            ggi;
        nFaces          40000;
        startFace       202080;
        shadowPatch     outsideSlider;
        zone            insideZone;
        bridgeOverlap   false;
    }
    outsideSlider
    {
        type            ggi;
        nFaces          40000;
        startFace       242080;
        shadowPatch     insideSlider;
        zone            outsideZone;
        bridgeOverlap   false;
    }
    defaultFaces
    {
        type            patch;
        nFaces          1600;
        startFace       282080;
    }
)
My ggi is comprised of 4 blocks, using the arc function, I formed into a ring. The outer ring is the outsideSlider. The inner is the insideSlider. The top and bottom are not incorporated into the ggi, therefore; the program calls from the default faces. Could the design of my ggi be the problem? Thank you for your help.
__________________
Debian Squeeze - OpenFOAM-2.1.x, Paraview-3.12.0
JulytoNovember 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
Ggi FabOr OpenFOAM 17 May 9, 2013 10:19
Difference between ggi and overlapGgi? GGI Tips and Tricks? philippose OpenFOAM Running, Solving & CFD 7 January 16, 2013 10:40
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
GGI in OpenFOAM-1.5-dev philippose OpenFOAM Running, Solving & CFD 14 November 13, 2011 15:55
Rotor-stator computation with Ggi interface (turbDyMFoam) Stylianos OpenFOAM Running, Solving & CFD 28 January 12, 2011 21:37


All times are GMT -4. The time now is 22:26.