[kcn]

Hi,

I used the following piece of code to change the boundary condition of a certain patch during each solver iteration.

label patchID = mesh.boundaryMesh().findPatchID("name of patch");

const polyPatch& cPatch = mesh.boundaryMesh()[patchID];


forAll (cPatch,faceI)
{

U.boundaryField()[patchID][faceI] == expression;

}

I need to define a source term which is only present at the cells that are attached to a particular patch. I hope to do this by first declaring a field that is every where zero and then during solution updating the values attached to the patch of interest, similar to updating the boundary face values in above code.

Can some one please tell me how to modify the code to access and change the values of cells attached to a boundary patch?

[fumiya]

Hi,

You can use faceCells() and might want to read my slide(http://www.slideshare.net/fumiyanoza...ogramming-tips).

Best regards,
Fumiya