CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

nonNewtonianIcoFoam: errors near boundaries

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

LinkBack Thread Tools Display Modes
Old   January 10, 2013, 10:26
Default nonNewtonianIcoFoam: errors near boundaries
Join Date: Jul 2012
Posts: 31
Rep Power: 7
kpax is on a distinguished road
hey folks,

i am currently using nonNewtonianIcoFoam to simulate the transient blood flow in an artery. BCs for the outlets are constant pressure, and for the inlet i prescribe a velocity profile. For the viscosity I use CrossPowerLaw, but the results with BirdCarreau are very similar.
Two things about the results are very strange:

1) When the inlet velocity is decreasing very rapidly, the viscosity field shows very large gradients at the boundary of the domain (see attached picture). This effect goes away as soon as the velocity is increasing again, but it clearly can't be right.
I thought it could have something to do with discretization, so I tried using the upwind and the QUICK schemes for div(phi,U) and laplacian(nu,U). This had no visible effect.
I also decreased the time step (although it was already very small before: 1e-4s) to 5e-5s. Again, this had no effect.

Apparently the reason for these large nu values are overly small gradients in the velocity field near the boundaries. How can these be avoided?

2) Even when the above effects do not appear, the range of viscosity values seems way to large to me: approximately 3e-6 to 5e-5 (or, multiplied by density: 0.003 Pas to 0.05 Pas). According to literature, this is usually not the case in arteries such as the one i'm modelling.
I realized that when I do steady-state simulations, i.e. constant inlet velocity, the values are in a much more realistic range after around 10 time steps.

I already did a couple of simulations with solvers based on icoFoam using the same mesh, BCs, and discretization schemes, and every thing worked fine. But maybe the mesh is not adequate near the boundaries for non-uniform viscosity?

any ideas/hints why the results are so strange?

thx in advance..
Attached Images
File Type: jpg asdf.jpg (21.8 KB, 12 views)
kpax is offline   Reply With Quote


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 write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Building OpenFOAM1.7.0 from source ata OpenFOAM Installation 42 May 14, 2012 20:48
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

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