[Vitor Monteiro]

I'm using multiphaseEulerFoam to simulate a fluidised bed reator. This is part of the code I'm using to quantify the components conversion of all finite volumes of my mesh to handle the mass balance calculations.
However, I would like to program into my multiphaseEulerFoam.C the mass fraction composition of the gas flux of the Outlet boundary, and implement the mass flow of each species into the mass balance calculations.
Could anyone help me do this? I don't know what functions to use or how to integrate the mass flux to the outlet face.

const scalarField& v = mesh.V();
forAll(fluid.multiComponentPhases(), multiComponentPhasei) 
 
    {
        phaseModel& phase = fluid.multiComponentPhases()[multiComponentPhasei]; 
        UPtrList<volScalarField>& Y = phase.YRef();
        const volScalarField& rho = phase.rho();
    
        forAll(Y, i)
        {
             const volScalarField& alpha = Y[i];                    
             alphaYrho[ii] =  gSum(alpha*phase*rho*v);  // mass of each species        

             const scalarField& outletPatchPhi = phi.boundaryField()[1];  // determines the Outlet patch
              outletPatchPhi = gSum((alpha*phase*rho)*phi.boundaryField()[1]);  // I don't understand this phi very well. My system has phi, phi.air, phi.sand and phi.particles fields
             Info << "    " << alpha.name() << " = " << outletPatchPhi << " kg/s" << endl;

         } 

   }

The red lines doesn't work, once I'm a new foam programmer, but it has the idea I was wondering about.
Thanks a lot!


I know I could do this at controlDict, but I just want to program it into the solver to make further calculations.

[Vitor Monteiro]

I have already made some progress with this at my multiphaseEulerFoam.C file:

const scalarField& phasex = phase.boundaryField()[outletPatchID]; 
const scalarField& Yip = alpha.boundaryField()[outletPatchID];   
const scalarField& rhoi = rho.boundaryField()[outletPatchID];  
const scalarField& outletPatchPhi = phi.boundaryField()[outletPatchID];
scalar outletMass = gSum(phasex*Yip*outletPatchPhi); // without rhoi just to match with controlDict tool
Info << "    " << alpha.name() << " = " << outletMass << " kg/s" << endl;

However, it seems like it doesn't matches with the controlDict results duirng simulation.
controlDict.C file I've set like this:

gasFlux
    {
        type            surfaceFieldValue;
        libs            ( "libfieldFunctionObjects.so" );
        writeControl    timeStep;
        log             true;
        writeFields     false;
        regionType      patch;
        name            outlet;
        operation       weightedAreaIntegrate;
        weightField     phi.air;
        fields            
        ( 
            H2O.air
            CO2.air
            CO.air
            H2.air
            O2.air
            N2.air
            CH4.air
        );
    }

Could anyone help?