CFD Online Discussion Forums - sonicFoam and reconstructPar

Hello all,

I am having some trouble using reconstructPar and sonicFoam. I have run a case in sonicFoam with an airfoil that changes pitch, and now I wish to use the results as the starting point for another case. I tried the mapFields utility, however there was an error on the k and omega boundary conditions in the mapped files. I traced the error to the reconstructPar utility. I ran the case in parallel, and in the directory for each of the time steps the wall functions for omega and k are the same as for my time 0 BCs: compressible::omegaWallFunction and compressible::kqRWallFunction, respectively. However, when the time steps are reconstructed, the BCs are changed to omegaWallFunction and kqRWallFunction. Therefore, I must edit the files manually before continuing to run the case. Does anyone know why that is, or how it can be fixed?

Problem 2: after fixing the files manually, the case diverges after one timestep with the famous error:

Code:

--> FOAM FATAL ERROR: 

Maximum number of iterations exceeded



    From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const

    in file /home/opencfd/OpenFOAM/OpenFOAM-2.0.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67.

I find this very odd, since it is the same case that I am trying to continue. Somehow, the reconstruct/decompose cycle makes it stop working. Does anyone know why this might be, or how it can be fixed?

Thanks in advance,
Dan