CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Simulation of drop with interFoam (https://www.cfd-online.com/Forums/openfoam-solving/134286-simulation-drop-interfoam.html)

abdessamad April 27, 2014 08:57

Simulation of drop with interFoam
 
1 Attachment(s)
Hi everybody!
i'm a beginner in openFoam. i have a problème in my simulation. i want to simulate a drop with a solver interFoam. but my drop moves to the right or left, but i want it remains in the middle.this image can explain more the situation.

massive_turbulence April 27, 2014 13:54

Quote:

Originally Posted by abdessamad (Post 488484)
Hi everybody!
i'm a beginner in openFoam. i have a problème in my simulation. i want to simulate a drop with a solver interFoam. but my drop moves to the right or left, but i want it remains in the middle.this image can explain more the situation.

I'm guessing the friction between the drop and the surface is very low. I'm also guessing that you have a file with transportProperties in your constant folder. You may need to check the phase (something that looks like this),

Code:

phase2
{
    transportModel  Newtonian;
    nu              nu [ 0 2 -1 0 0 0 0 ] 1.48e-05;
    rho            rho [ 1 -3 0 0 0 0 0 ] 1;
    CrossPowerLawCoeffs
    {
        nu0            nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
        nuInf          nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
        m              m [ 0 0 1 0 0 0 0 ] 1;
        n              n [ 0 0 0 0 0 0 0 ] 0;
    }

    BirdCarreauCoeffs
    {
        nu0            nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;
        nuInf          nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
        k              k [ 0 0 1 0 0 0 0 ] 99.6;
        n              n [ 0 0 0 0 0 0 0 ] 0.1003;
    }
}

you need to create a no-slip condition on the boundary.

http://www.princeton.edu/~asmits/Bic...tionflows.html

abdessamad April 27, 2014 19:53

Thank you Andrew Somorjai for your reply. i have the same file of transportpropreties like that:
/*--------------------------------*- 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 (water air);

water
{
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1e-06;
rho rho [ 1 -3 0 0 0 0 0 ] 1000;
CrossPowerLawCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
m m [ 0 0 1 0 0 0 0 ] 1;
n n [ 0 0 0 0 0 0 0 ] 0;
}

BirdCarreauCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
k k [ 0 0 1 0 0 0 0 ] 99.6;
n n [ 0 0 0 0 0 0 0 ] 0.1003;
}
}

air
{
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1.48e-05;
rho rho [ 1 -3 0 0 0 0 0 ] 1;
CrossPowerLawCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
m m [ 0 0 1 0 0 0 0 ] 1;
n n [ 0 0 0 0 0 0 0 ] 0;
}

BirdCarreauCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
k k [ 0 0 1 0 0 0 0 ] 99.6;
n n [ 0 0 0 0 0 0 0 ] 0.1003;
}
}

sigma sigma [ 1 0 -2 0 0 0 0 ] 0.07;


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


but i think the probleme is a "slip condition" but i dont knew how to add it in boundary condition. please if you have any ideas tell me about. thank u very much.

massive_turbulence April 27, 2014 20:21

Quote:

Originally Posted by abdessamad (Post 488544)
Thank you Andrew Somorjai for your reply. i have the same file of transportpropreties like that:
/*--------------------------------*- 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 (water air);

water
{
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1e-06;
rho rho [ 1 -3 0 0 0 0 0 ] 1000;
CrossPowerLawCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
m m [ 0 0 1 0 0 0 0 ] 1;
n n [ 0 0 0 0 0 0 0 ] 0;
}

BirdCarreauCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
k k [ 0 0 1 0 0 0 0 ] 99.6;
n n [ 0 0 0 0 0 0 0 ] 0.1003;
}
}

air
{
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1.48e-05;
rho rho [ 1 -3 0 0 0 0 0 ] 1;
CrossPowerLawCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
m m [ 0 0 1 0 0 0 0 ] 1;
n n [ 0 0 0 0 0 0 0 ] 0;
}

BirdCarreauCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
k k [ 0 0 1 0 0 0 0 ] 99.6;
n n [ 0 0 0 0 0 0 0 ] 0.1003;
}
}

sigma sigma [ 1 0 -2 0 0 0 0 ] 0.07;


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


but i think the probleme is a "slip condition" but i dont knew how to add it in boundary condition. please if you have any ideas tell me about. thank u very much.

According to this

http://www.cfd-online.com/Forums/ope...lawcoeffs.html

he says that the CrossPowerLawCoeffs and BirdCarreauCoeffs are for non-Newtonian viscosity models. So that might be one thing to think about. As for boundary conditions those would be in your 0 folder within the U file for the surface but I'm not sure what definition you would change for the type (if at all, maybe zeroGradient would work).

Here's a nice tutorial.

http://www.tfd.chalmers.se/~hani/kur...Hemida_VOF.pdf

abdessamad April 28, 2014 13:30

it does not work yet :confused:

massive_turbulence April 29, 2014 17:51

Quote:

Originally Posted by abdessamad (Post 488719)
it does not work yet :confused:

Have you tried fixedGradient for the surface? What about zeroGradient. It might be easier for someone to help you if you upload your project.

http://www.cfd-online.com/Forums/ope...ogradient.html

abdessamad April 29, 2014 19:16

this is my project
 
1 Attachment(s)
thanks Andrew for your time.this is my project .Attachment 30489

abdessamad April 29, 2014 19:28

1 Attachment(s)
I deleted the folder polymesh because it has a large capacity


All times are GMT -4. The time now is 15:11.