[Vader92]

Hi, I'm new to OpenFoam. I'm looking to develop an adjoint optimization solver in OpenFOAM. I have a mesh size of (20,20,1) .I added a new volScalarField vf which i defined in my 0 folder.

internalField nonuniform List<scalar>
400

I want to implement a coded boundary condition such that for each element on the boundary, if vf=1 , U=(5, 0 0) and if vf=0, U=(0,0,0). My implementation looks something like

inlet
{
type codedMixed;
refValue uniform (0 0 0);
refGradient uniform (0 0 0);
valueFraction uniform 1;

name inletBC;

code
{
const fvPatch& boundaryPatch = patch();
vectorField& field = *this;
vectorField vf= vf;

forAll(boundaryPatch, faceI)
{
field[faceI] = vector((1-vf[faceI])*0 + vf[faceI]*5,0,0);
}

}
}

I believe the bolded line is causing some errors as I don't think that's the correct way to get the volScalarField vf. I am unsure as to how I could access vf and use it in this BC. Any help will be appreciated.