CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

How should we define a multiphase interface

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 5, 2011, 09:24
Default How should we define a multiphase interface
  #1
New Member
 
Dijun Pan
Join Date: Jun 2011
Posts: 13
Rep Power: 6
zcesd47 is on a distinguished road
Hi, all

I am trying to define a inviscid water-air interface with a cos profile as the intial condition. The purpose is to model the sloshing behavior of water in a tank excited by sinusoidal external excitation. As the sloshing will eventually reach to the steady state (sloshing wave will has the same amplitude during the time history), we can defined a wave with a initial cos profile and start simulation.


I defined the cos interface with CEL and vof mehod, as following

height:0.03[m]*cos(2*pi[rad m^-1]/1.7680*x)+0.198[m]
water volum fraction: step((height-z)/1[m])

However, the interface profile defined as above is not smooth (like sawtooth) and there are some undefined mesh elements,


Can anybody give me some suggestion about how to define a smooth free surface profile, please, thank you very much.
Attached Images
File Type: jpg 1.jpg (54.7 KB, 42 views)
zcesd47 is offline   Reply With Quote

Old   July 5, 2011, 19:29
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,823
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
There should be any undefined elements. What elements are undefined?

You need to put a blurring function on the interface. This can be done with a linear blend over a few elements or more complex functions - but it will take a bit of thinking to develop the equation.

Alternately you can just start with the jagged interface and run it for a while to smooth out the jaggedness.
ghorrocks is offline   Reply With Quote

Old   July 6, 2011, 08:10
Default
  #3
New Member
 
Dijun Pan
Join Date: Jun 2011
Posts: 13
Rep Power: 6
zcesd47 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
There should be any undefined elements. What elements are undefined?

You need to put a blurring function on the interface. This can be done with a linear blend over a few elements or more complex functions - but it will take a bit of thinking to develop the equation.

Alternately you can just start with the jagged interface and run it for a while to smooth out the jaggedness.
Hi, Horrocks, thank you so much for you suggestions and help.

There are some cells in grey color (in the graph) when I check the water.volume fraction in post-processer, are those undefined cell? For other cases I did there are more 'grey cells'.

I will try to define to define the volume fraction of the cells on the interface as a hypobolic function,sorry, I am new to CFD, could you suggest me some name of the function may be useful, eg. function that may be applied for selecting cells on the interface. Great appriciated.

My purpose is to model the sloshing behavior in a fluid tank when it reach the steady state (constant wave crest level). So to check the linear wave behavior I define a initial water surface as above, with small amplitude, and set initial velocity to 0, pressure to hydrostatic pressure, but the wave amplitude still has some decay with 0 viscousity, so I thought the reason maybe the jaggedness. Is it possible to extract the volume fraction of the whole mesh from solution data after jaggedness are smoothed out, then put into mesh in CFX-pre and start another run?Thank you very much.


Sincerely

Dijun
zcesd47 is offline   Reply With Quote

Old   July 6, 2011, 19:20
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,823
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
No elements should be undefined. Have a close look at them and check they are OK before doing anything else.

The model will always have some decay due to numerical diffusion. You can reduce this with fine meshes and second order differencing schemes, but you cannot eliminate it.

Yes, you can use a previous result as the initial condition for a new run.
ghorrocks is offline   Reply With Quote

Old   July 7, 2011, 09:00
Default
  #5
New Member
 
Dijun Pan
Join Date: Jun 2011
Posts: 13
Rep Power: 6
zcesd47 is on a distinguished road
Hi Horrocks,

Thank you so much for your kind help. Could you give me some more help? Could you tell me the way of checking the volume fraction of the cells of the mesh,please.

I used the default second-order backward Euler scheme. But the wave amplitude actually increased by some reason in the first two periods, and decrease, and then increase again for two periods. As the water viscousity was set to zero, I expect the wave height should be constant, so I guess it may because of the initial condition. Do you feel it is wrong to set the initial velocity in all direction to zero? I am greatly appreciated for your help.

Sincerely

Dijun Pan
zcesd47 is offline   Reply With Quote

Old   July 7, 2011, 19:11
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,823
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Don't forget about the spatial differencing scheme. If you need to minimise diffusion then I recommend you use pure second order differencing (ie hybrid with 1 blend factor).
ghorrocks is offline   Reply With Quote

Old   July 8, 2011, 10:00
Default
  #7
New Member
 
Dijun Pan
Join Date: Jun 2011
Posts: 13
Rep Power: 6
zcesd47 is on a distinguished road
Thank you very much, really appreciated for your help.
zcesd47 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
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
REAL GAS UDF brian FLUENT 6 September 11, 2006 08:23
UDF FOR UNSTEADY TIME STEP mayur FLUENT 3 August 9, 2006 10:19
Multiphase Flow Problem Robi FLUENT 8 July 23, 2005 05:02
defining interface zone more than once Tomas FLUENT 2 November 26, 2004 07:05


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