[Geb1313]

Hello Foamers,
I want to solve a dam bream problem using interFoam. I divided the dam into certain compartments (say two). The requirement is two of them will not start flowing at the same time, say the 1st will start at t = 0 s and the 2nd will follow the 1st after 2 seconds.
Well, my first question is, is there any "easy way" to do it ?

My second question is I supposed to manipulate the solver using fvOptions by initiating the problem with setFields for the whole dam and constraining the cells inside the 2nd compartment until time = 2 s. Hence, I produced the following fvOptions dictionary:

velocitySource
{
    type          vectorCodedSource;
    active        yes;
    name          velocitySource;
    vectorCodedSourceCoeffs
    {
       selectionMode cellSet;
       cellSet box1;
       fields (U);
      
         codeInclude
            #{
            #};
        codeCorrect
        # {
        #};

        codeAddSup
        # {
        #};

        codeConstrain
        # {
            const  Time& time = mesh().time();            //get time 
            const  scalar startTime = 2.0;                // set value upto this time 
            cellSet  selectedCells(mesh_, cellSetName_);
            labelList  cells = selectedCells.toc ();
           if (time.value() <= startTime)
            {
            vector values(0, 0, 0); 
            eqn.setValues​(cells, values);
            }
        #};
  
        code
        # {
            $ codeInclude
            $ codeCorrect
            $ codeAddSup
            $ codeConstrain
        #};
    }
velocitySourceCoeffs
    {
    }
}

However, I got the following error:

 --> FOAM FATAL IO ERROR: (openfoam-2012)
[0] Entry 'selectionMode' not found in dictionary "/mnt/d/gebray/OpenFOAM/test/dambreak/constant/fvOptions.velocitySource.velocitySourceCoeffs"
[0]
[0]
[0] file: /mnt/d/gebray/OpenFOAM/test/dambreak/constant/fvOptions.velocitySource.velocitySourceCoeffs
[0]
[0]     From bool Foam::dictionary::readEntry(const Foam::word&, T&, Foam::keyType::option, bool) const [with T = Foam::word]
[0]     in file /home/pawan/OpenFOAM/OpenFOAM/OpenFOAM-v2012/src/OpenFOAM/lnInclude/dictionaryTemplates.C at line 334.

It is asking for keyword "SelectionMode" while it is there.
What did I miss?

Also, who is pawan in the file path

[olesen]

Originally Posted by Geb1313

Hello Foamers,
My second question is I supposed to manipulate the solver using fvOptions by initiating the problem with setFields for the whole dam and constraining the cells inside the 2nd compartment until time = 2 s. Hence, I produced the following fvOptions dictionary:

velocitySource
{
    type          vectorCodedSource;
    active        yes;
    name          velocitySource;
    vectorCodedSourceCoeffs
    {
       selectionMode cellSet;
       cellSet box1;
       fields (U);
      
         codeInclude
            #{
            #};
        codeCorrect
        # {
        #};

        codeAddSup
        # {
        #};

        codeConstrain
        # {
            const  Time& time = mesh().time();            //get time 
            const  scalar startTime = 2.0;                // set value upto this time 
            cellSet  selectedCells(mesh_, cellSetName_);
            labelList  cells = selectedCells.toc ();
           if (time.value() <= startTime)
            {
            vector values(0, 0, 0); 
            eqn.setValues​(cells, values);
            }
        #};
  
        code
        # {
            $ codeInclude
            $ codeCorrect
            $ codeAddSup
            $ codeConstrain
        #};
    }
velocitySourceCoeffs
    {
    }
}

You have fvOptions entries one level too deep.

 velocitySource
{
    type          vectorCodedSource;
    active        yes;
    name        velocitySource;


    selectionMode cellSet;
    cellSet box1;
    fields (U);



...  other entries    

}

[Geb1313]

Thank you olesen.

It is working now. I was mixing both openfom.com and openfoam.org manuals.
For openfoam.com (e.g OpenFOAM-v2012) your recommendation avoids the error.

[Geb1313]

Originally Posted by olesen

You have fvOptions entries one level too deep.

 velocitySource
{
    type          vectorCodedSource;
    active        yes;
    name        velocitySource;


    selectionMode cellSet;
    cellSet box1;
    fields (U);



...  other entries    

}

A quick question while I am at it. Is that possible to have multiple cellsets like:

cellSet         box1,box2;

So that I can access cells in both of the boxes at different time. I want to have two different times, t1 and t2 where t1 is when cells at box1 start to move and t2 is when cells inside box2 stat to move.


Gebray

[olesen]

Originally Posted by Geb1313

A quick question while I am at it. Is that possible to have multiple cellsets like:

cellSet         box1,box2;

So that I can access cells in both of the boxes at different time. I want to have two different times, t1 and t2 where t1 is when cells at box1 start to move and t2 is when cells inside box2 stat to move.

Gebray

That is a cool idea, but unfortunately not possible. At the moment you would have to use something like topoSet to create a cellSet box_1_2 (combining both).

No promises, but I think you should add this as an issue to develop.openfoam.com to keep your idea from being lost. Would give a chance to explore/discuss if your idea makes sense or how it could be added (or not).

[Geb1313]

Originally Posted by olesen

That is a cool idea, but unfortunately not possible. At the moment you would have to use something like topoSet to create a cellSet box_1_2 (combining both).

No promises, but I think you should add this as an issue to develop.openfoam.com to keep your idea from being lost. Would give a chance to explore/discuss if your idea makes sense or how it could be added (or not).

I somehow managed to do it. Instead of doing them in one, created box1,2, 3 using topSet and have separate dictionaries in fvoptions, like velocity source1, velocity source2, velocity source3, etc.

Gebray