View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0000389 | OpenFOAM | Bug | public | 2012-01-17 08:53 | 2012-01-20 14:46 |
Reporter | Assigned To | henry | |||
Priority | normal | Severity | minor | Reproducibility | always |
Status | closed | Resolution | fixed | ||
Platform | Linux | OS | Other | OS Version | (please specify) |
Summary | 0000389: inconsistent way of defining k for compressible Smagorinsky model | ||||
Description | Looking at the compressible Smagorinsky model i find two way how k gets evaluated. In the Smagorinsky.C in the updateSubgridScaleFields function: volScalarField k(sqr((-b + sqrt(sqr(b) + 4*a*c))/(2*a))); is used, whereas in the k() function (2.0*ck_/ce_)*sqr(delta())*magSqr(dev(symm(gradU))); is taken. | ||||
Tags | No tags attached. | ||||
|
Thanks for the bug report. Resolved in OpenFOAM-2.1.x by commit 09e2f6ca628f33fe8f5ccdd4a09f2492cc969924 |
|
Are there any reasons / references for the quadratic k equation ? In the literature (ie. Yoshizawa1985,Fureby1996) i find k only as k=c_I*Delta*||D||^2. thanks |
|
The quadratic form is a natural consequence of the effect of compressibility on the sub-grid stress. Take a look at the description in the header: B = 2/3*k*I - 2*nuSgs*dev(D) where D = symm(grad(U)); k from rho*D:B + ce*rho*k^3/2/delta = 0 muSgs = ck*rho*sqrt(k)*delta |
|
Ok i found my mistake. Thanks |
Date Modified | Username | Field | Change |
---|---|---|---|
2012-01-17 08:53 |
|
New Issue | |
2012-01-20 11:15 | henry | Note Added: 0000941 | |
2012-01-20 11:15 | henry | Status | new => resolved |
2012-01-20 11:15 | henry | Resolution | open => fixed |
2012-01-20 11:15 | henry | Assigned To | => henry |
2012-01-20 13:38 |
|
Note Added: 0000942 | |
2012-01-20 13:38 |
|
Status | resolved => feedback |
2012-01-20 13:38 |
|
Resolution | fixed => reopened |
2012-01-20 13:47 | henry | Status | feedback => acknowledged |
2012-01-20 13:49 | henry | Note Added: 0000943 | |
2012-01-20 13:49 | henry | Status | acknowledged => closed |
2012-01-20 13:49 | henry | Resolution | reopened => fixed |
2012-01-20 14:29 |
|
Note Added: 0000945 | |
2012-01-20 14:29 |
|
Status | closed => feedback |
2012-01-20 14:29 |
|
Resolution | fixed => reopened |
2012-01-20 14:46 | henry | Status | feedback => closed |
2012-01-20 14:46 | henry | Resolution | reopened => fixed |