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

VOF interface

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 1, 2005, 06:46
Default Aah, my mistake. I was wonderi
  #21
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
Aah, my mistake. I was wondering where g was hiding.

Still, the pressure does behave strangely toward the bottom of the outlet and the removal of gravity makes the problem as described above go away. I will have another look at the formulation when I have some time.
eugene is offline   Reply With Quote

Old   June 1, 2005, 06:55
Default Belay that, after a second loo
  #22
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
Belay that, after a second look I dont believe there is anything wrong at or near the outlet.
eugene is offline   Reply With Quote

Old   June 1, 2005, 07:07
Default Dear Mattijs, thanks for the
  #23
Member
 
Rattin Michele
Join Date: Mar 2009
Posts: 36
Blog Entries: 1
Rep Power: 8
michele is on a distinguished road
Dear Mattijs,
thanks for the help. I tried to compile the new zipUpMesh, but I experienced problems (the compilation of the old zipUpMesh executable was instead straightforward).
I obtained the following errors after issuing wmake:
Making dependency list for source file zipUpCells.C
Making dependency list for source file zipUpMesh.C

SOURCE_DIR=.
SOURCE=zipUpCells.C ; g++ -m64 -DlinuxAMD64 -Wall -W -Wno-unused-parameter -march=opteron -O3 -ffast-math -fno-gcse -DNoRepository -ftemplate-depth-30 -I/home/michele/OpenFOAM/OpenFOAM-1.1/src/OpenFOAM/lnInclude -IlnInclude -I. -fPIC -c $SOURCE -o Make/linuxAMD64Opt/zipUpCells.o
zipUpCells.C: In member function `void Foam::polyMesh::zipUpCells()':
zipUpCells.C:134: error: `WarningIn' undeclared (first use this function)
zipUpCells.C:134: error: (Each undeclared identifier is reported only once for each function it appears in.)
make: *** [Make/linuxAMD64Opt/zipUpCells.o] Error 1

However, if you are interested on a test mesh, the following link is that of the case under discussion:
www.lem3.it/tests/waves_050601/waves.tar.gz

Michele.
PS. Henry, I've just started a simulation with low order approximations, walls whose condition was set to 2 m/s (in rasInterFoam is not coded the slipWall condition, am I right?), and compression coefficient for gamma of 1.0. I will obtain results in about 3 hours...
michele is offline   Reply With Quote

Old   June 1, 2005, 07:29
Default Hi Michele, Sorry, didn't t
  #24
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Hi Michele,

Sorry, didn't try it under 1.1.

Replace in zupUpCells.C

WarningIn("XXXX")

with

Warning<< "XXXX"
mattijs is offline   Reply With Quote

Old   June 1, 2005, 08:38
Default Hi Mattijs, thank you. Now
  #25
Member
 
Rattin Michele
Join Date: Mar 2009
Posts: 36
Blog Entries: 1
Rep Power: 8
michele is on a distinguished road
Hi Mattijs,
thank you.
Now the utility zipUpMesh works correctly. The corrected mesh passes all the tests of the checkMesh.

Thanks,
Michele.
michele is offline   Reply With Quote

Old   June 1, 2005, 09:10
Default Dear Henry: >> Why do you t
  #26
Member
 
Rattin Michele
Join Date: Mar 2009
Posts: 36
Blog Entries: 1
Rep Power: 8
michele is on a distinguished road
Dear Henry:

>> Why do you think there is a problem with the oulet pressure BC? Remember we are solving for pd which is p - rho*g*h, i.e. is constant with depth at the outlet.

Excuse me, Henry, isn't pd = p - g*h ?
I expect for pd the dimensions [m^2/s^2], am I correct?
Another question about the pressure. I'm interested in computing the pressure around floating bodies (and integrate pressure and viscous/turbulent stress over surfaces in order to obtain overall forces).
It is rather complicated to obtain pressure: if I use the expression p = pd + gh, h is not a constant as it depends on the wave elevation above each point.
That is to say that
p = pd + int(-gamma*g, ds)
an integral to be computed along the g direction (gamma is zero in air and one in water, introduced in order to capture the surface).
But below a floating body it is impossible to integrate up to a free surface (this space is occupied by the floating body itself!).
Have you any suggestion?

Ragarding the simulation, I'm making the run with the parameters you suggested.
However my impression is the following: it seems to me that boundary conditions don't pose any problem on the solution: all the instability grows well after the flow enters the domain. It resembles more a turbulence generated instability...
By the way, I would like to get the eddy viscosity among the output results... is it automatically possible or does this require coding?
michele is offline   Reply With Quote

Old   June 1, 2005, 09:49
Default > Excuse me, Henry, isn't pd =
  #27
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
> Excuse me, Henry, isn't pd = p - g*h ?

No pd = p - rho*(g.h) and you can calculate p from pd using this expression. Notice the dimension of pd:

dimensions [1 -1 -2 0 0 0 0];

h is the the position vectors of the cell centres and g is the gavitational force vector.

You can get the eddy viscosity from the turbulence model and write it out.
henry is offline   Reply With Quote

Old   April 2, 2006, 09:24
Default Hello, I also have a simila
  #28
Member
 
nicolas
Join Date: Mar 2009
Location: Glasgow
Posts: 42
Rep Power: 8
nico765 is on a distinguished road
Hello,

I also have a similar problem at the inlet of my 3d domain. As shown in the following pictures (with isosurface at gamma=0.5) (t=0.03), soon after the start of the computation, a 'wave' is developing near the inlet, and gets bigger and bigger.

meshCheck is ok.

my boundaries
pd:
boundaryField
{
inlet
{
type zeroGradient;
}

inlet:002
{
type zeroGradient;
}

outlet
{
type zeroGradient;
}

bottom
{
type zeroGradient;
}

top
{
type totalPressure;
p0 uniform 0;
value uniform 0;
}


right
{
type zeroGradient;
}
hull
{
type zeroGradient;
}
sym
{
type zeroGradient;
}
}

U:
boundaryField
{
inlet
{
type fixedValue;
value uniform (0.5 0 0);
}
inlet:002
{
type fixedValue;
value uniform (2 0 0);
}
outlet
{
type zeroGradient;
}
top
{
type pressureInletOutletVelocity;
phi phi;
rho rho;
value uniform (0 0 0);
}
bottom
{
type fixedValue;
value uniform (0 0 0);
}
sym
{
type fixedValue;
value uniform (0 0 0);
}
right
{
type fixedValue;
value uniform (0 0 0);
}
hull
{
type fixedValue;
value uniform (0 0 0);
}
}




nico765 is offline   Reply With Quote

Old   April 2, 2006, 14:38
Default As I said before, OF is not re
  #29
liu
Senior Member
 
Xiaofeng Liu
Join Date: Mar 2009
Location: State College, PA, USA
Posts: 118
Rep Power: 8
liu is on a distinguished road
As I said before, OF is not ready for open channel yet. At least the inlet and outlet boundary conditions will cause problem.
__________________
Xiaofeng Liu, Ph.D., P.E.,
Assistant Professor
Department of Civil and Environmental Engineering
Penn State University
223B Sackett Building
University Park, PA 16802


Web: http://water.engr.psu.edu/liu/
liu is offline   Reply With Quote

Old   April 2, 2006, 20:56
Default Hi liu,nicolas i did numeri
  #30
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hi liu,nicolas

i did numerical simulations of a 2d hydrofoil in a channel. my preliminary results are 100% ok. i also had issues with the correct BC. but the VOF works for open channel flows. nicolas , please see my posts rasInterFoam - STRANGE RESULTS AT BOUNDARY . after i incorporated Hrv's suggestions my simulations went smooth

regards

kumar
kumar2 is offline   Reply With Quote

Old   April 3, 2006, 11:47
Default Hi, kumar, Which post are you
  #31
liu
Senior Member
 
Xiaofeng Liu
Join Date: Mar 2009
Location: State College, PA, USA
Posts: 118
Rep Power: 8
liu is on a distinguished road
Hi, kumar,
Which post are you refering to? I can't find it. Can you upload your case? I am so interested to see where I did wrong.

Thank you.
__________________
Xiaofeng Liu, Ph.D., P.E.,
Assistant Professor
Department of Civil and Environmental Engineering
Penn State University
223B Sackett Building
University Park, PA 16802


Web: http://water.engr.psu.edu/liu/
liu is offline   Reply With Quote

Old   April 3, 2006, 14:22
Default It's here http://www.cfd-onli
  #32
Member
 
nicolas
Join Date: Mar 2009
Location: Glasgow
Posts: 42
Rep Power: 8
nico765 is on a distinguished road
It's here
http://www.cfd-online.com/cgi-bin/Op...show.cgi?1/656
nico765 is offline   Reply With Quote

Old   April 3, 2006, 14:25
Default oops, actually it s this one
  #33
Member
 
nicolas
Join Date: Mar 2009
Location: Glasgow
Posts: 42
Rep Power: 8
nico765 is on a distinguished road
oops, actually it s this one

http://www.cfd-online.com/OpenFOAM_D...es/1/2014.html
nico765 is offline   Reply With Quote

Old   April 4, 2006, 04:07
Default I have checked the thread, but
  #34
Member
 
nicolas
Join Date: Mar 2009
Location: Glasgow
Posts: 42
Rep Power: 8
nico765 is on a distinguished road
I have checked the thread, but i do something similar in my case. My problem is at the inlet and not at the outlet.
nico765 is offline   Reply With Quote

Old   July 4, 2006, 10:06
Default Can anyone supply or point me
  #35
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
Can anyone supply or point me toward an implementation of outlet boundary conditions for free surface wave transmission?
eugene is offline   Reply With Quote

Old   November 5, 2007, 11:27
Default Hallo, I have the following
  #36
New Member
 
Karl-Heinz Leitz
Join Date: Mar 2009
Posts: 16
Rep Power: 8
khleitz is on a distinguished road
Hallo,
I have the following problem: I have a solid surrounded by air that moves from the left to the right with a constant velocity. However I get a very high pressure wave on the interface solid-air at the inlet and at the outlet. I have the feeling, that something with the boundary conditions is wrong. Can anybody give me a hit what the problem could be?
Regards,
Karl-Heinz
P.S. I would like to add a picture of my problem. Can anybody tell me how to do that?
khleitz is offline   Reply With Quote

Old   November 5, 2007, 11:35
Default Hallo, concerning my problem
  #37
New Member
 
Karl-Heinz Leitz
Join Date: Mar 2009
Posts: 16
Rep Power: 8
khleitz is on a distinguished road
Hallo,
concerning my problem with the strange pressure patterns on the interface solid air. Here is a picture of my current model.

I use a modified interFoam solver. However these pressure waves make my silulation unstable.
Best regards,
Karl-Heinz
khleitz is offline   Reply With Quote

Old   September 29, 2009, 18:31
Default
  #38
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
Hi all, I'm dealing with surface problems in interFoam, the post is:

Setting up interFoam for sloshing

Any comments are welcome. Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   January 7, 2012, 23:43
Default VOF Method
  #39
Member
 
ehsan
Join Date: Mar 2009
Posts: 92
Rep Power: 8
ehsan is on a distinguished road
Dear All

Could I ask you whether you could tell me which VOF technique is used in OpenFOAM? I mean there are different VOF models like Hirt and Nicols, Youngs, and so on. Which one is used in OpenFOAM?

Regards
Ehsan
ehsan is offline   Reply With Quote

Old   January 8, 2012, 10:45
Default
  #40
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
Dear ehsan
1) please ask your question just in one post, not in several posts!
2) look this paper:
Quote:
BRACKBILL, J. U., KOTHE, D. B. & ZEMACH, C. 1992. A Continuum Method for Modeling Surface Tension. Journal of Computational Physics, 100, 335-354.
this method reduces need for interface reconstruction, this method interface is smeared among 2 or 3 cells and it uses a high order schemes
and interface compression method to keep interface sharp
3) if you uses interFoam solver look at here:
Quote:
BERBEROVIĆ, E., HINSBERG, N. P. V., JAKIRLIĆ, S., ROISMAN, I. V. & TROPEA, C. 2009. Drop impact onto a liquid layer of finite thickness: Dynamics of the cavity evolution. The American Physical Society, 79, 036306 (15).
4)if you look for more, look at here:
Quote:
Rusche, H., Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions,
in Department of Mechanical Engineering. 2002, Imperial college of Science: University of London
nimasam 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
2d interface Rogerio Fernandes Brito CFX 2 July 26, 2008 12:25
interface rym FLUENT 2 January 16, 2008 05:55
About Interface BC Philip FLUENT 0 January 3, 2008 02:51
How to set BC at the interface K Tamemy FLUENT 0 June 5, 2006 10:45
interface braket FLUENT 1 November 15, 2005 02:58


All times are GMT -4. The time now is 05:53.