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

Problems with OpenFoaminterFoam and ACE

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 6, 2009, 10:08
Default Dear all. I'm modeling a dr
  #1
dingo1234
Guest
 
Posts: n/a
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
{
type zeroGradient;
}
##

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

pd
##
front
{
type symmetryPlane;
}

bottom
{
type symmetryPlane;
}

sciana
{
type zeroGradient;
}

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
  Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
UDF problems Winnie FLUENT 21 March 15, 2011 10:26
Problems with Fluent on simple 1D problems agg FLUENT 3 November 21, 2008 12:55
problems Rogerio Fernandes Brito CFX 1 May 5, 2008 21:08
MPI PROBLEMS gtg627e OpenFOAM Running, Solving & CFD 20 October 5, 2007 04:02
UDF problems Paolo Lampitella FLUENT 5 September 8, 2005 20:43


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