[MamboJambo]

Hello to all,


I would like to know in what conditions does the function correctBoundaryConditions() needs to be called?

[dlahaye]

Typically each type the equations and thus the matrix is updated within a S/Pimple loop.

[MamboJambo]

Originally Posted by dlahaye

Typically each type the equations and thus the matrix is updated within a S/Pimple loop.

Was not able to understand. Would you illustrate with some examples?
Are general rules of thumb to do this.


For instance, if you do:

        forAll(p, cellI)
        {
            p[cellI] = p[cellI] + 10*mesh.C()[cellI].x();
        }
        forAll(p.boundaryField(), patchI)
        {
            forAll(p.boundaryField()[patchI], faceI)
            p.boundaryFieldRef()[patchI][faceI] = p.boundaryField()[patchI][faceI] + 1010*mesh.Cf().boundaryField()[patchI][faceI].x(); 
        }

        U = fvc::grad(p)*dimensionedScalar("tmp", dimTime, 1.);

You should need p.correctBoundaryConditions() before the computation of U