[chrisb2244]

So it seems like the dnsFoam packaged with OF-2.2.2 doesn't support parallel running. Initially, it looks like the kMesh class is the one which creates some problems, but presumably these could be solved with not too much effort?

My question is essentially - is there some other (deeper) reason why dnsFoam does not work in parallel? I wanted to ask before altering lots of files and ending up with many slightly edited classes, just to find that the solver is for some (as now unknown to me) reason completely incompatible with parallel solving.

If the problem is only with the generation (and distribution?) of a kMesh, then could someone perhaps advise me on whether I should use some sort of

if (Pstream::master())
{
    // create a kMesh, and distribute to the other processes
    // possibly with some sort of scatter(..) function?
}

or alternatively, try and generate a kMesh for each subsection of a mesh, and only allow 1D decomposition (which I think should be ok?)

I am looking at this so that I can attempt to simulate turbulent pipe flow, and I'm just unsure as to the best way to go about it.

Any help or suggestions would be appreciated.