[Kurnibald]

Hello everyone,


I'm trying to set up a droplet impact case as an 2D axis symmetric case.
The basic case is running and working as long as i use interFoam with incompressibleVoF solver.


However as soon as I try to implement the alphaContactAngle boundary condition I get either the error stating that only contactAngle is allowed:
(using InterFoam incompressibleVoF solver)

--> FOAM FATAL IO ERROR: 
Unknown patchField type alphaContactAngle for patch substrate of type wall for field alpha.water

Valid patchField types are :
...
       )

Or if I run the case with foamrun using incompressibleMultiphaseVoF solver I get this Backtrace.

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

Create mesh for time = 0

Selecting solver incompressibleMultiphaseVoF
Selecting viscosity model constant
Selecting viscosity model constant
Generating stack trace...


Backtrace:
    ZN10StackTraceC1Ev [0x7ffda24a1805+0x25]
         module: E:\blueCFD-Core-2024\ThirdParty-12\platforms\mingw_w64Gcc122DPInt32\lib\libstack_trace.dll
    ZN4Foam5error10printStackERNS_7OstreamE [0x7ffd7bc1c0dc+0x23c]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libOpenFOAM.dll
    ZN4Foam6sigFpe13sigFpeHandlerEi [0x7ffd7bc1c774+0x34]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libOpenFOAM.dll
    (No symbol) [0x7ff6fcf92982]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\bin\foamRun.exe
    _C_specific_handler [0x7ffda6cd7ff8+0x98]
         module: C:\WINDOWS\System32\msvcrt.dll
    0_chkstk [0x7ffda7e928bf+0x11f]
         module: C:\WINDOWS\SYSTEM32\ntdll.dll
    RtlRaiseException [0x7ffda7e42554+0x484]
         module: C:\WINDOWS\SYSTEM32\ntdll.dll
    KiUserExceptionDispatcher [0x7ffda7e913ce+0x2e]
         module: C:\WINDOWS\SYSTEM32\ntdll.dll
    ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_ [0x7ffd7bbb88ee+0x1e]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libOpenFOAM.dll
    ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_RKS7_ [0x7ffd887aebe0+0x3f0]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libincompressibleMultiphaseVoF.dll
    ZN4Foam34incompressibleMultiphaseVoFMixture7correctEv [0x7ffd887472b9+0x209]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libincompressibleMultiphaseVoF.dll
    ZN4Foam34incompressibleMultiphaseVoFMixtureC1ERKNS_6fvMeshE [0x7ffd887478f5+0x365]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libincompressibleMultiphaseVoF.dll
    ZN4Foam7solvers27incompressibleMultiphaseVoFC1ERNS_6fvMeshE [0x7ffd887641b9+0x39]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libincompressibleMultiphaseVoF.dll
     ZN4Foam6solver27addfvMeshConstructorToTableINS_7solvers27incompressibleMultiphaseVoFEE3NewERNS_6fvMeshE  [0x7ffd887a3465+0x25]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libincompressibleMultiphaseVoF.dll
    ZN4Foam6solver3NewERKNS_4wordERNS_6fvMeshE [0x7ffd7a3ce5cb+0x1ab]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\lib\libfiniteVolume.dll
    (No symbol) [0x7ff6fcf9431e]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\bin\foamRun.exe
    (No symbol) [0x7ff6fcf912ee]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\bin\foamRun.exe
    (No symbol) [0x7ff6fcf91406]
         module: E:\blueCFD-Core-2024\OpenFOAM-12\platforms\mingw_w64Gcc122DPInt32Opt\bin\foamRun.exe
    BaseThreadInitThunk [0x7ffda7317374+0x14]
         module: C:\WINDOWS\System32\KERNEL32.DLL
    RtlUserThreadStart [0x7ffda7e3cc91+0x21]
         module: C:\WINDOWS\SYSTEM32\ntdll.dll

I also build a blockmesh file for a flat droplet impact in 2D planar.
I basically changed the damBreak4phase case to two phases and this is running fine as long as i do not implement the uTheta thetaA and TheatR

--> FOAM FATAL ERROR: 
The alphaContactAngle boundary conditions specified on patch leftWall for alpha.water and alpha.air contain inconsistent properties

    From function void Foam::correctContactAngle(const volScalarField&, const volScalarField&, const GeometricField<Vector<double>, fvPatchField, volMesh>::Boundary&, const dimensionedScalar&, GeometricField<Vector<double>, fvsPatchField, surfaceMesh>::Boundary&)
    in file correctContactAngle/correctContactAngle.C at line 92.

      FOAM exiting

So questions arise:

What does the Backtrace mean in case of using incompressibleMultiphaseVoF solver?

Can I use alphaContactAngle in a 2D axis symmetrical cases? ( I guess so since OF always uses a 3D Mesh and the paper i want to compare to does so as well)

Why can I not use alphaContactAngle in my axis symmetrical case and only get the option to choose contactAngle using InterFoam incompressibleVoF solver?
--> Answered myself: it is only implemetned with multiphase solvers see: https://doc.openfoam.com/2312/tools/...aContactAngle/

Why does openfoam complain about inconsistent properties in the flat case? (I tried a few different values adding up to 90°)

What does the value 1 or 0 do at the end of the alphaContactAnlge do?

    substrate
    {
        type            alphaContactAngle;
        contactAngleProperties
        {
            air             { theta0 90; uTheta 0.1; thetaA 60; thetaR 30;}
        }
        value           uniform 0;
        limit            alpha;
     }

I'm getting a bit frustrated, I've set up the case completely new multiple times. Basically it is a damBreak4phase case adapted to a axis symmetrical mesh with only one block which works as long as i don't implement the alphaContactAngle.


I use blueCFD-Core-2024-1 on Windows 10 Pro
which according to blueCFD is OpenFOAM 12



My 0/alpha files:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  12
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [];

internalField   uniform 0;

boundaryField
{
    atmosphere
    {
        type            inletOutlet;
        inletValue      $internalField;    //or uniform 0;
        value           $internalField; //or uniform 0; 
    }
    
    substrate
    {
        type            alphaContactAngle;
        contactAngleProperties
        {
            air             { theta0 90; uTheta 0.1; thetaA 60; thetaR 30;}
        }
        value           uniform 0;
        limit            alpha;
    }

//    substrate
//    {
//      type            contactAngle;
//      theta0             90;
//        UTheta             1;
//        thetaA             30;
//        thetaR             60;
//        limit            alpha;
//        value            $internalField;
//
//    }
    
//    substrate
//    {
//        type            contactAngle;
//        theta0             30;
//        limit            alpha;
//        value            $internalField;
//    }


//    substrate       
//    {
//        type            zeroGradient;
//    }

    front
    {
        type            wedge;
    }
    
    back
    {
        type            wedge;
    }
    
    axis
    {
        type    empty;
    }



}


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

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  12
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [];

internalField   uniform 0;

boundaryField
{

    atmosphere
    {
        type            inletOutlet;
        inletValue      $internalField;    //or uniform 1; 
        value           $internalField; //or uniform 1; 
    }
    
    substrate
    {
        type            alphaContactAngle;
        contactAngleProperties
        {
            water           { theta0 90; uTheta 0.1; thetaA 60; thetaR 30;}
        }
        value           uniform 0;
        limit            alpha;
    }

//    substrate
//    {
//      type            contactAngle;
//      theta0             90;
//        UTheta             1;
//        thetaA             30;
//        thetaR             60;
//        limit            alpha;
//        value            $internalField;
//    }

//    substrate
//    {
//        type            contactAngle;
//        theta0             value 90;
//        limit            alpha.water;
//        value            1;
//    }

//    substrate       
//    {
//        type            zeroGradient;
//    }

    front
    {
        type            wedge;
    }
    
    back
    {
        type            wedge;
    }
    
    axis
    {
        type    empty;
    }



}


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