Opened 3 years ago

Closed 2 years ago

Last modified 2 years ago

#2265 closed Defect (fixed)

Misleading minimum viscosity with Smagorinsky scheme ?

Reported by: jchanut Owned by: gm
Priority: low Milestone:
Component: LDF Version: release-4.0
Severity: minor Keywords:
Cc: Branch review:
MP ready?: Task progress:

Description (last modified by nicolasmartin)


From ldfdyn.F90, the Laplacian Smagorinsky scheme (formula 32) scales the viscosity minimum with the local velocity speed (U) such that:

A = MAX(Asmag, rn_minfac * U * dx / SQRT(4.) )

With the default namelist value rn_minfac=1., it is therefore assumed a Reynolds number (Re) that equals 2., a pretty low value that leads to values well above what gives the Smagorinsky scheme with a standard coefficient. One would have expected a Reynolds number that matches what is used with formula (31), i.e. Re=12 or… alternatively a default value rn_minfac=0.16.
Was this done on purpose ? Did I missed something ?

Commit History (2)


trunk: correct [12194], see #2265


trunk: minor bugfix in ldfdyn, see #2265

Change History (5)

comment:1 Changed 3 years ago by nicolasmartin

  • Description modified (diff)

comment:2 Changed 3 years ago by clevy

  • Owner changed from systeam to gm
  • Status changed from new to assigned

comment:3 Changed 2 years ago by smasson

In 12194:

trunk: minor bugfix in ldfdyn, see #2265

comment:4 Changed 2 years ago by smasson

  • Resolution set to fixed
  • Status changed from assigned to closed

comment:5 Changed 2 years ago by smasson

In 12201:

trunk: correct [12194], see #2265

Note: See TracTickets for help on using tickets.