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

Power-law viscosity model for blood

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   September 19, 2014, 05:21
Exclamation Power-law viscosity model for blood
  #1
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Hello friends,

Hope everything is ok.

I'm opening a new topic despite knowing that there are some other very similar. It happens that none clarifies me as I need...

So, I'm using solver simpleFoam to simulate blood in a 2D geometry as a newtonian and non-newtonian fluid and the only model that is giving me problems is power -law.

The power-law viscosity model has two parameters and, for blood, I found in the paper of Neofytou & Tsangaris (2006) the following values:
  • k = 14.6710−3 Pa.s^n
  • n = 0.7755

I know that OpenFOAM uses kinematic viscosity and when I analysed the code of power-law model, noticed that the consistency index, k, as the same units of nu. So, I simply divided the value of k by rho (1040 kg/m3), and got a value of 1.41110−5 m2/s for k.

In a first phase, I used this value in my simulation and it diverged rapidly.

In a second phase I have initialized all fields with potentialFoam and my simulation still doesn't converge.

I noticed, in some tutorials, that k has values in the order of 2500 and I have some questions:

  1. Is anything escaping me? Am I supposed to obtain a value of k of this magnitude for blood with the data that I have available?
  2. One thing I do not understand is the n exponent in the units of k. Can somebody explain me why is this exponent there? What does it mean?

I really need your help to proceed with my work. Sorry for my newbie question!

Cheers,

FS

Last edited by fs.chemech; September 19, 2014 at 22:00.
fs.chemech is offline   Reply With Quote

Old   September 19, 2014, 14:40
Default
  #2
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
People, can someone give me a help on this question, please?

It's really important to get my doubts solved!

FS
fs.chemech is offline   Reply With Quote

Old   September 21, 2014, 13:21
Default
  #3
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,215
Blog Entries: 35
Rep Power: 94
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings FS,

I quickly went to see your past posts and it seems I told you some months ago how important it is to increase the complexity in your test cases gradually, instead of jumping directly at the final problem

Personally I never use the Power-law viscosity model, therefore I'm not familiar with it. If you can create and/or provide a simple example case, preferably with the expect analytical solution, for example something like the very first "cavity" tutorial explained in the User Guide, it should become pretty clear as to what is what .

Best regards,
Bruno
__________________
___
I'll be at OFW11 in Portugal
wyldckat is offline   Reply With Quote

Old   September 23, 2014, 12:17
Default
  #4
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Hi Bruno!

First of all thank you for your answer! I remember well your advices!

However, I think now the problem is slightly different and I think you could help me...

I'm using a 2D T-shaped geometry (2 inlets and 1 outlet) and studying, with simpleFoam solver:

1) 3 levels of mesh refinement;
2) different schemes for divergence and pressure gradient;
3) Newtonian, Power-law and Bird-Carreau viscosity models.


In a first phase, I studied the coarsest mesh: my simulations successfully converged when applying both type of schemes and the three viscosity models.

For the intermediate and refined mesh, I only can attain convergence when applying newtonian viscosity model. I tried to manipulate relaxation factors for Power-law and Bird-Carreau but without success... Residuals for Ux and P stabilize far away from 1e-05!
I'm using potentialFoam to initialize the fields in all my non-Newtonian simulations...

Attached, there are several files for the intermediate refined mesh using a Power-law viscosity model, as well as a checkMesh file and my residuals evolution...

Can you see them and give me some hints, please?

Do you think that leaving simpleFoam and use nonNewtonianIcoFoam could be a solution?

FS

P.S.: as I said in my first message, I don't understand very well why there's an exponent "n" in the units of the consistency index "k"... In http://en.wikipedia.org/wiki/Power-law_fluid itit is possible to see that "K is the flow consistency index (SI units Pa•sn)"... Also, I don't know very well how this influences the value to input in powerLawCoeffs (transportProperties file)...


checkMesh.txt

p & U.txt

transportProperties.txt

controlDict fvScheme fvSolution.txt

Residuals.gif
fs.chemech is offline   Reply With Quote

Old   September 24, 2014, 09:06
Default
  #5
Member
 
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 3
itchy is on a distinguished road
Hi.

the things about k, nu and n. When I have a look into /src/transportModels/incompressible/viscosityModels/powerLaw/powerLaw.C I see this:

// * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * * //

Foam::tmp<Foam::volScalarField>
Foam::viscosityModels:owerLaw::calcNu() const
{
return max
(
nuMin_,
min
(
nuMax_,
k_*pow
(
max
(
dimensionedScalar("one", dimTime, 1.0)*strainRate(),
dimensionedScalar("VSMALL", dimless, VSMALL)
),
n_.value() - scalar(1.0)
)
)
);
}


k_: flow consistency index
n_:
flow behaviour index
strain-rate: strainRate()
nu_: kin. viscosity

k_ has the dimension of nu_. , because the strainRate() [1/sec] is multyplied by 1sec. ( "dimensionedScalar("one", dimTime, 1.0)"). The equation is like this:

nu_ = k*{strain-Rate()*sec.}^(n-1)

The flow behaviour index n is dimensionless. So everything is in the proper dimension.

Is this the answer of your question??

kind regards
Florian
itchy is offline   Reply With Quote

Old   September 24, 2014, 09:21
Default
  #6
Member
 
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 3
itchy is on a distinguished road
Hi,

-----------fvSchemes---------

Why you use this:

laplacianSchemes
{
default Gauss linear orthogonal;
}

snGradSchemes
{
default orthogonal;
}

instead of:

laplacianSchemes
{
default Gauss linear corrected;
}

snGradSchemes
{
default corrected;
}

????


Please post the complete case. Then we can see

Kind regards
Florian
itchy is offline   Reply With Quote

Old   September 24, 2014, 10:12
Default
  #7
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Florian, thank you very much for your answer!!

Regarding the units of "k", I understood your explanation! Thank you very much!

So, I'm using "orthogonal" and not "corrected" because my mesh is orthogonal. I saw (I don't remember where) that in these cases we can use this definition... Do you suggest me to change this?

By the complete case, you mean my mesh and all the other folders ("0", "constant" and "system")?

Cheers,

FS
fs.chemech is offline   Reply With Quote

Old   September 24, 2014, 10:31
Default
  #8
Member
 
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 3
itchy is on a distinguished road
Hi,

Ah I see in your checkMesh-file:

Mesh non-orthogonality Max: 0 average: 0

o.k. then everything is fine with that.

yes with complete case I mean "0", "constant","system". Mesh is in PolyMesh. I don't need the blockMesh file.

I will give it a shot at my PC. Then I will see what's going wrong with this case.

kind regards
Florian
itchy is offline   Reply With Quote

Old   September 24, 2014, 11:34
Default
  #9
Member
 
Nadish Saini
Join Date: Feb 2014
Location: Raleigh, North Carolina
Posts: 37
Rep Power: 4
90nash is on a distinguished road
Hi 'fs.chemech',

A comment on your original problem. I also faced a similar issue and could not figure out why the solution was diverging for non-newtonian law. Initalise your flow field with a newtonian model and then switch non-newtonian and it wont diverge. This solved the problem for me. Hope it helps!
90nash is offline   Reply With Quote

Old   September 24, 2014, 12:18
Default
  #10
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Florian,

my compressed case (zip) is much bigger (~800 KB) than the max size allowed by the Forum. How can I send it to you?




Nadish,

thank you for your suggestion. You change to non-Newtonian model while Newtonian simulation is running?
fs.chemech is offline   Reply With Quote

Old   September 24, 2014, 12:26
Default
  #11
Member
 
Nadish Saini
Join Date: Feb 2014
Location: Raleigh, North Carolina
Posts: 37
Rep Power: 4
90nash is on a distinguished road
In the controlDict file you may set the option rumTimeModifiable to yes and then switch during runtime. But i guess the better option is to run the computation for say 0.1 sec, halt it, change the power law model and start from 0.1 sec. it will give a good initial value to the solution and i am pretty sure it wont diverge (given you have defined appropriate BC's)
90nash is offline   Reply With Quote

Old   September 24, 2014, 13:02
Default
  #12
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Nadish,

I tried your suggestion... I know that with Newtonian model, my simulation converges in 290 iterations...

So, I ran it until the 260 iteration, stoped it, changed viscosity model to powerLaw and adapted startTime in controlDict file and resumed simulation again...

Unfortunately, residuals Ux and p began to rise and stabilized again near 1e-04...

Do you have any idea why this happens?

Thank you for your time!
fs.chemech is offline   Reply With Quote

Old   September 24, 2014, 13:08
Default
  #13
Member
 
Nadish Saini
Join Date: Feb 2014
Location: Raleigh, North Carolina
Posts: 37
Rep Power: 4
90nash is on a distinguished road
Hi,

You just need to initialize the values using Newtonian model. I don't think running upto 260 iterations is wise. Just run a few... Say 20... Iterations with Newtonian and then switch
90nash is offline   Reply With Quote

Old   September 24, 2014, 13:13
Default
  #14
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Nadish,

I've done it... Residuals continue to stabilize around 1e-04...
fs.chemech is offline   Reply With Quote

Old   September 24, 2014, 13:19
Default
  #15
Member
 
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 3
itchy is on a distinguished road
Hi,

do you have any experimental results for checking accuracy??? Pherhaps your results are good in case of experimental data.

Residuals are not good for checking accuracy. For example there can be some instationarities due to "fluctuations". (Keep in mind, laminar non-Newtonian flow is very similar to turbulent flow). Pherhaps you have transient effects in your simulation.

One hint:
If you want to have good results in case of residuals, then you can do a pseudo time-stepping. The only thing you have to do is to simulate your case unsteady (fvSchemes-ddtSchemes: Euler or backward). I would use backward for pseudo steady cases. You solve your case unsteady until there is no change in your residuals. Some people call this pseudo-transient continuation. This could help you a lot. Furthermore, with pseudo time-stepping, you do not need to init with newtonian flow.

Please let me know if this will work. In a lot of my cases this works really good.

kind regards
Florian
itchy is offline   Reply With Quote

Old   September 24, 2014, 13:36
Default
  #16
Member
 
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 3
itchy is on a distinguished road
Hi again,

if the fluctuations in your residuals will not decay, then you have transient effects in your flow. In this case your problem is unsteady and not steady. To get good avarage data, a steady simulation is not the best way. A better way is to track your field data over time and make an time-avarage (similar to LES or DNS). This can be very tricky, because the data are correlated over time (Autocorrelation). You need several integral time scale, until your results are good due to statistically accuracy.

kind regards
Florian
itchy is offline   Reply With Quote

Old   September 24, 2014, 13:40
Thumbs up
  #17
Member
 
Nadish Saini
Join Date: Feb 2014
Location: Raleigh, North Carolina
Posts: 37
Rep Power: 4
90nash is on a distinguished road
Quote:
Originally Posted by fs.chemech View Post
Nadish,

I've done it... Residuals continue to stabilize around 1e-04...

Glad to have helped out
90nash is offline   Reply With Quote

Old   September 24, 2014, 13:43
Default
  #18
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Hi Florian,

I don't have any experimental results to be based on...

Yes, one of my suspicions is also that I have transient effects affecting my results...
To solve my case as an unsteady case I have to change my solver to nonNewtonianIcoFoam isn't it? simpleFoam, doesn't run unsteady simulations, I think...

You said that residuals are not good for checking accuracy! What, in your opinion, should I be based on?

Cheers,

FS
fs.chemech is offline   Reply With Quote

Old   September 24, 2014, 13:52
Default
  #19
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Quote:
Originally Posted by itchy View Post

"A better way is to track your field data over time and make an time-avarage (similar to LES or DNS). This can be very tricky, because the data are correlated over time (Autocorrelation). You need several integral time scale, until your results are good due to statistically accuracy."

Florian, can you explain me this a little better?
Sorry for my newbie questions...
fs.chemech is offline   Reply With Quote

Old   September 24, 2014, 14:22
Default
  #20
New Member
 
Filippo
Join Date: May 2014
Posts: 27
Rep Power: 3
fs.chemech is on a distinguished road
Hi Nadish,

I wanted to say that I performed your suggestion. However, the simulation still does not converge ... Residuals continues to stabilize well above the desirable, around 1e-04!

Thank you very much for your help!
fs.chemech 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
Herschel-Bulkley non-Newtonian viscosity model has term with sign error pbryant OpenFOAM Bugs 5 June 18, 2013 23:53
Viscosity ratio in gamma-theta transition model based on k-w sst turb model Qiaol618 Main CFD Forum 8 June 9, 2012 06:43
Power Law for Non-Newtonian Viscosity mannobot FLUENT 1 April 23, 2010 09:40
Power Law Viscosity Model cpplabs OpenFOAM Running, Solving & CFD 1 February 13, 2008 09:09
Non-Newtonian power law fluids Tim Phoenics 0 September 17, 2003 16:08


All times are GMT -4. The time now is 15:05.