solver settings for reactingFoam
 

Another newbie's question, if I may impose on the forum users.

I am running reactingFoam on OpenFoam 1.5. The reactingFoam I downloaded was developed in v1.4 and the system keeps conplaining the input found in fvSolution. I did the following to get it to work but am quite clueless if what I did was even close to being right.

Any suggestions/comments will be grealy appreciated.

The original settings solvers choice and setting are:
solvers
{
rho ICCG 1e-06 0;
U BICCG 1e-06 0;
p ICCG 1e-09 0;
Yi BICCG 1e-06 0;
h BICCG 1e-06 0;
k BICCG 1e-06 0;
epsilon BICCG 1e-06 0;
}
and the one I put in to get it to work with v1.5 are:

solvers
{
rho ICCG
{
preconditioner DILU;
tolerance 1e-06;
relTol 0;
};
U PBiCG
{
preconditioner DILU;
tolerance 1e-06;
relTol 0;
};
p PCG
{
preconditioner DIC;
tolerance 1e-06;
relTol 0;
};
Yi PBiCG
{
preconditioner DILU;
tolerance 1e-06;
relTol 0;
};
h PBiCG
{
preconditioner DILU;
tolerance 1e-06;
relTol 0;
};
k PBiCG
{
preconditioner DILU;
tolerance 1e-06;
relTol 0;
};
epsilon PBiCG
{
preconditioner DILU;
tolerance 1e-06;
relTol 0;
};

Jerry,
Looks equivalent to me, unless I am missing something. Note that for rho in reacting foam, the fvSolution settings don't really matter, since you get a diagonal system from rhoEqn.H:

{
solve(fvm::ddt(rho) + fvc::div(phi));
}

-David