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

The viscosity field

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 20, 2001, 06:11
Default The viscosity field
  #1
Mehdi BEN ELHADJ
Guest
 
Posts: n/a
The density field can be calculated in SEM or in HOL as follows:

RHO = RHOG * (1-VFOL) + RHOL * VFOL (it's logical)

but

why ? the viscosity field is calculated as follows:

ENUL = ENULG * (1-VFOL) + ENULL * VFOL

  Reply With Quote

Old   May 21, 2001, 08:30
Default Re: The viscosity field
  #2
Kike
Guest
 
Posts: n/a
Hi, Medhi

In my oppinion it is a simple, but a very nice and strong question? You maybe are "touching" a delicate point inside those methods.

Could you be more specific about why you think it is improper to calculate the viscosity by this expression? In principle, you are averaging the fluids properties.

Kike
  Reply With Quote

Old   May 21, 2001, 16:11
Default Re: The viscosity field
  #3
Mehdi BEN ELHADJ
Guest
 
Posts: n/a
Hi, Enrique

If you take a cell and you want to calculate the density, you can write

Mt = Ml+Mg

where

Mt: total mass, Ml: mass of liquid and Mg: mass of gas.

so,

Mt/Vt = Ml/Vt+Mg/Vt

Vt: total volume of the cell.

and finally

Mt/Vt= (Ml/Vl)*(Vl/Vt) + (Mg/Vg)*(Vg/Vl)

Where

Vl: volume of liquid and Vg: volume of gaz

and we have:

RHO = RHOG * (1-VFOL) + RHOL * VFOL

but I don't know why we can calculate the viscosity in this form:

ENUL = ENULG * (1-VFOL) + ENULL * VFOL

  Reply With Quote

Old   June 12, 2001, 06:21
Default Re: The viscosity field
  #4
Kike
Guest
 
Posts: n/a
Dear Medhi

I am sorry to ansuwer you so late. Effectively there is no physical explanation for this "special average". It is a gross approximation to obtain a transport coefficient for momentum according to the presence of the two phases inside the cell.

In those cells you have the stress balance expresions (normal and tangential to the interface) and you can extract a different relation for the viscosity coefficients ENULG, ENULL and a total viscosity coefficient ENUL from the tangential relation. It is more consistent with the stress continuity at interfaces considering piecewise linear velocity profile in each phase. After some algebraic calculations you will have:

1/ENUL=VFOL*(1/ENULL)+((1-VFOL)*(1/ENULG))

This has been proposed for a professor of CALTEC but I don't remember his name. It was a friend of mine who gave me the hint.

Nevertheless, every body uses the relation you find into PHOENICS and it works for most cases (I mean with respect to experimental results). Try to implement the new relation and see what the differences are. Please tell me what you find after the change in your problem.

Regards

Kike
  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
How to access custom Scalar Field in viscosity model caw OpenFOAM Running, Solving & CFD 0 March 26, 2008 12:11
Putting submesh field values into field on parent mesh helmut OpenFOAM Running, Solving & CFD 2 June 20, 2006 07:31
Turbulent viscosity limited to viscosity ratio... Cyril FLUENT 8 February 7, 2006 21:34
The viscosity field Mehdi BEN ELHADJ Main CFD Forum 1 May 27, 2001 13:50
field by field decomposition vs. component-wise zhanglei Main CFD Forum 0 January 5, 2001 10:06


All times are GMT -4. The time now is 08:41.