CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Herschel Bulkley Papanastasiou model

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

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

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:
http://www.bsr.org.uk/rheology-revie...-Mitsoulis.pdf

Any suggestions regarding fixing this is appreciated.

Thanks,
Amith
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
  #2
New Member
 
Yaroslav
Join Date: Nov 2013
Posts: 3
Rep Power: 5
ignayars is on a distinguished road
Hi,

I think some initial velocity distribution should set up. Like that:
Code:
internalField   uniform (0 0 1);
It resolves problem with iterations. But simpleFoam speed convergence is very low.
residuals.png
Uz.png
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,
Yaroslav
ignayars is offline   Reply With Quote

Old   March 2, 2017, 10:24
Default
  #3
New Member
 
Térence
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.

Regards.
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
  #4
New Member
 
Yaroslav
Join Date: Nov 2013
Posts: 3
Rep Power: 5
ignayars is on a distinguished road
Hi!
Here my realization of Papanastasiou regularization of HerschelBulkley. Tested for OF4.1.

Use this post to implement the code.
Quote:
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:
Code:
/opt/openfoam4/src/transportModels/incompressible/viscosityModels/newModel
newModel.c
newModel.H
If the newModel includes:
addToRunTimeSelectionTable...

Then update the make file by appending the newModel to:
Code:
/opt/openfoam4/src/transportModels/incompressible/make/files
From the transportModels dir, run:
Code:
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:
Code:
libs ( newModel.so );
Permissions error...
Use su or sudo bash. Do not forget to source the openfoam toolbox:
Code:
 . /opt/openfoam4/etc/bashrc
With best regards,
Yaroslav
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

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
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.