
[Sponsors] 
October 23, 2012, 05:51 
pressure dependent viscosity

#1 
Senior Member
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 178
Rep Power: 6 
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 timestep, and at the ECCOMAS 2012 in Vienna it turned out that this is a general problem when solving incompressible NavierStokes 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. 

October 25, 2012, 11:31 

#2 
Senior Member
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA  San Diego, CA USA
Posts: 280
Rep Power: 9 
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.


April 10, 2013, 09:02 

#3 
Senior Member
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 178
Rep Power: 6 
It is not that trivial. See http://www.compstat2004.cuni.cz/~mal...okprehled.pdf
especially chapter 7.1. The recently published twophase 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... 

May 7, 2013, 11:37 
pressure dependent viscosity model of Pudasaini in OF

#4 
Senior Member
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 178
Rep Power: 6 
I implemented the actual MohrCoulomb plasticity formulation provided by Pudasaini (2012): 'A general twophase debris flow model' and tested the formulation with a single phase simulation of sand in a madeup 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 NavierStokes equations with a pressuredependent viscosity by Michael Renardy (1986) in 'Communications in partial differential equations' it sais in section 2:
"We conclude that, in contrast to the ordinary NavierStokes equations, the initial value problem is not always wellposed. 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 NavierStokes 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 experimentsized 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. 

September 26, 2013, 05:07 

#5 
Senior Member
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 178
Rep Power: 6 
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.


Thread Tools  
Display Modes  


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 