CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

MixerGgi with additional movement of the inner ring in a Taylor-Couette-System

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 13, 2012, 06:13
Default MixerGgi with additional movement of the inner ring in a Taylor-Couette-System
  #1
New Member
 
Matthias Neben
Join Date: Oct 2011
Location: Cottbus (Germany)
Posts: 28
Rep Power: 14
mneben is on a distinguished road
Hello Foamers,


I created a new motion solver based on the motion solver “mixerGgiFvMesh”. The result is a rotating inner ring which also moves horizontally or vertically.


Therefor I rewrote the last lines of the source code of mixerGgiFvMesh:


bool Foam::journalBearingMotionSolver::update() //instead of Foam::mixerGgiFvMesh::update()
{
//-points of the moving par:
pointField p = motionPtr_->newPoints();





//-moves the origin of the rotating part in the same way, like the inner ring moves:
if (time().timeIndex()==1)
{
# include "initMotion.H"
}
else
{
# include "calcMotion.H"
point& orgn=csPtr_().origin();
orgn=sIst_;
}
Info<<"origin"<<cs().origin()<<endl;


//-Rotational speed needs to be converted from rpm
p=csPtr_->globalPosition
(
csPtr_->localPosition(p)
+ vector(0, rpm_*360.0*time().deltaT().value()/60.0, 0)
*movingPointsMask()
);






fvMesh::movePoints(p);



//-Move points, returns volumes swept by faces in motion

return false;
}


Then I added a pointMotionU-file with the mesh-motion-boundary-conditions and everything seems to be fine for small cases on a single cpu-core. For bigger cases I wanted to use decomposePar, but every time I got an error message like this:


Processor 0
Number of cells = 452264
Number of faces shared with processor 1 = 387
Number of processor patches = 1
Number of processor faces = 387
Number of boundary faces = 941279
--> FOAM Warning :
From function dlLibraryTable:pen(const fileName& functionLibName)
in file db/dlLibraryTable/dlLibraryTable.C at line 86
could not load /home/matthias/OpenFOAM/matthias-1.6-ext/lib/libmyDynamicFvMesh.so: undefined symbol: _ZTIN4Foam13dynamicFvMeshE

Processor 1
Number of cells = 452264
Number of faces shared with processor 0 = 387
Number of processor patches = 1
Number of processor faces = 387
Number of boundary faces = 941317

Number of processor faces = 387
Max number of processor patches = 1
Max number of faces between processors = 387

Processor 0: field transfer
--> FOAM Warning :
From function dlLibraryTable:pen(const fileName& functionLibName)
in file db/dlLibraryTable/dlLibraryTable.C at line 86
could not load /home/matthias/OpenFOAM/matthias-1.6-ext/lib/libmyDynamicFvMesh.so: undefined symbol: _ZTIN4Foam13dynamicFvMeshE
Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider


--> FOAM FATAL ERROR:
size of field = 978086 is not the same as the size of mesh = 941670

From function DimensionedField<Type, GeoMesh>:imensionedField(const IOobject& io,const Mesh& mesh, const dimensionSet& dims, const Field<Type>& field)
in file /build/buildd/openfoam-1.6-ext-1.6.0/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/DimensionedField.C at line 72.

FOAM aborting



When I erase the pointMotionU-file then decomposePar runs till the end. With this fact and the text of the line “Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider” I come to the conclusion that the ggi-routines don't work with pointMotionU.
What can I do to solve this problem?


Awaiting for replies!

Regards
Matthias
mneben is offline   Reply With Quote

Reply

Tags
ggi, mesh motion, taylor-couette


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFX11 + Fortran compiler ? Mohan CFX 20 March 30, 2011 18:56
emag beta feature: charge density charlotte CFX 4 March 22, 2011 09:14
The mechanism system of wave movement to power catamaran Yogi_NA'03 ANSYS 2 December 8, 2009 07:21
plz rply urgent regrding vof model for my system garima chaudhary FLUENT 1 July 20, 2007 08:37
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 20:09


All times are GMT -4. The time now is 10:25.