[babakflame]

Dear Fellows

I want to redifine fixedGradient boundary condition for a specified patch inside my source code. i.e. whenever the code sees "fixedGradient" in 0/cPlus file (cPlus is a volScalarField and a dependent variable)

My snippet now looks like:

const polyPatchList& patches = mesh.boundaryMesh();
   forAll (patches,patchI){
     
      if(cPlus.boundaryField()[patchI].type() == 'fixedGradient'){


           fixedGradientPatchField& fgPatchField
         (
          refCast<fixedGradientPatchField>(cPlus.boundaryField()[patchI])
         );
          fgPatchField.gradient() = fvc::snGrad(cPlus) + cPlus * fvc::snGrad(Ue);
       }
    }

In above snippet, I want to tell the code that whenever it saw "fixedGradient" for a patch inside 0/cPlus file calculate the following relation at the boundary surface and assign it as boundary condition:

in which & are dependent scalar variables that should be calculated at the boundary patch from previous time step.

However the code compilation returns this error:

createFields.H:501:3: error: ‘fgPatchField’ was not declared in this scope
In file included from interFoamEHDIon.C:79:0:

My first question is How should I declare fgPatchField ? since I tried volScalarField or volVectorField but did not work. Moreover it does not accept fvc::snGrad. Should I try fvc::reconstruct(snGrad(cPlus))?


The second error is:

createFields.H:532:12: error: ‘fixedGradientFvPatchVectorField’ was not declared in this scope

createFields.H:533:74: error: no matching function for call to ‘refCast(Foam::fvPatchField<Foam::Vector<double> >&)’

By the way, My OpenFoam version is 2.1.x