[Rotidpor]

Hello,


I'm looking for a way to use the scalar value given by a surfaceFieldValue function as a fixedValue boundary condition after each iteration. I guess I can use codedFixedValue but I don't know how to write the piece of code that would do that.


As I don't know what to do, I'm trying to code a new BC that would be something like this :


#include "fvCFD.H"
#include "surfaceFields.H"

class fixedSurfaceScalarFieldValue : public fixedValueFvPatchScalarField
{
const surfaceScalarFieldValue& surfaceFieldValue_;

public:
fixedSurfaceScalarFieldValue(const fvPatch& patch, const DimensionedField<scalar, volMesh>& iF, const dictionary& dict)
:
fixedValueFvPatchScalarField(patch, iF),
surfaceFieldValue_(db().lookupObject<surfaceScalar FieldValue>(dict.lookup("functionObject")))
{
// empty
}

scalar value() const
{
return surfaceFieldValue_.value()(patch().Cf(), patch().index());
}
};

makePatchTypeField(fixedSurfaceScalarFieldValue);