[chefkoch89]

Dear community,

I'm working on the multiphase solver interFoam. I need to check, if the temperatures in the cells within a radius of 5 times the cell radius are in a special value range. After that I want to add their volumes either to the volScalarField named Vcons or to Vfeed.
Finally I need to calculate the ratio of these values.

Right now my code is:

forAll(mesh.cells(), i){

    float abstand;

    vector cellCenter = mesh.C()[i];
    labelList neighbours = mesh.cellCells()[i];

    forAll(neighbours, k){

        vector direction = mesh.C()[neighbours[k]] - cellCenter;
        abstand = mag(direction);
        
    }    


    for(int x=0; x < 5; x++){
            
        for(int y=0; y < 5; y++){

            for(int z=0; z < 5; z++){
        
                vector nachbar;
                nachbar[0] = cellCenter[0] + abstand*x;
                nachbar[1] = cellCenter[1] + abstand*y;
                nachbar[2] = cellCenter[2] + abstand*z;

                if(fluidNames[nachbar] == "fluid"){
                
                if( (Tcrit.value() <= T[nachbar]) && (T[nachbar] <= Trange.value()) ){
                    Vcons += mesh.V().field()[nachbar]);
                }
        
                else if( T[nachbar] > Trange.value() ){
                    Vfeed += mesh.V().field()[nachbar];
                
                }
    
                if(Vfeed[nachbar] > 0)
                fgeo = Vcons / Vfeed;

                }
    
            }        
        }
    }
}

My first problem is to check the temperatures of the neighbours of cell i and that the operand types are different. When compiling, I always get the error:

error: no match for ‘operator[]’ (operand types are ‘Foam::volScalarField {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ and ‘Foam::vector {aka Foam::Vector<double>}’)

The second difficulty is close to the boundary condition. I tried to solve this via the line

if(fluidNames[nachbar] == "fluid")

Do you got any proposals to solve these problems?

Thank you very much!!

Best regards

[santiagomarquezd]

Hi, cellCells() will give you only the first level neighbors of your cell so are probably not checking all the cells in the sphere you want. Please add in what lines are you having errors.

Regards!

[chefkoch89]

Ok, thank you for your response. I am now facing another problem with the neighbors of one cell: I need to get the velocities of the neighbors in each direction and store them in different new values.

Do you got an idea how to do that?

Thank you and best regards!

[santiagomarquezd]

Hi, if you have the neighbor cell numbers, you can do:

U.component(vector::X)[i]

Regards

[chefkoch89]

Hi, that is a good hint. Thank you!
Right now, I am struggling with the right definition of the values.

The picture is explaining what I want to do. Maybe you can help me out.

Thanks a lot!