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

interFoam: strange pressure with high fluid viscosity?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 1, 2010, 08:22
Default interFoam: strange pressure with high fluid viscosity?
  #1
New Member
 
Christian Kröner
Join Date: Nov 2009
Location: Bonn
Posts: 7
Rep Power: 7
ckroener is on a distinguished road
Hi,

I want to model a visco-plastic fluid (using OpenFOAM-1.6.x), therefore i need to handle high viscosities. When I take for example the damBreak tutorial and just increase the viscosity of the fluid to nu [ 0 2 -1 0 0 0 0 ] 1e3,
my pressure field gives very strange values (0.05s and 0.10s):




I have tried to use a finer grid, and to increade the number of nCorrectors, and to decrease the tolerance values for the pressure solver.

It seems to be a problem with the interface, if i change the viscosity of the second fluid to the same viscosity, my pressure field looks fine.

What should I change to get a stable pressure field?

Thanks,
ckroener
ckroener is offline   Reply With Quote

Old   April 1, 2010, 17:15
Default
  #2
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
What does the velocity field look like in these cases?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   April 6, 2010, 04:54
Default the velocity field
  #3
New Member
 
Christian Kröner
Join Date: Nov 2009
Location: Bonn
Posts: 7
Rep Power: 7
ckroener is on a distinguished road
This is how the velocity field looks like:

Time 0.05s, scale factor for vectors: 500


Time: 0.1s, scale factor 5:


I have uploaded the case here: http://geo.uni-bonn.de/members/kroen...roblem.tar.bz2

The problem seems to be in the air phase?

Tanks a lot!
ckroener
ckroener is offline   Reply With Quote

Old   April 7, 2010, 12:51
Default
  #4
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
As far as I know the CSF model for calculating the forces due to surface tension may have some problems when dealing with high density ratios.

In addition with high viscosity ratio you may run in even more problems.

Keeping in mind that the Navier-Stokes-Equations are solved with weighted average values for density \rho = \rho_1 \cdot \gamma + (1-\gamma)\rho_2 and vicosity \nu= \nu_1 \cdot \gamma + (1-\gamma)\nu_2 I'm picturing the situation in my head like the following:

The water column (\gamma=1, the visco-plastic fluid in your case) starts to move and generates high values (due to the high density and viscosity) for most of the terms in the Navier-Stokes-Equations.
These "high" forces act on the air-phase (\gamma=0) and since it incorporates very low (in comparison) values for density and viscosity has to maintain the force balance by creating high velocities.

Correct me if I'm wrong. I'm aware that this is not very usefull for you, but maybe it's good for some brainstorming ...
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Reply

Tags
interfoam

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 simulate exhaust of high pressure tank song CFX 3 January 11, 2009 17:54
My Revised "Time Vs Energy" Article For Review Abhi Main CFD Forum 2 July 9, 2002 09:08
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13
Hydrostatic pressure in 2-phase flow modeling (CFX4.2) HB &DS CFX 0 January 9, 2000 14:19


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