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

Bug in LES eddy viscosity GenEddyVis.C??-- compressible

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 12, 2013, 22:38
Default Bug in LES eddy viscosity GenEddyVis.C??-- compressible
  #1
Member
 
yijin Mao
Join Date: May 2010
Location: Columbia, MO
Posts: 62
Rep Power: 15
alundilong is on a distinguished road
Here is the function return effective stress tensor term in momentum equation from GenEddyVis.C (OpenFOAM-2.1.1).

tmp<fvVectorMatrix> GenEddyVisc::divDevRhoBeff(volVectorField& U) const
{
return
(
- fvm::laplacian(muEff(), U) - fvc::div(muEff()*dev2(T(fvc::grad(U))))
);
}

my question is
where is the term for the trace of the subgrid scale stresses?
It should looks like this or others..., anyway, k them is missing, and I didn't see anywhere account for it.

tmp<fvVectorMatrix> GenEddyVisc::divDevRhoBeff(volVectorField& U) const
{
return
(
- fvm::laplacian(muEff(), U) - fvc::div(muEff()*dev2(T(fvc::grad(U))))
+ 2.0/3.0*fvc::grad(rho()*k())
);
}

My derivation is attached, please give me some suggestions?
Attached Files
File Type: zip derivationLES.zip (26.0 KB, 7 views)
alundilong is offline   Reply With Quote

Old   January 17, 2013, 12:16
Default
  #2
Member
 
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 89
Rep Power: 16
gregor is on a distinguished road
Hi Edward,

as far as i know you can absorb the missing term into pressure by defining a new pressure as p*=p+p_SGS where p_SGS is -1/3*rho*k_SGS but following the argument of Erlebacher et. al 1991 (Toward the large-eddy simulation of compressible turbulent flows) p.159 the influence off the sgs pressure is neglectible compared to the thermodynamical pressure (for moderate Mach numbers).
gregor is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Boussinesq eddy viscosity assumption siw Main CFD Forum 2 February 17, 2011 13:30
Serious bug in LES interface fs82 OpenFOAM Bugs 21 November 16, 2009 09:15
Turbulent Eddy Dissipation in LES Patrick CFX 8 July 12, 2006 03:44
Help on the initial condition for compressible LES Hall Main CFD Forum 0 August 16, 2005 07:18
give me some advice about compressible LES. Bin Li Main CFD Forum 3 September 19, 2003 14:01


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