Hello everyone,

I’m studying the Poiseuille3D test case, trying to build a model very similar to the micron scale, and I’ve found myself asking some questions about the assignment of material numbers.

1. At the very beginning the material number 1 (internal domain) is set on an enlarged cylinder using an origin point of center0(-converter.getPhysDeltaX() * 0.2, radius, radius);

Why this 20% of displacement is added?

2. In the forced case this quantity is expanded by 3 cells on each side. Why exactly three?

a. The voxel are intended by cells and assigned by point values on the direction vertices?

3. In the non-forced case to describe the inlet (and similarly the outlet), two cylinders that are 4 cells thick are used (with origin[0] = -converter.getPhysDeltaX() * 2; and extend[0] = converter.getPhysDeltaX() * 2;).

Then, if I understood correctly, using the function superGeometry.rename(2, 3, 1, inflow), the material number is assigned by exploring only the two neighboring voxels in the chosen direction. The goal is to obtain a surface, so why is a cylinder that is 4 cells thick used instead of directly defining a surface?

4. Going back to the previous question (2a), does the assignment on neighboring voxels made by superGeometry.rename(2, 3, 1, inflow) explore points or cells?

If the assignment is made on the cells, when is it converted to point values?

I suspect that when using Paraview, cell values are somehow shifted to the points in the geometry visualization. I can’t find an exact match between what I see and what I expect in the code. If this is the case, I apologize in advance.

Thank you in advance for any clarification.

Alice