multiphaseInterfoam non-constant inlet

kaaja · February 20, 2018, 03:04

Edit2: Problem solved. See "edit" in the end.

Hi,

I run multiphaseInterfoam, and I have trouble with the inlet being non-constant (I want it to be constant.)

Here are my alpha-files

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.air;
}
// * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            alphaContactAngle;
        thetaProperties
        (
            ( freshWater air ) 90 0 0 0
            ( saltWater air ) 90 0 0 0
            ( freshWater saltWater ) 90 0 0 0
        );
        value           uniform 0;
    }

    outlet
    
    {
        type            alphaContactAngle;
        thetaProperties
        (
            ( freshWater air ) 90 0 0 0
            ( saltWater air ) 90 0 0 0
            ( freshWater saltWater ) 90 0 0 0
        );
        value           uniform 0;
    
    }

        atmosphere
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }
    

    barge
    {
        type            alphaContactAngle;
        thetaProperties
        (
            ( freshWater air ) 90 0 0 0
            ( saltWater air ) 90 0 0 0
            ( freshWater saltWater ) 90 0 0 0
        );
        value           uniform 0;
    }



}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.freshWater;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            variableHeightFlowRate;
        lowerBound      0;
        upperBound      1;
        value           $internalField;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    barge
    {
        type            zeroGradient;
    }
}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.saltWater;
}


dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            variableHeightFlowRate;
        lowerBound      0;
        upperBound      1;
        value           $internalField;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    barge
    {
        type            zeroGradient;
    }
}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alphas;
}
// * * * * * * * * * ** * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            variableHeightFlowRate;
        lowerBound      0;
        upperBound      1;
        value           $internalField;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    barge
    {
        type            zeroGradient;
    }
}

The above gives the following, wanted, distribtion of fluids for timestep 1

initialAlpha.png

However, after several time the above changes, also at the inlet.

I really don't understand the contactAngle functino used in alpha.air above. I have tried with the following alpha.air

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.air;
}
// * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            variableHeightFlowRate;
        lowerBound      0;
        upperBound      1;
        value           $internalField;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    barge
    {
        type            zeroGradient;
    }
}

With the alpha.air above I get the following error message

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 5.x-197d9d3bf20a
Exec   : multiphaseInterFoam -parallel
Date   : Feb 20 2018
Time   : 08:53:57
Host   : "..."
PID    : 27944
I/O    : uncollated
Case   :...
nProcs : 8
Slaves : 
7
(
"....27945"
"....27946"
"....27947"
"....27948"
"....27949"
"....27950"
"....27951"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
[1] #0[0] #0    Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[7] [6] #0  Foam::error::printStack(Foam::Ostream&)[3] #0  Foam::error::printStack(Foam::Ostream&)#0  Foam::error::printStack(Foam::Ostream&) at ??:?
[1] #1  Foam::sigFpe::sigHandler(int) at ??:?
[6] #1  Foam::sigFpe::sigHandler(int) at ??:?
[7] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[3] #1  Foam::sigFpe::sigHandler(int) at ??:?
[1] #2  ? at ??:?
[6] #2  ? at ??:?
[7] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[0] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[3] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[6] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[1] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[0] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[7] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[1] #5  Foam::multiphaseMixture::multiphaseMixture(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
[6] #5  Foam::multiphaseMixture::multiphaseMixture(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
[0] #5  Foam::multiphaseMixture::multiphaseMixture(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
[7] #5  Foam::multiphaseMixture::multiphaseMixture(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
[3] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[1] #6   at ??:?
[7] #6   at ??:?
[6] #6   at ??:?
[0] #6  ???? at ??:?
[3] #5  Foam::multiphaseMixture::multiphaseMixture(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[1] #7  __libc_start_main in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[6] #7  __libc_start_main in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[0] #7  __libc_start_main at ??:?
[3] #6   in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[7] #7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #8   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #8   in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #8  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #8  ?? in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[3] #7  __libc_start_main??[karl-HP-OMEN-Pro-15-Notebook:27945] *** Process received signal ***
 in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[...27945] Signal: Floating point exception (8)
[...:27945] Signal code:  (-6)
[...:27945] Failing at address: 0x3e800006d29
[...:27945] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7f5db443d140]
[...:27945] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f5db443d0bb]
[...:27945] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7f5db443d140]
[...:27945] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0x155 in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[...:27944] *** Process received signal ***
[...:27944] Signal: Floating point exception (8)
[...:27944] Signal code:  (-6)
[...:27944] Failing at address: 0x3e800006d28
[...:27944] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7f00cb829140]
[...:27944] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f00cb8290bb]
[...:27944] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7f00cb829140]
[...:27944] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0x)[0x7f5db598ae55]
[...:27945] [ 4] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_SA_+0x20e)[0x7f5db9e9e3de]
[...:27945] [ 5] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4Foam17multiphaseMixtureC2ERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEE+0x51d)[0x7f5db9e7172d]
[...:27945] [ 6] multiphaseInterFoam(+0x34329)[0x5556335d5329]
[...:27945] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf1)[0x7f5db44271c1]
[...:27945] [ 8] multiphaseInterFoam(+0x3a26a)[0x5556335db26a]
[...:27945] *** End of error message ***
155)[0x7f00ccd76e55]
[...:27944] [ 4] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_SA_+0x20e)[0x7f00d128a3de]
[...:27944] [ 5] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4Foam17multiphaseMixtureC2ERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEE+0x51d)[0x7f00d125d72d]
[...:27944] [ 6] multiphaseInterFoam(+0x34329)[0x558038293329]
[...:27944] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf1)[0x7f00cb8131c1]
[...:27944] [ 8] multiphaseInterFoam(+0x3a26a)[0x55803829926a]
[...:27944] *** End of error message ***
 in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
 in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[...:27950] *** Process received signal ***
[...:27950] Signal: Floating point exception (8)
[...:27950] Signal code:  (-6)
[...:27950] Failing at address: 0x3e800006d2e
[...:27950] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7feb936ac140]
[...:27950] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7feb936ac0bb]
[...:27950] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7feb936ac140]
[...:27950] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0x155)[0x7feb94bf9e55]
[...:27950] [ 4] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_SA_+0x20e)[0x7feb9910d3de]
[...:27950] [ 5] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4Foam17multiphaseMixtureC2ERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEE+0x51d)[0x7feb990e072d]
[...:27950] [ 6] multiphaseInterFoam(+0x34329)[0x562c0d3ad329]
[...:27950] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf1)[0x7feb936961c1]
[...:27950] [ 8] multiphaseInterFoam(+0x3a26a)[0x562c0d3b326a]
[...:27950] *** End of error message ***
[...:27951] *** Process received signal ***
[...27951] Signal: Floating point exception (8)
[...27951] Signal code:  (-6)
[...27951] Failing at address: 0x3e800006d2f
[...27951] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7f88852f6140]
[...27951] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f88852f60bb]
[...27951] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7f88852f6140]
[...27951] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0x155)[0x7f8886843e55]
[...27951] [ 4] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_SA_+0x20e)[0x7f888ad573de]
[...27951] [ 5] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4Foam17multiphaseMixtureC2ERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEE+0x51d)[0x7f888ad2a72d]
[...27951] [ 6] multiphaseInterFoam(+0x34329)[0x5566b2bd7329]
[...27951] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf1)[0x7f88852e01c1]
[...27951] [ 8] multiphaseInterFoam(+0x3a26a)[0x5566b2bdd26a]
[...27951] *** End of error message ***
 in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #8  ?[...27947] *** Process received signal ***
 in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
[...27947] Signal: Floating point exception (8)
[...27947] Signal code:  (-6)
[...27947] Failing at address: 0x3e800006d2b
[...27947] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7fa7776e8140]
[...27947] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7fa7776e80bb]
[...27947] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x37140)[0x7fa7776e8140]
[...27947] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0x155)[0x7fa778c35e55]
[...27947] [ 4] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_SA_+0x20e)[0x7fa77d1493de]
[...27947] [ 5] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libmultiphaseInterFoam.so(_ZN4Foam17multiphaseMixtureC2ERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEE+0x51d)[0x7fa77d11c72d]
[...27947] [ 6] multiphaseInterFoam(+0x34329)[0x55bd837b3329]
[...27947] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf1)[0x7fa7776d21c1]
[...27947] [ 8] multiphaseInterFoam(+0x3a26a)[0x55bd837b926a]
[...27947] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 6 with PID 0 on node ... exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

Does anybvody know how to modify the above files to make the distribution of alpha's at the inlet constant?

edit:
In the last version of alpa.air I changed from "internalField uniform 0" to "internalField uniform 1", as shown below, and not the inlet is constant. I think the internalFIeld value has to equal the default value set in setFieldsDict.

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    //- Set patchGroups for constraint patches
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            variableHeightFlowRate;
        lowerBound      0;
        upperBound      1;
        value           $internalField;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }

    barge
    {
        type            zeroGradient;
    }
}

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

wangchengan2003 · February 21, 2018, 20:13

Dear kaaja,

Have you solved this problem? I am working on the same problem like you. I want to set the three phases moving from left to right. I use groovybc, but it doesn't work.

Best regards,

Chengan

kaaja · February 22, 2018, 01:50

Chengan,

The Inlet became stable (constant alpha shares) when I did what I wrote in my "edit". So I did not use the contactAngle stuff for alpha.air, but used a alpha.air file that looked very similar to the other alpha-files. Also I made sure that the internalField was 1, not 0. That made it work. Hope that helps!

wangchengan2003 · February 22, 2018, 20:24

Dear kaaja,

Thank you for your reply. I still have a question. If we want to set the three phases moving from left to right, we still use the following code for all the alpha-files?

Quote:

inlet
{
type fixedValue;
value $internalField;
}

Best regards

Chengan

kaaja · February 23, 2018, 02:04

What I ended up doing is as you say, using the following on all inlets

Code:

inlet
    {
        type            fixedValue;
        value           $internalField;
    }

However, the internalField value is not the same in the alpha-files. I use the value "1" for the alpha that is set as the default phase in "setFieldsDict", and the value "0" for internalFIeld in the other alpha-files.

I do not understand this fully myself. It is maybe not necessary to choose internalFIeld "1" for the alpha that is set as default in "setFieldsDict". Maybe it is possible to have internalFIeld "1" for another alpa than the alpha that is set as default in "setFieldsDict". However, I think only one of the alpha-files can have internalfield equal to "1" at the inlet.

As I said, I do not understand this fully myself, so please let me know if you find out something more!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
multiphaseInterFoam vs multiphaseEulerFoam	rahulksoni	OpenFOAM	0	August 11, 2017 02:14
[snappyHexMesh] Normals abnormal? My snappyHexMesh mesh is crashing extrudeMesh :(	KarenRei	OpenFOAM Meshing & Mesh Conversion	9	October 5, 2016 22:50
How to keep the water level constant at inlet	Tanjina	FLUENT	10	January 21, 2016 10:29
How to set up the inlet boundary condition for a low pressure case?	beastieboys6	FLUENT	3	April 10, 2012 22:46
Inlet table in STAR-CD	Sachin	Siemens	1	March 26, 2008 10:22

February 21, 2018, 20:13		#2
wangchengan2003 Member chengan.wang Join Date: Jan 2016 Location: china Posts: 47 Rep Power: 10	Dear kaaja, Have you solved this problem? I am working on the same problem like you. I want to set the three phases moving from left to right. I use groovybc, but it doesn't work. Best regards, Chengan

February 22, 2018, 01:50		#3
kaaja Member Join Date: Apr 2017 Posts: 68 Rep Power: 9	Chengan, The Inlet became stable (constant alpha shares) when I did what I wrote in my "edit". So I did not use the contactAngle stuff for alpha.air, but used a alpha.air file that looked very similar to the other alpha-files. Also I made sure that the internalField was 1, not 0. That made it work. Hope that helps!

February 23, 2018, 02:04		#5
kaaja Member Join Date: Apr 2017 Posts: 68 Rep Power: 9	What I ended up doing is as you say, using the following on all inlets Code: inlet { type fixedValue; value $internalField; } However, the internalField value is not the same in the alpha-files. I use the value "1" for the alpha that is set as the default phase in "setFieldsDict", and the value "0" for internalFIeld in the other alpha-files. I do not understand this fully myself. It is maybe not necessary to choose internalFIeld "1" for the alpha that is set as default in "setFieldsDict". Maybe it is possible to have internalFIeld "1" for another alpa than the alpha that is set as default in "setFieldsDict". However, I think only one of the alpha-files can have internalfield equal to "1" at the inlet. As I said, I do not understand this fully myself, so please let me know if you find out something more!