|
[Sponsors] |
Report a bug in LienCubicKELowReSetWallDissipation.H |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 20, 2011, 00:37 |
Report a bug in LienCubicKELowReSetWallDissipation.H
|
#1 |
Member
|
This head file sets the value of epsilon in the near wall cell. In order to consider the possibility that one cell may contain more than one wall-boundary faces. It contains the average function as follow.
Code:
forAll(patches, patchi) { const fvPatch& curPatch = patches[patchi]; if (isA<wallFvPatch>(curPatch)) { forAll(curPatch, facei) { label faceCelli = curPatch.faceCells()[facei]; // For corner cells (with two boundary or more faces), // epsilon in the near-wall cell are calculated as an average cellBoundaryFaceCount[faceCelli]++; epsilon_[faceCelli] += Cmu75*pow(k_[faceCelli], 1.5) /( kappa_.value()*y_[faceCelli] *(1.0 - exp(-Aepsilon_.value()*yStar_[faceCelli])) ) *exp(-Amu_.value()*sqr(yStar_[faceCelli])); } } } // perform the averaging forAll(patches, patchi) { const fvPatch& curPatch = patches[patchi]; if (isA<wallFvPatch>(curPatch)) { forAll(curPatch, facei) { label faceCelli = curPatch.faceCells()[facei]; epsilon_[faceCelli] /= cellBoundaryFaceCount[faceCelli]; } } } Code:
forAll(patches, patchi) { const fvPatch& curPatch = patches[patchi]; if (isA<wallFvPatch>(curPatch)) { forAll(curPatch, facei) { label faceCelli = curPatch.faceCells()[facei]; // For corner cells (with two boundary or more faces), // epsilon in the near-wall cell are calculated as an average cellBoundaryFaceCount[faceCelli]++; label n = cellBoundaryFaceCount[faceCelli]; scalar eps_n = Cmu75*pow(k_[faceCelli], 1.5) /( kappa_.value()*y_[faceCelli] *(1.0 - exp(-Aepsilon_.value()*yStar_[faceCelli])) ) *exp(-Amu_.value()*sqr(yStar_[faceCelli])); // eps_ave = [(n-1)*eps_{n-1} + eps_n]/n epsilon_[faceCelli] = ((n-1)*epsilon_[faceCelli] + eps_n) / n;; } } } Jinbiao
__________________
Jinbiao |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Cyclic patch bug in version 141 patched | cosimobianchini | OpenFOAM Bugs | 10 | February 1, 2011 19:16 |
Please report this bug | egp | OpenFOAM Installation | 5 | December 8, 2006 12:56 |
bug report pro am | F.K. | Siemens | 3 | January 25, 2005 00:27 |
BÚG REPORT for CD | Piet | Siemens | 2 | September 20, 2004 02:45 |
BUG REPORT for CD | Piet | Siemens | 0 | August 18, 2004 04:24 |