problem with the ellipse curvededge
 

I want to use the ellipse curvededge and I write a simple test, the blockMeshDict is as follows

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

convertToMeters 1;

vertices
(
(1.2 0 -0.05)
(2 0 -0.05)
(2 0.42426 -0.05)
(0.84853 0.42426 -0.05)
(2 2 -0.05)
(0.84853 2 -0.05)
(0 2 -0.05)
(0 0.6 -0.05)

(1.2 0 0.05)
(2 0 0.05)
(2 0.42426 0.05)
(0.84853 0.42426 0.05)
(2 2 0.05)
(0.84853 2 0.05)
(0 2 0.05)
(0 0.6 0.05)
);

blocks
(
hex (0 1 2 3 8 9 10 11) (10 10 1) simpleGrading (1 1 1)
hex (3 2 4 5 11 10 12 13) (10 10 1) simpleGrading (1 1 1)
hex (7 3 5 6 15 11 13 14) (10 10 1) simpleGrading (1 1 1)
);

edges
(
ellipse 7 3 (0 0 -0.05) 0.86603
ellipse 3 0 (0 0 -0.05) 0.86603
ellipse 15 11 (0 0 0.05) 0.86603
ellipse 11 8 (0 0 0.05) 0.86603
);

patches
(

symmetryPlane down
(
(0 1 9 8)
)
symmetryPlane up
(
(6 5 13 14)
(5 4 12 13)
)
symmetryPlane ellipse
(
(7 3 11 15)
(3 0 8 11)
)
patch left
(
(6 7 15 14)
)
patch right
(
(4 2 10 12)
(2 1 9 10)
)

);

mergePatchPairs
(
);

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

and then I type blockMesh

it give some errors

Create time


Creating block mesh from
"/home/alvin/Linux/OpenFOAM/doctor/ellipse/constant/polyMesh/blockMeshDict"


Creating blockCorners

Creating curved edges


--> FOAM FATAL ERROR:
Root is not bracketed. f(x0) = -0.389993 f(x1) = -0.360002

From function Foam::scalar Foam::BisectionRoot<Func>::root
(
const scalar x0,
const scalar x1
) const
in file /home/alvin/OpenFOAM/OpenFOAM-1.6-ext/src/ODE/lnInclude/BisectionRoot.C at line 78.

FOAM aborting

Aborted

can someone help me ,
can someone give a case with the ellipse edge

problem with the ellipse curvededge
 

Hi,

Have you found a solution to this problem?

I am presently using Foam-Extend 3.1 and am struggling with creating an elliptic edge.

Sam

--------------------------
Creating curved edges


--> FOAM FATAL ERROR:
Root is not bracketed. f(x0) = -1.87498 f(x1) = -1.20396e-05

From function Foam::scalar Foam::BisectionRoot<Func>::root
(
const scalar x0,
const scalar x1
) const

blockMesh: Elliptic edge example
 

Does anyone have an example case of blockMeshDict with an elliptic edge?

Sam

Hi Sam,

I have just an example with an circular arc, but an elliptic edge should work in nearly the same way.

Is the simpleSpline like in the dokumentation (http://www.openfoam.org/docs/user/blockMesh.php) right for you?
You just have to set the right interpolation points.

greetings
Alex

Dear Sam,

I recently have the same question and find a way to build it. Are you still interested about this? This attached picture is my mesh.

Dear Chun

I would be greatfull if you could provide me the way you build the mesh with the ellipse

Regards
John