|
[Sponsors] |
January 30, 2011, 05:31 |
TurbDyMFoam-Ggi
|
#1 |
Member
Jason Eason
Join Date: Jan 2010
Location: Portage, Michigan
Posts: 45
Rep Power: 16 |
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
__________________
Debian Squeeze - OpenFOAM-2.1.x, Paraview-3.12.0 |
|
January 31, 2011, 00:45 |
|
#2 |
Senior Member
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 434
Rep Power: 22 |
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? |
|
January 31, 2011, 02:36 |
|
#3 |
Member
Jason Eason
Join Date: Jan 2010
Location: Portage, Michigan
Posts: 45
Rep Power: 16 |
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; } )
__________________
Debian Squeeze - OpenFOAM-2.1.x, Paraview-3.12.0 |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Floating point exception error | Alan | OpenFOAM Running, Solving & CFD | 11 | July 1, 2021 22:51 |
Ggi | FabOr | OpenFOAM | 17 | May 9, 2013 11:19 |
Difference between ggi and overlapGgi? GGI Tips and Tricks? | philippose | OpenFOAM Running, Solving & CFD | 7 | January 16, 2013 10:40 |
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 |