[MMRC]

Good day,

I wanted to implement some new features to actuationdisksource in Openfoam ESI 2012 and 2306 versions. I wanted to check if the scalar T computed was of the same value that the total amount applied in Usource[celli], so I added some new line codes to verify if there was a balance.

I checked it by adding to each method (Froude and variableScaling) something like:

     // (BJSB:Eq. 3.9)
    const scalar a = 1.0 - Cp/Ct;
    const scalar T = 2.0*rhoRef*diskArea_*magSqr(Uref & diskDir_)*a*(1 - a);
    scalar magSource = 0.0;
                scalar cellInCellSet = 0.0;
                scalar volumenCheck = 0.0; 
    for (const label celli : cells_)
    {
        Usource[celli] += ((cellsV[celli]/V())*T)*diskDir_;
        magSource += mag(Usource[celli] & diskDir_);
        cellInCellSet += 1.0;
        volumenCheck += (cellsV[celli]/V());
    }
                        /* DEBUG PART IV*/
                Info << "the variable cells_ is = " << cells() << endl;
                Info << "Volumen check via (cellsV[celli]/V()) " << volumenCheck << endl;
                Info << "Volumen check via V() " << V() << endl;
                Info << "Total cells check via cellInCellSet " << cellInCellSet << endl;
                Info << "Total thrust " << T << endl;
                Info << "Total Usource " << magSource << endl;

By running turbineSiting tutorial of ESI2306, the cellSets are of 484 cells but the counter 'cellInCellSet ' inside the for-loop finds 374, that is the same value seen when printing on screen cells(), that at the end is a ratio of 0.77 seen in volumenCheck. This also ends that magSource does not match with T



What could be the reason that original fvOption is not receiving the full list of cells selected in the cellZone/cellSet?


Or are there issues in my line codes?

[MMRC]

I'd say, issue solved! by adding this, magic happens:
reduce(magSource, sumOp<scalar>());
A new thing learned