[frobaux]

Hey foamers,
I'm currently trying to develop a solver based on chtMultiRegionFoam,


I would like to solve the equations from interFoam in the fluids zones,


for now I import fields following the logic of multiRegiondFoam, for example in createFluidFields.H:

PtrList<volVectorField> U(fluidRegions.size());
...
 U.set
    (
        i,
        new volScalarField
        (
            IOobject
            (
                "U",
                runTime.timeName(),
                fluidRegions[i],
                IOobject::MUST_READ,
                IOobject::AUTO_WRITE
            ),
            fluidRegions[i]
        )
    );

The problem is "setRefCell" needs pimple.dict()

setRefCell
(
    p[i],
    p_rgh[i],
    pimple.dict(),
    pRefCell[i],
    pRefValue[i]
);

And compilation returns an error as pimple is not defined in this scope.
when I comment the

#define NO_CONTROL

That part is ok to compile.

What is the use of NO_CONTROL, and why does it prevent pimple to be loaded?

Moreover, from what I know from C++, the #define should have two operands, STRING and REPLACEMENT_STRING right? so NO_CONTROL is replaced by nothing?


Thanks for your help
Fabien

[raumpolizei]

Hey,
Thanks for asking this question as I could also learn something valuable for myself .
I hope this part of the code will help you reason about the problem (from createControl.H)

 #if defined(NO_CONTROL)
 #elif defined(PISO_CONTROL)

     #include "createPisoControl.H"

 #elif defined(PIMPLE_CONTROL)

     #include "createPimpleControl.H"

 #elif defined(SIMPLE_CONTROL)

     #include "createSimpleControl.H"

 #endif

[frobaux]

Oh, thank you for answering,

This is strange that "CONTROL" has to do with Pimple (which is defined in fvSolutions)
At the end, and because I don't completely get It, I have decided to use the PIMPLE dict located on the region directory (because I don't have a general mesh to construct a pimpleControl instance)
The only thing in the "main" pimple (in system/fvSolution) will be the nOuterCorr, which, from what I understand, are common for all regions (at least in chtMultiRegionFoam). This number is imported manuallly without the use of a pimpleControl instance.

Only I have to do in every loop over the regions (which involve a pimple) an include of

 pimpleControl pimple(fluidRegions[i]);

if( pimple.dict().found( "nOuterCorrectors") )
{
    FatalError << "entry nOuterCorrectors is present in dictionnary of the subregion '" 
        << fluidRegions[i].name() << "', It should be defined in the 'global' PIMPLE dict (in constant/fvSolutions)" << exit(FatalError);
}

Another way to do it would be to instanciante List of pimpleControl, but I didn't managed to do it.



Anyway with those method just hope that only the main loop that I will write manually will use nOuterCorrectors.



Thanks again for your help!