[bou]

Dear,

I have developed a simple solver:

solve
        (
            S*fvm::ddt(H) == fvm::laplacian(K,H)
         );
U = -K*fvc::grad(H);

I precise that I am using groovyBc at the same time.

Everything works fine, however when I am setting K (volScalarField) using this from the solver :

 forAll(K,i){
     K[i] = 1.0;
}

OpenFoam computes only the first time step and just skip the rest without errors or warnings. I mean that for each time step it writes something like this in the log file as nothing changed in the BC:

DICPCG:  Solving for H, Initial residual = 2.64583e-08, Final residual = 2.64583e-08, No Iterations 0

I precise that I have tested the problem and this really comes form the internal setting of K, if I use the SetFieldDict instead and just comment the previous loop over K, everything works fine.

Thanks,

Bou