CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problems with OpenFoaminterFoam and ACE (http://www.cfd-online.com/Forums/openfoam-solving/58154-problems-openfoaminterfoam-ace.html)

 dingo1234 January 6, 2009 10:08

Dear all. I'm modeling a dr

Dear all.

I'm modeling a drop of the water drop.
First calculations were done using ACE program.
Now I have recalculate this problem using OpenFoam/interFoam.

Results between ACE and interFoam are different, despite fact that both configurations are similar.

All physical units are in mm, ms, g.

Configuration of my project:

##
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

// FoamX Case Dictionary.

FoamFile
{
version 2.0;
format ascii;

root "/opt/OpenFOAM/caelinux-1.4.1/run";
case "3d";
instance "constant/polyMesh";
local "";

class dictionary;
object blockMeshDict;
}

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

arguments "/opt/OpenFOAM/caelinux-1.4.1/run/3d" off;

convertToMeters 1;

vertices
(
(-7e-2 0 0)
( 0 0 0)
( 0 4e-2 0)
(-7e-2 4e-2 0)
(-7e-2 0 4e-2)
( 0 0 4e-2)
( 0 4e-2 4e-2)
(-7e-2 4e-2 4e-2)

);

blocks
(
hex (0 1 2 3 4 5 6 7) (35 20 20) simpleGrading (1 1 1)
);

edges
(
);

patches
(
symmetryPlane front
(
(0 1 2 3)
)
symmetryPlane bottom
(
(0 1 5 4)
)
wall sciana
(
(1 2 6 5)
)
patch atm
(
(4 5 6 7)
(3 7 6 2)
(0 3 7 4)
)
);

mergePatchPairs
(
);

// ************************************************** *********************** //
##

##
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

// FoamX Case Dictionary.

FoamFile
{
version 2.0;
format ascii;

root "/opt/OpenFOAM/caelinux-1.4.1/run";
case "3d";
instance "constant";
local "";

class dictionary;
object environmentalProperties;
}

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

g g [0 1 -2 0 0 0 0] (9.81e-3 0 0);

// ************************************************** *********************** //
##

##
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

// FoamX Case Dictionary.

FoamFile
{
version 2.0;
format ascii;

root "/opt/OpenFOAM/caelinux-1.4.1/run";
case "3d";
instance "constant";
local "";

class dictionary;
object transportProperties;
}

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

phase1
{
transportModel Newtonian;
rho rho [1 -3 0 0 0 0 0] 0.000997;
nu nu [0 2 -1 0 0 0 0] 1e-03;
}

phase2
{
transportModel Newtonian;
rho rho [1 -3 0 0 0 0 0] 1.1614e-6;
nu nu [0 2 -1 0 0 0 0] 1.589e-02;

}

sigma sigma [1 0 -2 0 0 0 0] 3e-4;

// ************************************************** *********************** //
##

##

Initial gamma configuration:
##
x0=-3e-2
y0=0.0
z0=0.0
r=2e-2
xmin=\$x0-\$r
xmax=\$x0+\$r
ymin=\$y0-\$r
ymax=\$y0+\$r
zmin=\$z0-\$r
zmax=\$z0+\$r

funkySetFields -time 0 -field gamma -keepPatches -expression "(sqr(pos().x-\$x0)+sqr(pos().y-\$y0) +sqr(pos().z-\$z0)) <sqr(\$r)>= \$y0 && pos().z >= \$z0 ? 1 : 0"
funkySetFields -time 0 -field U -keepPatches -expression "(sqr(pos().x-\$x0)+sqr(pos().y-\$y0) +sqr(pos().z-\$z0)) <sqr(\$r)>= \$y0 && pos().z >= \$z0 ? vector(6,0,0) : vector(0,0,0)"
##

Boundary conditions:
gamma
##
front
{
type symmetryPlane;
}
bottom
{
type symmetryPlane;
}
sciana
{
type constantGammaContactAngle;
theta0 90;
value uniform 0;
}
atm
{
}
##

U
##
front
{
type symmetryPlane;
}
bottom
{
type symmetryPlane;
}
sciana
{
type fixedValue;
value uniform (0 0 0);
}
atm
{
}
##

pd
##
front
{
type symmetryPlane;
}

bottom
{
type symmetryPlane;
}

sciana
{
}

atm
{
type fixedValue;
value uniform 0;
}
##

When I introduce ConvertToMeters 0.1, and rescale the radius of sphere by factor 10 I obtain correct physical behaviour.
Maybe some physical parameters is missing ?

Thank You.
Jarek

 All times are GMT -4. The time now is 12:28.