|
[Sponsors] |
July 2, 2007, 11:43 |
Hi foamers,
I would like to
|
#1 |
Member
|
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 Ergon Research s.r.l. Via Panciatichi, 92 50127 Florence - ITALY Tel: +39 055 0763716 Mob: +39 320 9460153 e-mail: cosimo.bianchini@ergonresearch.it URL: www.ergonresearch.it |
|
July 22, 2007, 23:02 |
Thank you Cosimo, your codes i
|
#2 |
New Member
Haibin Li
Join Date: Mar 2009
Posts: 8
Rep Power: 17 |
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? |
|
July 23, 2007, 05:11 |
Hi Haibin,
I hope you can g
|
#3 |
Member
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 17 |
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 |
|
July 24, 2007, 21:02 |
Thank you Luca, I think it's a
|
#4 |
New Member
Haibin Li
Join Date: Mar 2009
Posts: 8
Rep Power: 17 |
Thank you Luca, I think it's a very important feature, tons of cases need it. I will put my hands soon.
|
|
July 31, 2007, 08:26 |
Ciao Cosimo,
I tried to ins
|
#5 |
Member
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 62
Rep Power: 17 |
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 |
|
July 31, 2007, 09:03 |
Hi Thomas,
welcome in this
|
#6 |
Member
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 17 |
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 |
|
August 6, 2007, 12:28 |
Hi Luca,
thanks for the war
|
#7 |
Member
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 62
Rep Power: 17 |
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 |
|
August 6, 2007, 12:58 |
Hi Thomas,
just put the whol
|
#8 |
Member
|
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 Ergon Research s.r.l. Via Panciatichi, 92 50127 Florence - ITALY Tel: +39 055 0763716 Mob: +39 320 9460153 e-mail: cosimo.bianchini@ergonresearch.it URL: www.ergonresearch.it |
|
August 8, 2007, 09:07 |
Hi Cosimo,
I keep on gettin
|
#9 |
Member
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 62
Rep Power: 17 |
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 |
|
August 8, 2007, 10:00 |
Hi Thomas,
could it be you ad
|
#10 |
Member
|
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 Ergon Research s.r.l. Via Panciatichi, 92 50127 Florence - ITALY Tel: +39 055 0763716 Mob: +39 320 9460153 e-mail: cosimo.bianchini@ergonresearch.it URL: www.ergonresearch.it |
|
August 9, 2007, 05:59 |
Hi Cosimo,
that did the tri
|
#11 |
Member
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 62
Rep Power: 17 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 16:36 |
periodic b.c. | nick | FLUENT | 0 | September 1, 2006 10:46 |
Periodic BC | mourad | FLUENT | 0 | June 30, 2005 12:35 |
periodic | Ehsan Alborzi | FLUENT | 1 | February 17, 2004 04:29 |