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

nonNewtonianFluidFoam - no grad(nu) & grad(U)

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

Like Tree1Likes
  • 1 Post By JonW

Reply
 
LinkBack Thread Tools Display Modes
Old   November 29, 2012, 11:07
Default nonNewtonianFluidFoam - no grad(nu) & grad(U)
  #1
New Member
 
Janto
Join Date: Oct 2009
Location: Kassel
Posts: 8
Rep Power: 7
chilledkroete is on a distinguished road
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.
chilledkroete is offline   Reply With Quote

Old   November 29, 2012, 11:49
Default Answer: continuum mechanics
  #2
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 53
Rep Power: 6
JonW is on a distinguished road
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.
JonW is offline   Reply With Quote

Old   November 29, 2012, 12:16
Default
  #3
New Member
 
Janto
Join Date: Oct 2009
Location: Kassel
Posts: 8
Rep Power: 7
chilledkroete is on a distinguished road
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

Last edited by chilledkroete; November 29, 2012 at 12:36.
chilledkroete is offline   Reply With Quote

Old   November 29, 2012, 13:30
Default
  #4
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 53
Rep Power: 6
JonW is on a distinguished road
sorry, apparantly I was reading your post too fast,...

good luck,...
J.
JonW is offline   Reply With Quote

Old   August 26, 2013, 20:01
Default
  #5
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Excuse me, where can I download this book?

Theory and Problems of Continuum Mechanics
sharonyue is offline   Reply With Quote

Old   August 27, 2013, 12:00
Default
  #6
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 53
Rep Power: 6
JonW is on a distinguished road
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

Quote:
Originally Posted by sharonyue View Post
Excuse me, where can I download this book?

Theory and Problems of Continuum Mechanics
JonW is offline   Reply With Quote

Old   September 20, 2013, 19:50
Default
  #7
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by JonW View Post
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.
sharonyue is offline   Reply With Quote

Old   September 21, 2013, 13:15
Default
  #8
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 53
Rep Power: 6
JonW is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
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
JonW is offline   Reply With Quote

Old   September 21, 2013, 21:48
Default
  #9
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by JonW View Post
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 is offline   Reply With Quote

Old   November 4, 2013, 09:09
Default
  #10
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
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?
sharonyue is offline   Reply With Quote

Old   November 4, 2013, 11:35
Default
  #11
New Member
 
Janto
Join Date: Oct 2009
Location: Kassel
Posts: 8
Rep Power: 7
chilledkroete is on a distinguished road
Hi,

maybe this picture helps: gradv&gradeta.png

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

Regards,
Janto
chilledkroete is offline   Reply With Quote

Old   November 4, 2013, 11:54
Default
  #12
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by chilledkroete View Post
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]
Attached Images
File Type: jpg yourequation.jpg (4.7 KB, 15 views)
File Type: jpg myequation.jpg (28.7 KB, 17 views)
sharonyue is offline   Reply With Quote

Old   November 4, 2013, 14:00
Default
  #13
New Member
 
Janto
Join Date: Oct 2009
Location: Kassel
Posts: 8
Rep Power: 7
chilledkroete is on a distinguished road
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.
chilledkroete is offline   Reply With Quote

Old   November 4, 2013, 14:07
Default
  #14
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 53
Rep Power: 6
JonW is on a distinguished road
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
Attached Images
File Type: jpg gradU.jpg (96.2 KB, 30 views)
sharonyue likes this.
JonW is offline   Reply With Quote

Old   November 5, 2013, 00:43
Default
  #15
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
NEAT!!!
To Janto
Quote:
Originally Posted by chilledkroete View Post
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!
sharonyue 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
Problem to calculate grad(U) using swak4Foam Hugoles OpenFOAM Post-Processing 11 September 25, 2014 19:07
Isn't the trace of gradU the continuity equation? I need a quantity for elongation alexey7783 OpenFOAM 0 September 9, 2009 03:04
LES correct function luiz eduardo OpenFOAM Running, Solving & CFD 6 September 14, 2007 04:10
GradU at the faces sampaio OpenFOAM Running, Solving & CFD 3 May 11, 2005 15:43
GradU in Post sch CFX 5 September 4, 2003 16:56


All times are GMT -4. The time now is 00:07.