[CoolHand]

Foamers,

My case is a large mesh (~450e6) which causes blockMesh to fail (bad number). I got around that problem by starting with a coarse mesh and using refineMesh -parallel after decomposePar. I'm using timeVaryingMappedFixedValue at the inlet and when I got to run pisoFoam I get an error like this:

[3357]--> FOAM FATAL IO ERROR:
[3357] size 501 is not equal to the given value of 2004
[3357]
[3357] ... /processor3357/0/U.boundaryField.inlet at line 58 to 60.
[3357]     From Foam::Field<Type>::Field(const Foam::word&, const Foam::dictionary&, Foam::label) [with Type = Foam::Vector<double>; Foam::label = int]
[3357]     in file [3202]

...for each processor on the inlet. The numbers are different for each processor since I'm using scotch decomposition but they're always 1/4 ratio so clearly the problem is that the assignment that timeVaryingMappedFixedValue gets from blockMesh isn't updated by refineMesh.

I've tried several combinations like:

decomposePar  
refineMesh 
decomposePar -fields

decomposePar  
refineMesh 
decomposePar -copyZero

decomposePar -noZero
refineMesh 
decomposePar -copyZero

I've also tried

decomposePar
refineMesh
reconstructPar

...but reconstructPar gives a similar error. So far nothing has worked. Is there another tool or option that can update the decomposed inlet BCs to match the refinedMesh? I would also accept advice on how to get blockMesh to do the large mesh in the first place, although I'm skeptical that there is such a way since I'm running it on a large memory cluster node and it still fails. Thanks!