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

Stratified two-phase flow in a channel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 26, 2013, 16:35
Default Stratified two-phase flow in a channel
  #1
New Member
 
Join Date: Mar 2009
Posts: 16
Rep Power: 8
Sashankh is on a distinguished road
Hi,

I am trying to solve a liquid-liquid two phase laminar stratified flow in a channel (2D problem) for the velocity distribution. The liquids are immiscible and Newtonian. The interface is stationary. I am solving the Poisson equation in both domains. I impose the Neumann (stress balance) and Dirichlet (velocity balance) condition at the interface. I solve the problem as follows:

1. Initialize the velocity in both domains.
2. Solve domain A using the velocity condition (dirichlet) at the interface.
3. Determine the gradient of velocity at the interface.
4. Share the gradient with domain B.
5. Solve domain B using the gradient condition (Neumann) at the interface as the boundary condition.
6. Determine the velocity at the interface in domain B.
7. Share the velocity with domain A.
8. Solve domain A with the new velocity condition (dirichlet) at the interface.
9. Iterate until convergence.

However, I see that my velocity solution blows up if I increase the viscosity ratio between the two liquids. Is there a numerical stability criteria in solving such problems? I am using a finite-difference scheme.

Any help would be appreciated. Thanks.
Sashankh is offline   Reply With Quote

Old   March 29, 2013, 07:18
Default
  #2
Senior Member
 
kbaker's Avatar
 
Khalid Baker
Join Date: Mar 2009
Location: IRAQ
Posts: 118
Rep Power: 8
kbaker is on a distinguished road
Did you solve Navier Stokes Equations for steady or unsteady state case? what software you tried to solve the governing equations?

kbaker

Quote:
Originally Posted by Sashankh View Post
Hi,

I am trying to solve a liquid-liquid two phase laminar stratified flow in a channel (2D problem) for the velocity distribution. The liquids are immiscible and Newtonian. The interface is stationary. I am solving the Poisson equation in both domains. I impose the Neumann (stress balance) and Dirichlet (velocity balance) condition at the interface. I solve the problem as follows:

1. Initialize the velocity in both domains.
2. Solve domain A using the velocity condition (dirichlet) at the interface.
3. Determine the gradient of velocity at the interface.
4. Share the gradient with domain B.
5. Solve domain B using the gradient condition (Neumann) at the interface as the boundary condition.
6. Determine the velocity at the interface in domain B.
7. Share the velocity with domain A.
8. Solve domain A with the new velocity condition (dirichlet) at the interface.
9. Iterate until convergence.

However, I see that my velocity solution blows up if I increase the viscosity ratio between the two liquids. Is there a numerical stability criteria in solving such problems? I am using a finite-difference scheme.

Any help would be appreciated. Thanks.
kbaker is offline   Reply With Quote

Old   March 29, 2013, 11:09
Default
  #3
Member
 
Join Date: Jul 2011
Posts: 59
Rep Power: 6
rmh26 is on a distinguished road
I'm a little confused on the domain of your problem. Are you looking at flow down a channel (z direction) and solving for the z and y (height) components of the velocity and assuming there is no variation in the x direction?. You are also forcing the interface to be stationary. Is the normal component of the flow at the interface zero? This is the kinematic boundary condition. Also how are you initilializing the velocities?

I am currently working with a two phase flow solver for transient simulation and I do see a decrease in stability when using different viscosities. This comes from the fact that the "cross" terms in the divergence of the stress have to be handled explicitly.

Define your problem a little better and maybe we can find a solution.
rmh26 is offline   Reply With Quote

Old   March 29, 2013, 13:13
Default
  #4
New Member
 
Join Date: Mar 2009
Posts: 16
Rep Power: 8
Sashankh is on a distinguished road
I am solving the Poisson equation

grad.(grad(u)) = 1,

where u is the velocity. The interface conditions are equal velocity and equal shear stress.

I am currently using the Matlab pde toolbox to solve it. This toolbox uses FEM. I have also tried solving this problem using finite-difference in Matlab. Both times I get the same problem and the solution blows up when I use the algorithm listed in my earlier message.

Also, this is a 1D problem. My earlier post says it is 2D which is a typing error. Once I figure out the algorithm that will work for this simple problem I will extend it to my real problem in 2D.

I really appreciate your help. Thanks.
Sashankh is offline   Reply With Quote

Old   March 29, 2013, 13:29
Default
  #5
New Member
 
Join Date: Mar 2009
Posts: 16
Rep Power: 8
Sashankh is on a distinguished road
Hi,

I am sorry but there was an error in my previous post.

I am solving the 1D problem of channel flow (not 2D). The z axis points in the direction of the flow while the y axis is normal to the flow direction (height). Since the flow is fully developed the velocity w (z-component of velocity) does not vary in the z-direction but only in the y direction. Hence I am solving for the velocity w(z-component of velocity) as a function of only y. Hence it is a 1D problem.

There are two fluids flowing parallel to each other. I am solving the Poisson equation in both fluids as discussed in my above post. The interface is rigid and does not change shape. The normal component of velocity is zero. The fluids are immiscible.

Since we have two conditions at the interface (Dirichlet and Neumann) I use the algorithm mentioned in my previous post to solve each domain separately and then solve them iteratively.

The above problem has an analytical solution. I am trying to validate my numerical program with this simple problem before extending it to a more complicated domain.

Thanks for your help and please let me know if there is any other information you need with regards to this problem.
Sashankh is offline   Reply With Quote

Old   March 29, 2013, 14:46
Default
  #6
Senior Member
 
kbaker's Avatar
 
Khalid Baker
Join Date: Mar 2009
Location: IRAQ
Posts: 118
Rep Power: 8
kbaker is on a distinguished road
Hello Sashankh you mentioned that you solve multiphase flow problem? Usually there are three models to solve multiphase flow problems homogeneous model, in-homogeneous model and drift flux model? if you said that your fluids are immisible then you need to follow the homogeneous model but you said that you solve poisson eqn? what the relation of poison equation with the flow? Second thing did you consider the interface as flat?
kbaker is offline   Reply With Quote

Old   March 30, 2013, 13:46
Default
  #7
Member
 
Join Date: Jul 2011
Posts: 59
Rep Power: 6
rmh26 is on a distinguished road
I'm not sure of what exactly MATLAB is doing but assuming that the equation is being discretized and solved correctly then we can assume the issue is with your algorithm. Have you tried a under relaxing method yet. Everytime you solve for velocity you would update your solution as the some combination of your old solution and your new solution.

U(k-1) = solution from previous boundary condition
U* = solution from the current boundary conditions

U(k) = accept new solution

U(k) = (1-w)*U* + w*U(k-1)

if stability is an issue sometimes picking a w close to 1 will help.
rmh26 is offline   Reply With Quote

Old   March 30, 2013, 16:58
Default
  #8
New Member
 
Join Date: Mar 2009
Posts: 16
Rep Power: 8
Sashankh is on a distinguished road
Hi Khalid,

I am solving a Stokes flow problem. Hence my inertia terms drop out of the Navier Stokes equation leaving me with the Poisson equation to solve. I also consider the interface to be flat. Please let me know what you think. Thanks.
Sashankh is offline   Reply With Quote

Old   March 30, 2013, 17:28
Default
  #9
New Member
 
Join Date: Mar 2009
Posts: 16
Rep Power: 8
Sashankh is on a distinguished road
Hi,

I think you are correct about the stability problem with my algorithm. I am trying the under-relaxation method. I will keep you posted on the result.
Is there any other iterative algorithm (one in which I can solve each domain separately and then share node values or gradients at the interface) that you are aware of to solve such problems?
Thanks a lot for your help and I will let you know how the under-relaxation methods works out.
Sashankh is offline   Reply With Quote

Reply

Tags
finite-difference, stratified, two phase

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
Discontinuity at water level in stratified 2 phase flow kbaker CFX 24 June 14, 2012 07:37
[ICEM] Flow channel meshing problems StefanG ANSYS Meshing & Geometry 19 May 15, 2012 06:44
How to calculate phase flow rate? sangramroy FLUENT 0 January 11, 2012 14:02
gas-liquid stratified flow (two phase) farshad pargar FLUENT 6 July 22, 2009 14:46
compressible channel flow.. R.D.Prabhu Main CFD Forum 0 July 17, 1998 17:23


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