|
[Sponsors] |
October 23, 2024, 21:59 |
Weird physical observation in InterFOAM
|
#1 |
Senior Member
mohammad
Join Date: Sep 2015
Posts: 281
Rep Power: 12 |
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; } } 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; } } 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; } } 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; } } 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; } } 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); } } 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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |