[Kellis]

Foamers,

Good afternoon. I am in the process of modifying chtMultiRegionFoam (in OF 1806) to model an incompressible, two-phase fluid. I am doing this in several steps, and have already modified the solver for incompressible flow. Now, I am nearly done adding in the code from interFoam to allow for two phases to be solved, but am experiencing issues with the createAlphaFluxes.H file. When compiling, I get the following error:

./multiRegionVoF/createAlphaFluxes.H:47:59: error: array must be initialized with a brace-enclosed initializer
   alphaPhi10HeaderList[i].typeHeaderOk<surfaceScalarField>();

The issue is arising from my attempts to initialize the alpha fields for each fluid region. I have tried several variants of code, including:

// createAlphaFluxes.H

PtrList<IOobject> alphaPhi10HeaderList(fluidRegions.size());
PtrList<bool> alphaRestartList(fluidRegions.size());
PtrList<surfaceScalarField> alphaPhi10List(fluidRegions.size());
PtrList<surfaceScalarField> talphaPhi1Corr0List(fluidRegions.size());

forAll(fluidRegions, i)
 {

    alphaPhi10HeaderList.set
    (
        i,
        new IOobject
        (
            IOobject
            (
                IOobject::groupName("alphaPhi0", alpha1l[i].group()),
                runTime.timeName(),
                fluidRegions[i],
                IOobject::READ_IF_PRESENT,
                IOobject::AUTO_WRITE
            ),
            fluidRegions[i]
        )
    );
  
   // this line trips error
     const bool alphaRestartList[i] =
        alphaPhi10HeaderList[i].typeHeaderOk<surfaceScalarField>();

   // more stuff below here
     ...
  }

As well as:

...

    alphaRestartList.set
    (
        i,
        new const bool
        (
            "alphaRestart",
            alphaPhi10HeaderList[i].typeHeaderOk<surfaceScalarField>(true)
        )
    );

 ...

To be honest, I am not sure what this line of code is doing and have no real clue how to move forward other than blind trial and error. Any help diagnosing this issue would be greatly appreciated! I have attached my solver code as it stands currently for inspection. I have created a new multiRegionVoF folder within my solver directory holding the modified code for solving for alpha.

Thanks,
Kellis

[XinXin]

Hi,Kellis,have you solved this problem?

[Kellis]

Yes, I have had the code working for some time now. I'm not exactly sure what I changed to get past this issue, but below is the working code in the createAlphaFluxes.H file where I was having issues.

Best of luck,
Kellis

PtrList<IOobject> alphaPhi10HeaderList(fluidRegions.size());
PtrList<bool> alphaRestartList(fluidRegions.size());
PtrList<surfaceScalarField> alphaPhi10List(fluidRegions.size());

Info << "    Creating alpha fluxes\n" << endl;

forAll(fluidRegions, i)
{
    alphaPhi10HeaderList.set
    (
        i,
        new IOobject
        (
            IOobject
            (
                IOobject::groupName("alphaPhi0", alpha1l[i].group()),
                runTime.timeName(),
                fluidRegions[i],
                IOobject::READ_IF_PRESENT,
                IOobject::AUTO_WRITE
            ),
            fluidRegions[i]
        )
    );

    IOobject alphaRestartTmp = alphaPhi10HeaderList[i];

    const bool alphaRestartListTmp = alphaRestartTmp.typeHeaderOk<surfaceScalarField>(true);

    alphaRestartList.set
    (
        i,
        new bool(alphaRestartListTmp)
    );

    if (alphaRestartList[i])
    {
        Info << "Restarting alpha" << endl;
    }

    // MULES flux from previous time-step

    alphaPhi10List.set
    (
        i,
        new surfaceScalarField
        (
            "alphaPhi10Header",
            phiFluid[i]*fvc::interpolate(alpha1l[i])
        )
    );

    tmp<surfaceScalarField> talphaPhi1Corr0;

    Info << "    Done creating alpha fluxes\n" << endl;
}