CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Herschel Bulkley Papanastasiou model

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

LinkBack Thread Tools Display Modes
Old   April 25, 2016, 04:55
Default Herschel Bulkley Papanastasiou model
New Member
Amith Gubbi
Join Date: Apr 2016
Location: Sweden
Posts: 1
Rep Power: 0
Amith Gubbi is on a distinguished road

I am trying to simulate a herchel-bulkley model of blood through a 3D pipe and have implemented the Papanastasiou regularization. The model seems to work fine for nonNewtonianIcoFOAM(transient solver)but when i try to simulate it in a steady state solver(simpleFOAM), I get a floating point exception. I have attached the case setup with the log file and also the papanastasiou model that I have implemented.

The Papanastasiou regularization is according to the following paper:

Any suggestions regarding fixing this is appreciated.

Attached Files
File Type: gz papanastasiou.tar.gz (4.3 KB, 14 views)
Amith Gubbi is offline   Reply With Quote

Old   May 28, 2016, 07:29
Default Re:Herschel Bulkley Papanastasiou model Reply to Thread
New Member
Join Date: Nov 2013
Posts: 3
Rep Power: 5
ignayars is on a distinguished road

I think some initial velocity distribution should set up. Like that:
internalField   uniform (0 0 1);
It resolves problem with iterations. But simpleFoam speed convergence is very low.
I dont think that length of tube should be so big to find solution.
I hope you found solution, please share it

With best regards,
ignayars is offline   Reply With Quote

Old   March 2, 2017, 10:24
New Member
Join Date: Feb 2017
Location: Bordeaux, France.
Posts: 27
Rep Power: 2
TemC is on a distinguished road
Hi Amith Gubbi,

Can you please tell me how you managed to build AND use your new viscosity model?

I'am working with simpleFoam, and I have implemented a new viscosity model based on the Papanastasiou regularization. But each time I try to run a case using my new viscosity model, I got an error message suggesting that my new viscosity model is not valid (attachments).

Thanks in advance for your time, and have a nice week.

Attached Images
File Type: jpg Simulation_Failure1.jpg (123.1 KB, 4 views)
File Type: png Simulation_Failure2.PNG (96.9 KB, 3 views)

Last edited by TemC; March 3, 2017 at 04:31.
TemC is offline   Reply With Quote

Old   March 31, 2017, 02:24
Default Papanastasiou regularization
New Member
Join Date: Nov 2013
Posts: 3
Rep Power: 5
ignayars is on a distinguished road
Here my realization of Papanastasiou regularization of HerschelBulkley. Tested for OF4.1.

Use this post to implement the code.
Originally Posted by open View Post
I'm posting here since this is one of the most recent cries for help over "valid viscosityModels are" error message.

Put the new transport model into the source directory:
If the newModel includes:

Then update the make file by appending the newModel to:
From the transportModels dir, run:
wmake libso incompressible
Or, compile the individual library (wmake libso from the source folder) and then add the model to controlDict by appending the line:
libs ( );
Permissions error...
Use su or sudo bash. Do not forget to source the openfoam toolbox:
 . /opt/openfoam4/etc/bashrc
With best regards,
Attached Images
File Type: png U.png (99.5 KB, 1 views)
Attached Files
File Type: gz HerschelBulkleyPapanastasiou.tar.gz (2.2 KB, 2 views)
ignayars 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
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
Coupling RPI wall boiling model with population balance model in Fluent softice2006 Fluent Multiphase 0 February 1, 2015 16:55
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
Herschel Bulkley model Problem emad FLUENT 5 May 18, 2012 13:25
Regarding Implementation of Hershel Bulkley Model in Viscosity correlation Pathak FLOW-3D 0 August 24, 2009 06:34

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