CFD Online Discussion Forums - Permeability filling of container outlet condition

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Permeability filling of container outlet condition (https://www.cfd-online.com/Forums/openfoam-solving/59002-permeability-filling-container-outlet-condition.html)

Hi, I have 2 phases with in

Hi,

I have 2 phases with interFoam. I suppose they're both incompressible so I need an outlet in order to fill the container entirely with water. The shape is U-formed, so the problem is that the container will not be filled completely since now both phases are led out from outlet-section, *before* the container is entirely filled with water.

I want only to allow air to be led out, trapping the water inside the container. It's like applying permeability to the surrounding walls and I want to compare flow result with another program in which the permeability around the container is 250 cm^3/min. Perhaps it's not possible to apply a condition like this to Foam but any suggestions that will give comparable results are most welcome.

How do I do this? My outlet boundary condition so far is:

U:
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}

gamma:
outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

pd:
outlet
{
type fixedValue;
value uniform 0;
}

Well, question asked in anothe

Well, question asked in another way: What would you do to simulate filling of a container completely? Don't tell me it can't be done in Foam...

Well, a primitive suggestion.

Well, a primitive suggestion. Make some few short outlets in walls of a container. Walls of those outlets should be hydrophobic (consult Henry, whether you can set the contact angle to 180).
Dimensions of the outlet you may want to compute from capillary pressure.

Thanks Sergei, Henry, what

Thanks Sergei,

Henry, what are your comments on this problem?

Hi, I am also curious about t

Hi,
I am also curious about this question. Does anyone have any idea?
Thank you.

I am interested in this issue.

I am interested in this issue. Any developments?

Hi all, I would like to si

Hi all,

I would like to simulate the filling of a container, but don't know which boundary conditions to choose.

I've put inflow on the left side of my box, wall on the upper, lower and right side, empty on the front and back side (since i want to simulate a 2D domain).

I initialize the first cell on the left side to gamma = 1 and the other to gamma = 0.

My simulation fail to start, probably for the BC, How can I solve this issue? Does anyone have any suggestions?

Thanks in advance

Francesco

Hi Francesco, if you want t

Hi Francesco,

if you want to fill something with something which allready contains something, you need to let that something out or compress it really much.

br

/eric

Hi Eric, I've put an outlet

Hi Eric,

I've put an outlet at the right side of my box and now I succeed in filling it http://www.cfd-online.com/OpenFOAM_D...part/happy.gif.
Now the problem is that the outlet surface doesn't behave as a wall for the fluid I have let in (gamma = 1). I'd like this surface to behave as a filter letting pass only the fluid with gamma = 0. How could I reach this aim?

Thank you in advance

Francesco

Hmm... That would mean a BC

Hmm...

That would mean a BC which has zeroGradient for the volocity for gamma=0 and wall for gamma > 0. Could be implemented but doesn't exist to my knowledge.

Try drilling a hole in the upper wall to let the air out. Or change the upper wall BC to atmospheric.

Good luck

//Eric

Hi Eric, I cannot drill a h

Hi Eric,

I cannot drill a hole or put atmospheric BC since it is an injection molding process; this kind of BC would change the behaviour of the fluid. Any hint about how to implement my condition? Or about how to solve my problem?

Thanks in advance

Francesco

Hi Francesco! I think the s

Hi Francesco!

I think the solution to your problem would be to implement a boundary condition that is based on the mixed-condition and changes the valueFraction according to gamma. There are conditions in the sources that do similar things (have a look at inletOutlet) and use them as a template

Bernhard

Hi Bernhard, I could implem

Hi Bernhard,

I could implement this boundary condition as you suggest but I think this would not solve my problem since the flow is modified by the presence of the outflow and the result is not as I wish.

I've discuss about this problem also in the thread: "Injection of a fluid into a cavity"

You can see there my opinion about the whole problem. I appreciate any comments

Thanks in advance

Francesco

Hi Francesco! The question

Hi Francesco!

The question is: where does the air get out in "real life"? Through the inlet (as some backflow in the gamma=0 part you were describing above)? And what kind of condition are you using for the velocity there (propably fixedValue)?

Maybe a "gamma"-switched boundary conditon for the velocity (fixedValue for the parts with gamma=1, zeroGradient for the rest) would be the solution

Bernhard

PS: the "simple"=="no programming needs to be done" solution would be to split the inlet condition in two parts - a "lower" one with gamma=1, velocity fixedValue and a "upper" with gamma=0, velocity zeroGradient. But maybe I misunderstand your problem completely - could you provide us with a sketch?

Hi Bernhard, I'd like to us

Hi Bernhard,

I'd like to use OF to simulate the injection of a fluid into a cavity which is empty, this is what molding softwares actually do. The presence of air makes the simulation more expensive and complicated. I've chosen the interFoam solver for the interface tracking, but I've had outlet problems, this seems to be the wrong way.
Probably it could better to use a single fluid solver, but I don't know which solver could perform this tasks:

1) Computation of the interface of fluid
2) Computation of the flow behind the interface

some important issues are:
A) Domain change every time-step
B) Flow is compressible and termodynamical effects are to be considered

I think the most difficult problem is A. Any hints about which solver could be the starting point?

Thanks in advance

Francesco

Hi Bernhard, as I said in

Hi Bernhard,

as I said in a different thread. after many tries I think that interFoam could be a good start point, the problem of the outlet could be fixed with your suggestion. I'm trying now to realize this kind of "gamma-switched" BC.

Thank you for your help, I'll tell you if it works

Francesco

Hi all, I've tried to creat

Hi all,

I've tried to create the new boundary condition, but I've got some issues. I tell what I've done, so maybe someone could tell me which is my error.

I've copied the folder "outletInlet" which is in $FOAM_SRC/finiteVolume/fields/fvPatchFields/derived to $FOAM_SRC/, renamed it "filterFvPatchField" and I've replaced the word "outletInlet" with "filter" inside the files. Than I've typed

cp -r $WM_PROJECT_DIR/src/filter $WM_PROJECT_USER_DIR/applications
cd $WM_PROJECT_USER_DIR/applications/filter

In the Make folder I've changed files and options to:

files:

filterFvPatchField.C

LIB = $(FOAM_USER_LIBBIN)/libfilterFvPatchField

options:

EXE_INC = \
-I$(LIB_SRC)/triSurface/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/finiteVolume/lnInclude

LIB_LIBS = \
-ltriSurface \
-lmeshTools \
-lfiniteVolume \
-lfilterFvPatchField

When I use the command

wmake libso

It gives me errors. I think that there's something wrong in the "files" and/or "options" but I cannot find what. Is it what I'm doing correct?

Thanks in advance

Francesco

Ok, I think I've found the err

Ok, I think I've found the errors:

- in Make/files it should be filterFvPatchFields.C with the final "s".

- in Make/options I've add -L$(FOAM_USER_LIBBIN) to -lfilterFvPatchField

I've typed wmake libso and now it gave me the "up to date" message http://www.cfd-online.com/OpenFOAM_D...part/happy.gif

Hi all, now I'm having tro

Hi all,

now I'm having trouble with the boundary condition. I've copied outletInlet and I've tried to write valueFraction as a function of gamma, but when I compile I get some issues about gamma.

The file filterFvPatchField.C has

template<class>
filterFvPatchField<type>::filterFvPatchField
(
const fvPatch& p,
const DimensionedField<type,>& iF,
const volScalarField& gamma
)
:
mixedFvPatchField<type>(p, iF),
gamma_("gamma")
{
this->refValue() = *this;
this->refGrad() = pTraits<type>::zero;
this->valueFraction() = 0.0;
}

and the updateCoeffs is

template<class>
void filterFvPatchField<type>::updateCoeffs()
{
if (this->updated())
{
return;
}

this->valueFraction() = 1.0 - gamma_();

mixedFvPatchField<type>::updateCoeffs();
}

}

I know that there are some errors, but I'm new to C++ and it is hard to understand this kind of code for me.

Can someone explain me which are the errors and how to fix them?

Thanks in advance

Francesco

I think I've fixed some issues

I think I've fixed some issues but still now I got some errors

on line 131 I have:

const volScalarField& gammap = gamma().mesh().lookupObject<volscalarfield>
(
gammaName_,
reinterpret_cast<const>(NULL),
reinterpret_cast<const>(NULL)
);

and on line 138:

this->valueFraction() = 1.0 - gammap();

I got this errors:

In file included from filterFvPatchField.H:142,
from filterFvPatchFields.H:30,
from filterFvPatchFields.C:29:
/usr/include/bits/mathcalls.h: In member function 'virtual void Foam::filterFvPatchField<type>::updateCoeffs()':
/usr/include/bits/mathcalls.h:265: error: too few arguments to function 'double gamma(double)'
filterFvPatchField.C:131: error: at this point in file
filterFvPatchField.C:131: error: expected primary-expression before '>' token
filterFvPatchField.C:138: error: no match for call to '(const Foam::GeometricField<double,>) ()'

Can someone explain me how to fix them? Pls

Thanks in advance

Francesco

Can you attach your code?

Can you attach your code?

selective outlet BC

Hi Foamers,

I'm trying to redo the boundary condition that Francesco created.
I copied and renamed the inletOutletFvPatchField boundary conditions to vapMembraneFvPatchField and tried costumize them.
My .C is:

Code:

#include "vapMembraneFvPatchField.H"



// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //



template<class Type>

Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

( 

    const fvPatch& p,

    const DimensionedField<Type, volMesh>& iF,

    const volScalarField& alphaMatrix

)

:

    mixedFvPatchField<Type>(p, iF),

  // original: phiName_("phi")



    alphaMatrixName_("alphaMatrix")



{

    this->refValue() = pTraits<Type>::zero;

    this->refGrad() = pTraits<Type>::zero;

    this->valueFraction() = 0.0;

}



template<class Type>

Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

(

    const vapMembraneFvPatchField<Type>& ptf,

    const fvPatch& p,

    const DimensionedField<Type, volMesh>& iF,

    const fvPatchFieldMapper& mapper

)

:

    mixedFvPatchField<Type>(ptf, p, iF, mapper),

    phiName_(ptf.phiName_)

{}





template<class Type>

Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

(

    const fvPatch& p,

    const DimensionedField<Type, volMesh>& iF,

    const dictionary& dict

)

:

    mixedFvPatchField<Type>(p, iF),

    phiName_(dict.lookupOrDefault<word>("phi", "phi"))

{

    this->refValue() = Field<Type>("inletValue", dict, p.size());



    if (dict.found("value"))

    {

        fvPatchField<Type>::operator=

        (

            Field<Type>("value", dict, p.size())

        );

    }

    else

    {

        fvPatchField<Type>::operator=(this->refValue());

    }



    this->refGrad() = pTraits<Type>::zero;

    this->valueFraction() = 0.0;

}





template<class Type>

Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

(

    const vapMembraneFvPatchField<Type>& ptf

)

:

    mixedFvPatchField<Type>(ptf),

    phiName_(ptf.phiName_)

{}





template<class Type>

Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

(

    const vapMembraneFvPatchField<Type>& ptf,

    const DimensionedField<Type, volMesh>& iF

)

:

    mixedFvPatchField<Type>(ptf, iF),

    phiName_(ptf.phiName_)

{}





// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //



template<class Type>

void Foam::vapMembraneFvPatchField<Type>::updateCoeffs()

{

    if (this->updated())

    {

        return;

    }





    const Field<scalar>& alphaMatrix =

        this->patch().template lookupPatchField<surfaceScalarField, scalar>

        (

            alphaMatrixName_

        );



    this->valueFraction() = 1.0 - alphaMatrix();





    mixedFvPatchField<Type>::updateCoeffs();

}





template<class Type>

void Foam::vapMembraneFvPatchField<Type>::write(Ostream& os) const

{

    fvPatchField<Type>::write(os);

    if (phiName_ != "phi")

    {

        os.writeKeyword("phi") << phiName_ << token::END_STATEMENT << nl;

    }

    this->refValue().writeEntry("inletValue", os);

    this->writeEntry("value", os);

}





// * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //



template<class Type>

void Foam::vapMembraneFvPatchField<Type>::operator=

(

    const fvPatchField<Type>& ptf

)

{

    fvPatchField<Type>::operator=

    (

        this->valueFraction()*this->refValue()

        + (1 - this->valueFraction())*ptf

    );

}





// ************************************************************************* //

And my .H file:

Code:



#ifndef vapMembraneFvPatchField_H

#define vapMembraneFvPatchField_H



#include "mixedFvPatchField.H"



// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



namespace Foam

{



/*---------------------------------------------------------------------------*\

                   Class vapMembraneFvPatchField Declaration

\*---------------------------------------------------------------------------*/



template<class Type>

class vapMembraneFvPatchField

:

    public mixedFvPatchField<Type>

{



protected:



    // Protected data



        //- Name of flux field

        word phiName_;

            

            //- Name of the volume fraction field

            word alphaMatrixName_;

/*

        // by me, to avoid " class doesn not have any diled name "gamma_"

        word gamma_;

*/



public:



    //- Runtime type information

    TypeName("vapMembrane");





    // Constructors



        //- Construct from patch and internal field

        vapMembraneFvPatchField

        (

            const fvPatch&,

            const DimensionedField<Type, volMesh>&,

            const volScalarField& alphaMatrix



        );



        //- Construct from patch, internal field and dictionary

        vapMembraneFvPatchField

        (

            const fvPatch&,

            const DimensionedField<Type, volMesh>&,

            const dictionary&

        );



        //- Construct by mapping given vapMembraneFvPatchField onto a new patch

        vapMembraneFvPatchField

        (

            const vapMembraneFvPatchField<Type>&,

            const fvPatch&,

            const DimensionedField<Type, volMesh>&,

            const fvPatchFieldMapper&

        );



        //- Construct as copy

        vapMembraneFvPatchField

        (

            const vapMembraneFvPatchField<Type>&

        );



        //- Construct and return a clone

        virtual tmp<fvPatchField<Type> > clone() const

        {

            return tmp<fvPatchField<Type> >

            (

                new vapMembraneFvPatchField<Type>(*this)

            );

        }



        //- Construct as copy setting internal field reference

        vapMembraneFvPatchField

        (

            const vapMembraneFvPatchField<Type>&,

            const DimensionedField<Type, volMesh>&

        );



        //- Construct and return a clone setting internal field reference

        virtual tmp<fvPatchField<Type> > clone

        (

            const DimensionedField<Type, volMesh>& iF

        ) const

        {

            return tmp<fvPatchField<Type> >

            (

                new vapMembraneFvPatchField<Type>(*this, iF)

            );

        }





    // Member functions



        //- Update the coefficients associated with the patch field

        virtual void updateCoeffs();



        //- Write

        virtual void write(Ostream&) const;





    // Member operators



        virtual void operator=(const fvPatchField<Type>& pvf);

};





// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



} // End namespace Foam



// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



#ifdef NoRepository

#   include "vapMembraneFvPatchField.C"

#endif



// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



#endif



// ************************************************************************* //

When I try to do wclean; wmake I get the following error:

Code:

peter@peter-HP-Compaq-nw8440-RZ646EC-ABD:~/OpenFOAM/peter-2.3.0/mySolvers/vapMembrane$ wclean; wmake

wmakeLnInclude: linking include files to ./lnInclude

Making dependency list for source file vapMembrane/vapMembraneFvPatchFields.C

SOURCE=vapMembrane/vapMembraneFvPatchFields.C ;  g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam230/src/triSurface/lnInclude -I/opt/openfoam230/src/meshTools/lnInclude -I/opt/openfoam230/src/transportModels/twoPhaseMixture/lnInclude -I/opt/openfoam230/src/transportModels/incompressible/lnInclude -I/opt/openfoam230/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam230/src/finiteVolume/lnInclude -I/opt/openfoam230/src/meshTools/lnInclude -I/opt/openfoam230/src/fvOptions/lnInclude -I/opt/openfoam230/src/sampling/lnInclude -IlnInclude -I. -I/opt/openfoam230/src/OpenFOAM/lnInclude -I/opt/openfoam230/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linuxGccDPOpt/vapMembraneFvPatchFields.o

In file included from vapMembrane/vapMembraneFvPatchField.H:195:0,

                 from vapMembrane/vapMembraneFvPatchFields.H:29,

                 from vapMembrane/vapMembraneFvPatchFields.C:28:

vapMembrane/vapMembraneFvPatchField.C: In member function ‘virtual void Foam::vapMembraneFvPatchField<Type>::updateCoeffs()’:

vapMembrane/vapMembraneFvPatchField.C:133:47: error: no match for call to ‘(const Foam::Field<double>) ()’

     this->valueFraction() = 1.0 - alphaMatrix();

                                               ^

In file included from /opt/openfoam230/src/OpenFOAM/lnInclude/token.H:49:0,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/UListIO.C:28,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/UList.C:239,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/UList.H:473,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/List.H:43,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/fileNameList.H:42,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/OSspecific.H:39,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/regIOobject.H:43,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/GeometricField.H:42,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/GeometricScalarField.H:38,

                 from /opt/openfoam230/src/OpenFOAM/lnInclude/GeometricFields.H:34,

                 from /opt/openfoam230/src/finiteVolume/lnInclude/volFields.H:37,

                 from vapMembrane/vapMembraneFvPatchFields.C:26:

/opt/openfoam230/src/finiteVolume/lnInclude/fvPatchField.H: In instantiation of ‘static Foam::tmp<Foam::fvPatchField<Type> > Foam::fvPatchField<Type>::addpatchConstructorToTable<fvPatchFieldType>::New(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&) [with fvPatchFieldType = Foam::vapMembraneFvPatchField<double>; Type = double]’:

/opt/openfoam230/src/finiteVolume/lnInclude/fvPatchField.H:120:9:   required from ‘Foam::fvPatchField<Type>::addpatchConstructorToTable<fvPatchFieldType>::addpatchConstructorToTable(const Foam::word&) [with fvPatchFieldType = Foam::vapMembraneFvPatchField<double>; Type = double]’

vapMembrane/vapMembraneFvPatchFields.C:38:1:   required from here

/opt/openfoam230/src/OpenFOAM/lnInclude/runTimeSelectionTables.H:76:66: error: no matching function for call to ‘Foam::vapMembraneFvPatchField<double>::vapMembraneFvPatchField(const Foam::fvPatch&, const Foam::DimensionedField<double, Foam::volMesh>&)’

             return autoPtr< baseType >(new baseType##Type parList);           \

                                                                  ^

/opt/openfoam230/src/finiteVolume/lnInclude/fvPatchField.H:120:9: note: in expansion of macro ‘declareRunTimeSelectionTable’

         declareRunTimeSelectionTable

         ^

/opt/openfoam230/src/OpenFOAM/lnInclude/runTimeSelectionTables.H:76:66: note: candidates are:

             return autoPtr< baseType >(new baseType##Type parList);           \

                                                                  ^

/opt/openfoam230/src/finiteVolume/lnInclude/fvPatchField.H:120:9: note: in expansion of macro ‘declareRunTimeSelectionTable’

         declareRunTimeSelectionTable

         ^

In file included from vapMembrane/vapMembraneFvPatchField.H:195:0,

                 from vapMembrane/vapMembraneFvPatchFields.H:29,

                 from vapMembrane/vapMembraneFvPatchFields.C:28:

vapMembrane/vapMembraneFvPatchField.C:105:1: note: Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField(const Foam::vapMembraneFvPatchField<Type>&, const Foam::DimensionedField<Type, Foam::volMesh>&) [with Type = double]

 Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

 ^

vapMembrane/vapMembraneFvPatchField.C:105:1: note:   no known conversion for argument 1 from ‘const Foam::fvPatch’ to ‘const Foam::vapMembraneFvPatchField<double>&’

vapMembrane/vapMembraneFvPatchField.C:94:1: note: Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField(const Foam::vapMembraneFvPatchField<Type>&) [with Type = double]

 Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

 ^

vapMembrane/vapMembraneFvPatchField.C:94:1: note:   candidate expects 1 argument, 2 provided

vapMembrane/vapMembraneFvPatchField.C:50:1: note: Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField(const Foam::vapMembraneFvPatchField<Type>&, const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&, const Foam::fvPatchFieldMapper&) [with Type = double]

 Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

 ^

vapMembrane/vapMembraneFvPatchField.C:50:1: note:   candidate expects 4 arguments, 2 provided

vapMembrane/vapMembraneFvPatchField.C:64:1: note: Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&, const Foam::dictionary&) [with Type = double]

 Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

 ^

vapMembrane/vapMembraneFvPatchField.C:64:1: note:   candidate expects 3 arguments, 2 provided

vapMembrane/vapMembraneFvPatchField.C:31:1: note: Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&, const volScalarField&) [with Type = double; Foam::volScalarField = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>]

 Foam::vapMembraneFvPatchField<Type>::vapMembraneFvPatchField

 ^

vapMembrane/vapMembraneFvPatchField.C:31:1: note:   candidate expects 3 arguments, 2 provided

I just pasted the first of the error message.
It seems that alphaMatrix, which should contain the volume fraction of liquid polymer, is not accessible. I'm running openfoam 2.3.0. Please find in attachment my alpha.matrix file which defines the volume fractions of liquid at t=0.

Thank you in advance for your help.
Peter