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