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

dambreak switched properties boundary

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 28, 2013, 04:29
Default dambreak switched properties boundary
  #1
New Member
 
Join Date: Mar 2009
Posts: 11
Rep Power: 0
story is on a distinguished road
Hi everybody,

I was trying out the dambreak tutorial, but swithed the transport properties.
I wanted the air to rise out of the domain as a rising bubble. I have gone through several threads in this forum and also found the single bubble rising papers of Hysing et al. and Klostermann et al. However in these cases the rising bubble does not leave the domain.

I wanted to simulate my bubble leaving the top boundary. Therefor as an example I switched the properties in the dambreak tutorial. All sides are all wall except the top boundary. For U it is pressureInletOutletVelocity and p_rgh is the totalPressure. For alpha1 I have set to calculated. However, the simulation crashes once the air reaches the top boundary.
Does any of you have a different boundary condition?

For other boundary conditions the simulation runs fine, however I do get a free surface on top for zeroGradient, inletOutlet and outletInlet condition for the alpha. But I want my alpha1 to leave the top boundary.
Do you think I should change also the other boundary conditions for U and p_rgh?

Thank you all in advance.
story is offline   Reply With Quote

Old   October 28, 2013, 06:38
Default
  #2
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
your question is too vague, however choosing calculated B.C. for alpha1 is totally wrong.
you can simply use zeroGradient for U and alpha1 and fixedValue for pressure, then bubble would leave the domain.
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is online now   Reply With Quote

Old   October 28, 2013, 09:48
Default
  #3
New Member
 
Join Date: Mar 2009
Posts: 11
Rep Power: 0
story is on a distinguished road
Hi Nimasam. Thanks for the reply.

You're right. My question does seems vague.
As I stated I started with the dambreak laminar interfoam (openfoam 2.2.2) tutorial as an initial test for gas rising out of domain starting with a initial shape.

I just renamed phase1 tot phase2 and phase2 to phase1 in the transportProperties file. That is the only change I make besides the boundary conditions.

But If I set the BCs as you suggested to zerogradient for U and alpha on top and p on top to a fixedValue, my phase1 still does not leave the domain. I still get a free surface on top. I attached the initial and end state images of the simulation with your suggested BCs.
Got any extra pointers?

As per your suggestion the boundaries are:

In the original tutorial files:

0/alpha1:
atmosphere
{
type zeroGradient;
}

0/U:
atmosphere
{
type zeroGradient;
}

0/p_rgh:

atmosphere
{
type fixedValue;
value uniform 1;
}

Thanks again.

initial:
state0.0.jpg

end:
state1.0.jpg
story is offline   Reply With Quote

Old   October 28, 2013, 10:08
Default
  #4
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
consider for alpha1
Code:
type   inletOutlet;
inletValue   uniform 0;
value   uniform 0;
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is online now   Reply With Quote

Old   October 28, 2013, 10:22
Default
  #5
New Member
 
Join Date: Mar 2009
Posts: 11
Rep Power: 0
story is on a distinguished road
I have tried it also. But in the end, phase1 still does not exit my domain.
I attached the whole tutorial with the minor adjustment of myself. As you can see, I have set up the tutorial with no major changes. It is still strange that for zeroGradient of alpha1, it forms a free surface on top.
Do you think the OpenFOAM version matters?

damBreak.zip
story 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
Implementation of boundary conditions for FVM Tom Main CFD Forum 7 August 26, 2014 05:58
Low Mixing time Problem Mavier CFX 5 April 29, 2013 00:00
domain imbalance for enrgy equation happy CFX 14 September 6, 2012 01:54
Boundary conditions/flow properties Manish Main CFD Forum 0 May 10, 2007 06:16
Boundary conditions? Tom Main CFD Forum 0 November 5, 2002 02:54


All times are GMT -4. The time now is 18:52.