[Amith Gubbi]

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

[ignayars]

Hi,

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

[TemC]

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.

[ignayars]

Hi!
Here my realization of Papanastasiou regularization of HerschelBulkley. Tested for OF4.1.

Use this post to implement the code.

Originally Posted by open

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:

/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:

/opt/openfoam4/src/transportModels/incompressible/make/files

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 ( newModel.so );

Permissions error...
Use su or sudo bash. Do not forget to source the openfoam toolbox:

 . /opt/openfoam4/etc/bashrc

With best regards,
Yaroslav

[mostanad]

Originally Posted by ignayars

Hi!
Here my realization of Papanastasiou regularization of HerschelBulkley. Tested for OF4.1.

Use this post to implement the code.

With best regards,
Yaroslav

But you are still using as maximum viscosity? Is it really needed?