[Mahmoud_aboukhedr]

Dear Foamrs,

I started working with the MULES::implicitSolve in (interFoam) few month ago for some case depend restrictions.
The code used in the alphaEqu.H is

        MULES::implicitSolve
        (
            alpha1,
            phi,
            phiAlpha,
            1,
            0
        );

I found that the alpha (0 1) is increasing by time up to the limit of (1.1)..

In fact I had tried everthing to make it bounded to 0 and 1 (changing everthing, (mesh, schemes, BC) nothing works.

I followed the work posted in ;

http://www.cfd-online.com/Forums/openfoam-solving/115137-alpha1-grows-unbounded-interdymfoam.html

As I tried to add the source Su to the MULES in the alphaEqu.H and compile the following;

        volScalarField Su
        (
            IOobject
            (
                "Su",
                runTime.timeName(),
                mesh
            ),
            divU*alpha1
        );

        MULES::implicitSolve
        (
            alpha1,
            phi,
            phiAlpha,
            Su,
            1,
            0
        );

The code blow off and give me this error;

alphaEqn.H: In function ‘int main(int, char**)’:
alphaEqn.H:49:9: error: no matching function for call to ‘implicitSolve(Foam::volScalarField&, Foam::surfaceScalarField&, Foam::surfaceScalarField&, Foam::volScalarField&, int, int)’
alphaEqn.H:49:9: note: candidates are:
In file included from /home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULES.H:87:0,
                 from myInterFoamVOFSmooth.C:47:
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULESTemplates.C:59:6: note: template<class RhoType, class SpType, class SuType> void Foam::MULES::implicitSolve(const RhoType&, Foam::volScalarField&, const surfaceScalarField&, Foam::surfaceScalarField&, const SpType&, const SuType&, Foam::scalar, Foam::scalar)
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULESTemplates.C:59:6: note:   template argument deduction/substitution failed:
In file included from alphaEqnSubCycle.H:19:0,
                 from myInterFoamVOFSmooth.C:100:
alphaEqn.H:49:9: note:   cannot convert ‘phi’ (type ‘Foam::surfaceScalarField {aka Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>}’) to type ‘Foam::volScalarField& {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&}’
In file included from myInterFoamVOFSmooth.C:47:0:
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULES.H:70:6: note: void Foam::MULES::implicitSolve(Foam::volScalarField&, const surfaceScalarField&, Foam::surfaceScalarField&, Foam::scalar, Foam::scalar)
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULES.H:70:6: note:   candidate expects 5 arguments, 6 provided
In file included from alphaEqnSubCycle.H:27:0,
                 from myInterFoamVOFSmooth.C:100:
alphaEqn.H:49:9: error: no matching function for call to ‘implicitSolve(Foam::volScalarField&, Foam::surfaceScalarField&, Foam::surfaceScalarField&, Foam::volScalarField&, int, int)’
alphaEqn.H:49:9: note: candidates are:
In file included from /home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULES.H:87:0,
                 from myInterFoamVOFSmooth.C:47:
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULESTemplates.C:59:6: note: template<class RhoType, class SpType, class SuType> void Foam::MULES::implicitSolve(const RhoType&, Foam::volScalarField&, const surfaceScalarField&, Foam::surfaceScalarField&, const SpType&, const SuType&, Foam::scalar, Foam::scalar)
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULESTemplates.C:59:6: note:   template argument deduction/substitution failed:
In file included from alphaEqnSubCycle.H:27:0,
                 from myInterFoamVOFSmooth.C:100:
alphaEqn.H:49:9: note:   cannot convert ‘phi’ (type ‘Foam::surfaceScalarField {aka Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>}’) to type ‘Foam::volScalarField& {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&}’
In file included from myInterFoamVOFSmooth.C:47:0:
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULES.H:70:6: note: void Foam::MULES::implicitSolve(Foam::volScalarField&, const surfaceScalarField&, Foam::surfaceScalarField&, Foam::scalar, Foam::scalar)
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/IMULES.H:70:6: note:   candidate expects 5 arguments, 6 provided
In file included from myInterFoamVOFSmooth.C:68:0:
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable]
make: *** [Make/linux64GccDPOpt/myInterFoamVOFSmooth.o] Error 1

I will appreciate the help,

Mahmoud

[kmou]

Hey Mahmoud,
Did you make this work?
Camille

[rapierrz]

Hi mahmoud,

Do you try to change your maximum unboundness in fvSolution/alpha1?

it's work for me in interPhaseChangeFoam.

[kmou]

Originally Posted by rapierrz

Hi mahmoud,

Do you try to change your maximum unboundness in fvSolution/alpha1?

it's work for me in interPhaseChangeFoam.

Hi rapierrz,

In fvSolution/alpha1,
what is the flag you would use to set the maximum unboundedness?
I am not seeing any/I'm mising out on something.
Thanks
Camille

[Noémie]

Hi Mahmoud,

Do you solve the problem?
I'm also facing this problem with unbounded alpha1.

Hao