Problem horizontal separator -multiphaseEulerFoam

humber11 · November 12, 2022, 13:35

Hello

I'm trying to simulate a separator (water, oil and air) in the multiphaseEulerFoam, but I face a problem of divergence very fast, I believe it's a problem in the boundary conditions, but I don't find the error, if anyone can help me, I'd appreciate it.

Note: in multiphaseInterFoam the simulation worked.

Code:

/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      phaseProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

type basicMultiphaseSystem;

phases (water oil air);

water
{
    type            purePhaseModel;
    diameterModel   constant;
    constantCoeffs
    {
        d               1e-4;
    }

    residualAlpha 1e-3;
}

oil
{
    type            purePhaseModel;
    diameterModel   constant;
    constantCoeffs
    {
        d               1e-4;
    }

    residualAlpha 1e-3;
}

air
{
    type            purePhaseModel;
    diameterModel   constant;
    constantCoeffs
    {
        d              3e-3;
    }

    residualAlpha 1e-3;
}

blending
{
   default
    {
        type            none;
        minFullyContinuousAlpha.oil 0.7;
        minPartlyContinuousAlpha.oil 0.3;
        minFullyContinuousAlpha.water 0.7;
        minPartlyContinuousAlpha.water 0.3;
        minFullyContinuousAlpha.air 0.7;
        minPartlyContinuousAlpha.air 0.3;
    }
    
}

surfaceTension
(

    (air and water)
    {
        type            constant;
        sigma           0.07;
    }

    (air and oil)
    {
        type            constant;
        sigma           0.032;
    }

    (water and oil)
    {
        type            constant;
        sigma           0.03;
    }
    
);

interfaceCompression
(

);

aspectRatio
(

);

drag
(
/
);

virtualMass
(
/
);

heatTransfer
();

phaseTransfer
();

lift
();

wallLubrication
();

turbulentDispersion
();

// Minimum allowable pressure
pMin            10000;

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 1e5;

boundaryField
{
 inlet
    {
       type            fixedFluxPressure;
    	value           $internalField;
    }
    outletFluid1
    {
        type            totalPressure;
        p0              $internalField;
        U               U.water;
        phi             phi.water;
        rho             thermo:rho.water;
        value           $internalField;
    }
    outletFluid2
    {
        type            totalPressure;
        p0              $internalField;
        U               U.oil;
        phi             phi.oil;
        rho             thermo:rho.oil;
        value           $internalField;
    }
    walls
    {
         type            fixedFluxPressure;
        value           $internalField;
    }
    outlet_air
    {
        type            totalPressure;
        p0              $internalField;
        U               U.air;
        phi             phi.air;
        rho             thermo:rho.air;
        value           $internalField;
   }
}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions          [1 -1 -2 0 0 0 0];

internalField       uniform 1e5;

boundaryField
{	
	inlet
    {
        type            calculated;
        value           $internalField;
    }
    outletFluid1
    {
        type            calculated;
        value           $internalField;
    }
    outletFluid2
    {
        type            calculated;
        value           $internalField;
    }
    walls
    {
        type            calculated;
        value           $internalField;
    }
    outlet_air
    {
        type            calculated;
        value           $internalField;
    }
}

Code:

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

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    inlet
    {
    	type			fixedValue;
    	value           uniform 0.5;
    }
    outletFluid1
    {
        type            inletOutlet;
        phi             phi.water;
        inletValue      uniform 1;
        value           uniform 1;
    }
    outletFluid2
    {
        type            inletOutlet;
        phi             phi.water;
        inletValue      uniform 0;
        value           uniform 0;
    }
    walls
    {
        type            zeroGradient;
    }
        outlet_air
    {
        type            inletOutlet;
        phi             phi.water;
        inletValue      uniform 0;
        value           uniform 0;
    }
}

Code:

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

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
    	type			fixedValue;
    	value           uniform 0.5;
    }
    outletFluid1
    {
        type            inletOutlet;
        phi             phi.oil;
        inletValue      uniform 0;
        value           uniform 0;
    }
    outletFluid2
    {
        type            inletOutlet;
        phi             phi.oil;
        inletValue      uniform 1;
        value           uniform 1;
    }
    walls
    {
        type            zeroGradient;
    }
    outlet_air
    {
        type            inletOutlet;
        phi             phi.oil;
        inletValue      uniform 0;
        value           uniform 0;
    }
}

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
{
	inlet
    {
    	type			fixedValue;
    	value           uniform 0;
    }
    outletFluid1
    {
        type            inletOutlet;
        phi             phi.air;
        inletValue      uniform 0;
        value           uniform 0;
    }
    outletFluid2
    {
        type            inletOutlet;
        phi             phi.air;
        inletValue      uniform 0;
        value           uniform 0;
    }
    walls
    {
        type            zeroGradient;
    }
    outlet_air
    {
        type            inletOutlet;
        phi             phi.air;
        inletValue      uniform 1;
        value           uniform 1;
    }
}

Code:

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

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
	inlet
    {
        type           fixedValue;
        value           $internalField;
    }
    outletFluid1
    {
        type            pressureInletOutletVelocity;
        phi             phi.air;
        value           $internalField;
    }
    outletFluid2
    {
        type            pressureInletOutletVelocity;
        phi             phi.air;
        value           $internalField;
    }
    walls
    {
        type            noSlip;
    }
        outlet_air
    {
       type            pressureInletOutletVelocity;
        phi             phi.air;
        value           $internalField;
    }
}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U.oil;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
inlet
    {
        type            fixedValue;
        value           uniform (4 0 0);
    }
    outletFluid1
    {
        type            pressureInletOutletVelocity;
        phi             phi.oil;
        value           $internalField;
    }
    outletFluid2
    {
        type            pressureInletOutletVelocity;
        phi             phi.oil;
        value           $internalField;
    }
    walls
    {
        type            noSlip;
    }
       outlet_air
    {
       type            pressureInletOutletVelocity;
        phi             phi.oil;
        value           $internalField;
    }
}

Code:

{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
	inlet
    {
        type            fixedValue;
        value           uniform (5 0 0);
    }
    outletFluid1
    {
        type            pressureInletOutletVelocity;
        phi             phi.water;
        value           $internalField;
    }
    outletFluid2
    {
        type            pressureInletOutletVelocity;
        phi             phi.water;
        value           $internalField;
    }
    walls
    {
        type            noSlip;
    }
       outlet_air
    {
       type            pressureInletOutletVelocity;
        phi             phi.water;
        value           $internalField;
    }
}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "alpha.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 4;
    }

    p_rgh
    {
        solver          PBiCGStab;
        preconditioner  DIC;
        tolerance       1e-07;
        relTol          0;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "U.*"
    {
        solver          PBiCGStab;
        spreconditioner  DIC;
        tolerance       1e-5;
        relTol          0;
    }

    "e.*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-8;
        relTol          2;
      //  maxIter         0;
        minIter         0;
    }
    
}

PIMPLE
{
    nOuterCorrectors 1;
    nCorrectors         3;
    nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
    equations
    {
        ".*"            1;
    }
}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
	
    "div\(phi,alpha.*\)"    Gauss vanLeer;
    "div\(phir,alpha.*,alpha.*\)"   Gauss vanLeer;

    "div\(alphaRhoPhi.*,U.*\)"      Gauss limitedLinearV 1;
    "div\(phi.*,U.*\)"              Gauss limitedLinearV 1;

    "div\(alphaRhoPhi.*,(h|e).*\)"  Gauss limitedLinear 1;
    "div\(alphaRhoPhi.*,K.*\)"      Gauss limitedLinear 1;
    "div\(alphaPhi.*,p\)"           Gauss limitedLinear 1;

    "div\(\(\(\(alpha.*\*thermo:rho.*\)*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear uncorrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         uncorrected;
}

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     multiphaseEulerFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         20;

deltaT          0.001;

writeControl    adjustableRunTime;

writeInterval   0.1;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression uncompressed;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

adjustTimeStep  no;

maxCo           0.5;

maxDeltaT       1e-2;

I use BlueCFD, does anyone know the meaning of these errors?

1)ZN4Foam11phaseSystem5solveERKNS_7PtrListINS_14Ge ometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEERK NS1_INS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEEE E
2)ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tm pINS_14GeometricFieldIdT_T0_EEEERKS7_SA_

Code:

Backtrace:
        ZN10StackTraceC1Ev [0x626c1855+0x25]
                 module: C:\PROGRA~1\BLUECF~2\ThirdParty-8\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll
        ZN4Foam5error10printStackERNS_7OstreamE [0x6c30ae5a+0x23a]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
        ZN4Foam6sigFpe13sigFpeHandlerEi [0x6c30b693+0x33]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
        (No symbol) [0x409412]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\multiphaseEulerFoam.exe
        _C_specific_handler [0x7ffd3bf17ff8+0x98]
                 module: C:\WINDOWS\System32\msvcrt.dll
        0_chkstk [0x7ffd3c07241f+0x11f]
                 module: C:\WINDOWS\SYSTEM32\ntdll.dll
        RtlRaiseException [0x7ffd3c0214a4+0x434]
                 module: C:\WINDOWS\SYSTEM32\ntdll.dll
        KiUserExceptionDispatcher [0x7ffd3c070f4e+0x2e]
                 module: C:\WINDOWS\SYSTEM32\ntdll.dll
        ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_ [0x6c2bbf56+0x26]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
        ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_SA_ [0x61eff86f+0x25f]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libphaseSystem.dll
        ZN4Foam11phaseSystem5solveERKNS_7PtrListINS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEERKNS1_INS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEEEE [0x61df06d3+0x4b43]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libphaseSystem.dll
        (No symbol) [0x472431]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\multiphaseEulerFoam.exe
        (No symbol) [0x4013c1]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\multiphaseEulerFoam.exe
        (No symbol) [0x4014f6]
                 module: C:\PROGRA~1\BLUECF~2\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\multiphaseEulerFoam.exe
        BaseThreadInitThunk [0x7ffd3bb974b4+0x14]
                 module: C:\WINDOWS\System32\KERNEL32.DLL
        RtlUserThreadStart [0x7ffd3c0226a1+0x21]
                 module: C:\WINDOWS\SYSTEM32\ntdll.dll

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
area does not match neighbour by ... % -- possible face ordering problem	St.Pacholak	OpenFOAM	10	February 7, 2024 21:50
BuoyantBoussinesqSimpleFoam_Facing problem	Mondal131211	OpenFOAM Running, Solving & CFD	1	April 10, 2019 19:41
Mesh& steptime independant: conduction-convection problem	Fati1	Main CFD Forum	1	October 28, 2018 13:52
Adiabatic and Rotating wall (Convection problem)	ParodDav	CFX	5	April 29, 2007 19:13
@ horizontal silencer damping problem	vivek p patil	CFX	0	June 12, 2006 03:31