[charles4allme]

Hi,

Hopefully I can make this question easy to understand. I want to implement a boundary condition but I want to do it locally, meaning i would like to implement it on a cell-by-cell basis. I believe that I have to loop through the cells in order to implement the boundary condition. So far what I have is this:

forAll(patch(),facei)
{
   label celli = patch().faceCells()[facei];
   scalar mflux[celli] = {0};
        mflux[celli] += hm*(((avgPatch*WF)/(R*avgTs))*(log((avgYfg-1)/(XFl-1))));
    Info<<"mass flux is=" <<mflux<<nl<<endl;
    scalar avgU[celli] = {0};
    avgU[celli] = mflux[celli]/avgRhocomp;
    
Info<< "avgU="<< avgU<<nl<<endl;
   // tmp<vectorField> n = patch().nf();
   vectorField n = patch().nf();
 Info<< "n = " << n <<nl << endl;
Info<< "patchName"<< patch().name()<<nl<<endl;

    operator==(n*avgU[celli]);

    fixedValueFvPatchField<vector>::updateCoeffs();
    }

The values used to calculate mflux are scalar values that was defined before the loop. When i run the code I get value of 1 every iteration. I was wondering if anyone could advise me on what am doing wrong.

Thanks

[HPE]

scalar mflux[celli] = {0};
scalar avgU[celli] = {0};

scalar type is not a vector type, so you cant use [] operator.

{} initialization and = operator: you cant use both of them.

i dont think this can compile operator==(n*avgU[celli]);

[charles4allme]

Originally Posted by HPE

scalar mflux[celli] = {0};
scalar avgU[celli] = {0};

scalar type is not a vector type, so you cant use [] operator.

{} initialization and = operator: you cant use both of them.

i dont think this can compile operator==(n*avgU[celli]);

Hi HPE,

I am not very experienced with the C++ code but from your reply I have a couple of questions:

-First of all, I was able to compile the code like I said earlier but I get a value of 1 for mflux and avgU every iteration.
- What operator am I meant to use then because I believe what am defining is an array containing the number of cells in the patch.
- How do I initialize the array then.

Regards