CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Weird physical observation in InterFOAM

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 23, 2024, 21:59
Default Weird physical observation in InterFOAM
  #1
Senior Member
 
mohammad
Join Date: Sep 2015
Posts: 281
Rep Power: 12
mostanad is on a distinguished road
Hi foamers,

I have recently solved a simple case with three phases in openfoam and have seen a physical observation which cannot interpret.
The case is a vessel with open boundary condition at the top boundary.

The boundary conditions I set are as follows;


alpha.air:


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
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
{
    wall
    {
        type            alphaContactAngle;
        thetaProperties
        (
            ( slag air ) 90 0 0 0
            ( slag molten ) 90 0 0 0
            ( molten air ) 90 0 0 0

        );
        value           uniform 0;
    }
    inlet
    {
        type            alphaContactAngle;
        thetaProperties
        (
            ( slag air ) 90 0 0 0
            ( slag molten ) 90 0 0 0
            ( molten air ) 90 0 0 0
    
        );
        value           uniform 0;
    }
    outlet
    {
        type        fixedValue;
    value        uniform 1;
    } 
}
alpha.molten


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.molten;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    wall
    {
        type            zeroGradient;
    } 
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }
}
alpha.slag


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha.slag;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    wall
    {
        type            zeroGradient;
    }
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }
}
alphas


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alphas;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    wall
    {
        type            zeroGradient;
    }
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
    }

}
p_rgh


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{

    wall
    {
        type        fixedFluxPressure;
        value       uniform 0;
    }
    inlet
    {
        type        fixedFluxPressure;
        value       uniform 0;
    }
    outlet
    {
        type        fixedValue;
        value        uniform 0;
    }
}
U


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    outlet
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }
    wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}
transportProperties:


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

phases
(
     molten
     {
         transportModel Newtonian;
         nu nu   [ 0  2 -1  0 0 0 0 ] 1e-06;
         rho rho [ 1 -3  0  0 0 0 0 ] 7000;
                 Cp  Cp  [ 0  2 -2 -1 0 0 0 ] 440;
                 kf  kf  [ 1  1 -3 -1 0 0 0 ] 73; // sped up! 16.5*20;
                 D   D   [ 0  2 -1  0 0 0 0 ] 0.17e-6; // sped up! 0.85e-8*20;
                 Cs  Cs  [ 1 -3  0  0 0 0 0 ] 335;
     }
     slag
     {
         transportModel Newtonian;
         nu nu   [ 0  2 -1  0 0 0 0 ] 2.8e-05;  
         rho rho [ 1 -3  0  0 0 0 0 ] 2000; 
         Cp  Cp  [ 0  2 -2 -1 0 0 0 ] 440;     
         kf  kf  [ 1  1 -3 -1 0 0 0 ] 40; // sped up! 16.5*20; 
         D   D   [ 0  2 -1  0 0 0 0 ] 0.17e-6; // sped up! 0.85e-8*20;
         Cs  Cs  [ 1 -3  0  0 0 0 0 ] 335;
     }
     air
     {
         transportModel Newtonian;
         nu  nu  [ 0 2  -1  0 0 0 0 ] 1e-05;
         rho rho [ 1 -3  0  0 0 0 0 ] 1.25; 
         Cp  Cp  [ 0  2 -2 -1 0 0 0 ] 1000;
         kf  kf  [ 1  1 -3 -1 0 0 0 ] 0.025; // sped up! 0.1*20;
         D   D   [ 0  2 -1  0 0 0 0 ] 0.17e-6; // sped up! 0.85e-8*20; // dummy
         Cs  Cs  [ 1 -3  0  0 0 0 0 ] 0.00025;
     }
);

sigmas
(
    (slag molten) 0.07
    (air slag) 0.07
    (air molten) 0.07

);

Having these boundary conditions, once I set the fields for the phases, air, slag, and molten with suitable thickness like:


Attched Image 1 (below)

and the simulation goes well. But when I set the fields thinner like:

Attched Image 2 (below)

The phases on the bottom (slag in image) becomes disturbed;

Attached image 3 (below)


My question is why this only happens with thinner layers of bottom phases? Which value in properties or BC is making this thin layers case disturbed?

I hope someone can help me to generate more physical result.
Attached Images
File Type: png 1.png (8.3 KB, 7 views)
File Type: png 2.png (7.8 KB, 7 views)
File Type: png 3.png (16.0 KB, 9 views)
mostanad is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Gmsh does not ouput SU2 markers from physical surfaces stewels SU2 4 May 17, 2023 06:15
[Gmsh] Nasa Nozzle: Plot 3D -->gmsh -->fluent - no physical groups Nico89 OpenFOAM Meshing & Mesh Conversion 1 June 14, 2016 06:40
Mesh refinement with interFoam = solution that does no physical sense lyna OpenFOAM Running, Solving & CFD 2 January 19, 2016 12:05
Partition does not start on physical sector boundary!!?? immortality Lounge 1 December 8, 2013 16:04
Problem of InterFoam with LES SpalartAllmarasIDDES keepfit OpenFOAM 3 August 29, 2013 12:21


All times are GMT -4. The time now is 01:43.