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

2 Phase flow simulation in a T-junction microchannel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 13, 2013, 04:11
Exclamation 2 Phase flow simulation in a T-junction microchannel
  #1
New Member
 
Arun Vasudevan
Join Date: Feb 2013
Posts: 6
Rep Power: 4
prohackrav is on a distinguished road
Hi. I am doing a transient laminar flow simulation to get the slug flow of 2 phase immiscible liquids in a t junction micro channel using VOF model. Its a simple 2d geometry and surface tension modelling is enabled. Whenever my timestep advances I get an error stating that the global COurant number is greater than 250. I reduced the time step to values as low as 1e-5. Still the problem persists. I also reduced the UR factors and increased the iterations per time step. (Pressure-0.2, Density-0.5, Body forces-0.5 Momentum- 0.3 As seen in some forum). The velocity seems to be fine at time 0 but drops to zero as soon as the timestep advances. I have included the snapshot. Can somebody help me?


Attached Images
File Type: jpg Untitled.jpg (62.6 KB, 36 views)
prohackrav is offline   Reply With Quote

Old   November 13, 2013, 04:36
Default
  #2
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
Did you estimate the time step size necessary to achieve a Courant number below 250?
When you say "micro channel" and the velocity magnitude is around 2.5 m/s, I estimate a time step size in the order of 10^{-6} - 10^{-7} to meet this requirement.
flotus1 is offline   Reply With Quote

Old   November 14, 2013, 01:01
Default
  #3
New Member
 
Arun Vasudevan
Join Date: Feb 2013
Posts: 6
Rep Power: 4
prohackrav is on a distinguished road
Quote:
Originally Posted by flotus1 View Post
Did you estimate the time step size necessary to achieve a Courant number below 250?
When you say "micro channel" and the velocity magnitude is around 2.5 m/s, I estimate a time step size in the order of 10^{-6} - 10^{-7} to meet this requirement.
Thanks for the reply. Its 0.8mm microchannel t-junction. I have two opposite inlets one of kerosene and other water. Both with velocity ranging from 0.02-0.1 m/sec. So what would be the estimated time step size?
prohackrav is offline   Reply With Quote

Old   November 14, 2013, 04:20
Default
  #4
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
The Courant number is defined as

C = \frac{u \, \Delta t}{\Delta x} thus \Delta t = \frac{C \, \Delta x}{u}

where \Delta x is the size of a cell. Just to make sure, decrease the time step size further and see if the error persists.
Maybe the problem arises from the initial conditions. Did you try to initialize each half of the domain with the substance at the inlet?
flotus1 is offline   Reply With Quote

Old   November 21, 2013, 01:55
Default
  #5
New Member
 
Arun Vasudevan
Join Date: Feb 2013
Posts: 6
Rep Power: 4
prohackrav is on a distinguished road
Quote:
Originally Posted by flotus1 View Post
The Courant number is defined as

C = \frac{u \, \Delta t}{\Delta x} thus \Delta t = \frac{C \, \Delta x}{u}

where \Delta x is the size of a cell. Just to make sure, decrease the time step size further and see if the error persists.
Maybe the problem arises from the initial conditions. Did you try to initialize each half of the domain with the substance at the inlet?
Do you mean to patch both the inlets till the junction with their corresponding inlet fluids? No. I havent tried that. Will coarsening the mesh help for increasing the time step? I did that instead of reducing the time step further from 1e-5.
prohackrav is offline   Reply With Quote

Old   November 21, 2013, 09:47
Default
  #6
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
Would it be such a mayor inconvenience to reduce the time step size to lets say 10e-8s just to see if the error is really related to the time step size?
flotus1 is offline   Reply With Quote

Old   November 21, 2013, 17:00
Default
  #7
Senior Member
 
Join Date: Jan 2010
Location: Germany
Posts: 268
Rep Power: 8
Zaktatir is on a distinguished road
Could you please give more input about your Simulation?

-Reference density
-I think you are using a tracking method (explicit that's why CFL-dependency) Geo-REconstruct
-Are both phases well separated like a normal slug flow simulation or do they interpenetrate each other? Are small bubbles and liquid entrainment expected to occur or not?
-Which interpolation of Pressure are you using?
-Reference Location of the operating Pressure
-Operating Pressure..
Zaktatir is offline   Reply With Quote

Old   November 21, 2013, 17:53
Post "Monitoring bubble rising velocity in Fluent"
  #8
New Member
 
Ali Masoudi
Join Date: Nov 2013
Posts: 4
Rep Power: 3
amasoudi is on a distinguished road
Hi,
Does anybody know how can I monitor the rising velocity of single bubble in fluent?
amasoudi is offline   Reply With Quote

Reply

Tags
slug flow

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
Bubble plume simulation, multi phase flow Artvandelay Main CFD Forum 3 January 1, 2015 18:11
Mass flow rate of phase in post mat_cfd CFX 0 September 3, 2013 07:55
Convergence of jet flow simulation MiraLisa FLUENT 0 August 15, 2013 04:44
Two Phase immiscible Flow in a microchannel Flo FLUENT 0 October 23, 2008 11:39
low speed compressible two phase flow?? cat CFX 0 November 15, 2005 08:59


All times are GMT -4. The time now is 00:38.