[chilledkroete]

Hello folks,

does anybody know why in this solver no grad(nu) & grad (U) is taken into account in Ueqn?

In interFoam one finds
fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muEff, U)
- (fvc::grad(U) & fvc::grad(muEff))
);
due to the fact that muEff is variable in space.

In nonNewtonianIcoFoam
- (fvc::grad(U) & fvc::grad(muEff))
is ignored.

Thanks for replies.

[JonW]

Hello there

the additional term "- (fvc::grad(U) & fvc::grad(muEff))" is correct and is now present in nonNewtoninanIcoFoam in the OpenFOAM 2.1.x. You can also see this term in the interFoam solver.

If you are only used to standard fluid mechanics, you will be puzzled by this term. If you want to understand the physical reason for this term, you have to go into continuum mechanics. My favourite book about the subject is...

G. E. Mase. Schaums Outline Series: Theory and Problems of Continuum Mechanics. McGraw-Hill Inc., USA, 1970.


Hope this helps
J.

[chilledkroete]

Hello Jon,

thank you for your answer. I just wondered why this term was missing in the solver I used so far (v.1.6-ext).

I actually thought it was missing.

I recently found out that you reported this bug in April.

Thanks a lot,
janto

[JonW]

sorry, apparantly I was reading your post too fast,...

good luck,...
J.

[sharonyue]

Excuse me, where can I download this book?

Theory and Problems of Continuum Mechanics

[JonW]

You have to buy the book (I dont know of any download). But it is cheep.
See for example:

http://www.amazon.com/Schaums-Outlin...mechanics+mase

Originally Posted by sharonyue

Excuse me, where can I download this book?

Theory and Problems of Continuum Mechanics

[sharonyue]

Originally Posted by JonW

You have to buy the book (I dont know of any download). But it is cheep.
See for example:

http://www.amazon.com/Schaums-Outlin...mechanics+mase

Sorry, so late.

This book is not on sell in our country. Do you have another articles? Thanks.

[JonW]

Originally Posted by sharonyue

Sorry, so late.

This book is not on sell in our country. Do you have another articles? Thanks.

Ohh, that’s bad.
There is another book I use :
L. E. Malvern. Introduction to the Mechanics of Continuous Medium. Prentice-Hall, Inc., USA, 1969.

...but that book is much more expensive

The only thing I can suggest if you have to download is search by using the string: "continuum mechanics pdf"
hope this help

Cheers
Jon

[sharonyue]

Originally Posted by JonW

Ohh, that’s bad.
There is another book I use :
L. E. Malvern. Introduction to the Mechanics of Continuous Medium. Prentice-Hall, Inc., USA, 1969.

...but that book is much more expensive

The only thing I can suggest if you have to download is search by using the string: "continuum mechanics pdf"
hope this help

Cheers
Jon

Okay...Thanks, but this book is too long. My primary research is not this. So Jon, do you have any notes by yourself?

[sharonyue]

I still cant get out why there is (fvc::grad(U) & fvc::gradfluid.nu()). It exists in nonnewtonianIcoFoam and interFoam. Who can give me a hint pls?

[chilledkroete]

Hi,

maybe this picture helps: gradv&gradeta.png

The first term is equivalent with eta*laplace(v).

Regards,
Janto

[sharonyue]

Originally Posted by chilledkroete

Hi,

maybe this picture helps: Attachment 26578

The first term is equivalent with eta*laplace(v).

Regards,
Janto

Hi Janto,
You mean
[LaTeX Error: Syntax error]
Why shouldnt it be:
[LaTeX Error: Syntax error]

[chilledkroete]

Hello Forrest,

you would be right if the viscosity was constant all over the domain. However, in interFoam and nonNewtonianIcoFoam the viscosity is a field variable (have a look in createFields where the viscosity is defined) and depends therefore on the spatial coordinates. Owing to that the gradient of the viscosity does not vanish.

[JonW]

Hello Sharonyue

Here is something I did in haste (there might be some typo errors, like vector arrow symbol missing over the velocity U), but basically this is the approach.
You would not use array approach as you showed earlier, but rather directly work with dyads and vectors directly. Such is much easier. I hope you can use this document, but you need to know continuum mechanics,... there is no way around it, if you are going to do serious CFD work.

Cheers and good luck
Jon

[sharonyue]

NEAT!!!
To Janto

Originally Posted by chilledkroete

Owing to that the gradient of the viscosity does not vanish.

Yeah, I think this is where I was wrong.

To Jon:

Wonderful, It helps me a lot!!

Thank you both very much!