[Yujuan]

Hi all,


I want to make special treatment to the temperature on the wall boundary in the solver. And I wrote in the solver the following code:

forAll (T.boundaryField(),patchi)
{
if (mesh.boundary()[patchi].type() == "wall")
{
const fvPatchScalarField& Ts = T.boundaryField()[patchi];
forAll(Ts, facei)
Ts[facei] = specific formula;

}

}

The solver compiled without problem. However, it just crashes when running with error information about mpi:

I can't figure out the reason. What should I do to avoid this problem? Thanks a lot!