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

Periodic GGI

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 2, 2007, 10:43
Default Hi foamers, I would like to
  #1
Member
 
cosimo bianchini
Join Date: Mar 2009
Location: Florence, Tuscany, Italy
Posts: 88
Rep Power: 8
cosimobianchini is on a distinguished road
Send a message via Skype™ to cosimobianchini
Hi foamers,
I would like to know if someone can suggest any hint on how to proceed in implementing periodic GGI: that means cyclic boundaries that are not conformal.

I already implemented a beta version of such boundary conditions (unfortunately calculations explodes after some iteration and I don't understand why). NOTE: the code is working for conformal cyclic.

The idea I followed (inspired by ggiInterface) is to use patchToPatchInterpolation::faceInterpolate() to calculate a "mean" value associated with the reconstructed twin cell centroid.
I tried to reconstruct the reconstructed twin cell centroid in different ways (averaging centroid position on the shadow patch, reconstructing its position on the normal vector exiting the face centroid etc etc) but results are not changing that much.
I also tried to average the values on the shadow patch including more cells or excluding some neighbor cell in case the weight was too low and again no significant change.

I read somewhere here in the forum (thread named: "Computing a cell-centered scalarField from a vertex-centered scalarField") that inverse distance weighting could cause conservativeness problems, can someone be more specific? How can this problem be avoided / solved?

If someone want to have a look at the code I will appreciate any advice:
http://htc.de.unifi.it/docpub/ggiCyclicTools.tgz .
Here are the cases I have been running as tests:
http://htc.de.unifi.it/docpub/conf.tgz ,
http://htc.de.unifi.it/docpub/non-conf.tgz

Many thanks in advance
Cosimo
__________________
Cosimo Bianchini

Energy Engineering Department "S. Stecco"
University of Florence
Via di S.Marta, 3
50139 Florence - ITALY
Tel: +39 055 4796575
Fax: +39 055 4796342
Mob: +39 320 9460153
e-mail: cosimo.bianchini@htc.de.unifi.it
URL: www.htc.de.unifi.it
cosimobianchini is offline   Reply With Quote

Old   July 22, 2007, 22:02
Default Thank you Cosimo, your codes i
  #2
New Member
 
Haibin Li
Join Date: Mar 2009
Posts: 8
Rep Power: 8
haibinli is on a distinguished road
Thank you Cosimo, your codes is just what I want.
But I'm not very clear the line 265 in CalcPatchToPatchWeigths.C. How should I change it in a rotational Cyclic BC?
haibinli is offline   Reply With Quote

Old   July 23, 2007, 04:11
Default Hi Haibin, I hope you can g
  #3
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
Hi Haibin,

I hope you can give us a hand to develop this fondamental feature for OpenFOAM.
Only put the translation vector in this line that leads the cell center projection. In the case of rotational cyclic put the chord vector that defines your periodicity.

Regards

Luca
luca is offline   Reply With Quote

Old   July 24, 2007, 20:02
Default Thank you Luca, I think it's a
  #4
New Member
 
Haibin Li
Join Date: Mar 2009
Posts: 8
Rep Power: 8
haibinli is on a distinguished road
Thank you Luca, I think it's a very important feature, tons of cases need it. I will put my hands soon.
haibinli is offline   Reply With Quote

Old   July 31, 2007, 07:26
Default Ciao Cosimo, I tried to ins
  #5
Member
 
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 60
Rep Power: 8
anger is on a distinguished road
Ciao Cosimo,

I tried to install the files you linked above but got an error stating:
meshes/polyMesh/polyPatches/derivedPolyPatches/arbitraryCyclicPolyPatch/arbitrar yCyclicPolyPatch.C: In member function 'virtual void Foam::arbitraryCyclicPolyPatch::initUpdateTopology ()':
meshes/polyMesh/polyPatches/derivedPolyPatches/arbitraryCyclicPolyPatch/arbitrar yCyclicPolyPatch.C:364: error: 'initUpdateTopology' is not a member of 'Foam::polyPatch'
I assume that I put a file in the wrong directory so could you please give some more details how to install the files as given in the README file?

Best regards,
-Thomas
anger is offline   Reply With Quote

Old   July 31, 2007, 08:03
Default Hi Thomas, welcome in this
  #6
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
Hi Thomas,

welcome in this thread and thanks for contributing to PeriodicGGI developement.

You have to test the code in OpenFOAM version 1.3. It seems that you are using the latest version 1.4.

Please feel free to ask any question and to suggest every possible bug fix or change.

Keep in touch

Luca
luca is offline   Reply With Quote

Old   August 6, 2007, 11:28
Default Hi Luca, thanks for the war
  #7
Member
 
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 60
Rep Power: 8
anger is on a distinguished road
Hi Luca,

thanks for the warm welcome. However, I'm not sure whether I'm able to really contribute to the development of PeriodicGGI due to my lack of knowledge on this subject )-:
Anyway, I installed OF V1.3 and recopiled the vanilla source tree of $FOAM_SRC/OpenFOAM and $FOAM_SRC/finiteVolume which worked (almost) without problems. After applying the modifications found in the README file, I still get errors. So my question is: the README states:
"put arbitraryCyclicPolyPatch in $FOAM_SRC/OpenFOAM/meshes/polyMesh/polyPatches/derivedPolyPatches".
Does that mean that I should put the whole directory into this place, only the files found in the directory or anything else? I assume that I'm missing something, so could you please be more precise on the installation procedure?

Best regards,
-Thomas
anger is offline   Reply With Quote

Old   August 6, 2007, 11:58
Default Hi Thomas, just put the whol
  #8
Member
 
cosimo bianchini
Join Date: Mar 2009
Location: Florence, Tuscany, Italy
Posts: 88
Rep Power: 8
cosimobianchini is on a distinguished road
Send a message via Skype™ to cosimobianchini
Hi Thomas,
just put the whole directories in the places suggested by the README file.
Sorry for not having been clear enough.
Best regards,

Cosimo
__________________
Cosimo Bianchini

Energy Engineering Department "S. Stecco"
University of Florence
Via di S.Marta, 3
50139 Florence - ITALY
Tel: +39 055 4796575
Fax: +39 055 4796342
Mob: +39 320 9460153
e-mail: cosimo.bianchini@htc.de.unifi.it
URL: www.htc.de.unifi.it
cosimobianchini is offline   Reply With Quote

Old   August 8, 2007, 08:07
Default Hi Cosimo, I keep on gettin
  #9
Member
 
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 60
Rep Power: 8
anger is on a distinguished road
Hi Cosimo,

I keep on getting an errror:
SOURCE_DIR=interpolations/patchToPatchInterpolation
SOURCE=interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C ; g++ -m64 -DlinuxAMD64 -DDP -Wall -W -Wno-unused-parameter -Wold-style-cast -march=opteron -O3 -fno-gcse -DNoRepository -ftemplate-depth-30 -DWM_PROJECT_VERSION='"'1.3'"' -I/home/anger/OpenFOAM/OpenFOAM-1.3/src/zlib-1.2.1 -I/home/anger/OpenFOAM/OpenFOAM-1.3/src/OpenFOAM/lnInclude -IlnInclude -I. -fPIC -c $SOURCE -o Make/linuxAMD64DPOpt/CalcPatchToPatchWeights.o
interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C:41: error: redefinition of 'void Foam::PatchToPatchInterpolation<frompatch,>::calcP ointAddressing() const'
interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C:41: error: 'void Foam::PatchToPatchInterpolation<frompatch,>::calcP ointAddressing() const' previously declared here
interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C:235: error: redefinition of 'void Foam::PatchToPatchInterpolation<frompatch,>::calcF aceAddressing() const'
interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C:235: error: 'void Foam::PatchToPatchInterpolation<frompatch,>::calcF aceAddressing() const' previously declared here
make: *** [Make/linuxAMD64DPOpt/CalcPatchToPatchWeights.o] Fehler 1
Anyone can help me on this?

Best regards,
-Thomas
anger is offline   Reply With Quote

Old   August 8, 2007, 09:00
Default Hi Thomas, could it be you ad
  #10
Member
 
cosimo bianchini
Join Date: Mar 2009
Location: Florence, Tuscany, Italy
Posts: 88
Rep Power: 8
cosimobianchini is on a distinguished road
Send a message via Skype™ to cosimobianchini
Hi Thomas,
could it be you add the wrong files into the OpenFOAM/Make/files file?

It should look like:

interpolations = interpolations
interpolation = $(interpolations)/interpolation
$(interpolations)/patchToPatchInterpolation/PatchToPatchInterpolationName.C

Sorry but I got no other clue to think about.
Regards,
Cosimo
__________________
Cosimo Bianchini

Energy Engineering Department "S. Stecco"
University of Florence
Via di S.Marta, 3
50139 Florence - ITALY
Tel: +39 055 4796575
Fax: +39 055 4796342
Mob: +39 320 9460153
e-mail: cosimo.bianchini@htc.de.unifi.it
URL: www.htc.de.unifi.it
cosimobianchini is offline   Reply With Quote

Old   August 9, 2007, 04:59
Default Hi Cosimo, that did the tri
  #11
Member
 
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 60
Rep Power: 8
anger is on a distinguished road
Hi Cosimo,

that did the trick. I was putting the names of all files in the subdirectory of the unpacked archive ggiCyclicTools.tgz into the Make/files file, but obviuosly only these for which a .dep- file exists should be placed there.
I'm keen on getting your eample(s) runnning.

Best regards,
-Thomas
anger 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
Periodic BC Valerio FLUENT 7 January 7, 2008 08:11
Must periodic oscillation give periodic results? zonexo Main CFD Forum 6 May 13, 2007 15:36
periodic b.c. nick FLUENT 0 September 1, 2006 09:46
Periodic BC mourad FLUENT 0 June 30, 2005 11:35
periodic Ehsan Alborzi FLUENT 1 February 17, 2004 04:29


All times are GMT -4. The time now is 05:40.