[Ben UWIHANGANYE]

Hello Foamers,

i need your assistance, it may sound easy for you but it has challenged me as beginner.

I want to map a cylinder on my fluid domain created using blockMesh. The fluid domain has one block as you may see it bellow:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.1;

vertices
(
(200 -50 0)
(200 50 0)
(-25 50 0)
(-25 -50 0)

(200 -50 1)
(200 50 1)
(-25 50 1)
(-25 -50 1)
);

blocks
(
hex (0 1 2 3 4 5 6 7) (100 225 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(2 3 7 6)
);
}

outlet
{
type patch;
faces
(
(0 1 5 4)
);
}

topAndBottom
{
type patch;
faces
(
(1 2 6 5)
(0 3 7 4)
);
}

frontAndBack
{
type empty;
faces
(
(0 1 2 3)
(4 5 6 7)
);
}
);

mergePatchPairs
(
);

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

I have been reading online posts and thought that topoSetDict/ cylinderTocell can help to map that cylinder in that Cartesian grid.

The cylinder i want has: radius: 0.25; p1 (0 0 0), p2 (0 0 1).

May you help to know what i am missing in the attached commpressed case setup, i am not getting that cylinder after ./Allrun in paraview?

I will appreciate you support.

Best regard!

[lebc]

Hi Ben,

Your code seems to be Ok, maybe you forgot to enable the cellzone visualization in paraView, you can find where to click on this picture ("Include Zones"):

https://goo.gl/images/FVxs3p

I didn't run your case, so I may be wrong.

Are you trying to use it to simulate the flow over a cylinder or you have anything else in mind?

Best Regards,
Luis

[Ben UWIHANGANYE]

Yes, my ultimate goal is to simulate flow over that cylinder? (stationary and oscillating).
I thought that cylinder will be immersed not body fitted? Does it sound feasible?

Finally I got that cylinder in my paraview(see attached images), then is there away i can move that cylinder?(oscillating motion). which codes should i use?

Say if the solid region eta=1 and fluid region =0, is there a way i can alternate eta to move that cylinder or?

Just advise me how i can move that cylinder without dynamic mesh.


Thank you!

[lebc]

Hi Ben,

I'm not sure if this is the best approach for an oscillating cylinder... I've never worked on this kind of problems, so I can't rally help you on that, sorry!

Best Regards,
Luis