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

Foam-3.0-ext + MixingPlane Problems

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 18, 2014, 13:42
Default Foam-3.0-ext + MixingPlane Problems
  #1
Member
 
don lenardo
Join Date: Nov 2010
Posts: 75
Rep Power: 13
lentschi is on a distinguished road
Hello community,

I have 2 perfectly matching and non-rotating/moving blocks, which are coupled by a mixingPlane Interface.

During calculation I always get an Floating Point Exception error. I have tried many different settings (and copied the settings from a mixingPlane tutorial) but get this error again:

Create mesh for time = 0

Initializing the mixingPlane interpolator between master/shadow patches: STAT1/STAT2
Reading field p

Reading field U

Calculating potential flow
CG: Solving for p, Initial residual = 1, Final residual = 0.00992722, No Iterations 82
continuity error = 9.22837e-06
Floating point exception (core dumped)

My boundary file looks like this:

6
(
INLET
{
type patch;
nFaces 171;
startFace 27788;
}
STAT1
{
type mixingPlane;
nFaces 171;
startFace 27959;
shadowPatch STAT2;
zone STAT1ZONE;
coordinateSystem
{
type cylindrical;
name mixingCS;
origin (0 0 0);
e1 (0 0 1); //Direction of Streamlines in z-Direction
e3 (1 0 0);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis Z;
discretisation bothPatches;
}

}
WALL1
{
type wall;
nFaces 1624;
startFace 28130;
}
OUTLET
{
type patch;
nFaces 171;
startFace 29754;
}
STAT2
{
type mixingPlane;
nFaces 171;
startFace 29925;
shadowPatch STAT1;
zone STAT2ZONE;

}
WALL2
{
type wall;
nFaces 1624;
startFace 30096;
}
)

The mesh is ok, as can be seen in checkMesh:

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
*Number of regions: 2
The mesh has multiple regions which are not connected by any face.
<<Writing region information to "0/cellToRegion"

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Area [m^2] Surface topology
INLET 171 200 0.1 ok (non-closed singly connected)
STAT1 171 200 0.1 ok (non-closed singly connected)
WALL1 1624 1680 6.6 ok (non-closed singly connected)
OUTLET 171 200 0.1 ok (non-closed singly connected)
STAT2 171 200 0.1 ok (non-closed singly connected)
WALL2 1624 1680 6.6 ok (non-closed singly connected)

Checking geometry...
This is a 3-D mesh
Overall domain bounding box (-3.46945e-18 -2.77556e-17 0) (0.1 1 6)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Mesh (non-empty, non-wedge) dimensions 3
Boundary openness (5.95673e-17 -9.5346e-18 2.69053e-18) Threshold = 1e-06 OK.
Max cell openness = 1.96445e-16 OK.
Max aspect ratio = 13.0944 OK.
Minumum face area = 0.0001. Maximum face area = 0.0140221. Face area magnitudes OK.
Min volume = 2e-06. Max volume = 0.000160253. Total volume = 0.6. Cell volumes OK.
Mesh non-orthogonality Max: 0 average: 0 Threshold = 70
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.24444e-09 OK.

Mesh OK.


Has anyone an idea to solve this problem?

Thanks in advance.

Regards,
Markus

Last edited by lentschi; March 31, 2014 at 08:46.
lentschi is offline   Reply With Quote

Old   March 19, 2014, 04:02
Default
  #2
Member
 
don lenardo
Join Date: Nov 2010
Posts: 75
Rep Power: 13
lentschi is on a distinguished road
Problem solved!

Definition problem in the boundary file!
lentschi is offline   Reply With Quote

Old   March 22, 2014, 10:27
Default
  #3
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: India
Posts: 205
Rep Power: 16
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
hi don lenardo,
I am using MRFSimpleFoam, 1/8 th of the fan. So far I am using Conformal Mesh, but now I want to reduce the computational time by using AMI interface,

my questions are:

is it possible in MRFSimpleFoam? I have attached the fig, please have a look. (the fan is in the middle)

I have edited (manually, I didt use createPatch utility) as follows:

constant/polyMesh/boundary

AMI4
{
type cyclicAMI;
nFaces 13000;
startFace 40391000;
neighbourPatch AMI3;
bridgeOverlap false;
}
AMI3
{
type cyclicAMI;
nFaces 26000;
startFace 40404000;
neighbourPatch AMI4;
bridgeOverlap false;
}
AMI2
{
type cyclicAMI;
nFaces 23400;
startFace 40430000;
neighbourPatch AMI1;
bridgeOverlap false;
}
AMI1
{
type cyclicAMI;
nFaces 13000;
startFace 40453400;
neighbourPatch AMI2;
bridgeOverlap false;
}

when I run check Mesh, I am getting the warning message as you given in the first post.

Can you please help me solve the problem?

I am using OF-2.1.1

mixing plane, is it available in OF2.1.1?

can you please explain the procedure?




kind regards,
Sivakumar
Attached Images
File Type: jpg Fully HexMesh.jpg (85.8 KB, 36 views)

Last edited by sivakumar; April 1, 2014 at 04:21.
sivakumar is offline   Reply With Quote

Old   March 31, 2014, 08:40
Default
  #4
Member
 
don lenardo
Join Date: Nov 2010
Posts: 75
Rep Power: 13
lentschi is on a distinguished road
Hello Sivakumar,

Mixing Plane is not available in OF 2.1!
Re-check your boundary file - maybe the boundaries are incorrectly defined and AMI 1 is not the neighbour of AMI2 (and AMI3 to AMI4 respectively).

Additionally, you have to add the following lines to each AMI-boundary:

transform noOrdering;
matchTolerance 0.001;


Hopefully this will help!

//Markus
lentschi is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
[Commercial meshers] Several problems with the mesh conversion utility when converting the meshes from Gridgen su_junwei OpenFOAM Meshing & Mesh Conversion 2 July 26, 2008 23:58
[OpenFOAM] Problems using paraview 244 under open FOAM 141 or 13 sergio ParaView 0 March 13, 2008 13:52
[Commercial meshers] FOAM FATAL ERROR points deallocated hoerl OpenFOAM Meshing & Mesh Conversion 0 October 30, 2006 07:47
Installation problems shellbell1999 OpenFOAM Installation 9 April 6, 2006 13:29
FOAM licensed free to academic institutions Nabla Ltd. Main CFD Forum 0 November 18, 2002 10:24


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