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

[blockMesh] CyclicAMI boundaries Issues

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By zarox

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 24, 2020, 11:55
Default CyclicAMI boundaries Issues
  #1
New Member
 
Marco Sala
Join Date: Oct 2019
Posts: 3
Rep Power: 6
masal is on a distinguished road
Hello everyone, i'm a beginner in CFD and in OpenFoam in general. I want build a cylinder using CyclicAMI boundaries fot the inlet e for the outlet. I set up the blockMeshDict in this way:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

//- Half angle of wedge in degrees
halfAngle 2.5;

//- Radius of pipe [m]
radius 0.5;


radHalfAngle #calc "degToRad($halfAngle)";
x #calc "$radius*sin($radHalfAngle)";
x_ #calc "-1.0*$x";
y #calc "$radius*cos($radHalfAngle)";

//-Pipe lengt
L 15;

vertices
(
(0.0 0.0 0) //0
(0.0 0.0 $L)
(0.0 0.0 $L) //2
(0.0 0.0 0)

($x_ $y 0) //4
($x_ $y $L)
($x $y $L) //6
($x $y 0)

);

blocks
(
// cyc_inlet block
hex (0 3 7 4 1 2 6 5) (1 80 10) simpleGrading (1 0.4 1)
);

edges
(
arc 4 7 (0 $radius 0)
arc 5 6 (0 $radius $L)
);

boundary
(
cyc_inlet
{
type cyclicAMI;
neighbourPatch cyc_outlet;
faces
(
(0 4 7 3)
);
transform translational;
separationVector (0 0 15);
// matchTolerance 300;
}

cyc_outlet
{
type cyclicAMI;
neighbourPatch cyc_inlet;
faces
(
(1 2 6 5)
);
transform translational;
separationVector (0 0 -15);
// matchTolerance 300;
}

side1
{
type wedge;
faces
(
(0 1 5 4)
);
}

side2
{
type wedge;
faces
(
(7 6 2 3)
);

}

walls
{
type wall;
faces
(
(4 5 6 7)
// (3 2 1 0)
);
}
centerLine
{
type empty;
faces
(
(3 2 1 0)
);
}
);


// ************************************************** *********************** //


Running blockMesh I have the error:
--> FOAM FATAL ERROR:
The distance between the centre of patch cyc_inlet and the transformed centre of patch cyc_outlet is 30.
This is greater than the match tolerance of 1.04377e-05 for the patch.
Check that the patches are geometrically similar and that any transformations defined between them are correct
It might be possible to fix this problem by increasing the "matchTolerance" setting for this patch in the boundary



I set matchTolerance as 300 to resolve the error, and in this way I can compute the mesh, but when I'm tray to use decomposePar I obtain:


[AMI: Creating addressing and weights between 80 source faces and 80 target faces
--> FOAM Warning :
From function void Foam::AMIMethod::checkPatches() const
in file AMIInterpolation/AMIInterpolation/AMIMethod/AMIMethod/AMIMethod.C at line 65
Source and target patch bounding boxes are not similar
source box span : (0.0436194 0.499524 0)
target box span : (0.0436194 0.499524 0)
source box : (-0.0218097 0 0) (0.0218097 0.499524 0)
target box : (-0.0218097 0 30) (0.0218097 0.499524 30)
inflated target box : (-0.0468809 -0.0250712 29.9749) (0.0468809 0.524595 30.0251)


--> FOAM FATAL ERROR:
Unable to find initial target face

From function bool Foam::AMIMethod::initialise(Foam::labelListList&, Foam::scalarListList&, Foam::labelListList&, Foam::scalarListList&, Foam::label&, Foam::label&)
in file AMIInterpolation/AMIInterpolation/AMIMethod/AMIMethod/AMIMethod.C at line 134.

FOAM aborting



I'm sure that with OpenFoam4.x version the first solution works, I'm using OpenFoam8, it's a problem of syntax?
masal is offline   Reply With Quote

Old   February 11, 2022, 02:37
Default
  #2
Member
 
Emeline Noel
Join Date: Dec 2013
Location: Paris
Posts: 31
Rep Power: 12
zarox is on a distinguished road
Hello,


I have the same issue, do you find what we did wrong?
zarox is offline   Reply With Quote

Old   February 11, 2022, 03:02
Default SeparationVector in OF1712 and OF9
  #3
Member
 
Emeline Noel
Join Date: Dec 2013
Location: Paris
Posts: 31
Rep Power: 12
zarox is on a distinguished road
Hi,


With OF1712 the SeparationVector should be given from the patch to the neighbourPatch (for example minY patch to maxY neighbourPatch (0 51. 0.)).
But for the OF9, the distance should be given form the neighbourgPatch to the patch, so here (0 -51 0).



The exact same blockMeshDict and 0/ folder from OF1712 give the error in OF9:


The distance between the centre of patch minY and the transformed centre of patch maxY is 102.


The error vanish is the distance is given from the opposite reference.
vonedelmann likes this.
zarox is offline   Reply With Quote

Old   February 11, 2022, 03:05
Default
  #4
Member
 
Emeline Noel
Join Date: Dec 2013
Location: Paris
Posts: 31
Rep Power: 12
zarox is on a distinguished road
In addition,


It's properly working in parallel on serial with the above suggestion.
zarox is offline   Reply With Quote

Reply


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
Setting cyclicAMI Boundary conditions NARENDRAN OpenFOAM 0 July 8, 2020 10:14
Getting started with CyclicAMI Algebraist OpenFOAM 4 February 11, 2019 01:39
Setting cyclicAMI and cyclic Boundary conditions (ICEM Mesh to OpenFoam) bowen1024 OpenFOAM Pre-Processing 4 March 1, 2018 19:28
cyclicAMI issue for infinitely long circular cylinder doctorWho OpenFOAM Running, Solving & CFD 1 July 17, 2017 08:23
CyclicAMI with half cylinder/ Projection of boundaries OlliWer OpenFOAM Pre-Processing 0 December 10, 2013 09:34


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