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/)
-   -   TurbDyMFoam-Ggi (http://www.cfd-online.com/Forums/openfoam-solving/84447-turbdymfoam-ggi.html)

JulytoNovember January 30, 2011 05:31

TurbDyMFoam-Ggi
 
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?

cnsidero January 31, 2011 00:45

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?

JulytoNovember January 31, 2011 02:36

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.


All times are GMT -4. The time now is 23:21.