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

pressure dependent viscosity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 23, 2012, 05:51
Default pressure dependent viscosity
  #1
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
Dear FOAMers,

I tried implementing pressure dependent viscosity models for incompressible multiphase flow solvers.

First of all, the quick & dirty work around to access the pressure field from the viscosity model (for example in myPouliquenViscoPlastic.C in viscosityModels/myPouliquenViscoPlastic) without changing all related constructors and header files (inspired by http://freefoam.sourceforge.net/doc/...8C_source.html) was:

U_.db().lookupObject<volScalarField>("p_rgh")

The problem is that one gets pressure oscillations even with a very fine grid and small time-step, and at the ECCOMAS 2012 in Vienna it turned out that this is a general problem when solving incompressible Navier-Stokes eq. with pressure dependent viscosities. However I have heard that there are FOAMers that solved this problem so I would like to invite anybody to present his/her solutions here.

Last edited by vonboett; October 26, 2012 at 04:18.
vonboett is offline   Reply With Quote

Old   October 25, 2012, 11:31
Default
  #2
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 267
Rep Power: 8
kmooney is on a distinguished road
At first glance I would be hesitant to use an incompressible formulation for a pressure dependent viscosity. In this case, only the gradient of the pressure is used. The actual magnitudes of the pressure in p_rgh is arbitrary. if you're using a formulation like eta=f(p_rgh), I don't think you can count on accurate results.
kmooney is offline   Reply With Quote

Old   April 10, 2013, 09:02
Default
  #3
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
It is not that trivial. See http://www.compstat2004.cuni.cz/~mal...ok-prehled.pdf
especially chapter 7.1.
The recently published two-phase debris flow model of Pudasaini uses pressure dependent and shear viscosity for the granular phase by dividing the linearly pressure dependent part through the rate of deformation tensor, so I guess instabillity by pressure is equalized by shear rate and vorticity caused by local viscosity increase...
vonboett is offline   Reply With Quote

Old   May 7, 2013, 11:37
Default pressure dependent viscosity model of Pudasaini in OF
  #4
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
I implemented the actual Mohr-Coulomb plasticity formulation provided by Pudasaini (2012): 'A general two-phase debris flow model' and tested the formulation with a single phase simulation of sand in a made-up flume. The simulation runs fine. When I apply it to other cases (large rotating drum of Berkley, large scale flume experiments at the USGS flume), the simulation becomes unstable and it is not so obvious why. In 'Some remarks on the Navier-Stokes equations with a pressure-dependent viscosity by Michael Renardy (1986) in 'Communications in partial differential equations' it sais in section 2:
"We conclude that, in contrast to the ordinary Navier-Stokes equations, the initial value problem is not always well-posed. Only as long as the eigenvalues of the symmetric part of the velocity gradient are less than 1/(2 nu') [nu being the pressure dependent viscosity] can we expect to prove a local existence result. If this condition is violated, problems of nonexistence and nonuniqueness occur in the constant coefficient problem and hence can be expected in the full Navier-Stokes system"

Since I am not that familiar with math anymore, I am not shure what this means- could someone see if this could explain successful or failing simulations dependent on initial flow conditions?

He further goes into the "complementing condition" of the boundary with respect to the pressure dependent viscosity. I guess but only guess that the first statement means if viscosity (in this case due to pressure) varies strong spatially, solving incompressible NS might fail dependent on velocity gradient, and I further guess that solving incompressible NS might still be present if boundaries are close (due to "complementing condition" additionnally to ellipticity)
So laboratory experiment-sized simulations might run well, but local pressure inhomogenity for example due to impact to a small obstacle will fail in larger scale.

Anyone who is better informed about partial differential equations could confirm that?

Last edited by vonboett; May 8, 2013 at 05:59.
vonboett is offline   Reply With Quote

Old   September 26, 2013, 05:07
Default
  #5
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
by the way, note that nu' denotes derivative of viscosity with respect to pressure. Got a nice solver by now but still work on getting the stability criteria introduced in a good form.
vonboett 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
"Pressure Inlet" Boundary Setup Wijaya FLUENT 13 April 11, 2013 08:50
error message cuteapathy CFX 14 March 20, 2012 07:45
change pressure and viscosity of solid phase NAD Fluent UDF and Scheme Programming 0 January 9, 2012 01:01
Neumann pressure BC and velocity field Antech Main CFD Forum 0 April 25, 2006 02:15
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13


All times are GMT -4. The time now is 12:07.