[heather]

Hi,

I'm attempting to create a 2d axisymmetric (grid) with swirl case, for use with interFoam. I've created a simple geometry of a 5deg sector of a cylinder, one cell thick. The initial volume fraction is zero everywhere.

Boundary conditions:
- imagine a rectangle: left = inlet, right = pressure, top = wall, bottom = empty (axis)
- front and back boundaries declared as cyclic

Looking at the engineSwirl apps. I can see how to introduce the initial swirl (have added the swirl axis, point through which the axis passes, and a swirl component - defined as a frequency in Hz, to determine a volVectorField due to the swirl contribution). However - how do I include the swirl component in the equation set (for run-time calcs, not just initial U field setting)? - the swirl component is static wrt time.

Also - now that all cells have an associated constraint due to the boundary condition set, I'm getting the error message:

--> FOAM FATAL ERROR : Cannot find a cell not on a constraint boundary starting from cell 0

- Any ideas of how i can get around this - only route I can see is to include additional layers of cells?

Many thanks,

Andy