[amattosc]

Hi everyone,

I was hoping to get some help here about this solver.

I've tried to run the original tutorial "bubble2D_r0.75mm" and it worked fine, but I noticed that the option "cleanInterface" was set as "Yes" in "freeSurfaceProperties".

I've looked into the code and for what I understood with this option enabled, BubbleInterTrackFoam doesn't solve the surfactant's conservation equation (as expected).

So, I've tried to change this option to "No" with the add of "soluble yes" at the end of "freeSurfaceProperties". I've also added an initial field for the concentration of surfactant as a "volScalarField".

The result was the following error:

"--> FOAM FATAL IO ERROR:
cannot open file

file: /home/amanda/OpenFOAM/amanda-1.6-ext/run/tutorials/surfaceTracking/bubbleInterTrackFoam/bubble2D_r0.75mm-teste/1e-05/Cs at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 62.

FOAM exiting"

I've tried to locate where the solver tries to read Cs and found its declaration in the middle of "freeSurface" in the file "makeFreeSurfaceData.C":

"surfactConcPtr_ = new areaScalarField
(
IOobject
(
"Cs",
DB().timeName(),
mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
aMesh()
);"

But at this point the solver had already entered the time loop and was trying to read from a file that has not been created yet.

It might be a silly mistake because I'm a new user of OpenFoam, but could someone give me a hand?

Thanks in advance.

Amanda Mattos

[akidess]

Create the file yourself to specify the initial surfactant distribution.

[amattosc]

Originally Posted by akidess

Create the file yourself to specify the initial surfactant distribution.

Thanks for the answer but as I said in my previous post I've already created this file.

The problem is that the solver tries to read this concentration after it started the time loop so the concentration for time 0 exists but for the next time 1e-5 it doesnt.

Any more ideas?

[akidess]

Ah, yes. You said you created a volScalarField though, whereas what you need is an areaScalarField?

[amattosc]

Originally Posted by akidess

Ah, yes. You said you created a volScalarField though, whereas what you need is an areaScalarField?

I've tried to create a Cs file as an areaScalarField but it didn't work anyway.

[Chia]

Hi Amanda,

I have the same problem with bubbleInterTrackFoam and the field Cs.
Did you solve the problem then? If yes, can you tell me how?

Thank you!

Chiara

[amattosc]

Originally Posted by Chia

Hi Amanda,

I have the same problem with bubbleInterTrackFoam and the field Cs.
Did you solve the problem then? If yes, can you tell me how?

Thank you!

Chiara

Hello Chiara,

Try calling this function "interface.surfactantConcentration();" before line 22 in file "createSurfactantConcentrationField.H" as shown in the code below and then recompile the solver.

I hope it helps!

// Volume surfactant concentration

volScalarField* CPtr = NULL;

if(!interface.cleanInterface() && interface.surfactant().soluble())
{
    Info << "\nReading field C\n" << endl;
    CPtr = new volScalarField
    (
        IOobject
        (
            "C",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    );

    interface.surfactantConcentration();

    if(runTime.timeIndex() == 0)
    {
        const dimensionedScalar& CA =
            interface.surfactant().surfactBulkConc();

        dimensionedScalar CB("CB", CA.dimensions(), 0);

        (*CPtr) = fluidIndicator*(CA - CB) + CB;

        if
        (
            CPtr->boundaryField()[spacePatchID].type()
         == inletOutletFvPatchScalarField::typeName
        )
        {
            inletOutletFvPatchScalarField& spaceC =
                refCast<inletOutletFvPatchScalarField>
                (
                    CPtr->boundaryField()[spacePatchID]
                );

            spaceC.refValue() = CA.value();
        }
        else
        {
            FatalErrorIn(args.executable())
                << "Boundary condition for bulk surfactant concentration at "
                    << "space patch isn't and should be: "
                    << inletOutletFvPatchScalarField::typeName
                    << abort(FatalError);
        }

        CPtr->correctBoundaryConditions();

        CPtr->write();
    }
}

[Chia]

Hi Amanda,

sorry for the delay of my feedback. Your correction worked fine, thank you!

Best,

Chiara

[Johnzi]

Hi,
Could you please advise how to recompile the code?
Thanks

[ali_jan]

Hi Amanda
thank you for your advice
Thanks