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

Simulation of drop with interFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 27, 2014, 08:57
Default Simulation of drop with interFoam
  #1
New Member
 
abdessamad
Join Date: Apr 2014
Posts: 9
Rep Power: 3
abdessamad is on a distinguished road
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.
Attached Images
File Type: jpg ddropp.jpg (12.7 KB, 31 views)
abdessamad is offline   Reply With Quote

Old   April 27, 2014, 13:54
Default
  #2
Member
 
Andrew Somorjai
Join Date: May 2013
Posts: 77
Rep Power: 4
massive_turbulence is on a distinguished road
Quote:
Originally Posted by abdessamad View Post
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

Last edited by massive_turbulence; April 27, 2014 at 13:56. Reason: added details
massive_turbulence is offline   Reply With Quote

Old   April 27, 2014, 19:53
Default
  #3
New Member
 
abdessamad
Join Date: Apr 2014
Posts: 9
Rep Power: 3
abdessamad is on a distinguished road
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.
abdessamad is offline   Reply With Quote

Old   April 27, 2014, 20:21
Default
  #4
Member
 
Andrew Somorjai
Join Date: May 2013
Posts: 77
Rep Power: 4
massive_turbulence is on a distinguished road
Quote:
Originally Posted by abdessamad View Post
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

transportProperties: CrossPowerLawCoeffs...?

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
massive_turbulence is offline   Reply With Quote

Old   April 28, 2014, 13:30
Default
  #5
New Member
 
abdessamad
Join Date: Apr 2014
Posts: 9
Rep Power: 3
abdessamad is on a distinguished road
it does not work yet
abdessamad is offline   Reply With Quote

Old   April 29, 2014, 17:51
Default
  #6
Member
 
Andrew Somorjai
Join Date: May 2013
Posts: 77
Rep Power: 4
massive_turbulence is on a distinguished road
Quote:
Originally Posted by abdessamad View Post
it does not work yet
Have you tried fixedGradient for the surface? What about zeroGradient. It might be easier for someone to help you if you upload your project.

fixedGradient, zeroGradient?
massive_turbulence is offline   Reply With Quote

Old   April 29, 2014, 19:16
Default this is my project
  #7
New Member
 
abdessamad
Join Date: Apr 2014
Posts: 9
Rep Power: 3
abdessamad is on a distinguished road
thanks Andrew for your time.this is my project .project_drop.tar.gz
abdessamad is offline   Reply With Quote

Old   April 29, 2014, 19:28
Default
  #8
New Member
 
abdessamad
Join Date: Apr 2014
Posts: 9
Rep Power: 3
abdessamad is on a distinguished road
I deleted the folder polymesh because it has a large capacity
Attached Files
File Type: gz project_drop.tar.gz (2.8 KB, 12 views)
abdessamad is offline   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
interFoam simulation yields inconsistent results for alpha1 surface Ralinus OpenFOAM Running, Solving & CFD 8 January 13, 2014 09:54
Convergence of jet flow simulation MiraLisa FLUENT 0 August 15, 2013 04:44
Rigid body motion error when restarting a terminated FSI simulation lingdeer ANSYS 1 May 19, 2013 01:40
simulation of mass transfer of liquid drop in openfoam salman189 OpenFOAM Running, Solving & CFD 5 August 15, 2011 07:40
Interfoam... free surface simulation urgent lostin4ever Main CFD Forum 4 October 12, 2010 08:29


All times are GMT -4. The time now is 16:43.