 |
 |
 |
|
|
|
[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. |
|
|
|
|
 |
|
|
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 |
Linear Mode
