[tjliang]

Dear Foamers,

I wan to loop over all face values of a surfaceScalar field, and for some faces of them which could fulfill certain condition, I want to apply a value to its owner cell. Is it possible to do so in OpenFOAM? I have found that mesh().faceOwner() could possibly search for this, could anyone tell me a possible way to do it? Thanks.

Regards,

Peng

[bigphil]

Originally Posted by tjliang

Dear Foamers,

I wan to loop over all face values of a surfaceScalar field, and for some faces of them which could fulfill certain condition, I want to apply a value to its owner cell. Is it possible to do so in OpenFOAM? I have found that mesh().faceOwner() could possibly search for this, could anyone tell me a possible way to do it? Thanks.

Regards,

Peng

// Take a reference to the owner list: this gives the owner cell index for each internal face
// As a by the way, "faceOwner" is for all faces (internal and boundary) whereas "owner" is for just internal faces
const labelUList& owner = mesh.owner();

// Assuming you have some volScalarField called myField
// Take a reference to the internal field of myField, i.e. the raw list of cell values
scalarField& myFieldI = myField.primitiveFieldRef();

// Loop over all internal faces
forAll(owner, faceI)
{
    // Owner cell index
    const label ownCellID = owner[faceI];

    // Set value of owner cell in myField
    myFieldI[ownCellID] = 123456.789; // or what ever you want
}

// What about the boundary face values for myField? You may want to set something there too?
// In any case, you should update the processor boundaries:
myField.correctBoundaryConditions()