|
[Sponsors] |
[mesh manipulation] Cyclic BC how to obtain exact same meshes? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 12, 2014, 17:23 |
Cyclic BC how to obtain exact same meshes?
|
#1 |
Member
phil
Join Date: Oct 2013
Posts: 36
Rep Power: 12 |
Hi,
i want to simulate a steady state pipe flow the inlet and outlet are cyclic BC. I import my geometry through a .stl file. The problem is the mesh at the inlet and outlet are not from the same size and i cant make them cyclic with the "createPatch" utility. How can i achieve this? Hope i made my point Greetings & good night |
|
January 13, 2014, 00:11 |
|
#2 |
Member
Akshay Kumar
Join Date: Aug 2010
Location: India
Posts: 84
Rep Power: 15 |
Hi1
Use cyclicAMI for the two patches. This should do the job. |
|
January 13, 2014, 03:49 |
|
#3 |
Member
phil
Join Date: Oct 2013
Posts: 36
Rep Power: 12 |
Hey,
i've already tried to run it as cyclicAMI but the problem is the error: Code:
AMI: Creating addressing and weights between 142345 source faces and 142177 target faces --> FOAM Warning : From function AMIInterpolation<SourcePatch, TargetPatch>::checkPatches(const SourcePatch&, const TargetPatch&) in file lnInclude/AMIInterpolation.C at line 109 Source and target patch bounding boxes are not similar source box span : (0.000115329 0.0987398 0.152781) target box span : (8.22414e-05 0.0987483 0.152795) source box : (0.00671131 1.90699e-06 -0.00673927) (0.00682664 0.0987417 0.146042) target box : (0.0301843 1.448e-06 0.0167262) (0.0302665 0.0987498 0.169521) inflated target box : (0.0210879 -0.0090949 0.00762988) (0.0393629 0.107846 0.178617) Greetings |
|
January 13, 2014, 04:01 |
|
#4 |
Member
Akshay Kumar
Join Date: Aug 2010
Location: India
Posts: 84
Rep Power: 15 |
I see a warning not an error. Can you post the entire log?
|
|
January 13, 2014, 04:49 |
|
#5 |
Member
phil
Join Date: Oct 2013
Posts: 36
Rep Power: 12 |
Yep,
here it is. I create the cyclicAMI patches after my sHM parallel run with Code:
createPatch -overwrite Code:
t >createPatch -overwrite .... Build : 2.2.x-ae7a43cbbfe3 Exec : createPatch -overwrite Date : Jan 13 2014 Time : 10:46:56 Host : PID : 28596 Case : nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Reading createPatchDict Adding new patch inlet as patch 4 from { type cyclicAMI; neighbourPatch outlet; } Adding new patch outlet as patch 5 from { type cyclicAMI; neighbourPatch inlet; } Moving faces from patch eingang to patch 4 Moving faces from patch ausgang to patch 5 Doing topology modification to order faces. AMI: Creating addressing and weights between 142345 source faces and 142177 target faces --> FOAM Warning : From function AMIInterpolation<SourcePatch, TargetPatch>::checkPatches(const SourcePatch&, const TargetPatch&) in file lnInclude/AMIInterpolation.C at line 109 Source and target patch bounding boxes are not similar source box span : (0.000115329 0.0987398 0.152781) target box span : (8.22414e-05 0.0987483 0.152795) source box : (0.00671131 1.90699e-06 -0.00673927) (0.00682664 0.0987417 0.146042) target box : (0.0301843 1.448e-06 0.0167262) (0.0302665 0.0987498 0.169521) inflated target box : (0.0210879 -0.0090949 0.00762988) (0.0393629 0.107846 0.178617) --> FOAM FATAL ERROR: Unable to find initial target face From function void Foam::AMIInterpolation<SourcePatch, TargetPatch>::calcAddressing(const SourcePatch&, const TargetPatch&, label, label) in file lnInclude/AMIInterpolation.C at line 712. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::error::abort() in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 at cyclicAMIPolyPatch.C:0 #3 Foam::AMIInterpolation<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> >, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > >::update(Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&) in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libmeshTools.so" #4 Foam::AMIInterpolation<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> >, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > >::AMIInterpolation(Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::autoPtr<Foam::searchableSurface> const&, Foam::faceAreaIntersect::triangulationMode const&, bool) in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libmeshTools.so" #5 Foam::cyclicAMIPolyPatch::resetAMI() const in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libmeshTools.so" #6 Foam::polyBoundaryMesh::movePoints(Foam::Field<Foam::Vector<double> > const&) in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #7 Foam::polyMesh::movePoints(Foam::Field<Foam::Vector<double> > const&) in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #8 in "/sw/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/createPatch" #9 __libc_start_main in "/lib64/libc.so.6" #10 at /usr/src/packages/BUILD/glibc-2.11.2/csu/../sysdeps/x86_64/elf/start.S:116 Abgebrochen > Code:
FoamFile { version 2.0; format binary; class polyBoundaryMesh; location "constant/polyMesh"; object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 4 ( defaultFaces { type empty; inGroups 1(empty); nFaces 0; startFace 6220988; } eingang { type wall; nFaces 142345; startFace 6220988; } ausgang { type wall; nFaces 142177; startFace 6363333; } wand { type wall; nFaces 313335; startFace 6505510; } ) // ************************************************************************* // |
|
January 22, 2014, 05:20 |
|
#6 |
Member
Julian Langowski
Join Date: May 2011
Location: Bremen, Germany
Posts: 91
Rep Power: 14 |
How does your createPatchDict look?
Julian
__________________
πάντα ῥεῖ - Heraclitus |
|
January 22, 2014, 08:04 |
|
#7 | ||
Member
phil
Join Date: Oct 2013
Posts: 36
Rep Power: 12 |
Hey,
I'm 2 steps further. First of all i am able to use createPatch without errors. Here is the "createPatch" output and my createPatchDict file attached. file@bottom Quote:
But now i am hanging after calculating the velocity field (simpleFoam solver). Quote:
Have you any further suggestions? Greetings Phil ____________ UPDATE: createPatchDict: createPatchDict.txt Last edited by gelbebanane; January 22, 2014 at 10:45. Reason: file |
|||
January 22, 2014, 10:27 |
|
#8 |
Member
Chris L
Join Date: Sep 2012
Posts: 53
Rep Power: 13 |
You attachment "Attachment 28173" is a dead link can you repost?
Chris |
|
January 23, 2014, 08:18 |
|
#10 |
Member
Julian Langowski
Join Date: May 2011
Location: Bremen, Germany
Posts: 91
Rep Power: 14 |
Which solver did you use before? I am using GAMG and AMI in combination succesfully...
__________________
πάντα ῥεῖ - Heraclitus |
|
January 23, 2014, 11:02 |
|
#11 |
Member
phil
Join Date: Oct 2013
Posts: 36
Rep Power: 12 |
I have used the GAMG and now i switched to smoothSolver it is working but it needs over 1000 iterations to calculate. I use simpleFoam as CFD code what do you use with your AMI BC?
Greetings |
|
January 23, 2014, 12:24 |
|
#12 |
Member
Julian Langowski
Join Date: May 2011
Location: Bremen, Germany
Posts: 91
Rep Power: 14 |
I use MRFSimpleFoam.
__________________
πάντα ῥεῖ - Heraclitus |
|
January 29, 2014, 04:16 |
|
#13 |
Member
Akshay Kumar
Join Date: Aug 2010
Location: India
Posts: 84
Rep Power: 15 |
Hi Phil
Can you check if your 'mergeLevels' for the GAMG solver is >1 ? If yes then make it '1' and try running it . Akshay |
|
August 3, 2014, 23:28 |
cyclicAMI
|
#14 |
New Member
Mehid
Join Date: Nov 2013
Location: Tehran-Iran
Posts: 22
Rep Power: 0 |
Hi Phil,
I don't get how you passed the createPatch step..? you just said that you are 2 steps furthur. my problem is exactly there. after SHM when i want to change my BC's with "createPatch", i came to the error of [AMI: Creating addressing and weights between 98 source faces and 99 target faces --> FOAM Warning : From function AMIInterpolation<SourcePatch, TargetPatch>::checkPatches(const SourcePatch&, const TargetPatch&) in file lnInclude/AMIInterpolation.C at line 111 Source and target patch bounding boxes are not similar source box span : (105.417 80.3913 20.3297) target box span : (104.946 81.1113 20.031) source box : (-105 -79.9299 4.69853) (0.417155 0.461361 25.0282) target box : (-104.448 77.2405 5) (0.497791 158.352 25.031) inflated target box : (-111.155 70.5334 -1.70707) (7.20486 165.059 31.7381) --> FOAM FATAL ERROR: Unable to find initial target face] what's the procedure..? |
|
August 6, 2014, 05:14 |
|
#15 |
New Member
Mehid
Join Date: Nov 2013
Location: Tehran-Iran
Posts: 22
Rep Power: 0 |
yep..,
the problem is solved. I got it myslfe. Thanks anyway. |
|
September 12, 2014, 09:50 |
|
#16 | |
New Member
Mehid
Join Date: Nov 2013
Location: Tehran-Iran
Posts: 22
Rep Power: 0 |
Hey Foamers,
I'm being stopped. after SHM & createPatch for cyclicAMI faces when running simpleFOAM, it crashes out with the following error: Quote:
I'm still stuck in my progress. share your Tips if possible for help. Thanks in Advance, best regards, Metti. |
||
December 29, 2014, 15:24 |
|
#17 |
New Member
Ferdinand Leinbach
Join Date: Nov 2014
Posts: 9
Rep Power: 11 |
||
December 30, 2014, 08:32 |
|
#18 |
New Member
Mehid
Join Date: Nov 2013
Location: Tehran-Iran
Posts: 22
Rep Power: 0 |
dear Ferdinand,
I know nothing of your case details but ... what I did, was to get along with the Phill's procedure BUT the final point that he didn't mentioned was to change the option "writePrecision" in controlDict from 6 to a higher level like 10 or ... . by this, I finally became able to change my patches' type to cyclicAMI without any error, BUT even after that till now I wasn't able to RUN since it crashes out and I haven't found out why ,yet. Best, |
|
January 24, 2015, 09:50 |
|
#19 |
New Member
Ferdinand Leinbach
Join Date: Nov 2014
Posts: 9
Rep Power: 11 |
Dear Metti,
i got mine to work and not crash! What I did was (I generated my mesh in Pointwise): Set the boundaries I wanted to be cyclic to patches. Used the createPatch utility to transform these to cyclic boundary conditions. Then renumbered the mesh with renumberMesh. Voila, it stopped complaining for me! |
|
January 25, 2015, 15:52 |
|
#20 |
New Member
Mehid
Join Date: Nov 2013
Location: Tehran-Iran
Posts: 22
Rep Power: 0 |
Hey effi,
So, u suggest to mesh in another software? (From what i can undrestand from u're post )The way that u are mentioning, is to import the mesh from other tools (and i'm aware of that and of course i did it before and made my case RUN) But the problem was to deal with meshing for complex curved geometries including periodic patches in OF. specifically with SHM. Correct me if i misundrestood that. Best. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Setting up cyclic boundary condition | KateEisenhower | OpenFOAM Pre-Processing | 6 | January 11, 2017 17:17 |
LES in OF 3.0.1, no Output for nuSgs | Mirage | OpenFOAM Programming & Development | 1 | October 4, 2016 18:00 |
Cyclic Boundary Condition Errors? | nyflyer | OpenFOAM Running, Solving & CFD | 2 | April 26, 2016 14:14 |
Problem with cyclic patches | sven | OpenFOAM Running, Solving & CFD | 0 | December 5, 2011 14:27 |
Quadratic interpolation to obtain exact solution of Poiseuille flow | kar | OpenFOAM Running, Solving & CFD | 0 | March 25, 2008 03:31 |