[leroyv]

Dear foamers,

I am trying to set patch values of a given field equal to the corresponding internal field values. This shouldn't be too complicated, but it seems to be not as trivial as I would have expected.

Here is what I wrote (let's name our volScalarField 'myField'):

const fvPatchList & patches = mesh.boundary();
forAll(patches, iPatch)
{
    UList<scalar> myList(myField.boundaryField()[iPatch].patchInternalField());

    forAll(myField.boundaryField()[iPatch], iCell)
    {   
        myField.boundaryField()[iPatch][iCell] = myList[iCell];
    }
}

I end up with all patch values set as expected, except a few ones. I guess I must have done something wrong. Any ideas?

[alexeym]

Hi,

what's wrong with fixedInternalValue patch type? It does the thing this way:

template<class Type>
void Foam::fixedInternalValueFvPatchField<Type>::manipulateMatrix
(
    fvMatrix<Type>& matrix
)
{
    // Apply the patch internal field as a constraint in the matrix
    matrix.setValues(this->patch().faceCells(), this->patchInternalField());
}