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

GGI bug??

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 25, 2010, 00:39
Default GGI bug??
  #1
Senior Member
 
Pavan
Join Date: May 2009
Location: Melbourne
Posts: 101
Rep Power: 8
rieuk is on a distinguished road
Hey guys, a couple of things firstly,

does moveDynamicMesh work with GGI? Not working for me. Output is:

Code:
[pab565@vayu1 pitchEllipse4]$ moveDynamicMesh
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5-dev                               |
|   \\  /    A nd           | Revision: 1669                                  |
|    \\/     M anipulation  | Web:      http://www.OpenFOAM.org               |
\*---------------------------------------------------------------------------*/
Exec   : moveDynamicMesh
Date   : Apr 25 2010
Time   : 14:31:01
Host   : vayu1
PID    : 6875
Case   : /short/g75/Pavan_FYP/pitchEllipse4
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 0 1)
    rpm   : 1
--> 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 = 0.002
Clear VTK directory
Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 2.22045e-16 average: 7.51536e-17
  Largest master weighting factor correction: 2.66454e-15 average: 1.34081e-16

volume continuity errors : volume = 2.62147e-05, max error = 6.07437e-12, sum local = 1.78474e-15, global = -1.30073e-17
ExecutionTime = 0.2 s  ClockTime = 0 s




    request for surfaceScalarField phi from objectRegistry region0 failed
    available objects of type surfaceScalarField are

2
(
weightingFactors
meshPhi
)


    From function objectRegistry::lookupObject<Type>(const word&) const
    in file /home/565/pab565/OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 142.

FOAM aborting

Abort
Secondly, I'm trying to use sinusoidal rotation and to this end I add this line in the bool Foam::mixerGgiFvMesh::update() function of the mixerGgiFvMesh.C file in the mixerGgiFvMesh folder:

Code:
rpm_=10/9*pi*cos(pi/6*time().value());
Upon compiling, this works fine and rotates my mesh sinusoidally as expected - HOWEVER when I change this line to:

Code:
rpm_=5/9*pi*cos(pi/6*time().value());
and compile, the code solves and everything, says it's doing the GGI interpolation and all - but the mesh does not move at all. What the???
rieuk is offline   Reply With Quote

Old   April 25, 2010, 05:24
Default
  #2
Senior Member
 
Henrik Rusche
Join Date: Mar 2009
Location: Braunschweig, Niedersachsen, Germany
Posts: 275
Rep Power: 9
henrik is on a distinguished road
Dear rieuk,

a) Try switching off the ggiCheck function object.

b) You may be a victim of implicit conversion. Try:

rpm_=5.0/9.0*pi*cos(pi/6.0*time().value());
Best Regards,

Henrik
henrik is offline   Reply With Quote

Old   April 25, 2010, 10:47
Default
  #3
Senior Member
 
Pavan
Join Date: May 2009
Location: Melbourne
Posts: 101
Rep Power: 8
rieuk is on a distinguished road
Solved and solved. Cheers, Henrik! And sorry about crying wolf in the thread topic
rieuk 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
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
Performance of GGI case in parallel hannes OpenFOAM Running, Solving & CFD 26 August 3, 2011 03:07
BCs for parallel ggi NickG OpenFOAM 1 March 15, 2010 12:12
CFX GGI Interface Error (non-overlapping) surge519 CFX 1 August 3, 2009 18:54


All times are GMT -4. The time now is 02:03.