[MdoNascimento]

Hello,

I am trying to work with a wedge geometry with a solver based on interFoam. One phase has the shape of a sphere in the beginning.

So my geometry is very similar to this one:
https://openfoamwiki.net/index.php/M...C4.2C5.2C2_.3F

When I was working with two- and three-dimensional meshes I used to initialize my alpha1 phase with the cylinderToCell resp. sphereToCell functions in the setFields file.

So, now that I have changed to the wedge shape, those functions don't seem to result in a reasonable cylindric or spherical shape.

My mesh:

convertToMeters 1.0e-2;

vertices
(
(0 -0.5 0) // V00
(0.49952 -0.5 0.02181) // V01
(0.49952 0.5 0.02181) // V02

(0 0.5 0) // V03
(0.49952 -0.5 -0.02181) //
(0.49952 0.5 -0.02181) //
);

blocks
(
hex (0 4 1 0 3 5 2 3) (350 700 1) simpleGrading (1 1 1)
);

patch top
(
(3 5 2 3)
)

patch bottom
(
(0 4 1 0)
)

patch walls
(
(1 4 5 2)
)

empty achse
(
(0 3 3 0)
)

wedge wedge1
(
(0 1 2 3)
)

wedge wedge2
(
(0 3 5 4)


My setFields

cylinderToCell
{
p1 (0 0 -0.001);
p2 (0 0 0.001);
radius 0.001;

fieldValues
(
volScalarFieldValue alpha1 1

[MdoNascimento]

The result:

It looks like the cylinder was rotated 90°, but I am sure it is not. The cylinder's axys should point in z-direction.

[MdoNascimento]

I installed swak4foam to ues funkySetFields:

funkySetFields -field alpha1 -expression 0 -time 0 -keepPatches -condition "pow(pos().x,2) + pow(pos().y,2) < pow(0.001,2)"

SAME result ...

any ideas?