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

divDevReff?

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

Like Tree1Likes
  • 1 Post By sharonyue

Reply
 
LinkBack Thread Tools Display Modes
Old   September 2, 2013, 23:08
Default divDevReff?
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 669
Rep Power: 8
sharonyue is on a distinguished road
Hi guys,

In KEpsilon.C we have:

Code:
tmp<fvVectorMatrix> kEpsilon::divDevReff(volVectorField& U) const
{
    return
    (
      - fvm::laplacian(nuEff(), U)
      - fvc::div(nuEff()*dev(T(fvc::grad(U))))
    );
}
According to my equation, see the image attached.

Should it be:
Code:
return
    (
      - fvm::laplacian(nuEff(), U)
      - fvc::div(nuEff()*T(fvc::grad(U)))
    );
?
Attached Images
File Type: jpg 1.jpg (8.4 KB, 84 views)
vbnhfylbh likes this.
sharonyue is offline   Reply With Quote

Old   September 5, 2013, 04:06
Default The math behind
  #2
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
Take a look at the equations I attached.
Attached Images
File Type: png someMath.png (23.3 KB, 131 views)
GerhardHolzinger is offline   Reply With Quote

Old   September 5, 2013, 09:15
Default
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 669
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by GerhardHolzinger View Post
Take a look at the equations I attached.
Long time no see Gernard, Excellent!!!
Thanks so much
sharonyue is offline   Reply With Quote

Old   September 8, 2013, 23:30
Default
  #4
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 669
Rep Power: 8
sharonyue is on a distinguished road
Woo, but if these two terms are the same? Why should we use the dev instead of neglecting it?
sharonyue is offline   Reply With Quote

Old   September 9, 2013, 02:28
Default
  #5
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
It two terms are equal you are neglecting nothing.

There may be numerical or other reasons for doing so.
GerhardHolzinger 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
Usinng divDevReff vs div(R) + laplacian(nu,U) mchurchf OpenFOAM 5 April 8, 2015 08:29
Calculating divDevReff jposunz OpenFOAM Running, Solving & CFD 13 November 6, 2014 04:31
turbulent energy in UEqn / divDevReff Andreas.Herwig OpenFOAM Running, Solving & CFD 0 April 12, 2011 10:50
About devReff() and divDevReff() in kEpsilon.C Zheng.Zhi OpenFOAM 0 August 24, 2009 02:05


All times are GMT -4. The time now is 23:33.