divDevReff? -- CFD Online Discussion Forums

September 2, 2013, 23:08

divDevReff?

sharonyue

Senior Member

Dongyue Li

Join Date: Jun 2012

Location: Beijing, China

Posts: 839

Rep Power: 17

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

1.jpg (8.4 KB, 190 views)

vbnhfylbh likes this.

« Previous Thread | Next Thread »

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 Off Pingbacks are On Refbacks are On Forum Rules

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Calculating divDevReff	jposunz	OpenFOAM Running, Solving & CFD	82	June 12, 2022 14:01
Usinng divDevReff vs div(R) + laplacian(nu,U)	mchurchf	OpenFOAM	5	April 8, 2015 08:29
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 11:28.